Identification of fibril-like tertiary contacts in soluble monomeric α-synuclein

Structures prediction and replica exchange molecular dynamics simulations of α-synuclein: A case study for intrinsically disordered proteins

In recent years, a variety of three-dimensional structure prediction tools, including AlphaFold2, AlphaFold3, I-TASSER, C-I-TASSER, Phyre2, ESMFold, and RoseTTAFold, have been employed in the investigation of intrinsically disordered proteins. However, a comprehensive validation of these tools specifically for intrinsically disordered proteins has yet to be conducted. In this study, we utilize AlphaFold2, AlphaFold3, I-TASSER, C-I-TASSER, Phyre2, ESMFold, and RoseTTAFold to predict the structure of a model intrinsically disordered α-synuclein protein. Additionally, extensive replica exchange molecular dynamics simulations of the intrinsically disordered protein are conducted. The resulting structures from both structure prediction tools and replica exchange molecular dynamics simulations are analyzed for radius of gyration, secondary and tertiary structure properties, as well as Cα and Hα chemical shift values. A comparison of the obtained results with experimental data reveals that replica exchange molecular dynamics simulations provide results in excellent agreement with experimental observations. However, none of the structure prediction tools utilized in this study can fully capture the structural characteristics of the model intrinsically disordered protein. This study shows that a cluster of ensembles are required for intrinsically disordered proteins. Artificial-intelligence based structure prediction tools such as AlphaFold3 and C-I-TASSER could benefit from stochastic sampling or Monte Carlo simulations for generating an ensemble of structures for intrinsically disordered proteins.

E46K Mutation of α-Synuclein Preorganizes the Intramolecular Interactions Crucial for Aggregation

Aggregation of α-synuclein is central to the pathogenesis of Parkinson's disease. The most toxic familial mutation E46K accelerates the aggregation process by an unknown mechanism. Herein, we provide a clue by investigating the influence of E46K on monomeric α-synuclein and its relation to aggregation with molecular dynamics simulations. The E46K mutation suppresses β-sheet structures in the N-terminus while promoting those at the key fibrillization region named NACore. Even though WT and E46K monomers share conserved intramolecular interactions with fibrils, E46K abolishes intramolecular contacts within the N-terminus which are present in the WT monomer but absent in fibrils. Network analysis identifies residues 36-53 as the interaction core of the WT monomer. Upon mutation, residues 36-46 are expelled to water due to aggravated electrostatic repulsion in the 43KTKK46 segment. Instead, NACore (residues 68-78) becomes the interaction hub and connects preceding residues 47-56 and the C-terminus. Consequently, residues 47-95 which belong to the fibril core form more compact β-sheets. Overall, the interaction network of E46K is more like fibrils than WT, stabilizing the fibril-like conformations. Our work provides mechanistic insights into the faster aggregation of the E46K mutant. It implies a close link between monomeric conformations and fibrils, which would spur the development of therapeutic strategies.

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.

Modulation of α-synuclein phase separation by biomolecules

Liquid-liquid phase separation (LLPS) is currently recognized as a common mechanism involved in the regulation of a number of cellular functions. On the other hand, aberrant phase separation has been linked to the biogenesis of several neurodegenerative disorders since many proteins that undergo LLPS are also found in pathological aggregates. The formation of mixed protein coacervates may constitute a risk factor in overlapping neuropathologies, such as Parkinson's (PD) and Alzheimer's (AD) diseases. In this work, we evaluated the homotypic and heterotypic phase behaviour of the PD-related protein α-synuclein (AS) in the presence of the biologically relevant molecules ATP, polyamines, and the AD-related protein Tau. We found that AS exhibits a low propensity to form homotypic liquid droplets, yet phase separates into liquid-like or solid-like phases depending on the interacting biomolecule. We further demonstrated the synergistic droplet formation of AS and Tau providing support for a mechanism in which mixed condensates might contribute to the biogenesis of AS/Tau pathologies.

Concerted enhanced-sampling simulations to elucidate the helix-fibril transition pathway of intrinsically disordered α-Synuclein

Fibril formation of α-synuclein is linked with Parkinson's disease. The intrinsically disordered nature of α-syn provides extensive conformational plasticity and becomes difficult to characterize its transition pathway from native monomeric to disease-associated fibril form. We implemented different simulation methods such as steered dynamics-umbrella sampling, and replica-exchange and conventional MD simulations to access various conformational states of α-syn. Nineteen distinct intermediate structures were identified by free energy landscape and cluster analysis. They were then sorted based on secondary structure and solvent exposure of fibril-core residues to illustrate the fibril dissociation pathway. The analysis showed that following the initial dissociation of the polypeptide chain from the fibril, α-syn might form either compact-conformations by long-range interactions or extended-conformations stabilized by local interactions. This leads α-syn to adapt two different pathways. The secondary structure, solvation, contact distance, interaction energies and backbone dihedrals of thirty-two selected residues were analyzed for all the 19 intermediates. The results suggested that formation of β-turns, reorganization of salt bridges, and dihedral changes in the hydrophobic regions are the major driving forces for helix-fibril transition. Structural features of the intermediates also correlated with the earlier experimental and computational studies. The study provides critical information on the fibrillation pathway of α-syn.

Molecular Insights into the Misfolding and Dimerization Dynamics of the Full-Length α-Synuclein from Atomistic Discrete Molecular Dynamics Simulations

The misfolding and pathological aggregation of α-synuclein forming insoluble amyloid deposits is associated with Parkinson's disease, the second most common neurodegenerative disease in the world population. Characterizing the self-assembly mechanism of α-synuclein is critical for discovering treatments against synucleinopathies. The intrinsically disordered property, high degrees of freedom, and macroscopic timescales of conformational conversion make its characterization extremely challenging in vitro and in silico. Here, we systematically investigated the dynamics of monomer misfolding and dimerization of the full-length α-synuclein using atomistic discrete molecular dynamics simulations. Our results suggested that both α-synuclein monomers and dimers mainly adopted unstructured formations with partial helices around the N-terminus (residues 8-32) and various β-sheets spanning the residues 35-56 (N-terminal tail) and residues 61-95 (NAC region). The C-terminus mostly assumed an unstructured formation wrapping around the lateral surface and the elongation edge of the β-sheet core formed by an N-terminal tail and NAC regions. Dimerization enhanced the β-sheet formation along with a decrease in the unstructured content. The inter-peptide β-sheets were mainly formed by the N-terminal tail and NACore (residues 68-78) regions, suggesting that these two regions played critical roles in the amyloid aggregation of α-synuclein. Interactions of the C-terminus with the N-terminal tail and the NAC region were significantly suppressed in the α-synuclein dimer, indicating that the interaction of the C-terminus with the N-terminal tail and NAC regions could prevent α-synuclein aggregation. These results on the structural ensembles and early aggregation dynamics of α-synuclein will help understand the nucleation and fibrillization of α-synuclein.

Aggregation of a Parkinson's Disease-Related Peptide: When Does Urea Weaken Hydrophobic Interactions?

While the exact cause of neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease is not completely understood, compelling evidence implicates the aggregation of specific proteins and peptides. Co-solvents can provide molecular insight into protein aggregation mechanisms and the chemical nature of potential aggregation inhibitors. Here, we study, through molecular simulations, the hydration and binding free energies of an amphiphilic peptide from the nonamyloid-β component (NAC), a key aggregation-prone domain of α-synuclein, in water and an 8 M aqueous urea solution. Isoleucine, glycine, and serine peptides of the same length are also studied to unravel the role of urea in the hydration and aggregation of hydrophobic and hydrophilic domains. A strong impact of urea in hindering the aggregation of the NAC subdomain is observed. A slightly weaker aggregation inhibition is observed for the Gly and Ser peptides, whereas a much lower aggregation inhibitory activity is found for the Ile peptide, seemingly contrasting with urea's protein unfolding mechanism. This behavior is shown to derive from a lower profusion of urea next to the hydrophobic side chains and the backbone of the Ile's peptide in the dimeric form. As a consequence, β-sheets, formed upon aggregation, remain nearly intact. Hydrophilic neighbor groups in the amphiphilic NAC subdomain, however, are shown to anchor enough urea to weaken hydrophobic interactions and disrupt β-sheet structures. Our results indicate that urea's activity is potentiated in amphiphilic domains and that potential drugs could disrupt hydrophobic β-sheet-rich regions while not binding primarily to hydrophobic amino acids.

Effect of an Amyloidogenic SARS-COV-2 Protein Fragment on α-Synuclein Monomers and Fibrils

Aggregates of α-synuclein are thought to be the disease-causing agent in Parkinson's disease. Various case studies have hinted at a correlation between COVID-19 and the onset of Parkinson's disease. For this reason, we use molecular dynamics simulations to study whether amyloidogenic regions in SARS-COV-2 proteins can initiate and modulate aggregation of α-synuclein. As an example, we choose the nine-residue fragment SFYVYSRVK (SK9), located on the C-terminal of the envelope protein of SARS-COV-2. We probe how the presence of SK9 affects the conformational ensemble of α-synuclein monomers and the stability of two resolved fibril polymorphs. We find that the viral protein fragment SK9 may alter α-synuclein amyloid formation by shifting the ensemble toward aggregation-prone and preferentially rod-like fibril seeding conformations. However, SK9 has only a small effect on the stability of pre-existing or newly formed fibrils. A potential mechanism and key residues for potential virus-induced amyloid formation are described.

Effect of an amyloidogenic SARS-COV-2 protein fragment on α-synuclein monomers and fibrils

Using molecular dynamic simulations we study whether amyloidogenic regions in viral proteins can initiate and modulate formation of α-synuclein aggregates, thought to be the disease-causing agent in Parkinson's Disease. As an example we choose the nine-residue fragment SFYVYSRVK (SK9), located on the C-terminal of the Envelope protein of SARS-COV-2. We probe how the presence of SK9 affects the conformational ensemble of α-synuclein monomers and the stability of two resolved fibril polymorphs. We find that the viral protein fragment SK9 may alter α-synuclein amyloid formation by shifting the ensemble toward aggregation-prone and preferentially rod-like fibril seeding conformations. However, SK9 has only little effect of the stability of pre-existing or newly-formed fibrils.

Importance of Solvent in Guiding the Conformational Properties of an Intrinsically Disordered Peptide

Aggregated form of α-synuclein in the brain has been found to be the major component of Lewy bodies that are hallmarks of Parkinson's disease (PD), the second most devastating neurodegenerative disorder. We have carried out room-temperature all-atom molecular dynamics (MD) simulations of an ensemble of widely different α-synuclein_1-95 peptide monomer conformations in aqueous solution. Attempts have been made to obtain a generic understanding of the local conformational motions of different repeat unit segments, namely R1-R7, of the peptide and the correlated properties of the solvent at the interface. The analyses revealed relatively greater rigidity of the hydrophobic R6 unit as compared to the other repeat units of the peptide. Besides, water molecules around R6 have been found to be less structured and weakly interacting with the peptide. These are important observations as the R6 unit with reduced conformational motions can act as the nucleation site for the aggregation process, while less structured weakly interacting water around it can become displaced easily, thereby facilitating the hydrophobic collapse of the peptide monomers and their association during the nucleation phase at higher concentrations. In addition, we demonstrated presence of doubly coordinated highly ordered as well as triply coordinated relatively disordered water molecules at the interface. We believe that while the ordered water molecules can favor water-mediated interactions between different peptide monomers, the randomly ordered ones on the other hand are likely to be expelled easily from the interface, thereby facilitating direct peptide-peptide interactions during the aggregation process.

Insights into the molecular mechanism of amyloid filament formation: Segmental folding of α-synuclein on lipid membranes

Recent advances in the structural biology of disease-relevant α-synuclein fibrils have revealed a variety of structures, yet little is known about the process of fibril aggregate formation. Characterization of intermediate species that form during aggregation is crucial; however, this has proven very challenging because of their transient nature, heterogeneity, and low population. Here, we investigate the aggregation of α-synuclein bound to negatively charged phospholipid small unilamellar vesicles. Through a combination of kinetic and structural studies, we identify key time points in the aggregation process that enable targeted isolation of prefibrillar and early fibrillar intermediates. By using solid-state nuclear magnetic resonance, we show the gradual buildup of structural features in an α-synuclein fibril filament, revealing a segmental folding process. We identify distinct membrane-binding domains in α-synuclein aggregates, and the combined data are used to present a comprehensive mechanism of the folding of α-synuclein on lipid membranes.

Atomic Resolution Map of Hierarchical Self-Assembly for an Amyloidogenic Protein Probed through Thermal 15N-R2 Correlation Matrices

Soluble oligomers formed by amyloidogenic intrinsically disordered proteins are some of the most cytotoxic species linked to neurodegeneration. Due to the transient and heterogeneous nature of such oligomeric intermediates, the underlying self-association events often remain elusive. NMR relaxation measurements sensitive to zero-frequency spectral densities (J(0)), such as the ¹⁵N - R₂ rates, are ideally suited to map sites of self-association at atomic resolution without the need of exogenous labels. Such experiments exploit the dynamic exchange between NMR visible monomers and slowly tumbling oligomers. However,¹⁵N - R₂ rates are also sensitive to intrinsic monomer dynamics, and it is often difficult to discern these contributions from those arising from exchange with oligomers. Another challenge pertains to defining a hierarchy of self-association. Here, using the archetypical amyloidogenic protein alpha synuclein (αS), we show that the temperature-dependence of ¹⁵N - R₂ effectively identifies self-association sites with reduced bias from internal dynamics. The key signature of the residues involved in self-association is a nonlinear temperature-dependence of ¹⁵N - R₂ with a positive ΔR₂/ΔT slope. These two hallmarks are systematically probed through a thermal R₂ correlation matrix, from which the network of residues involved in self-association as well as the hierarchy of αS self-association sites is extracted through agglomerative clustering. We find that aggregation is initiated by residues within the NAC region that is solvent inaccessible in αS fibrils and eventually extends to the N-terminal segment harboring familial PD mutations. These hierarchical self-association maps help dissect the essential drivers of oligomerization and reveal how amyloid inhibitors affect oligomer formation.

β-Turn exchanges in the α-synuclein segment 44-TKEG-47 reveal high sequence fidelity requirements of amyloid fibril elongation

The folding of turns and β-hairpins has been implicated in amyloid formation, with diverse potential consequences such as promotion or inhibition of fibril nucleation, fibril elongation, or off-pathway oligomer formation. In the Parkinson's disease-associated protein α-synuclein (αS), a β-hairpin comprised of residues 36-56 was detected in complex with an engineered binding protein, with a turn formed by the αS sequence segment 44-TKEG-47. Molecular dynamics simulations revealed extensive populations of transient β-hairpin conformations in this region in free, monomeric αS. Here, we investigated potential effects of turn formation on αS fibril formation by studying the aggregation kinetics of an extensive set of αS variants with between two and four amino acid exchanges in the 44-TKEG-47 segment. The exchanges were chosen to specifically promote formation of β1-, β1'-, or β2'-turns. All variants assembled into amyloid fibrils, with increased β1'- or β2'-turn propensity associated with faster aggregation and increased β1-turn propensity with slower aggregation compared to wild-type (WT) αS. Atomic force microscopy demonstrated that β-turn exchanges altered fibril morphology. In cross-elongation experiments, the turn variants showed a low ability to elongate WT fibril seeds, and, vice versa, WT monomer did not efficiently elongate turn variant fibril seeds. This demonstrates that sequence identity in the turn region is crucial for efficient αS fibril elongation. Elongation experiments of WT fibril seeds in the presence of both WT and turn variant monomers suggest that the turn variants can bind and block WT fibril ends to different degrees, but cannot efficiently convert into the WT fibril structure. Our results indicate that modifications in the 44-TKEG-47 segment strongly affect amyloid assembly by driving αS into alternative fibril morphologies, whose elongation requires high sequence fidelity.

Mounting evidence of FKBP12 implication in neurodegeneration

Intrinsically disordered proteins, such as tau or α-synuclein, have long been associated with a dysfunctional role in neurodegenerative diseases. In Alzheimer's and Parkinson's' diseases, these proteins, sharing a common chemical-physical pattern with alternating hydrophobic and hydrophilic domains rich in prolines, abnormally aggregate in tangles in the brain leading to progressive loss of neurons. In this review, we present an overview linking the studies on the implication of the peptidyl-prolyl isomerase domain of immunophilins, and notably FKBP12, to a variety of neurodegenerative diseases, focusing on the molecular origin of such a role. The involvement of FKBP12 dysregulation in the aberrant aggregation of disordered proteins pinpoints this protein as a possible therapeutic target and, at the same time, as a predictive biomarker for early diagnosis in neurodegeneration, calling for the development of reliable, fast and cost-effective detection methods in body fluids for community-based screening campaigns.

Extent of N-terminus exposure of monomeric alpha-synuclein determines its aggregation propensity

As an intrinsically disordered protein, monomeric alpha-synuclein (aSyn) occupies a large conformational space. Certain conformations lead to aggregation prone and non-aggregation prone intermediates, but identifying these within the dynamic ensemble of monomeric conformations is difficult. Herein, we used the biologically relevant calcium ion to investigate the conformation of monomeric aSyn in relation to its aggregation propensity. We observe that the more exposed the N-terminus and the beginning of the NAC region of aSyn are, the more aggregation prone monomeric aSyn conformations become. Solvent exposure of the N-terminus of aSyn occurs upon release of C-terminus interactions when calcium binds, but the level of exposure and aSyn's aggregation propensity is sequence and post translational modification dependent. Identifying aggregation prone conformations of monomeric aSyn and the environmental conditions they form under will allow us to design new therapeutics targeted to the monomeric protein.

A short motif in the N-terminal region of α-synuclein is critical for both aggregation and function

Aggregation of human α-synuclein (αSyn) is linked to Parkinson’s disease (PD) pathology. The central region of the αSyn sequence contains the non-amyloid β-component (NAC) crucial for aggregation. However, how NAC flanking regions modulate αSyn aggregation remains unclear. Using bioinformatics, mutation, and NMR we identify a 7-residue sequence, named P1 (residues 36-42), that controls αSyn aggregation. Deletion or substitution of this ‘master-controller’ prevents aggregation at pH 7.5 in vitro. At lower pH, P1 synergises with a sequence containing the PreNAC region (P2, residues 45-57) to prevent aggregation. Deleting P1 (ΔP1) or both P1 and P2 (ΔΔ) also prevents age-dependent αSyn aggregation and toxicity in C. elegans models and prevents αSyn-mediated vesicle fusion by altering the conformational properties of the protein when lipid-bound. The results highlight the importance of a master-controller sequence motif that controls both αSyn aggregation and function- a region that could be targeted to prevent aggregation in disease.

Blocking the FKBP12 induced dendrimeric burst in aberrant aggregation of α-synuclein by using the ElteN378 synthetic inhibitor

α-Synuclein (α-syn), a disordered cytoplasmatic protein, plays a fundamental role in the pathogenesis of Parkinson’s disease (PD). Here, we have shown, using photophysical measurements, that addition of FKBP12 to α-syn solutions, dramatically accelerates protein aggregation, leading to an explosion of dendritic structures revealed by fluorescence and phase-contrast microscopy. We have further demonstrated that this aberrant α-syn aggregation can be blocked using a recently discovered non-immunosuppressive synthetic inhibitor of FKBP12, ElteN378. The role of FKBP12 and of ElteN378 in the α-syn aggregation mechanism has been elucidated using molecular dynamics simulations based on an effective coarse-grained model. The reported data not only reveal a new potent synthetic drug as a candidate for early stage treatment of α-syn dependent neurodegenerations but also pave the way to a deeper understanding of the mechanism of action of FKBP12 on α-syn oligomeric aggregation, a topic which is still controversial.

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

Modulation of α-Synuclein Aggregation by Cytochrome c Binding and Hetero-dityrosine Adduct Formation

The aggregation of α-synuclein (α-Syn) has been implicated strongly in Parkinson's disease (PD). The intrinsically disordered nature of α-Syn makes this protein prone to self-association or heteroassociation with another protein or lipid. While conformational fluctuation and free radical chemistry have been shown to play important roles in its ability toward self- and heteroassociation, any systematic understanding of their contributions is missing. Here, we report an in vitro investigation of the interaction between α-Syn and cytochrome c in the oxidized (cyt c III) and reduced forms (cyt c II), in which cyt c III was found to induce a large compaction of α-Syn and inhibit the aggregation by favoring a hetero-dityrosine bond formation. In contrast, the presence of cyt c II did not result in any compaction and its presence was found to facilitate α-Syn aggregation. The variation in the charge distribution of the surface residues of cyt c III and cyt c II is expected to play a decisive role in their interaction with α-Syn.

The Cellular Environment Affects Monomeric α-Synuclein Structure

The presynaptic protein α-synuclein (aSyn) is an 'intrinsically disordered protein' that is highly dynamic in conformation. Transient intramolecular interactions between its charged N and C termini, and between its hydrophobic region and the C terminus, prevent self-association. These interactions inhibit the formation of insoluble inclusions, which are the pathological hallmark of Parkinson's disease and many other synucleinopathies. This review discusses how these intramolecular interactions are influenced by the specific environment aSyn is in. We discuss how charge, pH, calcium, and salt affect the physiological structure of monomeric aSyn, and how they may favour the formation of toxic structures. The more we understand the dynamic conformations of aSyn, the better we can design desperately needed therapeutics to prevent disease progression.

Oleuropein aglycone stabilizes the monomeric α-synuclein and favours the growth of non-toxic aggregates

α-synuclein plays a key role in the pathogenesis of Parkinson’s disease (PD); its deposits are found as amyloid fibrils in Lewy bodies and Lewy neurites, the histopathological hallmarks of PD. Amyloid fibrillation is a progressive polymerization path starting from peptide/protein misfolding and proceeding through the transient growth of oligomeric intermediates widely considered as the most toxic species. Consequently, a promising approach of intervention against PD might be preventing α-synuclein build-up, misfolding and aggregation. A possible strategy involves the use of small molecules able to slow down the aggregation process or to alter oligomer conformation favouring the growth of non-pathogenic species. Here, we show that oleuropein aglycone (OleA), the main olive oil polyphenol, exhibits anti-amyloidogenic power in vitro by interacting with, and stabilizing, α-synuclein monomers thus hampering the growth of on-pathway oligomers and favouring the growth of stable and harmless aggregates with no tendency to evolve into other cytotoxic amyloids. We investigated the molecular basis of such interference by both biophysical techniques and limited proteolysis; aggregate morphology was monitored by electron microscopy. We also found that OleA reduces the cytotoxicity of α-synuclein aggregates by hindering their binding to cell membrane components and preventing the resulting oxidative damage to cells.

Perturbation in Long-Range Contacts Modulates the Kinetics of Amyloid Formation in α-Synuclein Familial Mutants

The characteristic cross-β-sheet-rich amyloid fibril formation by intrinsically disordered α-synuclein proteins is one of the pathological hallmarks of Parkinson's disease. Although unstructured in solution, the presence of autoinhibitory long-range contacts in monomeric form prevents protein aggregation. Out of the various factors that affect the rate of amyloid formation, familial mutations play an important role in α-synuclein aggregation. Even though these mutations are believed to form an aggregation-prone intermediate by perturbing these contacts, the correlation between perturbation and rate of fibril formation is not very straightforward. A combination of solution and solid-state NMR in conjunction with other biophysical methods has been used to identify the underlying mechanism behind the anomaly in the rate of aggregation for the novel mutants H50Q (fast aggregating) and G51D (slow aggregating). Perturbation of long-range contacts at the mutation sites and C-termini in all of the six familial mutants of α-synuclein during the diseased condition (acidic pH) was observed. These contacts get rearranged at physiological pH resulting in the shielding of mutation sites. Additional contacts at the mutation site in a slow aggregating mutant could be the reason for slower aggregation. Indeed, these contacts provide more rigidity to the monomeric G51D. Nonetheless, these mutations did not alter the overall secondary structure. The differential pattern of the long-range contacts at the monomeric level resulted in the perturbation of the fibrillar-core region, which was evident in the solid-state NMR spectra. Our results provide valuable insights in understanding the effect of long-range contacts on the aggregation of α-synuclein and its mutants.

Discrete Molecular Dynamics Approach to the Study of Disordered and Aggregating Proteins

We present a refinement of the Coarse Grained PACSAB force field for Discrete Molecular Dynamics (DMD) simulations of proteins in aqueous conditions. As the original version, the refined method provides good representation of the structure and dynamics of folded proteins but provides much better representations of a variety of unfolded proteins, including some very large, impossible to analyze by atomistic simulation methods. The PACSAB/DMD method also reproduces accurately aggregation properties, providing good pictures of the structural ensembles of proteins showing a folded core and an intrinsically disordered region. The combination of accuracy and speed makes the method presented here a good alternative for the exploration of unstructured protein systems.

Phosphorylation induces distinct alpha-synuclein strain formation

Synucleinopathies are a group of neurodegenerative diseases associated with alpha-synuclein (α-Syn) aggregation. Recently, increasing evidence has demonstrated the existence of different structural characteristics or 'strains' of α-Syn, supporting the concept that synucleinopathies share several common features with prion diseases and possibly explaining how a single protein results in different clinical phenotypes within synucleinopathies. In earlier studies, the different strains were generated through the regulation of solution conditions, temperature, or repetitive seeded fibrillization in vitro. Here, we synthesize homogeneous α-Syn phosphorylated at serine 129 (pS129 α-Syn), which is highly associated with the pathological changes, and demonstrate that phosphorylation at Ser129 induces α-Syn to form a distinct strain with different structures, propagation properties, and higher cytotoxicity compared with the wild-type α-Syn. The results are the first demonstration that post-translational modification of α-Syn can induce different strain formation, offering a new mechanism for strain formation.

The neural chaperone proSAAS blocks α-synuclein fibrillation and neurotoxicity

Emerging evidence strongly suggests that chaperone proteins are cytoprotective in neurodegenerative proteinopathies involving protein aggregation; for example, in the accumulation of aggregated α-synuclein into the Lewy bodies present in Parkinson's disease. Of the various chaperones known to be associated with neurodegenerative disease, the small secretory chaperone known as proSAAS (named after four residues in the amino terminal region) has many attractive properties. We show here that proSAAS, widely expressed in neurons throughout the brain, is associated with aggregated synuclein deposits in the substantia nigra of patients with Parkinson's disease. Recombinant proSAAS potently inhibits the fibrillation of α-synuclein in an in vitro assay; residues 158-180, containing a largely conserved element, are critical to this bioactivity. ProSAAS also exhibits a neuroprotective function; proSAAS-encoding lentivirus blocks α-synuclein-induced cytotoxicity in primary cultures of nigral dopaminergic neurons, and recombinant proSAAS blocks α-synuclein-induced cytotoxicity in SH-SY5Y cells. Four independent proteomics studies have previously identified proSAAS as a potential cerebrospinal fluid biomarker in various neurodegenerative diseases. Coupled with prior work showing that proSAAS blocks β-amyloid aggregation into fibrils, this study supports the idea that neuronal proSAAS plays an important role in proteostatic processes. ProSAAS thus represents a possible therapeutic target in neurodegenerative disease.

X-ray Crystallographic Structure of Oligomers Formed by a Toxic β-Hairpin Derived from α-Synuclein: Trimers and Higher-Order Oligomers

Oligomeric assemblies of the protein α-synuclein are thought to cause neurodegeneration in Parkinson’s disease and related synucleinopathies. Characterization of α-synuclein oligomers at high resolution is an outstanding challenge in the field of structural biology. The absence of high-resolution structures of oligomers formed by α-synuclein impedes understanding the synucleinopathies at the molecular level. This paper reports the X-ray crystallographic structure of oligomers formed by a peptide derived from residues 36–55 of α-synuclein. The peptide 1a adopts a β-hairpin structure, which assembles in a hierarchical fashion. Three β-hairpins assemble to form a triangular trimer. Three copies of the triangular trimer assemble to form a basket-shaped nonamer. Two nonamers pack to form an octadecamer. Molecular modeling suggests that full-length α-synuclein may also be able to assemble in this fashion. Circular dichroism spectroscopy demonstrates that peptide 1a interacts with anionic lipid bilayer membranes, like oligomers of full-length α-synuclein. LDH and MTT assays demonstrate that peptide 1a is toxic toward SH-SY5Y cells. Comparison of peptide 1a to homologues suggests that this toxicity results from nonspecific interactions with the cell membrane. The oligomers formed by peptide 1a are fundamentally different than the proposed models of the fibrils formed by α-synuclein and suggest that α-Syn_36–55, rather than the NAC, may nucleate oligomer formation.

Transient β-hairpin formation in α-synuclein monomer revealed by coarse-grained molecular dynamics simulation

Parkinson's disease, originating from the intrinsically disordered peptide α-synuclein, is a common neurodegenerative disorder that affects more than 5% of the population above age 85. It remains unclear how α-synuclein monomers undergo conformational changes leading to aggregation and formation of fibrils characteristic for the disease. In the present study, we perform molecular dynamics simulations (over 180 μs in aggregated time) using a hybrid-resolution model, Proteins with Atomic details in Coarse-grained Environment (PACE), to characterize in atomic detail structural ensembles of wild type and mutant monomeric α-synuclein in aqueous solution. The simulations reproduce structural properties of α-synuclein characterized in experiments, such as secondary structure content, long-range contacts, chemical shifts, and (3)J(HNHCα )-coupling constants. Most notably, the simulations reveal that a short fragment encompassing region 38-53, adjacent to the non-amyloid-β component region, exhibits a high probability of forming a β-hairpin; this fragment, when isolated from the remainder of α-synuclein, fluctuates frequently into its β-hairpin conformation. Two disease-prone mutations, namely, A30P and A53T, significantly accelerate the formation of a β-hairpin in the stated fragment. We conclude that the formation of a β-hairpin in region 38-53 is a key event during α-synuclein aggregation. We predict further that the G47V mutation impedes the formation of a turn in the β-hairpin and slows down β-hairpin formation, thereby retarding α-synuclein aggregation.

Information flow and protein dynamics: the interplay between nuclear magnetic resonance spectroscopy and molecular dynamics simulations

Proteins participate in information pathways in cells, both as links in the chain of signals, and as the ultimate effectors. Upon ligand binding, proteins undergo conformation and motion changes, which can be sensed by the following link in the chain of information. Nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations represent powerful tools for examining the time-dependent function of biological molecules. The recent advances in NMR and the availability of faster computers have opened the door to more detailed analyses of structure, dynamics, and interactions. Here we briefly describe the recent applications that allow NMR spectroscopy and MD simulations to offer unique insight into the basic motions that underlie information transfer within and between cells.

Water dispersion interactions strongly influence simulated structural properties of disordered protein states

Many proteins can be partially or completely disordered under physiological conditions. Structural characterization of these disordered states using experimental methods can be challenging, since they are composed of a structurally heterogeneous ensemble of conformations rather than a single dominant conformation. Molecular dynamics (MD) simulations should in principle provide an ideal tool for elucidating the composition and behavior of disordered states at an atomic level of detail. Unfortunately, MD simulations using current physics-based models tend to produce disordered-state ensembles that are structurally too compact relative to experiments. We find that the water models typically used in MD simulations significantly underestimate London dispersion interactions, and speculate that this may be a possible reason for these erroneous results. To test this hypothesis, we create a new water model, TIP4P-D, that approximately corrects for these deficiencies in modeling water dispersion interactions while maintaining compatibility with existing physics-based models. We show that simulations of solvated proteins using this new water model typically result in disordered states that are substantially more expanded and in better agreement with experiment. These results represent a significant step toward extending the range of applicability of MD simulations to include the study of (partially or fully) disordered protein states.

Molecular details of α-synuclein membrane association revealed by neutrons and photons

α-Synuclein (α-syn) is an abundant neuronal protein associated with Parkinson's disease that is disordered in solution, but it exists in equilibrium between a bent-helix and an elongated-helix on negatively charged membranes. Here, neutron reflectometry (NR) and fluorescence spectroscopy were employed to uncover molecular details of the interaction between α-syn and two anionic lipids, phosphatidic acid (PA) and phosphatidylserine (PS). Both NR and site-specific Trp measurements indicate that penetration depth of α-syn is similar for either PA- or PS-containing membranes (∼9-11 Å from bilayer center) even though there is a preference for α-syn binding to PA. However, closer examination of the individual Trp quenching profiles by brominated lipids reveals differences into local membrane interactions especially at position 39 where conformational heterogeneity was observed. The data also indicate that while W94 penetrates the bilayer as deeply as W4, W94 resides in a more polar surrounding. Taken together, we suggest the N- and C-terminal regions near positions 4 and 94 are anchored to the membrane, while the putative linker spanning residue 39 samples multiple conformations, which are sensitive to the chemical nature of the membrane surface. This flexibility may enable α-syn to bind diverse biomembranes in vivo.

Quantitative biophysical characterization of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are broadly defined as protein regions that do not cooperatively fold into a spatially or temporally stable structure. Recent research strongly supports the hypothesis that a conserved functional role for structural disorder renders IDPs uniquely capable of functioning in biological processes such as cellular signaling and transcription. Recently, the frequency of application of rigorous mechanistic biochemistry and quantitative biophysics to disordered systems has increased dramatically. For example, the launch of the Protein Ensemble Database (pE-DB) demonstrates that the potential now exists to refine models for the native state structure of IDPs using experimental data. However, rigorous assessment of which observables place the strongest and least biased constraints on those ensembles is now needed. Most importantly, the past few years have seen strong growth in the number of biochemical and biophysical studies attempting to connect structural disorder with function. From the perspective of equilibrium thermodynamics, there is a clear need to assess the relative significance of hydrophobic versus electrostatic forces in IDP interactions, if it is possible to generalize at all. Finally, kinetic mechanisms that invoke conformational selection and/or induced fit are often used to characterize coupled IDP folding and binding, although application of these models is typically built upon thermodynamic observations. Recently, the reaction rates and kinetic mechanisms of more intrinsically disordered systems have been tested through rigorous kinetic experiments. Motivated by these exciting advances, here we provide a review and prospectus for the quantitative study of IDP structure, thermodynamics, and kinetics.

Structural insights into amyloid oligomers of the Parkinson disease-related protein α-synuclein

The presence of intraneuronal deposits mainly formed by amyloid fibrils of the presynaptic protein α-synuclein (AS) is a hallmark of Parkinson disease. Currently, neurotoxicity is attributed to prefibrillar oligomeric species rather than the insoluble aggregates, although their mechanisms of toxicity remain elusive. Structural details of the supramolecular organization of AS oligomers are critically needed to decipher the structure-toxicity relationship underlying their pathogenicity. In this study, we employed site-specific fluorescence to get a deeper insight into the internal architecture of AS oligomeric intermediates. We demonstrate that AS oligomers are ordered assemblies possessing a well defined pattern of intermolecular contacts. Some of these contacts involve regions that form the β-sheet core in the fibrillar state, although their spatial arrangement may differ in the two aggregated forms. However, even though the two termini are excluded from the fibrillar core, they are engaged in a number of intermolecular interactions within the oligomer. Therefore, substantial structural remodeling of early oligomeric interactions is essential for fibril growth. The intermolecular contacts identified in AS oligomers can serve as targets for the rational design of anti-amyloid compounds directed at preventing oligomeric interactions/reorganizations.

A maximum entropy approach to the study of residue-specific backbone angle distributions in α-synuclein, an intrinsically disordered protein

α-Synuclein is an intrinsically disordered protein of 140 residues that switches to an α-helical conformation upon binding phospholipid membranes. We characterize its residue-specific backbone structure in free solution with a novel maximum entropy procedure that integrates an extensive set of NMR data. These data include intraresidue and sequential H(N) − H(α) and H(N) − H(N) NOEs, values for (3) JHNHα, (1) JHαCα, (2) JCαN, and (1) JCαN, as well as chemical shifts of (15)N, (13)C(α), and (13)C' nuclei, which are sensitive to backbone torsion angles. Distributions of these torsion angles were identified that yield best agreement to the experimental data, while using an entropy term to minimize the deviation from statistical distributions seen in a large protein coil library. Results indicate that although at the individual residue level considerable deviations from the coil library distribution are seen, on average the fitted distributions agree fairly well with this library, yielding a moderate population (20-30%) of the PPII region and a somewhat higher population of the potentially aggregation-prone β region (20-40%) than seen in the database. A generally lower population of the αR region (10-20%) is found. Analysis of (1)H − (1)H NOE data required consideration of the considerable backbone diffusion anisotropy of a disordered protein.

Digested disorder: Quarterly intrinsic disorder digest (July-August-September, 2013)

The current literature on intrinsically disordered proteins grows fast. To keep interested readers up to speed with this literature, we continue a "Digested Disorder" project and represent a new issue of reader's digest of the research papers and reviews on intrinsically disordered proteins. The only 2 criteria for inclusion in this digest are the publication date (a paper should be published within the covered time frame) and topic (a paper should be dedicated to any aspect of protein intrinsic disorder). The current digest issue covers papers published during the third quarter of 2013; i.e., during the period of June, July, and September of 2013. Similar to previous issues, the papers are grouped hierarchically by topics they cover, and for each of the included paper a short description is given on its major findings.