Determination of the free energy landscape of alpha-synuclein using spin label nuclear magnetic resonance measurements

Alpha-synuclein phosphorylation induces amyloid conversion via enhanced electrostatic bridging: Insights from molecular modeling of the full-length protein

Fibril formation from alpha-synuclein is a key point in Parkinson's disease, multiple system atrophy, and other synucleinopathies. The mechanism of the amyloid-like conversion followed by the formation of pre-fibrillar soluble oligomers and fibrils is not completely clear; furthermore, it is unclear how the Parkinson's disease-related point mutations located in the pre-NAC region enhance fibrillation. In the present paper, atomistic replica exchange molecular dynamics simulations of the full-length alpha-synuclein and its two mutants, A53T and E46K, elucidated amyloid conversion intermediates. Both mutants demonstrated an enhanced tendency for the conversion but in different manners; the main intermediate conformations populated in the WT alpha-synuclein conformational ensemble disappeared due to mutations, indicating a different conversion pathway. Analysis of the preferable beta-hairpin positions and intermediate conformations seems to reflect a tendency to form a particular amyloid fibril polymorph. A strong elevation of amyloid transformation level was shown also for Ser129-phosphorylated alpha-synuclein. Altered intermediate conformations, the most preferable beta-hairpin positions in the NAC region, and prevalent salt bridges propose the formation of so-called polymorph 2 or even a novel type of fibrils. A better understanding of the detailed mechanism of the amyloid conversion sheds light on the effect of Lewy body-related phosphorylation and might help in the development of new therapeutics for synucleinopathies.

Early-Onset Parkinson Mutation Remodels Monomer-Fibril Interactions to Allosterically Amplify Synuclein's Amyloid Cascade

Alpha synuclein (αS) aggregates are the main component of Lewy bodies (LBs) associated with Parkinson's disease (PD). A longstanding question about αS and PD pertains to the autosomal dominant E46K αS mutant, which leads to the early onset of PD and LB dementias. The E46K mutation not only promotes αS aggregation but also stabilizes αS monomers in "closed" conformers, which are compact and aggregation-incompetent. Hence, the mechanism of action of the E46K mutation is currently unclear. Here, we show that αS monomers harboring the E46K mutation exhibit more extensive interactions with fibrils compared to those of WT. Such monomer-fibril interactions are sufficient to allosterically drive transitions of αS monomers from closed to open conformations, enabling αS aggregation. We also show that E46K promotes head-to-tail monomer-monomer interactions in early self-association events. This multipronged mechanism provides a new framework to explain how the E46K mutation and possibly other αS variants trigger early-onset PD.

The role of solvation on the conformational landscape of α-synuclein

Native ion mobility mass spectrometry has been used extensively to characterize ensembles of intrinsically disordered protein (IDP) conformers, but the extent to which the gaseous measurements provide realistic pictures of the solution conformations for such flexible proteins remains unclear. Therefore, we systematically studied the relationship between the solution and gaseous structural ensembles by measuring electrospray charge state and collision cross section (CCS) distributions for cationic and anionic forms of α-synuclein (αSN), an anionic protein in solution, as well as directly probed gas phase residue to residue distances via ion/ion reactions between gaseous α-synuclein cations and disulfonic acid linkers that form strong electrostatic bonds. We also combined results from in-solution protein crosslinking identified from native tandem mass spectrometry (MS/MS) with an initial αSN ensemble generated computationally by IDPConformerGenerator to generate an experimentally restrained solution ensemble of αSN. CCS distributions were directly calculated for the solution ensembles determined by NMR and compared to predicted gaseous conformers. While charge state and collision cross section distributions are useful for qualitatively describing the relative structural dynamics of proteins and major conformational changes induced by changes to solution states, the predicted and measured gas phase conformers include subpopulations that are significantly different than those expected from completely "freezing" solution conformations and preserving them in the gas phase. However, insights were gained on the various roles of solvent in stabilizing various conformers for extremely dynamic proteins like α-synuclein.

Conformational and Interaction Landscape of Histone H4 Tails in Nucleosomes Probed by Paramagnetic NMR Spectroscopy

The fundamental repeat unit of chromatin, the nucleosome, consists of approximately 147 base pairs of double-stranded DNA and a histone protein octamer containing two copies each of histones H2A, H2B, H3, and H4. Each histone possesses a dynamically disordered N-terminal tail domain, and it is well-established that the tails of histones H3 and H4 play key roles in chromatin compaction and regulation. Here we investigate the conformational ensemble and interactions of the H4 tail in nucleosomes by means of solution NMR measurements of paramagnetic relaxation enhancements (PREs) in recombinant samples reconstituted with 15N-enriched H4 and nitroxide spin-label tagged H3. The experimental PREs, which report on the proximities of individual H4 tail residues to the different H3 spin-label sites, are interpreted by using microsecond time-scale molecular dynamics simulations of the nucleosome core particle. Collectively, these data enable improved localization of histone H4 tails in nucleosomes and support the notion that H4 tails engage in a fuzzy complex interaction with nucleosomal DNA.

Structural characterization of an intrinsically disordered protein complex using integrated small-angle neutron scattering and computing

Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure-function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Lastly, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.

A New Tool to Study the Binding Behavior of Intrinsically Disordered Proteins

Understanding the binding behavior and conformational dynamics of intrinsically disordered proteins (IDPs) is crucial for unraveling their regulatory roles in biological processes. However, their lack of stable 3D structures poses challenges for analysis. To address this, we propose an algorithm that explores IDP binding behavior with protein complexes by extracting topological and geometric features from the protein surface model. Our algorithm identifies a geometrically favorable binding pose for the IDP and plans a feasible trajectory to evaluate its transition to the docking position. We focus on IDPs from Homo sapiens and Mus-musculus, investigating their interaction with the Plasmodium falciparum (PF) pathogen associated with malaria-related deaths. We compare our algorithm with HawkDock and HDOCK docking tools for quantitative (computation time) and qualitative (binding affinity) measures. Our results indicated that our method outperformed the compared methods in computation performance and binding affinity in experimental conformations.

The Action of Chemical Denaturants: From Globular to Intrinsically Disordered Proteins

Proteins perform their many functions by adopting either a minimal number of strictly similar conformations, the native state, or a vast ensemble of highly flexible conformations. In both cases, their structural features are highly influenced by the chemical environment. Even though a plethora of experimental studies have demonstrated the impact of chemical denaturants on protein structure, the molecular mechanism underlying their action is still debated. In the present review, after a brief recapitulation of the main experimental data on protein denaturants, we survey both classical and more recent interpretations of the molecular basis of their action. In particular, we highlight the differences and similarities of the impact that denaturants have on different structural classes of proteins, i.e., globular, intrinsically disordered (IDP), and amyloid-like assemblies. Particular attention has been given to the IDPs, as recent studies are unraveling their fundamental importance in many physiological processes. The role that computation techniques are expected to play in the near future is illustrated.

Learning to evolve structural ensembles of unfolded and disordered proteins using experimental solution data

The structural characterization of proteins with a disorder requires a computational approach backed by experiments to model their diverse and dynamic structural ensembles. The selection of conformational ensembles consistent with solution experiments of disordered proteins highly depends on the initial pool of conformers, with currently available tools limited by conformational sampling. We have developed a Generative Recurrent Neural Network (GRNN) that uses supervised learning to bias the probability distributions of torsions to take advantage of experimental data types such as nuclear magnetic resonance J-couplings, nuclear Overhauser effects, and paramagnetic resonance enhancements. We show that updating the generative model parameters according to the reward feedback on the basis of the agreement between experimental data and probabilistic selection of torsions from learned distributions provides an alternative to existing approaches that simply reweight conformers of a static structural pool for disordered proteins. Instead, the biased GRNN, DynamICE, learns to physically change the conformations of the underlying pool of the disordered protein to those that better agree with experiments.

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.

Inhibition of Amyloid Nucleation by Steric Hindrance

Despite recent experiments and simulations suggesting that small-molecule inhibitors and some post-translational modifications (e.g., glycosylation and ubiquitination) can suppress the pathogenic aggregation of proteins due to steric hindrance, the effect of steric hindrance on amyloid formation has not been systematically studied. Based on Monte Carlo simulations using a coarse-grained model for amyloidogenic proteins and a hard sphere acting as steric hindrance, we investigated how steric hindrance on proteins could affect amyloid formation, particularly two steps of primary nucleation, namely, oligomerization and conformational conversion into a β-sheet-enriched nucleus. We found that steric spheres played an inhibitory role in oligomerization with the effect proportional to the sphere radius R_S, which we attributed to the decline in the nonspecific attractions between proteins. During the second step, small steric spheres facilitated the conformational conversion of proteins while large ones suppressed the conversion. The overall steric effect on amyloid nucleation was inhibitory regardless of R_S. As R_S increased, oligomeric assemblies changed from amorphous into sheet-like, structurally ordered species, reminiscent of the structure of amyloid fibrils. The oligomers with large R_S were off-pathway with their ordered structures induced by the competition between steric hindrance and nonspecific attractions of soluble proteins. Interestingly, the equimolar mixture of proteins with and without steric hindrance amplified the sterically inhibitory effect by increasing the energy barrier of protein's conformational conversion. The physical mechanisms and biological implications of the above results are discussed. Our findings improve the current understanding of how nature regulates protein aggregation and amyloid formation by steric hindrance.

Low-resolution description of the conformational space for intrinsically disordered proteins

Intrinsically disordered proteins (IDP) are at the center of numerous biological processes, and attract consequently extreme interest in structural biology. Numerous approaches have been developed for generating sets of IDP conformations verifying a given set of experimental measurements. We propose here to perform a systematic enumeration of protein conformations, carried out using the TAiBP approach based on distance geometry. This enumeration was performed on two proteins, Sic1 and pSic1, corresponding to unphosphorylated and phosphorylated states of an IDP. The relative populations of the obtained conformations were then obtained by fitting SAXS curves as well as Ramachandran probability maps, the original finite mixture approach RamaMix being developed for this second task. The similarity between profiles of local gyration radii provides to a certain extent a converged view of the Sic1 and pSic1 conformational space. Profiles and populations are thus proposed for describing IDP conformations. Different variations of the resulting gyration radius between phosphorylated and unphosphorylated states are observed, depending on the set of enumerated conformations as well as on the methods used for obtaining the populations.

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.

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.

Dynamic Heterogeneity at the Interface of an Intrinsically Disordered Peptide

It is believed that water around an intrinsically disordered protein or peptide (IDP) in an aqueous environment plays an important role in guiding its conformational properties and aggregation behavior. However, despite its importance, only a handful of studies exploring the correlation between the conformational motions of an IDP and the microscopic properties of water at its surface are reported. Attempts have been made in this work to study the dynamic properties of water present in the vicinity of α-synuclein, an IDP associated with Parkinson's disease (PD). Room temperature molecular dynamics (MD) simulations of eight α-synuclein_1-95 peptides with a wide range of initial conformations have been carried out in aqueous media. The calculations revealed that due to solid-like caging motions, the translational and rotational mobility of water molecules near the surfaces of the peptide repeat unit segments R1 to R7 are significantly restricted. A small degree of dynamic heterogeneity in the hydration environment around the repeat units has been observed with water near the hydrophobic R6 unit exhibiting relatively more restricted diffusivity. The time scales involving the overall structural relaxations of peptide-water and water-water hydrogen bonds near the peptide have been found to be correlated with the time scale of diffusion of the interfacial water molecules. We believe that the relatively more hindered dynamic environment near R6 can help create water-mediated contacts centered around R6 between peptide monomers at a higher concentration, thereby enhancing the early stages of peptide aggregation.

Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation

Alzheimer’s disease is understood to be caused by amyloid fibrils and oligomers formed by aggregated amyloid-β (Aβ) peptides. This review article presents molecular dynamics (MD) simulation studies of Aβ peptides and Aβ fragments on their aggregation, aggregation inhibition, amyloid fibril conformations in equilibrium, and disruption of the amyloid fibril by ultrasonic wave and infrared laser irradiation. In the aggregation of Aβ, a β-hairpin structure promotes the formation of intermolecular β-sheet structures. Aβ peptides tend to exist at hydrophilic/hydrophobic interfaces and form more β-hairpin structures than in bulk water. These facts are the reasons why the aggregation is accelerated at the interface. We also explain how polyphenols, which are attracting attention as aggregation inhibitors of Aβ peptides, interact with Aβ. An MD simulation study of the Aβ amyloid fibrils in equilibrium is also presented: the Aβ amyloid fibril has a different structure at one end from that at the other end. The amyloid fibrils can be destroyed by ultrasonic wave and infrared laser irradiation. The molecular mechanisms of these amyloid fibril disruptions are also explained, particularly focusing on the function of water molecules. Finally, we discuss the prospects for developing treatments for Alzheimer’s disease using MD simulations.

Generating Ensembles of Dynamic Misfolding Proteins

The early stages of protein misfolding and aggregation involve disordered and partially folded protein conformers that contain a high degree of dynamic disorder. These dynamic species may undergo large-scale intra-molecular motions of intrinsically disordered protein (IDP) precursors, or flexible, low affinity inter-molecular binding in oligomeric assemblies. In both cases, generating atomic level visualization of the interconverting species that captures the conformations explored and their physico-chemical properties remains hugely challenging. How specific sub-ensembles of conformers that are on-pathway to aggregation into amyloid can be identified from their aggregation-resilient counterparts within these large heterogenous pools of rapidly moving molecules represents an additional level of complexity. Here, we describe current experimental and computational approaches designed to capture the dynamic nature of the early stages of protein misfolding and aggregation, and discuss potential challenges in describing these species because of the ensemble averaging of experimental restraints that arise from motions on the millisecond timescale. We give a perspective of how machine learning methods can be used to extract aggregation-relevant sub-ensembles and provide two examples of such an approach in which specific interactions of defined species within the dynamic ensembles of α-synuclein (αSyn) and β2-microgloblulin (β2m) can be captured and investigated.

Consequences of variability in α-synuclein fibril structure on strain biology

Synucleinopathies are a group of clinically and neuropathologically distinct protein misfolding diseases caused by unique α-synuclein conformations, or strains. While multiple atomic resolution cryo-electron microscopy structures of α-synuclein fibrils are now deposited in Protein Data Bank, significant gaps in the biological consequences arising from each conformation have yet to be unraveled. Mutations in the α-synuclein gene (SNCA), cofactors, and the solvation environment contribute to the formation and maintenance of each disease-causing strain. This review highlights the impact of each of these factors on α-synuclein misfolding and discusses the implications of the resulting structural variability on therapeutic development.

Replica permutation with solute tempering for molecular dynamics simulation and its application to the dimerization of amyloid-β fragments

We propose the replica permutation with solute tempering (RPST) by combining the replica-permutation method (RPM) and the replica exchange with solute tempering (REST). Temperature permutations are performed among more than two replicas in RPM, whereas temperature exchanges are performed between two replicas in the replica-exchange method (REM). The temperature transition in RPM occurs more efficiently than in REM. In REST, only the temperatures of the solute region, the solute temperatures, are exchanged to reduce the number of replicas compared to REM. Therefore, RPST is expected to be an improved method taking advantage of these methods. For comparison, we applied RPST, REST, RPM, and REM to two amyloid-β(16–22) peptides in explicit water. We calculated the transition ratio and the number of tunneling events in the temperature space and the number of dimerization events of amyloid-β(16–22) peptides. The results indicate that, in RPST, the number of replicas necessary for frequent random walks in the temperature and conformational spaces is reduced compared to the other three methods. In addition, we focused on the dimerization process of amyloid-β(16–22) peptides. The RPST simulation with a relatively small number of replicas shows that the two amyloid-β(16–22) peptides form the intermolecular antiparallel β-bridges due to the hydrophilic side-chain contact between Lys and Glu and hydrophobic side-chain contact between Leu, Val, and Phe, which stabilizes the dimer of the peptides.

Exploration of Protein Aggregations in Parkinson's Disease Through Computational Approaches and Big Data Analytics

Protein aggregation has been implicated in numerous neurodegenerative disorders whose etiologies are poorly understood, and for which there are no effective treatments. Here we show that the computational approaches may help us to better understand the basics of Parkinson's disease (PD). The high-resolution structural, dynamical, and mechanistic insights delivered by computational studies of protein aggregation have a unique potential to enable the rational manipulation of oligomer formation. Additionally, big data and machine learning methods may provide valuable insights to better understand the nature of proteins involved in PD and their aggregative behavior for the betterment of PD treatment.

Characterization of Amyloidogenic Peptide Aggregability in Helical Subspace

Prototypical amyloidogenic peptides amyloid-β (Aβ) and α-synuclein (αS) can undergo helix-helix associations via partially folded helical conformers, which may influence pathological progression to Alzheimer's (AD) and Parkinson's disease (PD), respectively. At the other extreme, stable folded helical conformers have been reported to resist self-assembly and amyloid formation. Experimental characterisation of such disparities in aggregation profiles due to subtle differences in peptide stabilities is precluded by the conformational heterogeneity of helical subspace. The diverse physical models used in molecular simulations allow sampling distinct regions of the phase space and are extensive in capturing the ensemble of rich helical subspace. Robust and powerful computational predictive methods utilizing network theory and free energy mapping can model the origin of helical population shifts in amyloidogenic peptides, which highlight their inherent aggregability. In this chapter, we discuss computational models, methods, design rules, and strategies to identify the driving force behind helical self-assembly and the molecular origin of aggregation resistance in helical intermediates of Aβ42 and αS. By extensive multiscale mapping of intrapeptide interactions, we show that the computational models can capture features that are otherwise imperceptible to experiments. Our models predict that targeting terminal residues may allow modulation and control of initial pathogenic aggregability of amyloidogenic peptides.

Missense Mutations Modify the Conformational Ensemble of the α-Synuclein Monomer Which Exhibits a Two-Phase Characteristic

α-Synuclein is an intrinsically disordered protein occurring in different conformations and prone to aggregate in β-sheet structures, which are the hallmark of the Parkinson disease. Missense mutations are associated with familial forms of this neuropathy. How these single amino-acid substitutions modify the conformations of wild-type α-synuclein is unclear. Here, using coarse-grained molecular dynamics simulations, we sampled the conformational space of the wild type and mutants (A30P, A53P, and E46K) of α-synuclein monomers for an effective time scale of 29.7 ms. To characterize the structures, we developed an algorithm, CUTABI (CUrvature and Torsion based of Alpha-helix and Beta-sheet Identification), to identify residues in the α-helix and β-sheet from C^α-coordinates. CUTABI was built from the results of the analysis of 14,652 selected protein structures using the Dictionary of Secondary Structure of Proteins (DSSP) algorithm. DSSP results are reproduced with 93% of success for 10 times lower computational cost. A two-dimensional probability density map of α-synuclein as a function of the number of residues in the α-helix and β-sheet is computed for wild-type and mutated proteins from molecular dynamics trajectories. The density of conformational states reveals a two-phase characteristic with a homogeneous phase (state B, β-sheets) and a heterogeneous phase (state HB, mixture of α-helices and β-sheets). The B state represents 40% of the conformations for the wild-type, A30P, and E46K and only 25% for A53T. The density of conformational states of the B state for A53T and A30P mutants differs from the wild-type one. In addition, the mutant A53T has a larger propensity to form helices than the others. These findings indicate that the equilibrium between the different conformations of the α-synuclein monomer is modified by the missense mutations in a subtle way. The α-helix and β-sheet contents are promising order parameters for intrinsically disordered proteins, whereas other structural properties such as average gyration radius, R_g, or probability distribution of R_g cannot discriminate significantly the conformational ensembles of the wild type and mutants. When separated in states B and HB, the distributions of R_g are more significantly different, indicating that global structural parameters alone are insufficient to characterize the conformational ensembles of the α-synuclein monomer.

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.

Tackling Stereochemistry in Drug Molecules with Vibrational Optical Activity

Chirality plays a crucial role in drug discovery and development. As a result, a significant number of commercially available drugs are structurally dissymmetric and enantiomerically pure. The determination of the exact 3D structure of drug candidates is, consequently, of paramount importance for the pharmaceutical industry in different stages of the discovery pipeline. Traditionally the assignment of the absolute configuration of druggable molecules has been carried out by means of X-ray crystallography. Nevertheless, not all molecules are suitable for single-crystal growing. Additionally, valuable information about the conformational dynamics of drug candidates is lost in the solid state. As an alternative, vibrational optical activity (VOA) methods have emerged as powerful tools to assess the stereochemistry of drug molecules directly in solution. These methods include vibrational circular dichroism (VCD) and Raman optical activity (ROA). Despite their potential, VCD and ROA are still unheard of to many organic and medicinal chemists. Therefore, the present review aims at highlighting the recent use of VOA methods for the assignment of the absolute configuration of chiral small-molecule drugs, as well as for the structural analysis of biologics of pharmaceutical interest. A brief introduction on VCD and ROA theory and the best experimental practices for using these methods will be provided along with selected representative examples over the last five years. As VCD and ROA are commonly used in combination with quantum calculations, some guidelines will also be presented for the reliable simulation of chiroptical spectra. Special attention will be paid to the complementarity of VCD and ROA to unambiguously assess the stereochemical properties of pharmaceuticals.

Transient Tertiary Structures of Disordered Dynein Intermediate Chain Regulate its Interactions with Multiple Partners

The N-terminal domain of dynein intermediate chain (N-IC) is central to the cytoplasmic dynein 'cargo attachment subcomplex' and regulation of motor activity. It is a prototypical intrinsically disordered protein (IDP), serving as a primarily disordered polybivalent molecular scaffold for numerous binding partners, including three dimeric dynein light chains and coiled coil domains of dynein partners dynactin p150^Glued and NudE. At the very N-terminus, a 40 amino acid single alpha helix (SAH) forms the major binding site for both p150^Glued and NudE, while a shorter nascent helix (H2) separated from SAH by a disordered linker, is necessary for tight binding to dynactin p150^Glued but not to NudE. Here we demonstrate that transient tertiary interactions in this highly dynamic protein underlie the differences in its interactions with p150^Glued and NudE. NMR paramagnetic relaxation enhancement experiments and restrained molecular dynamics simulations identify interactions between the two non-contiguous SAH and H2 helical regions, the extent of which correlates with the length and stability of H2, showing clearly that tertiary and secondary structure formation are coupled in IDPs. These interactions are significantly attenuated when N-IC is bound to NudE, suggesting that NudE binding shifts the conformational ensemble to one that is more extended and with less structure in H2. While the intrinsic disorder and flexibility in N-IC modulate its ability to serve as a binding platform for numerous partners, deviations of this protein from random-coil behavior provide a process for regulating these binding interactions and potentially the dynein motor.

α-Synuclein: An All-Inclusive Trip Around its Structure, Influencing Factors and Applied Techniques

Alpha-synuclein (αSyn) is a highly expressed and conserved protein, typically found in the presynaptic terminals of neurons. The misfolding and aggregation of αSyn into amyloid fibrils is a pathogenic hallmark of several neurodegenerative diseases called synucleinopathies, such as Parkinson’s disease. Since αSyn is an Intrinsically Disordered Protein, the characterization of its structure remains very challenging. Moreover, the mechanisms by which the structural conversion of monomeric αSyn into oligomers and finally into fibrils takes place is still far to be completely understood. Over the years, various studies have provided insights into the possible pathways that αSyn could follow to misfold and acquire oligomeric and fibrillar forms. In addition, it has been observed that αSyn structure can be influenced by different parameters, such as mutations in its sequence, the biological environment (e.g., lipids, endogenous small molecules and proteins), the interaction with exogenous compounds (e.g., drugs, diet components, heavy metals). Herein, we review the structural features of αSyn (wild-type and disease-mutated) that have been elucidated up to present by both experimental and computational techniques in different environmental and biological conditions. We believe that this gathering of current knowledge will further facilitate studies on αSyn, helping the planning of future experiments on the interactions of this protein with targeting molecules especially taking into consideration the environmental conditions.

Protein Phase Separation Arising from Intrinsic Disorder: First-Principles to Bespoke Applications

The phase separation of biomolecules has become the focus of intense research in the past decade, with a growing body of research implicating this phenomenon in essentially all biological functions, including but not limited to homeostasis, stress responses, gene regulation, cell differentiation, and disease. Excellent reviews have been published previously on the underlying physical basis of liquid-liquid phase separation (LLPS) of biological molecules (Nat. Phys. 2015, 11, 899-904) and LLPS as it occurs natively in physiology and disease (Science 2017, 357, eaaf4382; Biochemistry 2018, 57, 2479-2487; Chem. Rev. 2014, 114, 6844-6879). Here, we review how the theoretical physical basis of LLPS has been used to better understand the behavior of biomolecules that undergo LLPS in natural systems and how this understanding has also led to the development of novel synthetic systems that exhibit biomolecular phase separation, and technologies that exploit these phenomena. In part 1 of this Review, we explore the theory behind the phase separation of biomolecules and synthetic macromolecules and introduce a few notable phase-separating biomolecules. In part 2, we cover experimental and computational methods used to study phase-separating proteins and how these techniques have uncovered the mechanisms underlying phase separation in physiology and disease. Finally, in part 3, we cover the development and applications of engineered phase-separating polypeptides, ranging from control of their self-assembly to create defined supramolecular architectures to reprogramming biological processes using engineered IDPs that exhibit LLPS.

Role of Water Molecules and Helix Structure Stabilization in the Laser-Induced Disruption of Amyloid Fibrils Observed by Nonequilibrium Molecular Dynamics Simulations

Water plays a crucial role in the formation and destruction of biomolecular structures. The mechanism for destroying biomolecular structures was thought to be an active breaking of hydrogen bonds by water molecules. However, using nonequilibrium molecular dynamics simulations, in which an amyloid-β amyloid fibril was destroyed via infrared free-electron laser (IR-FEL) irradiation, we discovered a new mechanism, in which water molecules disrupt protein aggregates. The intermolecular hydrogen bonds formed by C═O and N-H in the fibril are broken at each pulse of laser irradiation. These bonds spontaneously re-form after the irradiation in many cases. However, when a water molecule happens to enter the gap between C═O and N-H, it inhibits the re-formation of the hydrogen bonds. Such sites become defects in the regularly aligned hydrogen bonds, from which all hydrogen bonds in the intermolecular β-sheet are broken as the fraying spreads. This role of water molecules is entirely different from other known mechanisms. This new mechanism can explain the recent experiments showing that the amyloid fibrils are not destroyed by laser irradiation under dry conditions. Additionally, we found that helix structures form more after the amyloid disruption; this is because the resonance frequency is different in a helix structure. Our findings provide a theoretical basis for the application of IR-FEL to the future treatment of amyloidosis.

NMR unveils an N-terminal interaction interface on acetylated-α-synuclein monomers for recruitment to fibrils

Amyloid fibril formation of α-synuclein (αS) is associated with multiple neurodegenerative diseases, including Parkinson's disease (PD). Growing evidence suggests that progression of PD is linked to cell-to-cell propagation of αS fibrils, which leads to seeding of endogenous intrinsically disordered monomer via templated elongation and secondary nucleation. A molecular understanding of the seeding mechanism and driving interactions is crucial to inhibit progression of amyloid formation. Here, using relaxation-based solution NMR experiments designed to probe large complexes, we probe weak interactions of intrinsically disordered acetylated-αS (Ac-αS) monomers with seeding-competent Ac-αS fibrils and seeding-incompetent off-pathway oligomers to identify Ac-αS monomer residues at the binding interface. Under conditions that favor fibril elongation, we determine that the first 11 N-terminal residues on the monomer form a common binding site for both fibrils and off-pathway oligomers. Additionally, the presence of off-pathway oligomers within a fibril seeding environment suppresses seeded amyloid formation, as observed through thioflavin-T fluorescence experiments. This highlights that off-pathway αS oligomers can act as an auto-inhibitor against αS fibril elongation. Based on these data taken together with previous results, we propose a model in which Ac-αS monomer recruitment to the fibril is driven by interactions between the intrinsically disordered monomer N terminus and the intrinsically disordered flanking regions (IDR) on the fibril surface. We suggest that this monomer recruitment may play a role in the elongation of amyloid fibrils and highlight the potential of the IDRs of the fibril as important therapeutic targets against seeded amyloid formation.

Membrane Interactions of α-Synuclein Revealed by Multiscale Molecular Dynamics Simulations, Markov State Models, and NMR

α-Synuclein (αS) is a presynaptic protein that binds to cell membranes and is linked to Parkinson's disease (PD). Binding of αS to membranes is a likely first step in the molecular pathophysiology of PD. The αS molecule can adopt multiple conformations, being largely disordered in water, adopting a β-sheet conformation when present in amyloid fibrils, and forming a dynamic multiplicity of α-helical conformations when bound to lipid bilayers and related membrane-mimetic surfaces. Multiscale molecular dynamics simulations in conjunction with nuclear magnetic resonance (NMR) and cross-linking mass spectrometry (XLMS) measurements are used to explore the interactions of αS with an anionic lipid bilayer. The simulations and NMR measurements together reveal a break in the helical structure of the central non-amyloid-β component (NAC) region of αS in the vicinity of residues 65-70, which may facilitate subsequent oligomer formation. Coarse-grained simulations of αS starting from the structure of αS when bound to a detergent micelle reveal the overall pattern of protein contacts to anionic lipid bilayers, while subsequent all-atom simulations provide details of conformational changes upon membrane binding. In particular, simulations and NMR data for liposome-bound αS indicate incipient β-strand formation in the NAC region, which is supported by intramolecular contacts seen via XLMS and simulations. Markov state models based on the all-atom simulations suggest a mechanism of conformational change of membrane-bound αS via a dynamic helix break in the region of residue 65 in the NAC region. The emergent dynamic model of membrane-interacting αS advances our understanding of the mechanism of PD, potentially aiding the design of novel therapeutic approaches.

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.

A lowly populated, transient β-sheet structure in monomeric Aβ1-42 identified by multinuclear NMR of chemical denaturation

Chemical denaturation is a well-established approach for probing the equilibrium between folded and unfolded states of proteins. We demonstrate applicability of this method to the detection of a small population of a transiently folded structural element in a system that is often considered to be intrinsically fully disordered. The ¹H^N, ¹⁵N, ¹³C^α, and ¹³C' chemical shifts of Aβ^1-40 and Aβ^1-42 peptides and their M35-oxidized variants were monitored as a function of urea concentration and compared to analogous urea titrations of synthetic pentapeptides of homologous sequence. Fitting of the chemical shift titrations yields a 10 ± 1% population for a structured element at the C-terminus of Aβ^1-42 that folds with a cooperativity of m = 0.06 kcal/mol·M. The fit also yields the chemical shifts of the folded state and, using a database search, for Aβ^1-42 these shifts identified an antiparallel intramolecular β-sheet for residues I32-A42, linked by a type I' β-turn at G37 and G38. The structure is destabilized by oxidation of M35. Paramagnetic relaxation rates and two previously reported weak, medium-range NOE interactions are consistent with this transient β-sheet. Introduction of the requisite A42C mutation and tagging with MTSL resulted in a small stabilization of this β-sheet. Chemical shift analysis suggests a C-terminal β-sheet may be present in Aβ^1-40 too, but the turn type at G37 is not type I'. The approach to derive Transient Structure from chemical Denaturation by NMR (TSD-NMR), demonstrated here for Aβ peptides, provides a sensitive tool for identifying the presence of lowly populated, transiently ordered elements in proteins that are considered to be intrinsically disordered, and permits extraction of structural data for such elements.

Folding Free-Energy Landscape of α-Synuclein (35-97) Via Replica Exchange Molecular Dynamics

The misfolding and aggregation of α-synuclein (α-syn) in Lewy bodies are implicated in the pathogenesis of various neurodegenerative disorders, such as Parkinson's disease and dementia. The formation of α-syn fibrils is a complex process, involving various intermediates and oligomeric forms. These intermediates establish at an early stage of aggregation and subsequently lead to fibrillation. Determining which conformations are accessible to monomeric α-syn and especially, as shown in a recent work, to the central amino acids from residue 35 to residue 97 (63 residues) is thus crucial to understand the formation of these oligomers. Here, we carry out extensive replica exchange molecular dynamics (total time-18 μs) with an all-atom model and explicit solvent to characterize the free-energy landscape of human α-syn (residue 35 to residue 97). The simulation results lead us to identify two free-energy basins. Clustering analysis for the deepest free-energy minimum reveals a compact structure, with a secondary structure predominantly α-helix, while the shallower minimum corresponds to an elongated conformation, also predominantly α-helix. Furthermore, at physiological temperature, we find that conformational rearrangements happen via helix breaks due to the presence of glycine. We also show that the most likely conformations are characterized by the α-helix structure rather than the β-hairpin structure (for residue 38 to residue 53), in contrast with prior simulation studies using coarse-grained models or an implicit solvent. For higher temperatures, we observe a shift in secondary structure with a decrease in the population of α-helix in favor of random coils, β-bend, and β-turns.

Amyloid-Like Peptide Aggregates

Self-assembly of proteins and peptides into the amyloid fold is a widespread phenomenon in the natural world. The structural hallmark of self-assembly into amyloid fibrillar assemblies is the cross-beta motif, which conveys distinct morphological and mechanical properties. The amyloid fibril formation has contrasting results depending on the organism, in the sense that it can bestow an organism with the advantages of mechanical strength and improved functionality or, on the contrary, could give rise to pathological states. In this chapter we review the existing information on amyloid-like peptide aggregates, which could either be derived from protein sequences, but also could be rationally or de novo designed in order to self-assemble into amyloid fibrils under physiological conditions. Moreover, the development of self-assembled fibrillar biomaterials that are tailored for the desired properties towards applications in biomedical or environmental areas is extensively analyzed. We also review computational studies predicting the amyloid propensity of the natural amino acid sequences and the structure of amyloids, as well as designing novel functional amyloid materials.

Capturing the Conformational Ensemble of the Mixed Folded Polyglutamine Protein Ataxin-3

Ataxin-3 is a deubiquitinase involved in protein quality control and other essential cellular functions. It preferentially interacts with polyubiquitin chains of four or more units attached to proteins delivered to the ubiquitin-proteasome system. Ataxin-3 is composed of an N-terminal Josephin domain and a flexible C terminus that contains two or three ubiquitin-interacting motifs (UIMs) and a polyglutamine tract, which, when expanded beyond a threshold, leads to protein aggregation and misfolding and causes spinocerebellar ataxia type 3. The high-resolution structure of the Josephin domain is available, but the structural and dynamical heterogeneity of ataxin-3 has so far hindered the structural description of the full-length protein. Here, we characterize non-expanded and expanded variants of ataxin-3 in terms of conformational ensembles adopted by the proteins in solution by jointly using experimental data from nuclear magnetic resonance and small-angle X-ray scattering with coarse-grained simulations. Our results pave the way to a molecular understanding of polyubiquitin recognition.

The Role of Lipids in the Initiation of α-Synuclein Misfolding

The aggregation of α-synuclein (α-syn) is inseparably connected to Parkinson’s disease (PD). It is now well-established that certain forms of α-syn aggregates, oligomers and fibrils, can exert neurotoxicity in synucleinopathies. With the exception of rare familial forms, the vast majority of PD cases are idiopathic. Understanding the earliest molecular mechanisms that cause initial α-syn misfolding could help to explain why PD affects only some individuals and others not. Factors that chaperone the transition of α-syn’s physiological to pathological function are of particular interest, since they offer opportunities for intervention. The relationship between α-syn and lipids represents one of those factors. Membrane interaction is crucial for normal cellular function, but lipids also induce the aggregation of α-syn, causing cell toxicity. Also, disease-causing or risk-factor mutations in genes related to lipid metabolism like PLA2G6, SCARB2 or GBA1 highlight the close connection between PD and lipids. Despite the clear link, the ambivalent interaction has not been studied sufficiently so far. In this review, we address how α-syn interacts with lipids and how they can act as key factor for orchestrating toxic conversion of α-syn. Furthermore, we will discuss a scenario in which initial α-syn aggregation is determined by shifts in lipid/α-syn ratio as well as by dyshomeostasis of membrane bound/unbound state of α-syn.

Dynamics of oligomer populations formed during the aggregation of Alzheimer's Aβ42 peptide

Oligomeric species populated during the aggregation of the Aβ42 peptide have been identified as potent cytotoxins linked to Alzheimer's disease, but the fundamental molecular pathways that control their dynamics have yet to be elucidated. By developing a general approach that combines theory, experiment and simulation, we reveal, in molecular detail, the mechanisms of Aβ42 oligomer dynamics during amyloid fibril formation. Even though all mature amyloid fibrils must originate as oligomers, we found that most Aβ42 oligomers dissociate into their monomeric precursors without forming new fibrils. Only a minority of oligomers converts into fibrillar structures. Moreover, the heterogeneous ensemble of oligomeric species interconverts on timescales comparable to those of aggregation. Our results identify fundamentally new steps that could be targeted by therapeutic interventions designed to combat protein misfolding diseases.

Avidity within the N-terminal anchor drives α-synuclein membrane interaction and insertion

In the brain, α-synuclein (aSN) partitions between free unbound cytosolic and membrane bound forms modulating both its physiological and pathological role and complicating its study due to structural heterogeneity. Here, we use an interdisciplinary, synergistic approach to characterize the properties of aSN:lipid mixtures, isolated aSN:lipid co-structures, and aSN in mammalian cells. Enabled by the isolation of the membrane-bound state, we show that within the previously described N-terminal membrane anchor, membrane interaction relies both on an N-terminal tail (NTT) head group layer insertion of 14 residues and a folded-upon-binding helix at the membrane surface. Both binding events must be present; if, for example, the NTT insertion is lost, the membrane affinity of aSN is severely compromised and formation of aSN:lipid co-structures hampered. In mammalian cells, compromised cooperativity results in lowered membrane association. Thus, avidity within the N-terminal anchor couples N-terminal insertion and helical surface binding, which is crucial for aSN membrane interaction and cellular localization, and may affect membrane fusion.

Unravelling the effect of N(ε)-(carboxyethyl)lysine on the conformation, dynamics and aggregation propensity of α-synuclein

α-Synuclein (αS) aggregation is a hallmark in several neurodegenerative diseases. Among them, Parkinson's disease is highlighted, characterized by the intraneuronal deposition of Lewy bodies (LBs) which causes the loss of dopaminergic neurons. αS is the main component of LBs and in them, it usually contains post-translational modifications. One of them is the formation of advanced glycation end-products (mainly CEL and MOLD) arising from its reaction with methylglyoxal. Despite its biological relevance, there are no data available proving the effect of glycation on the conformation of αS, nor on its aggregation mechanism. This has been hampered by the formation of a heterogeneous set of compounds that precluded conformational studies. To overcome this issue, we have here produced αS homogeneously glycated with CEL. Its use, together with different biophysical techniques and molecular dynamics simulations, allowed us to study for the first time the effect of glycation on the conformation of a protein. CEL extended the conformation of the N-terminal domain as a result of the loss of transient N-/C-terminal long-range contacts while increasing the heterogeneity of the conformational population. CEL also inhibited the αS aggregation, but it was not able to disassemble preexisting amyloid fibrils, thus proving that CEL found on LBs must be formed in a later event after aggregation.

We study the effect of an advanced glycation end product (N(ε)-(carboxyethyl)lysine), found on the Lewy bodies of people suffering from Parkinson’s disease, on the conformational and aggregation features of alpha-synuclein.

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.

Entropy and Information within Intrinsically Disordered Protein Regions

Bioinformatics and biophysical studies of intrinsically disordered proteins and regions (IDRs) note the high entropy at individual sequence positions and in conformations sampled in solution. This prevents application of the canonical sequence-structure-function paradigm to IDRs and motivates the development of new methods to extract information from IDR sequences. We argue that the information in IDR sequences cannot be fully revealed through positional conservation, which largely measures stable structural contacts and interaction motifs. Instead, considerations of evolutionary conservation of molecular features can reveal the full extent of information in IDRs. Experimental quantification of the large conformational entropy of IDRs is challenging but can be approximated through the extent of conformational sampling measured by a combination of NMR spectroscopy and lower-resolution structural biology techniques, which can be further interpreted with simulations. Conformational entropy and other biophysical features can be modulated by post-translational modifications that provide functional advantages to IDRs by tuning their energy landscapes and enabling a variety of functional interactions and modes of regulation. The diverse mosaic of functional states of IDRs and their conformational features within complexes demands novel metrics of information, which will reflect the complicated sequence-conformational ensemble-function relationship of IDRs.

Using a FRET Library with Multiple Probe Pairs To Drive Monte Carlo Simulations of α-Synuclein

We describe a strategy for experimentally-constraining computational simulations of intrinsically disordered proteins (IDPs), using α-synuclein, an IDP with a central role in Parkinson's disease pathology, as an example. Previously, data from single-molecule Förster Resonance Energy Transfer (FRET) experiments have been effectively utilized to generate experimentally constrained computational models of IDPs. However, the fluorophores required for single-molecule FRET experiments are not amenable to the study of short-range (<30 Å) interactions. Using ensemble FRET measurements allows one to acquire data from probes with multiple distance ranges, which can be used to constrain Monte Carlo simulations in PyRosetta. To appropriately employ ensemble FRET data as constraints, we optimized the shape and weight of constraining potentials to afford ensembles of structures that are consistent with experimental data. We also used this approach to examine the structure of α-synuclein in the presence of the compacting osmolyte trimethylamine-N-oxide. Despite significant compaction imparted by 2 M trimethylamine-N-oxide, the underlying ensemble of α-synuclein remains largely disordered and capable of aggregation, also in agreement with experimental data. These proof-of-concept experiments demonstrate that our modeling protocol enables one to efficiently generate experimentally constrained models of IDPs that incorporate atomic-scale detail, allowing one to study an IDP under a variety of conditions.

Femtosecond Hydration Map of Intrinsically Disordered α-Synuclein

Protein hydration water plays a fundamentally important role in protein folding, binding, assembly, and function. Little is known about the hydration water in intrinsically disordered proteins that challenge the conventional sequence-structure-function paradigm. Here, by combining experiments and simulations, we show the existence of dynamical heterogeneity of hydration water in an intrinsically disordered presynaptic protein, namely α-synuclein, implicated in Parkinson's disease. We took advantage of nonoccurrence of cysteine in the sequence and incorporated a number of cysteine residues at the N-terminal segment, the central amyloidogenic nonamyloid-β component (NAC) domain, and the C-terminal end of α-synuclein. We then labeled these cysteine variants using environment-sensitive thiol-active fluorophore and monitored the solvation dynamics using femtosecond time-resolved fluorescence. The site-specific femtosecond time-resolved experiments allowed us to construct the hydration map of α-synuclein. Our results show the presence of three dynamically distinct types of water: bulk, hydration, and confined water. The amyloidogenic NAC domain contains dynamically restrained water molecules that are strikingly different from the water molecules present in the other two domains. Atomistic molecular dynamics simulations revealed longer residence times for water molecules near the NAC domain and supported our experimental observations. Additionally, our simulations allowed us to decipher the molecular origin of the dynamical heterogeneity of water in α-synuclein. These simulations captured the quasi-bound water molecules within the NAC domain originating from a complex interplay between the local chain compaction and the sequence composition. Our findings from this synergistic experimental simulation approach suggest longer trapping of interfacial water molecules near the amyloidogenic hotspot that triggers the pathological conversion into amyloids via chain sequestration, chain desolvation, and entropic liberation of ordered water molecules.

Conformational ensemble of native α-synuclein in solution as determined by short-distance crosslinking constraint-guided discrete molecular dynamics simulations

Combining structural proteomics experimental data with computational methods is a powerful tool for protein structure prediction. Here, we apply a recently-developed approach for de novo protein structure determination based on the incorporation of short-distance crosslinking data as constraints in discrete molecular dynamics simulations (CL-DMD) for the determination of conformational ensemble of the intrinsically disordered protein α-synuclein in the solution. The predicted structures were in agreement with hydrogen-deuterium exchange, circular dichroism, surface modification, and long-distance crosslinking data. We found that α-synuclein is present in solution as an ensemble of rather compact globular conformations with distinct topology and inter-residue contacts, which is well-represented by movements of the large loops and formation of few transient secondary structure elements. Non-amyloid component and C-terminal regions were consistently found to contain β-structure elements and hairpins.

As the population ages, neurodegenerative diseases such as Parkinson’s disease will become an increasing problem in many countries. Aggregation of the protein α-synuclein is the primary cause of Parkinson’s disease, but there is still a dearth of structural information pertaining to the native, non-aggregating form of this protein. A better understanding the structural state of the native protein may prove useful for the design of new therapeutics to combat this disease. In order to obtain more structural information on this protein, we have recently modelled the native α-synuclein protein. These models were generated using a novel approach which combines protein crosslinking and discrete molecular dynamics simulations. We have found that the α-synuclein protein can adopt several shapes, all with a similar topology, resembling a three fingered closed claw. A region of the protein important for aggregation was found to be protected from the surrounding biological environment in these conformations, and the stabilization of these structures may be a fruitful avenue for future drug research into mitigating the cause and effect of Parkinson’s disease.

Dynamic Studies on Intrinsically Disordered Regions of Two Paralogous Transcription Factors Reveal Rigid Segments with Important Biological Functions

Long stretches of intrinsically disordered regions (IDRs) are abundantly present in eukaryotic transcription factors. Although their biological significance is well appreciated, the underlying structural and dynamic mechanisms of their function are still not clear. Using solution NMR spectroscopy, we have studied the structural and dynamic features of two paralogous HOX transcription factors, SCR and DFD, from Drosophila. Both proteins have a conserved DNA-binding homeodomain and a long stretch of functionally important IDR. Using NMR dynamics, we determined flexibility of each residue in these proteins. The flexibility of the residues in the disordered region is not uniform. In both proteins, the IDRs have short stretches of consecutive residues with relatively less flexibility, that is, higher rigidity. We show that one such rigid segment is specifically recognized by another co-transcription factor, thus highlighting the importance of these rigid segments in IDR-mediated protein-protein interactions. Using molecular dynamics simulation, we further show that the rigid segments sample less conformations compared to the rest of the residues in the disordered region. The restrained conformational sampling of these rigid residues should lower the loss in conformational entropy during their interactions with binding partners resulting in sequence specific binding. This work provides experimental evidence of a "rigid-segment" model of IDRs, where functionally important rigid segments are connected by highly flexible linkers. Furthermore, a comparative study of IDRs in paralogous proteins reveals that in spite of low-sequence conservation, the rigid and flexible segments are sequentially maintained to preserve related functions and regulations of these proteins.

In-Cell NMR Spectroscopy of Intrinsically Disordered Proteins

This review summarizes the results of in-cell Nuclear Magnetic Resonance, NMR, spectroscopic investigations of the eukaryotic and prokaryotic intrinsically disordered proteins, IDPs: α-synuclein, prokaryotic ubiquitin-like protein, Pup, tubulin-related neuronal protein, Tau, phenylalanyl-glycyl-repeat-rich nucleoporins, FG Nups, and the negative regulator of flagellin synthesis, FlgM. The results show that the cellular behavior of IDPs may differ significantly from that observed in the test tube.

AWSEM-IDP: A Coarse-Grained Force Field for Intrinsically Disordered Proteins

The associative memory, water-mediated, structure and energy model (AWSEM) has been successfully used to study protein folding, binding, and aggregation problems. In this work, we introduce AWSEM-IDP, a new AWSEM branch for simulating intrinsically disordered proteins (IDPs), where the weights of the potentials determining secondary structure formation have been finely tuned, and a novel potential is introduced that helps to precisely control both the average extent of protein chain collapse and the chain's fluctuations in size. AWSEM-IDP can efficiently sample large conformational spaces, while retaining sufficient molecular accuracy to realistically model proteins. We applied this new model to two IDPs, demonstrating that AWSEM-IDP can reasonably well reproduce higher-resolution reference data, thus providing the foundation for a transferable IDP force field. Finally, we used thermodynamic perturbation theory to show that, in general, the conformational ensembles of IDPs are highly sensitive to fine-tuning of force field parameters.

Characterization of intrinsically disordered proteins and their dynamic complexes: From in vitro to cell-like environments

Over the last two decades, it has become increasingly clear that a large fraction of the human proteome is intrinsically disordered or contains disordered segments of significant length. These intrinsically disordered proteins (IDPs) play important regulatory roles throughout biology, underlining the importance of understanding their conformational behavior and interaction mechanisms at the molecular level. Here we review recent progress in the NMR characterization of the structure and dynamics of IDPs in various functional states and environments. We describe the complementarity of different NMR parameters for quantifying the conformational propensities of IDPs in their isolated and phosphorylated states, and we discuss the challenges associated with obtaining structural models of dynamic protein-protein complexes involving IDPs. In addition, we review recent progress in understanding the conformational behavior of IDPs in cell-like environments such as in the presence of crowding agents, in membrane-less organelles and in the complex environment of the human cell.