Understanding the dynamics of monomeric, dimeric, and tetrameric α-synuclein structures in water

Differentially heterogeneous hydration environment of the familial mutants of α-synuclein

The behavior of hydration water around familial Parkinson's disease linked mutants of α-synuclein may be linked to the early-onset of Parkinson's disease. For the first time, this study compares the local structure and dynamics of hydration water around different segments of some of the natural mutants of α-synuclein, i.e., E46K, G51D, A30P, and A53E, with that of the wild-type protein through explicit water MD simulations. The results show that the C-terminal segments of the fast aggregating mutants such as E46K and A30P are less exposed to water, while those of the slow aggregating ones such as A53E and G51D are more exposed to water relative to that of the wild-type protein. In addition, the water molecules are found to be more ordered around the C-terminal segment of the A53E and G51D mutants as compared to the wild-type protein. This is due to an increase in the overall charge of α-syn upon A53E and G51D mutations. The translational and rotational motions of water molecules in the hydration shell of the C-terminal segment of A53E and G51D mutants are found to be faster relative to that of the wild-type protein. This study validates the differential hydration environment around the C-terminal segment for the causative and protective mutants of α-synuclein.

Wild-Type α-Synuclein Structure and Aggregation: A Comprehensive Coarse-Grained and All-Atom Molecular Dynamics Study

α-Synuclein (α-syn) is a 140 amino acid intrinsically disordered protein (IDP) and the primary component of cytotoxic oligomers implicated in the etiology of Parkinson's disease (PD). While IDPs lack a stable three-dimensional structure, they sample a heterogeneous ensemble of conformations that can, in principle, be assessed through molecular dynamics simulations. However, describing the structure and aggregation of large IDPs is challenging due to force field (FF) accuracy and sampling limitations. To cope with the latter, coarse-grained (CG) FFs emerge as a potential alternative at the expense of atomic detail loss. Whereas CG models can accurately describe the structure of the monomer, less is known about aggregation. The latter is key for assessing aggregation pathways and designing aggregation inhibitor drugs. Herein, we investigate the structure and dynamics of α-syn using different resolution CG (Martini3 and Sirah2) and all-atom (Amber99sb and Charmm36m) FFs to gain insight into the differences and resemblances between these models. The dependence of the magnitude of protein-water interactions and the putative need for enhanced sampling (replica exchange) methods in CG simulations are analyzed to distinguish between force field accuracy and sampling limitations. The stability of the CG models of an α-syn fibril was also investigated. Additionally, α-syn aggregation was studied through umbrella sampling for the CG models and CG/all-atom models for an 11-mer peptide (NACore) from an amyloidogenic domain of α-syn. Our results show that despite the α-syn structures of Martini3 and Sirah2 with enhanced protein-water interactions being similar, major differences exist concerning aggregation. The Martini3 fibril is not stable, and the binding free energy of α-syn and NACore is positive, opposite to Sirah2. Sirah2 peptides in a zwitterionic form, in turn, display termini interactions that are too strong, resulting in end-to-end orientation. Sirah2, with enhanced protein-water interactions and neutral termini, provides, however, a peptide aggregation free energy profile similar to that found with all-atom models. Overall, we find that Sirah2 with enhanced protein-water interactions is suitable for studying protein-protein and protein-drug aggregation.

Structural and dynamic insights into α-synuclein dimer conformations

Parkinson disease is associated with the aggregation of the protein α-synuclein. While α-synuclein can exist in multiple oligomeric states, the dimer has been a subject of extensive debates. Here, using an array of biophysical approaches, we demonstrate that α-synuclein in vitro exhibits primarily a monomer-dimer equilibrium in nanomolar concentrations and up to a few micromolars. We then use spatial information from hetero-isotopic cross-linking mass spectrometry experiments as restrains in discrete molecular dynamics simulations to obtain the ensemble structure of dimeric species. Out of eight structural sub-populations of dimers, we identify one that is compact, stable, abundant, and exhibits partially exposed β-sheet structures. This compact dimer is the only one where the hydroxyls of tyrosine 39 are in proximity that may promote dityrosine covalent linkage upon hydroxyl radicalization, which is implicated in α-synuclein amyloid fibrils. We propose that this α-synuclein dimer features etiological relevance to Parkinson disease.

Essential collective dynamics analysis reveals nonlocal interactions of alpha-synuclein38-95 monomers with fibrillar seeds

Alpha-synuclein (α-syn) is a cytoplasmic protein commonly found in the nervous system. In solution, α-syn adopts disordered unfolded conformations, although it can also form α-helices upon binding to lipid membranes. Under conditions that are not yet fully understood, α-syn can misfold and aggregate, giving rise to β-sheet rich amyloid fibrils, which then tend to accumulate in degenerating neurons. This leads to Parkinson's disease (PD) and several other conditions collectively termed synucleinopathies. Development of disease-modifying treatments requires detailed understanding of structure and dynamics of α-syn's misfolded aggregates. We have employed 1000 ns long, all-atom molecular dynamics simulations to investigate the interaction of monomeric α-syn_38-95 fragments, which contain the most important amyloidogenic regions, with preformed fibrillar seeds composed of staggered, β-sheet rich α-syn chains of matching length. The simulations indicate that α-syn_38-95 monomers tend to form aggregates with the fibrillar seeds, although we have not observed alignment of the monomeric chains with β-strands of the fibril. To analyze the stability of these aggregates, we have employed the essential collective dynamics method, which allows making accurate assessment of dynamical coupling across individual atoms in macromolecules and supramolecular complexes. The analysis revealed extensive dynamical coupling across initially monomeric α-syn chains and the fibrillar seeds including distal regions thereof that did not contact the monomer directly. We have discussed structural origins of these long-range interactions, their impacts for the stability of α-syn aggregates, and potential implications for the development of anti-PD treatments.

Insights into the Interactions that Trigger the Primary Nucleation of Polymorphic α-Synuclein Dimers

Parkinson's disease is associated with the accumulation of α-synuclein (AS) aggregates that include polymorphic AS oligomers and polymorphic fibrils. There have been advances in solving the polymorphic state of AS fibrils, both by experimental techniques and molecular modeling tools. Yet, the polymorphic AS oligomers are now considered as the neurotoxic species, thus current and future studies making efforts to solve their structures at the molecular level. Importantly, it is crucial to explore the specific interactions between AS monomers within the dimer that stabilize the dimer and yield nucleation. Herein, we present a first work that probes at the molecular level the specific interactions between monomers in polymorphic AS dimers are derived from AS fibrils by applying molecular modeling tools. Our work reveals that both N-terminal and the non-amyloidogenic component domains play a role in the dimerization of all polymorphic AS dimers. In addition, helices along the N-terminal of AS monomers impede the contacts between AS monomers, thus preventing the nucleation or the dimerization of AS. This work provides insights into several mechanisms of the production of polymorphic AS dimers. Thus, the findings obtained in this work may assist in developing new therapeutic strategies for inhibiting the formation of the early-stage neurotoxic AS dimers.

Engineering of a protein probe with multiple inputs and multiple outputs for evaluation of alpha synuclein aggregation states

The aggregation of α-synuclein (αS) into oligomers and fibrils is implicated in the pathology of Parkinson's Disease (PD). While a molecular probe for rapid and comprehensive evaluation of αS aggregation states is critical for a better understanding of PD pathology, identification of therapeutic candidates, and the development of early diagnostic strategies, no such probe has yet to be developed. A structurally flexible αS variant, PG65, was previously developed as a target binding-driven, conformation-switching molecular probe for rapid αS oligomer detection. Though informative, detection using PG65 provides no comprehensive assessment of the αS aggregation states. In the present study, we report engineering of a molecular probe, PG65-MIMO (a PG65 variant with Multiple-Inputs and Multiple-Outputs), that rapidly (within 2 hr) produces comprehensive information on αS aggregation states. PG65-MIMO generates distinct fluorescence responses to the three major αS conformers (monomers, oligomers, and fibrils). PG65-MIMO also displays unique fluorescent signals for αS oligomers, depending on the tris(2-carboxyethyl)phosphine (TCEP) concentration. Our results suggest that the TCEP dependent signaling of PG65-MIMO may be associated with its conformational states. Overall, our study illustrates engineering of an αS variant to create a molecular probe for handling multiple inputs and multiple outputs, addressing the technological gap in αS detection.

The structural heterogeneity of α-synuclein is governed by several distinct subpopulations with interconversion times slower than milliseconds

α-Synuclein plays an important role in synaptic functions by interacting with synaptic vesicle membrane, while its oligomers and fibrils are associated with several neurodegenerative diseases. The specific monomer structures that promote its membrane binding and self-association remain elusive due to its transient nature as an intrinsically disordered protein. Here, we use inter-dye distance distributions from bulk time-resolved Förster resonance energy transfer as restraints in discrete molecular dynamics simulations to map the conformational space of the α-synuclein monomer. We further confirm the generated conformational ensemble in orthogonal experiments utilizing far-UV circular dichroism and cross-linking mass spectrometry. Single-molecule protein-induced fluorescence enhancement measurements show that within this conformational ensemble, some of the conformations of α-synuclein are surprisingly stable, exhibiting conformational transitions slower than milliseconds. Our comprehensive analysis of the conformational ensemble reveals essential structural properties and potential conformations that promote its various functions in membrane interaction or oligomer and fibril formation.

Aggregation of Aβ40/42 chains in the presence of cyclic neuropeptides investigated by molecular dynamics simulations

Alzheimer’s disease is associated with the formation of toxic aggregates of amyloid beta (Aβ) peptides. Despite tremendous efforts, our understanding of the molecular mechanisms of aggregation, as well as cofactors that might influence it, remains incomplete. The small cyclic neuropeptide somatostatin-14 (SST₁₄) was recently found to be the most selectively enriched protein in human frontal lobe extracts that binds Aβ₄₂ aggregates. Furthermore, SST₁₄’s presence was also found to promote the formation of toxic Aβ₄₂ oligomers in vitro. In order to elucidate how SST₁₄ influences the onset of Aβ oligomerization, we performed all-atom molecular dynamics simulations of model mixtures of Aβ₄₂ or Aβ₄₀ peptides with SST₁₄ molecules and analyzed the structure and dynamics of early-stage aggregates. For comparison we also analyzed the aggregation of Aβ₄₂ in the presence of arginine vasopressin (AVP), a different cyclic neuropeptide. We observed the formation of self-assembled aggregates containing the Aβ chains and small cyclic peptides in all mixtures of Aβ₄₂–SST₁₄, Aβ₄₂–AVP, and Aβ₄₀–SST₁₄. The Aβ₄₂–SST₁₄ mixtures were found to develop compact, dynamically stable, but small aggregates with the highest exposure of hydrophobic residues to the solvent. Differences in the morphology and dynamics of aggregates that comprise SST₁₄ or AVP appear to reflect distinct (1) regions of the Aβ chains they interact with; (2) propensities to engage in hydrogen bonds with Aβ peptides; and (3) solvent exposures of hydrophilic and hydrophobic groups. The presence of SST₁₄ was found to impede aggregation in the Aβ₄₂–SST₁₄ system despite a high hydrophobicity, producing a stronger “sticky surface” effect in the aggregates at the onset of Aβ₄₂–SST₁₄ oligomerization.

Improper folding of proteins causes disorders known as protein misfolding diseases. Under normal conditions most proteins adopt particular folds, which allow them functioning properly. However, for reasons that are not yet fully understood, proteins may misfold and aggregate, forming deposits known as amyloid fibrils, which accumulate in the brain or other tissues. This process affects functioning of the nervous system, gradually causing loss of cognitive abilities. Alzheimer’s disease is one of the most common diseases from this group. A better understanding of the aggregation of peptides implicated in Alzheimer’s disease, known as amyloid beta (Aβ) peptides, may facilitate the development of treatments that ameliorate or prevent the disease. We use detailed molecular dynamics simulations to investigate the influence of somatostatin-14 (SST₁₄), a cyclic neuropeptide that might be involved in the amyloidogenic aggregation of Aβ, on molecular processes occurring at the onset of Aβ aggregation. Results of these simulations explain how the presence of SST₁₄ might alter pathways of aggregation of Aβ, shedding light upon the possible role of extrinsic factors in the aggregation at a molecular level.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Paradoxical Effect of Trehalose on the Aggregation of α-Synuclein: Expedites Onset of Aggregation yet Reduces Fibril Load

Aggregation of α-synuclein is closely connected to the pathology of Parkinson's disease. The phenomenon involves multiple steps, commenced by partial misfolding and eventually leading to mature amyloid fibril formation. Trehalose, a widely accepted osmolyte, has been shown previously to inhibit aggregation of various globular proteins owing to its ability to prevent the initial unfolding of protein. In this study, we have examined if it behaves in a similar fashion with intrinsically disordered protein α-synuclein and possesses the potential to act as therapeutic agent against Parkinson's disease. It was observed experimentally that samples coincubated with trehalose fibrillate faster compared to the case in its absence. Molecular dynamics simulations suggested that this initial acceleration is manifestation of trehalose's tendency to perturb the conformational transitions between different conformers of monomeric protein. It stabilizes the aggregation prone "extended" conformer of α-synuclein, by binding to its exposed acidic residues of the C terminus. It also favors the β-rich oligomers once formed. Interestingly, the total fibrils formed are still promisingly less since it accelerates the competing pathway toward formation of amorphous aggregates.

Conformational Ensembles of α-Synuclein Derived Peptide with Different Osmolytes from Temperature Replica Exchange Sampling

Intrinsically disordered proteins (IDP) are a class of proteins that do not have a stable three-dimensional structure and can adopt a range of conformations playing various vital functional role. Alpha-synuclein is one such IDP which can aggregate into toxic protofibrils and has been associated largely with Parkinson's disease (PD) along with other neurodegenerative diseases. Osmolytes are small organic compounds that can alter the environment around the proteins by acting as denaturants or protectants for the proteins. In the present study, we have conducted a series of replica exchange molecular dynamics simulations to explore the role of osmolytes, urea which is a denaturant and TMAO (trimethylamine N-oxide), a protecting osmolyte, in aggregation and conformations of the synuclein peptide. We observed that both the osmolytes have significantly distinct impacts on the peptide and led to transitions of the conformations of the peptide from one state to other. Our findings highlighted that urea attenuated peptide aggregation and resulted in the formation of extended peptide structures whereas TMAO led to compact and folded forms of the peptide.