Molecular Basis of Small-Molecule Binding to α-Synuclein

Simulating the anti-aggregative effect of fasudil in early dimerisation process of α-synuclein

The aggregation of the protein α-synuclein into amyloid deposits is associated with multiple neurological disorders, including Parkinson's disease. Soluble amyloid oligomers are reported to exhibit higher toxicity than insoluble amyloid fibrils, with dimers being the smallest toxic oligomer. Small molecule drugs, such as fasudil, have shown potential in targeting α-synuclein aggregation and reducing its toxicity. In this study, we use atomistic molecular dynamics simulations to demonstrate how fasudil affects the earliest stage of aggregation, namely dimerization. Our results show that the presence of fasudil reduces the propensity for intermolecular contact formation between protein chains. Consistent with previous reports, our analysis confirms that fasudil predominantly interacts with the negatively charged C-terminal region of α-synuclein. However, we also observe transient interactions with residues in the charged N-terminal and hydrophobic NAC regions. Our simulations indicate that while fasudil prominently interacts with the C-terminal region, it is the transient interactions with residues in the N-terminal and NAC regions that effectively block the formation of intermolecular contacts between protein chains and prevent early dimerization of this disordered protein.

Exploitation of the nitro- and/or 4-Trifluoromethyl-decorated phenyl fragment to develop small inhibitors of Alpha-Syn fibril aggregation

Here, we report new 2-nitro and/or 4-trifluoromethylphenyl-based small molecules developed as inhibitors of alpha-Syn fibril formation. The set of eighteen compounds was inspired by well-known alpha-Syn aggregation modulators retrieved from literature. The preliminary biochemical data suggested that the two molecules out of eighteen compounds exerted activity comparable to that of reference compound SynuClean-D (SC-D, 5-nitro-6-(3-nitrophenyl)-2-oxo-4-(trifluoromethyl)-1H-pyridine-3-carbonitrile), according to Thioflavin T kinetics. Pharmacophore modelling deciphered the main structural requirements for alpha-Syn aggregation modulators. Moreover, docking and molecular dynamics simulations depicted the binding mode with the targeted alpha-Syn fibrils. The structural data of these new potential α-Syn binders might furnish additional information for understanding the mechanism of action of the ligands that specifically target the NAC domain as theranostic agents for α-synucleopathies.

Biomolecular simulations at the exascale: From drug design to organelles and beyond

The rapid advancement in computational power available for research offers to bring not only quantitative improvements, but also qualitative changes in the field of biomolecular simulation. Here, we review the state of biomolecular dynamics simulations at the threshold to exascale resources becoming available. Both developments in parallel and distributed computing will be discussed, providing a perspective on the state of the art of both. A main focus will be on obtaining binding and conformational free energies, with an outlook to macromolecular complexes and (sub)cellular assemblies.

Exploring high hydrostatic pressure effects on anthocyanin binding to serum albumin and food-derived transferrins

This study investigated the effects of high hydrostatic pressure (HHP) on the binding interactions of cyanindin-3-O-glucoside (C3G) to bovine serum albumin, human serum albumin (HSA), bovine lactoferrin, and ovotransferrin. Fluorescence quenching revealed that HHP reduced C3G-binding affinity to HSA, while having a largely unaffected role for the other proteins. Notably, pretreating HSA at 500 MPa significantly increased its dissociation constant with C3G from 24.7 to 34.3 μM. Spectroscopic techniques suggested that HSA underwent relatively pronounced tertiary structural alterations after HHP treatments. The C3G-HSA binding mechanisms under pressure were further analyzed through molecular dynamics simulation. The localized structural changes in HSA under pressure might weaken its interaction with C3G, particularly polar interactions such as hydrogen bonds and electrostatic forces, consequently leading to a decreased binding affinity. Overall, the importance of pressure-induced structural alterations in proteins influencing their binding with anthocyanins was highlighted, contributing to optimizing HHP processing for anthocyanin-based products.

Polymers for Disrupting Protein-Protein Interactions: Where Are We and Where Should We Be?

Protein-protein interactions (PPIs) are central to the cellular signaling and regulatory networks that underlie many physiological and pathophysiological processes. It is challenging to target PPIs using traditional small molecule or peptide-based approaches due to the frequent lack of well-defined binding pockets at the large and flat PPI interfaces. Synthetic polymers offer an opportunity to circumvent these challenges by providing unparalleled flexibility in tuning their physiochemical properties to achieve the desired binding properties. In this review, we summarize the current state of the field pertaining to polymer-protein interactions in solution, highlighting various polyelectrolyte systems, their tunable parameters, and their characterization. We provide an outlook on how these architectures can be improved by incorporating sequence control, foldability, and machine learning to mimic proteins at every structural level. Advances in these directions will enable the design of more specific protein-binding polymers and provide an effective strategy for targeting dynamic proteins, such as intrinsically disordered proteins.

Investigating the Mechanisms of Antibody Binding to Alpha-Synuclein for the Treatment of Parkinson's Disease

Parkinson's disease (PD) is an idiopathic neurodegenerative disorder with the second-highest prevalence rate behind Alzheimer's disease. The pathophysiological hallmarks of PD are both degeneration of dopaminergic neurons in the substantia nigra pars compacta and the inclusion of misfolded α-synuclein (α-syn) aggregates known as Lewy bodies. Despite decades of research for potential PD treatments, none have been developed, and developing new therapeutic agents is a time-consuming and expensive process. Computational methods can be used to investigate the properties of drug candidates currently undergoing clinical trials to determine their theoretical efficiency at targeting α-syn. Monoclonal antibodies (mAbs) are biological drugs with high specificity, and Prasinezumab (PRX002) is an mAb currently in Phase II, which targets the C-terminus (AA 118-126) of α-syn. We utilized BioLuminate and PyMol for the structure prediction and preparation of the fragment antigen-binding (Fab) region of PRX002 and 34 different conformations of α-syn. Protein-protein docking simulations were performed using PIPER, and 3 of the docking poses were selected based on the best fit. Molecular dynamics simulations were conducted on the docked protein structures in triplicate for 1000 ns, and hydrogen bonds and electrostatic and hydrophobic interactions were analyzed using MDAnalysis to determine which residues were interacting and how often. Hydrogen bonds were shown to form frequently between the HCDR2 region of PRX002 and α-syn. Free energy was calculated to determine the binding affinity. The predicted binding affinity shows a strong antibody-antigen attraction between PRX002 and α-syn. RMSD was calculated to determine the conformational change of these regions throughout the simulation. The mAb's developability was determined using computational screening methods. Our results demonstrate the efficiency and developability of this therapeutic agent.

Exploring the molecular mechanism of ginseng against anthracycline-induced cardiotoxicity based on network pharmacology, molecular docking and molecular dynamics simulation

BACKGROUND

Previous clinical and basic studies have revealed that ginseng might have cardioprotective properties against anthracycline-induced cardiotoxicity (AIC). However, the underlying mechanism of ginseng action against AIC remains insufficiently understood. The aim of this study was to explore the related targets and pathways of ginseng against AIC using network pharmacology, molecular docking, cellular thermal shift assay (CETSA) and molecular dynamics (MD) simulations.

RESULTS

Fourteen drug-disease common targets were identified. Enrichment analysis showed that the AGE-RAGE in diabetic complications, fluid shear stress and atherosclerosis, and TNF signaling pathway were potentially involved in the action of ginseng against AIC. Molecular docking demonstrated that the core components including Kaempferol, beta-Sitosterol, and Fumarine had notable binding activity with the three core targets CCNA2, STAT1, and ICAM1. Furthermore, the stable complex of STAT1 and Kaempferol with favorable affinity was further confirmed by CETSA and MD simulation.

CONCLUSIONS

This study suggested that ginseng might exert their protective effects against AIC through the derived effector compounds beta-Sitosterol, Kaempferol and Fumarine by targeting CCNA2, STAT1, and ICAM1, and modulating AGE-RAGE in diabetic complications, fluid shear stress and atherosclerosis, and TNF signaling pathways.

Predicting Conformational Ensembles of Intrinsically Disordered Proteins: From Molecular Dynamics to Machine Learning

Intrinsically disordered proteins and regions (IDP/IDRs) are ubiquitous across all domains of life. Characterized by a lack of a stable tertiary structure, IDP/IDRs populate a diverse set of transiently formed structural states that can promiscuously adapt upon binding with specific interaction partners and/or certain alterations in environmental conditions. This malleability is foundational for their role as tunable interaction hubs in core cellular processes such as signaling, transcription, and translation. Tracing the conformational ensemble of an IDP/IDR and its perturbation in response to regulatory cues is thus paramount for illuminating its function. However, the conformational heterogeneity of IDP/IDRs poses several challenges. Here, we review experimental and computational methods devised to disentangle the conformational landscape of IDP/IDRs, highlighting recent computational advances that permit proteome-wide scans of IDP/IDRs conformations. We briefly evaluate selected computational methods using the disordered N-terminal of the human copper transporter 1 as a test case and outline further challenges in IDP/IDRs ensemble prediction.

Predicting the Sequence-Dependent Backbone Dynamics of Intrinsically Disordered Proteins

How the sequences of intrinsically disordered proteins (IDPs) code for functions is still an enigma. Dynamics, in particular residue-specific dynamics, holds crucial clues. Enormous efforts have been spent to characterize residue-specific dynamics of IDPs, mainly through NMR spin relaxation experiments. Here we present a sequence-based method, SeqDYN, for predicting residue-specific backbone dynamics of IDPs. SeqDYN employs a mathematical model with 21 parameters: one is a correlation length and 20 are the contributions of the amino acids to slow dynamics. Training on a set of 45 IDPs reveals aromatic, Arg, and long-branched aliphatic amino acids as the most active in slow dynamics whereas Gly and short polar amino acids as the least active. SeqDYN predictions not only provide an accurate and insightful characterization of sequence-dependent IDP dynamics but may also serve as indicators in a host of biophysical processes, including the propensities of IDP sequences to undergo phase separation.

Amyloid Accelerator Polyphosphate Implicated as the Mystery Density in α-Synuclein Fibrils

Aberrant aggregation of α-Synuclein is the pathological hallmark of a set of neurodegenerative diseases termed synucleinopathies. Recent advances in cryo-electron microscopy have led to the structural determination of the first synucleinopathy-derived α-Synuclein fibrils, which contain a non-proteinaceous, "mystery density" at the core of the protofilaments, hypothesized to be highly negatively charged. Guided by previous studies that demonstrated that polyphosphate (polyP), a universally conserved polyanion, significantly accelerates α-Synuclein fibril formation, we conducted blind docking and molecular dynamics simulation experiments to model the polyP binding site in α-Synuclein fibrils. Here we demonstrate that our models uniformly place polyP into the lysine-rich pocket, which coordinates the mystery density in patient-derived fibrils. Subsequent in vitro studies and experiments in cells revealed that substitution of the two critical lysine residues K43 and K45 leads to a loss of all previously reported effects of polyP binding on α-Synuclein, including stimulation of fibril formation, change in filament conformation and stability as well as alleviation of cytotoxicity. In summary, our study demonstrates that polyP fits the unknown electron density present in in vivo α-Synuclein fibrils and suggests that polyP exerts its functions by neutralizing charge repulsion between neighboring lysine residues.

In situ single-molecule investigations of the impacts of biochemical perturbations on conformational intermediates of monomeric α-synuclein

α-Synuclein aggregation is a common trait in synucleinopathies, including Parkinson's disease. Being an unstructured protein, α-synuclein exists in several distinct conformational intermediates, contributing to both its function and pathogenesis. However, the regulation of these monomer conformations by biochemical factors and potential drugs has remained elusive. In this study, we devised an in situ single-molecule manipulation approach to pinpoint kinetically stable conformational intermediates of monomeric α-synuclein and explore the effects of various biochemical factors and drugs. We uncovered a partially folded conformation located in the non-amyloid-β component (NAC) region of monomeric α-synuclein, which is regulated by a preNAC region. This conformational intermediate is sensitive to biochemical perturbations and small-molecule drugs that influencing α-synuclein's aggregation tendency. Our findings reveal that this partially folded intermediate may play a role in α-synuclein aggregation, offering fresh perspectives for potential treatments aimed at the initial stage of higher-order α-synuclein aggregation. The single-molecule approach developed here can be broadly applied to the study of disease-related intrinsically disordered proteins.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

Picosecond Dynamics of a Small Molecule in Its Bound State with an Intrinsically Disordered Protein

Intrinsically disordered proteins (IDPs) are highly dynamic biomolecules that rapidly interconvert among many structural conformations. These dynamic biomolecules are involved in cancers, neurodegeneration, cardiovascular illnesses, and viral infections. Despite their enormous therapeutic potential, IDPs have generally been considered undruggable because of their lack of classical long-lived binding pockets for small molecules. Currently, only a few instances are known where small molecules have been observed to interact with IDPs, and this situation is further exacerbated by the limited sensitivity of experimental techniques to detect such binding events. Here, using experimental nuclear magnetic resonance (NMR) spectroscopy 19F transverse spin-relaxation measurements, we discovered that a small molecule, 5-fluoroindole, interacts with the disordered domains of non-structural protein 5A from hepatitis C virus with a Kd of 260 ± 110 μM. Our analysis also allowed us to determine the rotational correlation times (τc) for the free and bound states of 5-fluoroindole. In the free state, we observed a rotational correlation time of 27.0 ± 1.3 ps, whereas in the bound state, τc only increased to 46 ± 10 ps. Our findings imply that it is possible for small molecules to engage with IDPs in exceptionally dynamic ways, in sharp contrast to the rigid binding modes typically exhibited when small molecules bind to well-defined binding pockets within structured proteins.

Multi-eGO: Model Improvements toward the Study of Complex Self-Assembly Processes

Structure-based models have been instrumental in simulating protein folding and suggesting hypotheses about the mechanisms involved. Nowadays, at least for fast-folding proteins, folding can be simulated in explicit solvent using classical molecular dynamics. However, other self-assembly processes, such as protein aggregation, are still far from being accessible. Recently, we proposed that a hybrid multistate structure-based model, multi-eGO, could help to bridge the gap toward the simulation of out-of-equilibrium, concentration-dependent self-assembly processes. Here, we further improve the model and show how multi-eGO can effectively and accurately learn the conformational ensemble of the amyloid β42 intrinsically disordered peptide, reproduce the well-established folding mechanism of the B1 immunoglobulin-binding domain of streptococcal protein G, and reproduce the aggregation as a function of the concentration of the transthyretin 105-115 amyloidogenic peptide. We envision that by learning from the dynamics of a few minima, multi-eGO can become a platform for simulating processes inaccessible to other simulation techniques.

1,4-Diurea- and 1,4-Dithiourea-Substituted Aromatic Derivatives Selectively Inhibit α-Synuclein Oligomer Formation In Vitro

Parkinson's disease (PD) is the second most common neurodegenerative disease, affecting the elderly population worldwide. In PD, the misfolding of α-synuclein (α-syn) results in the formation of inclusions referred to as Lewy bodies (LB) in midbrain neurons of the substantia nigra and other specific brain localizations, which is associated with neurodegeneration. There are no approved strategies to reduce the formation of LB in the neurons of patients with PD. Our drug discovery program focuses on the synthesis of urea and thiourea compounds coupled with aminoindole moieties to abrogate α-syn aggregation and to slow down the progression of PD. We synthesized several urea and thiourea analogues with a central 1,4-phenyl diurea/thiourea linkage and evaluated their effectiveness in reducing α-syn aggregation with a special focus on the selective inhibition of oligomer formation among other proteins. We utilized biophysical methods such as thioflavin T (ThT) fluorescence assays, transmission electron microscopy (TEM), photoinduced cross-linking of unmodified proteins (PICUP), as well as M17D intracellular inclusion cell-based assays to evaluate the antiaggregation properties and cellular protection of our best compounds. Our results identified compound 1 as the best compound in reducing α-syn fibril formation via ThT assays. The antioligomer formation of compound 1 was subsequently superseded by compound 2. Both compounds selectively curtailed the oligomer formation of α-syn but not tau 4R isoforms (0N4R, 2N4R) or p-tau (isoform 1N4R). Compounds 1 and 2 failed to abrogate tau 0N3R fibril formation by ThT and atomic force microscopy. Compound 2 was best at reducing the formation of recombinant α-syn fibrils by TEM. In contrast to compound 2, compound 1 reduced the formation of α-syn inclusions in M17D neuroblastoma cells in a dose-dependent manner. Compound 1 may provide molecular scaffolds for the optimization of symmetric molecules for its α-syn antiaggregation activity with potential therapeutic applications and development of small molecules in PD.

Sensitized Emission Imaging Allows Nanoscale Surface Polarity Mapping of α-Synuclein Amyloid Fibrils

When misfolded, α-Synuclein (α-Syn), a natively disordered protein, aggregates to form amyloid fibrils responsible for the neurodegeneration observed in Parkinson's disease. Structural studies revealed distinct molecular packing of α-Syn in different fibril polymorphs and variations of interprotofilament connections in the fibrillar architecture. Fibril polymorphs have been hypothesized to exhibit diverse surface polarities depending on the folding state of the protein during aggregation; however, the spatial variation of surface polarity in amyloid fibrils remains unexplored. To map the local polarity (or hydrophobicity) along α-Syn fibrils, we visualized the spectral characteristics of two dyes with distinct polarities-hydrophilic Thioflavin T (ThT) and hydrophobic Nile red (NR)─when both are bound to α-Syn fibrils. Dual-channel fluorescence imaging reveals uneven partitioning of ThT and NR along individual fibrils, implying that relatively more polar/hydrophobic patches are spread over a few hundred nanometers. Remarkably, spectrally resolved sensitized emission imaging of α-Syn fibrils provides unambiguous evidence of energy transfer from ThT to NR, implying that dyes of dissimilar polarity are in close proximity. Furthermore, spatially resolved fluorescence spectroscopy of the solvatochromic probe NR allowed us to quantitatively map the range and variation of the polarity parameter ET30 along individual fibrils. Our results suggest the existence of interlaced polar and nonpolar nanoscale domains throughout the fibrils; however, the relative populations of these patches vary considerably over larger length scales likely due to heterogeneous packing of α-Syn during fibrilization and dissimilar exposed polarities of polymorphic segments. The employed method may provide a foundation for imaging modalities of other similar structurally unresolved systems with diverse hydrophobic-hydrophilic topology.

Biomolecular NMR spectroscopy in the era of artificial intelligence

Biomolecular nuclear magnetic resonance (NMR) spectroscopy and artificial intelligence (AI) have a burgeoning synergy. Deep learning-based structural predictors have forever changed structural biology, yet these tools currently face limitations in accurately characterizing protein dynamics, allostery, and conformational heterogeneity. We begin by highlighting the unique abilities of biomolecular NMR spectroscopy to complement AI-based structural predictions toward addressing these knowledge gaps. We then highlight the direct integration of deep learning approaches into biomolecular NMR methods. AI-based tools can dramatically improve the acquisition and analysis of NMR spectra, enhancing the accuracy and reliability of NMR measurements, thus streamlining experimental processes. Additionally, deep learning enables the development of novel types of NMR experiments that were previously unattainable, expanding the scope and potential of biomolecular NMR spectroscopy. Ultimately, a combination of AI and NMR promises to further revolutionize structural biology on several levels, advance our understanding of complex biomolecular systems, and accelerate drug discovery efforts.

Effect of flavonoids on the destabilization of α-synuclein fibrils and their conversion to amorphous aggregate: A molecular dynamics simulation and experimental study

The second most prevalent neurodegenerative disease, Parkinson's disease (PD), is caused by the accumulation and deposition of fibrillar aggregates of the α-Syn into the Lewy bodies. To create a potent pharmacological candidate to destabilize the preformed α-Syn fibril, it is important to understand the precise molecular mechanism underlying the destabilization of the α-Syn fibril. Through molecular dynamics simulations and experiments, we have examined the molecular mechanisms causing the destabilization and suppression of a newly discovered α-Syn fibril with a Greek-key-like shape and an aggregation prone state (APS) of α-Syn in the presence and absence of various Flvs. According to MD simulation and experimental evidence, morin, quercetin, and myricetin are the Flvs, most capable of destabilizing the fibrils and converting them into amorphous aggregates. Compared to galangin and kaempferol, they have more hydroxyl groups and form more hydrogen bonds with fibrils.The processes by which morin and myricetin prevent new fibril production from APS and destabilize the fibrils are different. According to linear interaction energy analysis, van der Waals interaction predominates with morin, and electrostatic interaction dominates with myricetin. Our MD simulation and experimental findings provide mechanistic insights into how Flvs destabilize α-Syn fibrils and change their morphology, opening the door to developing structure-based drugs for treating Parkinson's disease.

Exploring the Curvature-Dependence of Boron Nitride Nanoparticles on the Inhibition of hIAPP Aggregation

Nanoparticles, particularly carbon nanoparticles, have gathered significant interest in the field of anti-aggregation research. However, due to their cytotoxicity, the exploration of biocompatible nanoparticles has become a new frontier in the quest for drugs against human amyloid diseases. The application of non-cytotoxic and biocompatible boron nitride (BN) nanoparticles against amyloid aggregation has been probed to tackle this issue. BN nanoparticles displayed inhibitory activity against the aggregation of Aβ and α-syn peptides. In this work, the effect of BN nanoparticles on the dimerization of hIAPP, which is associated with the pathogenesis of type 2 diabetes, is studied. BN nanoparticles prevent the misfolding of hIAPP into β-sheet-rich aggregates. On varying the curvature, the nanoparticles display variation in the interaction preference with hIAPP. Interestingly, as the hydrophobicity of the nanoparticles increases from (5,5) BN nanotube to BN nanosheet, the interaction propensity shifts from N-terminal to the amyloid prone C-terminal of hIAPP. The hydrophobic and aromatic stacking interactions are a contributing factor toward the binding between hIAPP and BN. Due to this, the flat surface of the nanosheet shows better interaction potential toward hIAPP, compared to the nanotubes. Further, the nanoparticles can also disassemble preformed hIAPP fibrils, and the effect is more pronounced for (5,5) nanotube and the nanosheet. This study provides insight into the inhibitory mechanism of hIAPP aggregation by boron nitride nanoparticles and also an understanding of the significance of the curvature of nanoparticles in their interaction with amyloid peptides, which is valuable for the design of antiamyloid drugs.

Cross-Linked α-Synuclein as Inhibitor of Amyloid Formation

The aggregation and amyloid formation of α-synuclein is associated with Parkinson's disease and other synucleinopathies. In its native, monomeric form α-synuclein is an intrinsically disordered protein represented by highly dynamic conformational ensembles. Inhibition of α-synuclein aggregation using small molecules, peptides, or proteins has been at the center of interest in recent years. Our aim was to explore the effects of cross-linking on the structure and aggregation/amyloid formation properties of α-synuclein. Comparative analysis of available high-resolution amyloid structures and representative structural models and MD trajectory of monomeric α-synuclein revealed that potential cross-links in the monomeric protein are mostly incompatible with the amyloid forms and thus might inhibit fibrillation. Monomeric α-synuclein has been intramolecularly chemically cross-linked under various conditions using different cross-linkers. We determined the location of cross-links and their frequency using mass spectrometry and found that most of them cannot be realized in the amyloid structures. The inhibitory potential of cross-linked proteins has been experimentally investigated using various methods, including thioflavin-T fluorescence and transmission electron microscopy. We found that conformational constraints applied by cross-linking fully blocked α-synuclein amyloid formation. Moreover, DTSSP-cross-linked molecules exhibited an inhibitory effect on the aggregation of unmodified α-synuclein as well.

Origin of stronger binding of ionic pair (IP) inhibitor to Aβ42 than the equimolar neutral counterparts: synergy mechanism of IP in disrupting Aβ42 protofibril and inhibiting Aβ42 aggregation under two pH conditions

Fibrous aggregates of beta-amyloid (Aβ) is a hallmark of Alzheimer's disease (AD). Several major strategies of drugs or inhibitors, including neutral molecules, positive or negative ions, and dual-inhibitor, are used to inhibit the misfolding or aggregation of Aβ42, among which a kind of dual-inhibitor composed of a pair of positive and negative ions is emerging as the most powerful candidate. This knowledge lacks the origin of the strong inhibitory effect and synergy mechanisms blocking the development and application of such inhibitors. To this end, we employed 1 : 1 ionic pairs (IP) of oppositely charged benzothiazole molecules (+)BAM1-EG6 (Pos) and (-)BAM1-EG6 (Neg) as well as equimolar neutral BAM1-EG6 (Neu) counterpart at two pH conditions (5.5 and 7.0) to bind Aβ42 targets, Aβ42 monomer (AβM), soluble pentamer (AβP), and pentameric protofibril (AβF) models, respectively, corresponding to the products of three toxic Aβ42 development pathways, lag, exponential and fibrillation phases. Simulated results illustrated the details of the inhibitory mechanisms of IP and Neu for the AβY (Y = M, P, or F) in the three different phases, characterizing the roles of Pos and Neg of IP as well as their charged, hydrophobic groups and linker playing in the synergistic interaction, and elucidated a previously unknown molecular mechanism governing the IP-Aβ42 interaction. Most importantly, we first revealed the origin of the stronger binding of IP inhibitors to Aβ42 than that of the equimolar neutral counterparts, observing a perplexing phenomenon that the physiological condition (pH = 7.0) than the acidic one (pH = 5.5) is more favorable to the enhancement of IP binding, and finally disclosed that solvation is responsible to the enhancement because at pH 7.0, AβP and AβF act as anionic membranes, where solvation plays a critical role in the chemoelectromechanics. The result not only provides a new dimension in dual-inhibitor/drug design and development but also a new perspective for uncovering charged protein disaggregation under IP-like inhibitors.

Clustering Heterogeneous Conformational Ensembles of Intrinsically Disordered Proteins with t-Distributed Stochastic Neighbor Embedding

Intrinsically disordered proteins (IDPs) populate a range of conformations that are best described by a heterogeneous ensemble. Grouping an IDP ensemble into "structurally similar" clusters for visualization, interpretation, and analysis purposes is a much-desired but formidable task, as the conformational space of IDPs is inherently high-dimensional and reduction techniques often result in ambiguous classifications. Here, we employ the t-distributed stochastic neighbor embedding (t-SNE) technique to generate homogeneous clusters of IDP conformations from the full heterogeneous ensemble. We illustrate the utility of t-SNE by clustering conformations of two disordered proteins, Aβ42, and α-synuclein, in their APO states and when bound to small molecule ligands. Our results shed light on ordered substates within disordered ensembles and provide structural and mechanistic insights into binding modes that confer specificity and affinity in IDP ligand binding. t-SNE projections preserve the local neighborhood information, provide interpretable visualizations of the conformational heterogeneity within each ensemble, and enable the quantification of cluster populations and their relative shifts upon ligand binding. Our approach provides a new framework for detailed investigations of the thermodynamics and kinetics of IDP ligand binding and will aid rational drug design for IDPs.

Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model

A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long-time scale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We do not, however, find evidence for the existence of canonical conformational selection or induced fit binding pathways. We observe four kinetically separated native-like bound states that interconvert on time scales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.

Thermodynamic and kinetic approaches for drug discovery to target protein misfolding and aggregation

INTRODUCTION

Protein misfolding diseases, including Alzheimer's and Parkinson's diseases, are characterized by the aberrant aggregation of proteins. These conditions are still largely untreatable, despite having a major impact on our healthcare systems and societies.

AREAS COVERED

We describe drug discovery strategies to target protein misfolding and aggregation. We compare thermodynamic approaches, which are based on the stabilization of the native states of proteins, with kinetic approaches, which are based on the slowing down of the aggregation process. This comparison is carried out in terms of the current knowledge of the process of protein misfolding and aggregation, the mechanisms of disease and the therapeutic targets.

EXPERT OPINION

There is an unmet need for disease-modifying treatments that target protein misfolding and aggregation for the over 50 human disorders known to be associated with this phenomenon. With the approval of the first drugs that can prevent misfolding or inhibit aggregation, future efforts will be focused on the discovery of effective compounds with these mechanisms of action for a wide range of conditions.

Understanding the pH-dependent interaction of anthocyanin with two food-derived transferrins

The potential binding of cyanidin-3-O-glucoside (C3G) to bovine lactoferrin (BLF) and ovotransferrin (OTF) at pH 3, 5, and 7 was investigated for the first time. Multiple spectroscopic techniques demonstrated pH-dependent alterations in the conformational characteristics of BLF and OTF upon complexation with C3G. Fluorescence quenching assays showed that their highest binding affinity was at pH 7. Hydrophobic interactions and hydrogen bonds were found to be crucial in molecular dynamics simulations but with significantly lower probabilities of formation at pH 3 (p < 0.05). At pH 7, electrostatic attraction can occur for the negatively charged forms of C3G, and the well-maintained native structures of BLF and OTF may be favorable for stabilizing the C3G binding sites. This study sheds light on the stronger interaction of C3G with BLF/OTF at pH 7, which may have implications for future applications such as anthocyanin stabilization or the development of functional food ingredients.

Conformational Plasticity in α-Synuclein and How Crowded Environment Modulates It

A 140-residue intrinsically disordered protein (IDP), α-synuclein (αS), is known to adopt conformations that are vastly plastic and susceptible to environmental cues and crowders. However, the inherently heterogeneous nature of αS has precluded a clear demarcation of its monomeric precursor between aggregation-prone and functionally relevant aggregation-resistant states and how a crowded environment could modulate their mutual dynamic equilibrium. Here, we identify an optimal set of distinct metastable states of αS in aqueous media by dissecting a 73 μs-long molecular dynamics ensemble via building a comprehensive Markov state model (MSM). Notably, the most populated metastable state corroborates with the dimension obtained from PRE-NMR studies of αS monomer, and it undergoes kinetic transition at diverse time scales with a weakly populated random-coil-like ensemble and a globular protein-like state. However, subjecting αS to a crowded environment results in a nonmonotonic compaction of these metastable conformations, thereby skewing the ensemble by either introducing new tertiary contacts or by reinforcing the innate contacts. The early stage of dimerization process is found to be considerably expedited in the presence of crowders, albeit promoting nonspecific interactions. Together with this, using an extensively sampled ensemble of αS, this exposition demonstrates that crowded environments can potentially modulate the conformational preferences of IDP that can either promote or inhibit aggregation events.

Development of Small Molecules Targeting α-Synuclein Aggregation: A Promising Strategy to Treat Parkinson's Disease

Parkinson's disease, the second most common neurodegenerative disorder worldwide, is characterized by the accumulation of protein deposits in the dopaminergic neurons. These deposits are primarily composed of aggregated forms of α-Synuclein (α-Syn). Despite the extensive research on this disease, only symptomatic treatments are currently available. However, in recent years, several compounds, mainly of an aromatic character, targeting α-Syn self-assembly and amyloid formation have been identified. These compounds, discovered by different approaches, are chemically diverse and exhibit a plethora of mechanisms of action. This work aims to provide a historical overview of the physiopathology and molecular aspects associated with Parkinson's disease and the current trends in small compound development to target α-Syn aggregation. Although these molecules are still under development, they constitute an important step toward discovering effective anti-aggregational therapies for Parkinson's disease.

Structural ensembles of disordered proteins from hierarchical chain growth and simulation

Disordered proteins and nucleic acids play key roles in cellular function and disease. Here, we review recent advances in the computational exploration of the conformational dynamics of flexible biomolecules. While atomistic molecular dynamics (MD) simulation has seen a lot of improvement in recent years, large-scale computing resources and careful validation are required to simulate full-length disordered biopolymers in solution. As a computationally efficient alternative, hierarchical chain growth (HCG) combines pre-sampled chain fragments in a statistically reproducible manner into ensembles of full-length atomically detailed biomolecular structures. Experimental data can be integrated during and after chain assembly. Applications to the neurodegeneration-linked proteins α-synuclein, tau, and TDP-43, including as condensate, illustrate the use of HCG. We conclude by highlighting the emerging connections to AI-based structural modeling including AlphaFold2.

Functional and Pathological Effects of α-Synuclein on Synaptic SNARE Complexes

α-Synuclein is an abundant protein at the neuronal synapse that has been implicated in Parkinson's disease for over 25 years and characterizes the hallmark pathology of a group of neurodegenerative diseases now known as the synucleinopathies. Physiologically, α-synuclein exists in an equilibrium between a synaptic vesicle membrane-bound α-helical multimer and a cytosolic largely unstructured monomer. Through its membrane-bound state, α-synuclein functions in neurotransmitter release by modulating several steps in the synaptic vesicle cycle, including synaptic vesicle clustering and docking, SNARE complex assembly, and homeostasis of synaptic vesicle pools. These functions have been ascribed to α-synuclein's interactions with the synaptic vesicle SNARE protein VAMP2/synaptobrevin-2, the synaptic vesicle-attached synapsins, and the synaptic vesicle membrane itself. How α-synuclein affects these processes, and whether disease is due to loss-of-function or gain-of-toxic-function of α-synuclein remains unclear. In this review, we provide an in-depth summary of the existing literature, discuss possible reasons for the discrepancies in the field, and propose a working model that reconciles the findings in the literature.

Controlling amyloid formation of intrinsically disordered proteins and peptides: slowing down or speeding up?

The pathological assembly of intrinsically disordered proteins/peptides (IDPs) into amyloid fibrils is associated with a range of human pathologies, including neurodegeneration, metabolic diseases and systemic amyloidosis. These debilitating disorders affect hundreds of millions of people worldwide, and the number of people affected is increasing sharply. However, the discovery of therapeutic agents has been immensely challenging largely because of (i) the diverse number of aggregation pathways and the multi-conformational and transient nature of the related proteins or peptides and (ii) the under-development of experimental pipelines for the identification of disease-modifying molecules and their mode-of-action. Here, we describe current approaches used in the search for small-molecule modulators able to control or arrest amyloid formation commencing from IDPs and review recently reported accelerators and inhibitors of amyloid formation for this class of proteins. We compare their targets, mode-of-action and effects on amyloid-associated cytotoxicity. Recent successes in the control of IDP-associated amyloid formation using small molecules highlight exciting possibilities for future intervention in protein-misfolding diseases, despite the challenges of targeting these highly dynamic precursors of amyloid assembly.

Chemical tools for study and modulation of biomolecular phase transitions

Biomolecular phase transitions play an important role in organizing cellular processes in space and time. Methods and tools for studying these transitions, and the intrinsically disordered proteins (IDPs) that often drive them, are typically less developed than tools for studying their folded protein counterparts. In this perspective, we assess the current landscape of chemical tools for studying IDPs, with a specific focus on protein liquid-liquid phase separation (LLPS). We highlight methodologies that enable imaging and spectroscopic studies of these systems, including site-specific labeling with small molecules and the diverse range of capabilities offered by inteins and protein semisynthesis. We discuss strategies for introducing post-translational modifications that are central to IDP and LLPS function and regulation. We also investigate the nascent field of noncovalent small-molecule modulators of LLPS. We hope that this review of the state-of-the-art in chemical tools for interrogating IDPs and LLPS, along with an associated perspective on areas of unmet need, can serve as a valuable and timely resource for these rapidly expanding fields of study.

Unveiling the Metal-Dependent Aggregation Properties of the C-terminal Region of Amyloidogenic Intrinsically Disordered Protein Isoforms DPF3b and DPF3a

Double-PHD fingers 3 (DPF3) is a BAF-associated human epigenetic regulator, which is increasingly recognised as a major contributor to various pathological contexts, such as cardiac defects, cancer, and neurodegenerative diseases. Recently, we unveiled that its two isoforms (DPF3b and DPF3a) are amyloidogenic intrinsically disordered proteins. DPF3 isoforms differ from their C-terminal region (C-TERb and C-TERa), containing zinc fingers and disordered domains. Herein, we investigated the disorder aggregation properties of C-TER isoforms. In agreement with the predictions, spectroscopy highlighted a lack of a highly ordered structure, especially for C-TERa. Over a few days, both C-TERs were shown to spontaneously assemble into similar antiparallel and parallel β-sheet-rich fibrils. Altered metal homeostasis being a neurodegeneration hallmark, we also assessed the influence of divalent metal cations, namely Cu²⁺, Mg²⁺, Ni²⁺, and Zn²⁺, on the C-TER aggregation pathway. Circular dichroism revealed that metal binding does not impair the formation of β-sheets, though metal-specific tertiary structure modifications were observed. Through intrinsic and extrinsic fluorescence, we found that metal cations differently affect C-TERb and C-TERa. Cu²⁺ and Ni²⁺ have a strong inhibitory effect on the aggregation of both isoforms, whereas Mg²⁺ impedes C-TERb fibrillation and, on the contrary, enhances that of C-TERa. Upon Zn²⁺ binding, C-TERb aggregation is also hindered, and the amyloid autofluorescence of C-TERa is remarkably red-shifted. Using electron microscopy, we confirmed that the metal-induced spectral changes are related to the morphological diversity of the aggregates. While metal-treated C-TERb formed breakable and fragmented filaments, C-TERa fibrils retained their flexibility and packing properties in the presence of Mg²⁺ and Zn²⁺ cations.

Small molecules targeting the disordered transactivation domain of the androgen receptor induce the formation of collapsed helical states

Intrinsically disordered proteins, which do not adopt well-defined structures under physiological conditions, are implicated in many human diseases. Small molecules that target the disordered transactivation domain of the androgen receptor have entered human trials for the treatment of castration-resistant prostate cancer (CRPC), but no structural or mechanistic rationale exists to explain their inhibition mechanisms or relative potencies. Here, we utilize all-atom molecular dynamics computer simulations to elucidate atomically detailed binding mechanisms of the compounds EPI-002 and EPI-7170 to the androgen receptor. Our simulations reveal that both compounds bind at the interface of two transiently helical regions and induce the formation of partially folded collapsed helical states. We find that EPI-7170 binds androgen receptor more tightly than EPI-002 and we identify a network of intermolecular interactions that drives higher affinity binding. Our results suggest strategies for developing more potent androgen receptor inhibitors and general strategies for disordered protein drug design.

REMD Simulations of Full-Length Alpha-Synuclein Together with Ligands Reveal Binding Region and Effect on Amyloid Conversion

Alpha-synuclein is a key protein involved in the development and progression of Parkinson’s disease and other synucleinopathies. The intrinsically disordered nature of alpha-synuclein hinders the computational screening of new drug candidates for the treatment of these neurodegenerative diseases. In the present work, replica exchange molecular dynamics simulations of the full-length alpha-synuclein together with low-molecular ligands were utilized to predict the binding site and effect on the amyloid-like conversion of the protein. This approach enabled an accurate prediction of the binding sites for three tested compounds (fasudil, phthalocyanine tetrasulfonate, and spermine), giving good agreement with data from experiments by other groups. Lots of information about the binding and protein conformational ensemble enabled the suggestion of a putative effect of the ligands on the amyloid-like conversion of alpha-synuclein and the mechanism of anti- and pro-amyloid activity of the tested compounds. Therefore, this approach looks promising for testing new drug candidates for binding with alpha-synuclein or other intrinsically disordered proteins and at the same time the estimation of the effect on protein behavior, including amyloid-like aggregation.

Threonine Cavities Are Targetable Motifs That Control Alpha-Synuclein Fibril Growth

Recent high-resolution structures of alpha-synuclein (aSyn) fibrils offer promise for rational approaches to drug discovery for Parkinson's disease and Lewy body dementia. Harnessing the first such structures, we previously used molecular dynamics and free energy calculations to suggest that threonines 72 and 75─which line water-filled cavities within the fibril stacks─may be of central importance in stabilizing fibrils. Here, we used experimental mutagenesis of both wild-type and A53T aSyn to show that both threonine residues play important but surprisingly disparate roles in fibril nucleation and elongation. The T72A mutant, but not T75A, resulted in a large increase in the extent of fibrillization during primary nucleation, leading us to posit that T72 acts as a "brake" on run-away aggregation. An expanded set of simulations of five recent high-resolution fibril structures suggests that confinement of cavity waters around T72 correlates with this finding. In contrast, the T75A mutation led to a modest decrease in the extent of fibrillization. Furthermore, both T72A and T75A completely blocked the initial fibril elongation in seeded fibrillization. To test whether these threonine-lined cavities are druggable targets, we used computational docking to identify potential small-molecule binders. We show that the top-scoring hit, aprepitant, strongly promotes fibril growth while specifically interacting with aSyn fibrils and not monomer, and we offer speculation as to how such compounds could be used therapeutically.

Baicalein exhibits differential effects and mechanisms towards disruption of α-synuclein fibrils with different polymorphs

Parkinson's disease (PD) is the second most common neurodegenerative diseases with no cure yet and its major hallmark is α-synuclein fibrillary aggregates. The crucial role of α-synuclein aggregation in PD makes it an attractive target for potential disease-modifying therapies. Disaggregation of α-synuclein fibrils is considered as one of the promising therapeutic strategies to treat PD. The wild type (WT) and mutant α-synuclein fibrils exhibit different polymorphs and provide therapeutic targets for PD. Recent experiments reported that a flavonoid baicalein can disrupt WT α-synuclein fibrils. However, the underlying disruptive mechanism remains largely elusive, and whether BAC is capable of disrupting mutant α-synuclein fibrils is also unknown. Herein, we performed microsecond molecular dynamics simulations on cryo-EM-determined WT and two familial PD-associated mutant (E46K and H50Q) α-synuclein fibrils with and without baicalein. We find that baicalein destructs WT fibril by disrupting E46-K80 salt-bridge and β-sheets, and by remodeling the inter-protofilament interface. And baicalein can also damage E46K and H50Q mutant fibrils, but to different extents and via different mechanisms. The E46K fibril disruption is initiated from E61-K80 salt-bridge and N-terminal β-sheet, while the H50Q fibril disruption starts from the inter-protofilament interface and N-terminal β-sheet. These results reveal that disruptive effects and modes of baicalein on α-synuclein fibrils are polymorphism-dependent. This study suggests that baicalein may be a potential drug candidate to disrupt both WT and E46K/H50Q mutant α-synuclein fibrils and alleviate the pathological process of PD.

Multi-eGO: An in silico lens to look into protein aggregation kinetics at atomic resolution

Protein aggregation into amyloid fibrils is the archetype of aberrant biomolecular self-assembly processes, with more than 50 associated diseases that are mostly uncurable. Understanding aggregation mechanisms is thus of fundamental importance and goes in parallel with the structural characterization of the transient oligomers formed during the process. Oligomers have been proven elusive to high-resolution structural techniques, while the large sizes and long time scales, typical of aggregation processes, have limited the use of computational methods to date. To surmount these limitations, we here present multi-eGO, an atomistic, hybrid structure-based model which, leveraging the knowledge of monomers conformational dynamics and of fibril structures, efficiently captures the essential structural and kinetics aspects of protein aggregation. Multi-eGO molecular dynamics simulations can describe the aggregation kinetics of thousands of monomers. The concentration dependence of the simulated kinetics, as well as the structural features of the resulting fibrils, are in qualitative agreement with in vitro experiments carried out on an amyloidogenic peptide from Transthyretin, a protein responsible for one of the most common cardiac amyloidoses. Multi-eGO simulations allow the formation of primary nuclei in a sea of transient lower-order oligomers to be observed over time and at atomic resolution, following their growth and the subsequent secondary nucleation events, until the maturation of multiple fibrils is achieved. Multi-eGO, combined with the many experimental techniques deployed to study protein aggregation, can provide the structural basis needed to advance the design of molecules targeting amyloidogenic diseases.

Enhanced-Sampling Simulations for the Estimation of Ligand Binding Kinetics: Current Status and Perspective

The dissociation rate (k off) associated with ligand unbinding events from proteins is a parameter of fundamental importance in drug design. Here we review recent major advancements in molecular simulation methodologies for the prediction of k off. Next, we discuss the impact of the potential energy function models on the accuracy of calculated k off values. Finally, we provide a perspective from high-performance computing and machine learning which might help improve such predictions.

Energetics of a protein disorder-order transition in small molecule recognition

Many proteins recognise other proteins via mechanisms that involve the folding of intrinsically disordered regions upon complex formation. Here we investigate how the selectivity of a drug-like small molecule arises from its modulation of a protein disorder-to-order transition. Binding of the compound AM-7209 has been reported to confer order upon an intrinsically disordered ‘lid’ region of the oncoprotein MDM2. Calorimetric measurements revealed that truncation of the lid region of MDM2 increases the apparent dissociation constant of AM-7209 250-fold. By contrast, lid truncation has little effect on the binding of the ligand Nutlin-3a. Insights into these differential binding energetics were obtained via a complete thermodynamic analysis that featured adaptive absolute alchemical free energy of binding calculations with enhanced-sampling molecular dynamics simulations. The simulations reveal that in apo MDM2 the ordered lid state is energetically disfavoured. AM-7209, but not Nutlin-3a, shows a significant energetic preference for ordered lid conformations, thus shifting the balance towards ordering of the lid in the AM-7209/MDM2 complex. The methodology reported herein should facilitate broader targeting of intrinsically disordered regions in medicinal chemistry.

Molecular simulations and biophysical measurements elucidate why the ligand AM-7209 orders a disordered region of the protein MDM2 on binding. This work expands strategies available to medicinal chemists for targeting disordered proteins.

Tuning the rate of aggregation of hIAPP into amyloid using small-molecule modulators of assembly

Human islet amyloid polypeptide (hIAPP) self-assembles into amyloid fibrils which deposit in pancreatic islets of type 2 diabetes (T2D) patients. Here, we applied chemical kinetics to study the mechanism of amyloid assembly of wild-type hIAPP and its more amyloidogenic natural variant S20G. We show that the aggregation of both peptides involves primary nucleation, secondary nucleation and elongation. We also report the discovery of two structurally distinct small-molecule modulators of hIAPP assembly, one delaying the aggregation of wt hIAPP, but not S20G; while the other enhances the rate of aggregation of both variants at substoichiometric concentrations. Investigation into the inhibition mechanism(s) using chemical kinetics, native mass spectrometry, fluorescence titration, SPR and NMR revealed that the inhibitor retards primary nucleation, secondary nucleation and elongation, by binding peptide monomers. By contrast, the accelerator predominantly interacts with species formed in the lag phase. These compounds represent useful chemical tools to study hIAPP aggregation and may serve as promising starting-points for the development of therapeutics for T2D.