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

⍺-Synuclein levels in Parkinson's disease - Cell types and forms that contribute to pathogenesis

Parkinson's disease (PD) has two main pathological hallmarks, the loss of nigral dopamine neurons and the proteinaceous aggregations of ⍺-synuclein (⍺Syn) in neuronal Lewy pathology. These two co-existing features suggest a causative association between ⍺Syn aggregation and the underpinning mechanism of neuronal degeneration in PD. Both increased levels and post-translational modifications of ⍺Syn can contribute to the formation of pathological aggregations of ⍺Syn in neurons. Recent studies have shown that the protein is also expressed by multiple types of non-neuronal cells in the brain and peripheral tissues, suggesting additional roles of the protein and potential diversity in non-neuronal pathogenic triggers. It is important to determine (1) the threshold levels triggering ⍺Syn to convert from a biological to a pathologic form in different brain cells in PD; (2) the dominant form of pathologic ⍺Syn and the associated post-translational modification of the protein in each cell type involved in PD; and (3) the cell type associated biological processes impacted by pathologic ⍺Syn in PD. This review integrates these aspects and speculates on potential pathological mechanisms and their impact on neuronal and non-neuronal ⍺Syn in the brains of patients with PD.

Dissecting diazirine photo-reaction mechanism for protein residue-specific cross-linking and distance mapping

While photo-cross-linking (PXL) with alkyl diazirines can provide stringent distance restraints and offer insights into protein structures, unambiguous identification of cross-linked residues hinders data interpretation to the same level that has been achieved with chemical cross-linking (CXL). We address this challenge by developing an in-line system with systematic modulation of light intensity and irradiation time, which allows for a quantitative evaluation of diazirine photolysis and photo-reaction mechanism. Our results reveal a two-step pathway with mainly sequential generation of diazo and carbene intermediates. Diazo intermediate preferentially targets buried polar residues, many of which are inaccessible with known CXL probes for their limited reactivity. Moreover, we demonstrate that tuning light intensity and duration enhances selectivity towards polar residues by biasing diazo-mediated cross-linking reactions over carbene ones. This mechanistic dissection unlocks the full potential of PXL, paving the way for accurate distance mapping against protein structures and ultimately, unveiling protein dynamic behaviors.

Structural Elucidation of Ubiquitin via Gas-Phase Ion/Ion Cross-Linking Reactions Using Sodium-Cationized Reagents Coupled with Infrared Multiphoton Dissociation

Accurate structural determination of proteins is critical to understanding their biological functions and the impact of structural disruption on disease progression. Gas-phase cross-linking mass spectrometry (XL-MS) via ion/ion reactions between multiply charged protein cations and singly charged cross-linker anions has previously been developed to obtain low-resolution structural information on proteins. This method significantly shortens experimental time relative to conventional solution-phase XL-MS but has several technical limitations: (1) the singly deprotonated N-hydroxysulfosuccinimide (sulfo-NHS)-based cross-linker anions are restricted to attachment at neutral amine groups of basic amino acid residues and (2) analyzing terminal cross-linked fragment ions is insufficient to unambiguously localize sites of linker attachment. Herein, we demonstrate enhanced structural information for alcohol-denatured A-state ubiquitin obtained from an alternative gas-phase XL-MS approach. Briefly, singly sodiated ethylene glycol bis(sulfosuccinimidyl succinate) (sulfo-EGS) cross-linker anions enable covalent cross-linking at both ammonium and amine groups. Additionally, covalently modified internal fragment ions, along with terminal b-/y-type counterparts, improve the determination of linker attachment sites. Molecular dynamics simulations validate experimentally obtained gas-phase conformations of denatured ubiquitin. This method has identified four cross-linking sites across 8+ ubiquitin, including two new sites in the N-terminal region of the protein that were originally inaccessible in prior gas-phase XL approaches. The two N-terminal cross-linking sites suggest that the N-terminal half of ubiquitin is more compact in gas-phase conformations. By comparison, the two C-terminal linker sites indicate the signature transformation of this region of the protein from a native to a denatured conformation. Overall, the results suggest that the solution-phase secondary structures of the A-state ubiquitin are conserved in the gas phase. This method also provides sufficient sensitivity to differentiate between two gas-phase conformers of the same charge state with subtle structural variations.

The Power of Molecular Dynamics Simulations and Their Applications to Discover Cysteine Protease Inhibitors

A large family of enzymes with the function of hydrolyzing peptide bonds, called peptidases or cysteine proteases (CPs), are divided into three categories according to the peptide chain involved. CPs catalyze the hydrolysis of amide, ester, thiol ester, and thioester peptide bonds. They can be divided into several groups, such as papain-like (CA), viral chymotrypsin-like CPs (CB), papainlike endopeptidases of RNA viruses (CC), legumain-type caspases (CD), and showing active residues of His, Glu/Asp, Gln, Cys (CE). The catalytic mechanism of CPs is the essential cysteine residue present in the active site. These mechanisms are often studied through computational methods that provide new information about the catalytic mechanism and identify inhibitors. The role of computational methods during drug design and development stages is increasing. Methods in Computer-Aided Drug Design (CADD) accelerate the discovery process, increase the chances of selecting more promising molecules for experimental studies, and can identify critical mechanisms involved in the pathophysiology and molecular pathways of action. Molecular dynamics (MD) simulations are essential in any drug discovery program due to their high capacity for simulating a physiological environment capable of unveiling significant inhibition mechanisms of new compounds against target proteins, especially CPs. Here, a brief approach will be shown on MD simulations and how the studies were applied to identify inhibitors or critical information against cysteine protease from several microorganisms, such as Trypanosoma cruzi (cruzain), Trypanosoma brucei (rhodesain), Plasmodium spp. (falcipain), and SARS-CoV-2 (Mpro). We hope the readers will gain new insights and use our study as a guide for potential compound identifications using MD simulations.

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.

Advanced Techniques for Detecting Protein Misfolding and Aggregation in Cellular Environments

Protein misfolding and aggregation, a key contributor to the progression of numerous neurodegenerative diseases, results in functional deficiencies and the creation of harmful intermediates. Detailed visualization of this misfolding process is of paramount importance for improving our understanding of disease mechanisms and for the development of potential therapeutic strategies. While in vitro studies using purified proteins have been instrumental in delivering significant insights into protein misfolding, the behavior of these proteins in the complex milieu of living cells often diverges significantly from such simplified environments. Biomedical imaging performed in cell provides cellular-level information with high physiological and pathological relevance, often surpassing the depth of information attainable through in vitro methods. This review highlights a variety of methodologies used to scrutinize protein misfolding within biological systems. This includes optical-based methods, strategies leaning on mass spectrometry, in-cell nuclear magnetic resonance, and cryo-electron microscopy. Recent advancements in these techniques have notably deepened our understanding of protein misfolding processes and the features of the resulting misfolded species within living cells. The progression in these fields promises to catalyze further breakthroughs in our comprehension of neurodegenerative disease mechanisms and potential therapeutic interventions.

Saturation mutagenesis of α-synuclein reveals monomer fold that modulates aggregation

α-Synuclein (aSyn) aggregation underlies neurodegenerative synucleinopathies. aSyn seeds are proposed to replicate and propagate neuronal pathology like prions. Seeding of aSyn can be recapitulated in cellular systems of aSyn aggregation; however, the mechanism of aSyn seeding and its regulation are not well understood. We developed an mEos-based aSyn seeding assay and performed saturation mutagenesis to identify with single-residue resolution positive and negative regulators of aSyn aggregation. We not only found the core regions that govern aSyn aggregation but also identified mutants outside of the core that enhance aggregation. We identified local structure within the N terminus of aSyn that hinders the fibrillization propensity of its aggregation-prone core. Based on the screen, we designed a minimal aSyn fragment that shows a ~4-fold enhancement in seeding activity and enabled discrimination of synucleinopathies. Our study expands the basic knowledge of aSyn aggregation and advances the design of cellular systems of aSyn aggregation to diagnose synucleinopathies based on protein conformation.

Exploring Intra- and Inter-Regional Interactions in the IDP α-Synuclein Using smFRET and MD Simulations

Theoretical concepts from polymer physics are often used to describe intrinsically disordered proteins (IDPs). However, amino acid interactions within and between regions of the protein can lead to deviations from typical polymer scaling behavior and even to short-lived secondary structures. To investigate the key interactions in the dynamic IDP α-synuclein (αS) at the amino acid level, we conducted single-molecule fluorescence resonance energy transfer (smFRET) experiments and coarse-grained molecular dynamics (CG-MD) simulations. We find excellent agreement between experiments and simulations. Our results show that a physiological salt solution is a good solvent for αS and that the protein is highly dynamic throughout its entire chain, with local intra- and inter-regional interactions leading to deviations from global scaling. Specifically, we observe expansion in the C-terminal region, compaction in the NAC region, and a slightly smaller distance between the C- and N-termini than expected. Our simulations indicate that the compaction in the NAC region results from hydrophobic aliphatic contacts, mostly between valine and alanine residues, and cation-π interactions between lysine and tyrosine. In addition, hydrogen bonds also seem to contribute to the compaction of the NAC region. The expansion of the C-terminal region is due to intraregional electrostatic repulsion and increased chain stiffness from several prolines. Overall, our study demonstrates the effectiveness of combining smFRET experiments with CG-MD simulations to investigate the key interactions in highly dynamic IDPs at the amino acid level.

Big versus small: The impact of aggregate size in disease

Protein aggregation results in an array of different size soluble oligomers and larger insoluble fibrils. Insoluble fibrils were originally thought to cause neuronal cell deaths in neurodegenerative diseases due to their prevalence in tissue samples and disease models. Despite recent studies demonstrating the toxicity associated with soluble oligomers, many therapeutic strategies still focus on fibrils or consider all types of aggregates as one group. Oligomers and fibrils require different modeling and therapeutic strategies, targeting the toxic species is crucial for successful study and therapeutic development. Here, we review the role of different-size aggregates in disease, and how factors contributing to aggregation (mutations, metals, post-translational modifications, and lipid interactions) may promote oligomers opposed to fibrils. We review two different computational modeling strategies (molecular dynamics and kinetic modeling) and how they are used to model both oligomers and fibrils. Finally, we outline the current therapeutic strategies targeting aggregating proteins and their strengths and weaknesses for targeting oligomers versus fibrils. Altogether, we aim to highlight the importance of distinguishing the difference between oligomers and fibrils and determining which species is toxic when modeling and creating therapeutics for protein aggregation in disease.

Dissecting the Self-assembly Dynamics of Imperfect Repeats in α-Synuclein

The pathological aggregation of α-synuclein (αS) into amyloid fibrils is the hallmark of Parkinson's disease (PD). The self-assembly and membrane interactions of αS are mainly governed by the seven imperfect 11-residue repeats of the XKTKEGVXXXX motif around residues 1-95. However, the particular role of each repeat in αS fibrillization remains unclear. To answer this question, we studied the aggregation dynamics of each repeat with up to 10 peptides in silico by conducting multiple independent micro-second atomistic discrete molecular dynamics simulations. Our simulations revealed that only repeats R3 and R6 readily self-assembled into β-sheet-rich oligomers, while the other repeats remained as unstructured monomers with weak self-assembly and β-sheet propensities. The self-assembly process of R3 featured frequent conformational changes with β-sheet formation mainly in the non-conserved hydrophobic tail, whereas R6 spontaneously self-assembled into extended and stable cross-β structures. These results of seven repeats are consistent with their structures and organization in recently solved αS fibrils. As the primary amyloidogenic core, R6 was buried inside the central cross-β core of all αS fibrils, attracting the hydrophobic tails of adjacent R4, R5, and R7 repeats forming β-sheets around R6 in the core. Further away from R6 in the sequence but with a moderate amyloid aggregation propensity, the R3 tail could serve as a secondary amyloidogenic core and form independent β-sheets in the fibril. Overall, our results demonstrate the critical role of R3 and R6 repeats in αS amyloid aggregation and suggest their potential as targets for the peptide-based and small-molecule amyloid inhibitors.

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.

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.

Taking Charge: Metal Ions Accelerate Amyloid Aggregation in Sequence Variants of α-Synuclein

Αlpha-synuclein (αS) is an intrinsically disordered protein which exhibits a high degree of conformational heterogeneity. In vivo, αS experiences various environments which cause adaptation of its structural ensemble. Divalent metal ions are prominent in synaptic terminals where αS is located and are thought to bind to the αS C-terminal region. Herein, we used native nanoelectrospray ionization ion mobility-mass spectrometry to investigate changes in the charge state distribution and collision cross sections of wild-type N-terminally acetylated (NTA) αS, along with a deletion variant (ΔΔNTA) which inhibits amyloid formation and a C-terminal truncated variant (119NTA) which increases the rate of amyloid formation. We also examine the effect of the addition of divalent metal ions, Ca²⁺, Mn²⁺, and Zn²⁺, and correlate the conformational properties of the αS monomer with the ability to aggregate into amyloid, measured using Thioflavin T fluorescence and negative stain transmission electron microscopy. We find a correlation between the population of species with a low collision cross section and accelerated amyloid assembly kinetics, with the presence of metal ions resulting in protein compaction and causing ΔΔ to regain its ability to form an amyloid. The results portray how the αS conformational ensemble is governed by specific intramolecular interactions that influence its amyloidogenic behavior.

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.

Sequence-Dependent Backbone Dynamics of Intrinsically Disordered Proteins

For intrinsically disordered proteins (IDPs), a pressing question is how sequence codes for function. Dynamics serves as a crucial link, reminiscent of the role of structure in sequence-function relations of structured proteins. To define general rules governing sequence-dependent backbone dynamics, we carried out long molecular dynamics simulations of eight IDPs. Blocks of residues exhibiting large amplitudes in slow dynamics are rigidified by local inter-residue interactions or secondary structures. A long region or an entire IDP can be slowed down by long-range contacts or secondary-structure packing. On the other hand, glycines promote fast dynamics and either demarcate rigid blocks or facilitate multiple modes of local and long-range inter-residue interactions. The sequence-dependent backbone dynamics endows IDPs with versatile response to binding partners, with some blocks recalcitrant while others readily adapting to intermolecular interactions.

Molecular insights into the oligomerization dynamics and conformations of amyloidogenic and non-amyloidogenic amylin from discrete molecular dynamics simulations

The amyloid aggregation of human islet amyloid polypeptide (hIAPP) is associated with pancreatic β-cell death in type 2 diabetes. The S20G substitution of hIAPP (hIAPP(S20G)), found in Japanese and Chinese people, is more amyloidogenic and cytotoxic than wild-type hIAPP. Rat amylin (rIAPP) does not have aggregation propensity or cytotoxicity. Mounting evidence suggests that soluble low-molecular-weight amyloid oligomers formed during early aggregation are more cytotoxic than mature fibrils. The self-assembly dynamics and oligomeric conformations remain unknown because the oligomers are heterogeneous and transient. The molecular mechanism of sequence-variation rendering dramatically different aggregation propensity and cytotoxicity is also elusive. Here, we investigate the oligomerization dynamics and conformations of amyloidogenic hIAPP, hIAPP(S20G), and non-amyloidogenic rIAPP using atomistic discrete molecular dynamics (DMD) simulations. Our simulation results demonstrated that all three monomeric amylin peptides mainly adopted an unstructured formation with partial dynamical helices near the N-terminus. Relatively transient β-hairpins were more abundant in hIAPP and hIAPP(S20G) than in rIAPP. The S20G-substituting mutant of hIAPP altered the turn region of the β-hairpin motif, resulting in more hydrophobic residue-pairwise contacts within the β-hairpin. Oligomerization dynamic investigation revealed that all three peptides spontaneously accumulated into helix-populated oligomers. The conformational conversion to form β-sheet-rich oligomers was only observed in hIAPP and hIAPP(S20G). The population of high-β-sheet-content oligomers was enhanced by S20G substitution. Interestingly, both hIAPP and hIAPP(S20G) could form β-barrel formations, and the β-barrel propensity of hIAPP(S20G) was three times larger than that of hIAPP. No β-sheet-rich or β-barrel formations were observed in rIAPP. Our direct observation of the correlation between β-barrel oligomer formation and cytotoxicity suggests that β-barrels might play a critically important role in the cytotoxicity of amyloidosis.

Exploring Structural Flexibility and Stability of α-Synuclein by the Landau-Ginzburg-Wilson Approach

α-Synuclein (αS) is the principal protein component of the Lewy body and Lewy neurite deposits that are found in the brains of the victims of one of the most prevalent neurodegenerative disorders, Parkinson's disease. αS can be qualified as a chameleon protein because of the large number of different conformations that it is able to adopt: it is disordered under physiological conditions in solution, in equilibrium with a minor α-helical tetrameric form in the cytoplasm, and is α-helical when bound to a cell membrane. Also, in vitro, αS forms polymorphic amyloid fibrils with unique arrangements of cross-β-sheet motifs. Therefore, it is of interest to elucidate the origins of the structural flexibility of αS and what makes αS stable in different conformations. We address these questions here by analyzing the experimental structures of the micelle-bound, tetrameric, and fibrillar αS in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger equation. It is illustrated that without molecular dynamics simulations the kinks are capable of identifying the key residues causing structural flexibility of αS. Also, the stability of the experimental structures of αS is investigated by simulating heating/cooling trajectories using the Glauber algorithm. The findings are consistent with experiments.

Single residue modulators of amyloid formation in the N-terminal P1-region of α-synuclein

Alpha-synuclein (αSyn) is a protein involved in neurodegenerative disorders including Parkinson’s disease. Amyloid formation of αSyn can be modulated by the ‘P1 region’ (residues 36-42). Here, mutational studies of P1 reveal that Y39A and S42A extend the lag-phase of αSyn amyloid formation in vitro and rescue amyloid-associated cytotoxicity in C. elegans. Additionally, L38I αSyn forms amyloid fibrils more rapidly than WT, L38A has no effect, but L38M does not form amyloid fibrils in vitro and protects from proteotoxicity. Swapping the sequence of the two residues that differ in the P1 region of the paralogue γSyn to those of αSyn did not enhance fibril formation for γSyn. Peptide binding experiments using NMR showed that P1 synergises with residues in the NAC and C-terminal regions to initiate aggregation. The remarkable specificity of the interactions that control αSyn amyloid formation, identifies this region as a potential target for therapeutics, despite their weak and transient nature.

The authors of this work characterize the effect of amino acid substitution on α-synuclein (α-Syn) aggregation. Residues 38 and 42 (in addition to 39) within the P1 region of α-Syn affect amyloid formation. The effect of substitution at position 38 is dependent on the amino-acid introduced, suggesting that specific interactions control α -Syn aggregation.

Molecular dynamics simulations of an α-synuclein NAC domain fragment with a ff14IDPSFF IDP-specific force field suggest β-sheet intermediate states of fibrillation

For the discovery of treatments against synucleinopathies, it is necessary to unravel and fully understand the mechanism of fibrillation of proteins involved. Among them, α-synuclein (αS) plays a key role in the development of these diseases through its aggregation into oligomers found in Lewy bodies. However, its structural disorder as an intrinsically disordered protein (IDP) makes its characterization by experimental techniques arduously difficult. Atomistic simulations aim to provide insights into this blank canvas and, fortunately, some studies have already suggested promising mechanisms. Still, it is urgent to consider the IDP features in simulations, so recently a lot of force fields designed to deal with IDPs have been developed. In this study, we have carried out a total of 12 μs simulations of an αS core fragment using a popular ff14SB AMBER force field and the ff14IDPSFF variation that includes a grid-based energy correction map (CMAP) method. The predicted chemical shifts from the simulations and those measured from the αS protein in the NMR solution indicate that ff14IDPSFF reproduces the experimental data more accurately. Moreover, structural analysis exhibits opposite trends between secondary structure propensities. The ff14SB force field preserves the α-helices found in the micelle-bound αS structure, which is used as an initial conformation, while ff14IDPSFF stands out with increased structural disorder and the formation of β-sheets, which suggests that the IDP-specific force field can capture more suitable conformations representing the possible intermediate states of the fibrillation process.

Wide-Line NMR Melting Diagrams, Their Thermodynamic Interpretation, and Secondary Structure Predictions for A30P and E46K α-Synuclein

Parkinson's disease is thought to be caused by aggregation of the intrinsically disordered protein, α-synuclein. Two amyloidogenic variants, A30P, and E46K familial mutants were investigated by wide-line ¹H NMR spectrometry as a completion of our earlier work on wild-type and A53T α-synuclein (Bokor M. et al. WT and A53T α-synuclein systems: melting diagram and its new interpretation. Int. J. Mol. Sci.2020, 21, 3997.). A monolayer of mobile water molecules hydrates A30P α-synuclein at the lowest potential barriers (temperatures), while E46K α-synuclein has here third as much mobile hydration, insufficient for functionality. According to wide-line ¹H NMR results and secondary structure predictions, E46K α-synuclein is more compact than the A30P variant and they are more compact than the wild type (WT) and A53T variants. Linear hydration vs potential barrier sections of A30P α-synuclein shows one and E46K shows two slopes. The different slopes of the latter between potential barriers E _a,1 and E _a,2 reflect a change in water-protein interactions. The 31-32% homogeneous potential barrier distribution of the protein-water bonds refers to a non-negligible amount of secondary structures in all four α-synuclein variants. The secondary structures detected by wide-line ¹H NMR are solvent-exposed α-helices, which are predicted by secondary structure models. β-sheets are only minor components of the protein structures as three- and eight-state predicted secondary structures are dominated by α-helices and coils.

Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena

Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.

Computational Models for the Study of Protein Aggregation

Protein aggregation has been studied by many groups around the world for many years because it can be the cause of a number of neurodegenerative diseases that have no effective treatment. Obtaining the structure of related fibrils and toxic oligomers, as well as describing the pathways and main factors that govern the self-organization process, is of paramount importance, but it is also very difficult. To solve this problem, experimental and computational methods are often combined to get the most out of each method. The effectiveness of the computational approach largely depends on the construction of a reasonable molecular model. Here we discussed different versions of the four most popular all-atom force fields AMBER, CHARMM, GROMOS, and OPLS, which have been developed for folded and intrinsically disordered proteins, or both. Continuous and discrete coarse-grained models, which were mainly used to study the kinetics of aggregation, are also summarized.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Protein Chemistry Combined with Mass Spectrometry for Protein Structure Determination

The advent of soft-ionization mass spectrometry for biomolecules has opened up new possibilities for the structural analysis of proteins. Combining protein chemistry methods with modern mass spectrometry has led to the emergence of the distinct field of structural proteomics. Multiple protein chemistry approaches, such as surface modification, limited proteolysis, hydrogen-deuterium exchange, and cross-linking, provide diverse and often orthogonal structural information on the protein systems studied. Combining experimental data from these various structural proteomics techniques provides a more comprehensive examination of the protein structure and increases confidence in the ultimate findings. Here, we review various types of experimental data from structural proteomics approaches with an emphasis on the use of multiple complementary mass spectrometric approaches to provide experimental constraints for the solving of protein structures.

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.

The characterization of protein interactions - what, how and how much?

Protein interactions underlie most molecular events in biology. Many methods have been developed to identify protein partners, to measure the affinity with which these biomolecules interact and to characterise the structures of the complexes. Each approach has its own advantages and limitations, and it can be difficult for the newcomer to determine which methodology would best suit their system. This review provides an overview of many of the techniques most widely used to identify protein partners, assess stoichiometry and binding affinity, and determine low-resolution models for complexes. Key methods covered include: yeast two-hybrid analysis, affinity purification mass spectrometry and proximity labelling to identify partners; size-exclusion chromatography, scattering methods, native mass spectrometry and analytical ultracentrifugation to estimate stoichiometry; isothermal titration calorimetry, biosensors and fluorometric methods (including microscale thermophoresis, anisotropy/polarisation, resonance energy transfer, AlphaScreen, and differential scanning fluorimetry) to measure binding affinity; and crosslinking and hydrogen-deuterium exchange mass spectrometry to probe the structure of complexes.

Advances in Hydrogen/Deuterium Exchange Mass Spectrometry and the Pursuit of Challenging Biological Systems

Solution-phase hydrogen/deuterium exchange (HDX) coupled to mass spectrometry (MS) is a widespread tool for structural analysis across academia and the biopharmaceutical industry. By monitoring the exchangeability of backbone amide protons, HDX-MS can reveal information about higher-order structure and dynamics throughout a protein, can track protein folding pathways, map interaction sites, and assess conformational states of protein samples. The combination of the versatility of the hydrogen/deuterium exchange reaction with the sensitivity of mass spectrometry has enabled the study of extremely challenging protein systems, some of which cannot be suitably studied using other techniques. Improvements over the past three decades have continually increased throughput, robustness, and expanded the limits of what is feasible for HDX-MS investigations. To provide an overview for researchers seeking to utilize and derive the most from HDX-MS for protein structural analysis, we summarize the fundamental principles, basic methodology, strengths and weaknesses, and the established applications of HDX-MS while highlighting new developments and applications.

Structure of prion β-oligomers as determined by short-distance crosslinking constraint-guided discrete molecular dynamics simulations

The conversion of the native monomeric cellular prion protein (PrP^C ) into an aggregated pathological β-oligomeric form (PrP^β ) and an infectious form (PrP^Sc ) is the central element in the development of prion diseases. The structure of the aggregates and the molecular mechanisms of the conformational changes involved in the conversion are still unknown. We applied mass spectrometry combined with chemical crosslinking, hydrogen/deuterium exchange, limited proteolysis, and surface modification for the differential characterization of the native and the urea+acid-converted prion β-oligomer structures to obtain insights into the mechanisms of conversion and aggregation. For the determination of the structure of the monomer and the dimer unit of the β-oligomer, we applied a recently-developed approach for de novo protein structure determination which is based on the incorporation of zero-length and short-distance crosslinking data as intra- and inter-protein constraints in discrete molecular dynamics simulations (CL-DMD). Based on all of the structural-proteomics experimental data and the computationally predicted structures of the monomer units, we propose the potential mode of assembly of the β-oligomer. The proposed β-oligomer assembly provides a clue on the β-sheet nucleation site, and how template-based conversion of the native prion molecule occurs, growth of the prion aggregates, and maturation into fibrils may occur.

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.

α-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.

FRET-based dynamic structural biology: Challenges, perspectives and an appeal for open-science practices

Single-molecule FRET (smFRET) has become a mainstream technique for studying biomolecular structural dynamics. The rapid and wide adoption of smFRET experiments by an ever-increasing number of groups has generated significant progress in sample preparation, measurement procedures, data analysis, algorithms and documentation. Several labs that employ smFRET approaches have joined forces to inform the smFRET community about streamlining how to perform experiments and analyze results for obtaining quantitative information on biomolecular structure and dynamics. The recent efforts include blind tests to assess the accuracy and the precision of smFRET experiments among different labs using various procedures. These multi-lab studies have led to the development of smFRET procedures and documentation, which are important when submitting entries into the archiving system for integrative structure models, PDB-Dev. This position paper describes the current 'state of the art' from different perspectives, points to unresolved methodological issues for quantitative structural studies, provides a set of 'soft recommendations' about which an emerging consensus exists, and lists openly available resources for newcomers and seasoned practitioners. To make further progress, we strongly encourage 'open science' practices.

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.

Reinforcement learning to boost molecular docking upon protein conformational ensemble

Intrinsically disordered proteins (IDPs) are widely involved in human diseases and thus are attractive therapeutic targets. In practice, however, it is computationally prohibitive to dock large ligand libraries to thousands and tens of thousands of conformations. Here, we propose a reversible upper confidence bound (UCB) algorithm for the virtual screening of IDPs to address the influence of the conformation ensemble. The docking process is dynamically arranged so that attempts are focused near the boundary to separate top ligands from the bulk accurately. It is demonstrated in the example of transcription factor c-Myc that the average docking number per ligand can be greatly reduced while the performance is merely slightly affected. This study suggests that reinforcement learning is highly efficient in solving the bottleneck of virtual screening due to the conformation ensemble in the rational drug design of IDPs.

Unfolding bovine α-lactalbumin with T-jump: Characterizing disordered intermediates via time-resolved x-ray solution scattering and molecular dynamics simulations

The protein folding process often proceeds through partially folded transient states. Therefore, a structural understanding of these disordered states is crucial for developing mechanistic models of the folding process. Characterization of unfolded states remains challenging due to their disordered nature, and incorporating multiple methods is necessary. Combining the time-resolved x-ray solution scattering (TRXSS) signal with molecular dynamics (MD), we are able to characterize transient partially folded states of bovine α-lactalbumin, a model system widely used for investigation of molten globule states, during its unfolding triggered by a temperature jump. We track the unfolding process between 20 µs and 70 ms and demonstrate that it passes through three distinct kinetic states. The scattering signals associated with these transient species are then analyzed with TRXSS constrained MD simulations to produce protein structures that are compatible with the input signals. Without utilizing any experimentally extracted kinetic information, the constrained MD simulation successfully drove the protein to an intermediate molten globule state; signals for two later disordered states are refined to terminal unfolded states. From our examination of the structural characteristics of these disordered states, we discuss the implications disordered states have on the folding process, especially on the folding pathway. Finally, we discuss the potential applications and limitations of this method.

Structural Proteomics Methods to Interrogate the Conformations and Dynamics of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) and regions of intrinsic disorder (IDRs) are abundant in proteomes and are essential for many biological processes. Thus, they are often implicated in disease mechanisms, including neurodegeneration and cancer. The flexible nature of IDPs and IDRs provides many advantages, including (but not limited to) overcoming steric restrictions in binding, facilitating posttranslational modifications, and achieving high binding specificity with low affinity. IDPs adopt a heterogeneous structural ensemble, in contrast to typical folded proteins, making it challenging to interrogate their structure using conventional tools. Structural mass spectrometry (MS) methods are playing an increasingly important role in characterizing the structure and function of IDPs and IDRs, enabled by advances in the design of instrumentation and the development of new workflows, including in native MS, ion mobility MS, top-down MS, hydrogen-deuterium exchange MS, crosslinking MS, and covalent labeling. Here, we describe the advantages of these methods that make them ideal to study IDPs and highlight recent applications where these tools have underpinned new insights into IDP structure and function that would be difficult to elucidate using other methods.

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.

Direct Observation of β-Barrel Intermediates in the Self-Assembly of Toxic SOD128-38 and Absence in Nontoxic Glycine Mutants

Soluble low-molecular-weight oligomers formed during the early stage of amyloid aggregation are considered the major toxic species in amyloidosis. The structure-function relationship between oligomeric assemblies and the cytotoxicity in amyloid diseases are still elusive due to the heterogeneous and transient nature of these aggregation intermediates. To uncover the structural characteristics of toxic oligomeric intermediates, we compared the self-assembly dynamics and structures of SOD1_28-38, a cytotoxic fragment of the superoxide dismutase 1 (SOD1) associated with the amyotrophic lateral sclerosis, with its two nontoxic mutants G33V and G33W using molecular dynamics simulations. Single-point glycine substitutions in SOD1_28-38 have been reported to abolish the amyloid toxicity. Our simulation results showed that the toxic SOD1_28-38 and its nontoxic mutants followed different aggregation pathways featuring distinct aggregation intermediates. Specifically, wild-type SOD1_28-38 initially self-assembled into random-coil-rich oligomers, among which fibrillar aggregates composed of well-defined curved single-layer β-sheets were nucleated via coil-to-sheet conversions and the formation of β-barrels as intermediates. In contrast, the nontoxic G33V/G33W mutants readily assembled into small β-sheet-rich oligomers and then coagulated with each other into cross-β fibrils formed by two-layer β-sheets without forming β-barrels as the intermediates. The direct observation of β-barrel oligomers during the assembly of toxic SOD1_28-38 fragments but not the nontoxic glycine-substitution mutants strongly supports β-barrels as the toxic oligomers in amyloidosis, probably via interactions with the cell membrane and forming amyloid pores. With well-defined structures, the β-barrel might serve as a novel therapeutic target against amyloid-related diseases.

Unraveling the NaCl Concentration Effect on the First Stages of α-Synuclein Aggregation

Intraneuronal aggregation of the intrinsically disordered protein α-synuclein is at the core of Parkinson's disease and related neurodegenerative disorders. Several reports show that the concentration of salts in the medium heavily affects its aggregation rate and fibril morphology, but a characterization of the individual monomeric conformations underlying these effects is still lacking. In this work, we have applied our α-synuclein-optimized coarse-grained molecular dynamics approach to decipher the structural features of the protein monomer under a range of NaCl concentrations (0.0-1.0 M). The results show that key intramolecular contacts between the terminal domains are lost at intermediate concentrations (leading to extended conformations likely to fibrillate), but recovered at high concentrations (leading to compact conformations likely to evolve toward amorphous aggregates). The pattern of direct interactions of the terminal α-synuclein domains with Na⁺ and Cl^- ions plays a key role in explaining this effect. Our results are consistent with a recent study reporting a fibrillation enhancement at moderate NaCl concentrations but an inhibition at higher concentrations. The present work will contribute to improving our understanding of the structural features of monomeric α-synuclein, determining its NaCl-induced fibrillation propensity and the molecular basis of synucleinopathies, necessary for the future development of disease-halting therapies.

αB-Crystallin Chaperone Inhibits Aβ Aggregation by Capping the β-Sheet-Rich Oligomers and Fibrils

Inhibiting the cytotoxicity of amyloid aggregation by endogenous proteins is a promising strategy against degenerative amyloid diseases due to their intrinsically high biocompatibility and low immunogenicity. In this study, we investigated the inhibition mechanism of the structured core region of αB-crystallin (αBC) against Aβ fibrillization using discrete molecular dynamics simulations. Our computational results recapitulated the experimentally observed Aβ binding sites in αBC and suggested that αBC could bind to various Aβ aggregate species during the aggregation process-including monomers, dimers, and likely other high molecular weight oligomers, protofibrils, and fibrils-by capping the exposed β-sheet elongation surfaces. Thus, the nucleation of Aβ oligomers into fibrils and the fibril growth could be inhibited. Mechanistic insights obtained from our systematic computational studies may aid in the development of novel therapeutic strategies to modulate the aggregation of pathological, amyloidogenic protein in degenerative diseases.

Emerging Role of Mass Spectrometry-Based Structural Proteomics in Elucidating Intrinsic Disorder in Proteins

Inherent disorder is an integral part of all proteomes, represented as fully or partially unfolded proteins. The lack of order in intrinsically disordered proteins (IDPs) results in an incredibly flexible, floppy, and heterogeneous ensemble, contrary to the well-structured and unique organization of folded proteins. Despite such unusual demeanor, IDPs are crucial for numerous cellular processes and are increasingly being associated with disease-causing pathologies. These warrant more intensive investigation of this atypical class of protein. Traditional biophysical tools, however, fall short of analyzing IDPs, thus making their structure-function characterization challenging. Mass spectrometry (MS) in recent years has evolved as a valuable tool for elucidating the unusual conformational facets of IDPs. In this review, the features of advanced MS techniques such as Hydrogen-deuterium exchange (HDX)-MS, native MS, limited proteolysis (LiP)-MS, chemical cross-linking (XL)-MS, and Fast photochemical oxidation of proteins (FPOP)-MS are briefly discussed. Recent MS studies on IDPs and the unique advantages/shortfalls associated with the above methods while evaluating structural proteomics of IDPs, are illustrated. Eventually the future scope of the MS methods in further decoding the unexplored landscapes of IDPs is presented.

Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations

Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.

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.

Experimentally-driven protein structure modeling

Revolutions in natural and exact sciences started at the dawn of last century have led to the explosion of theoretical, experimental, and computational approaches to determine structures of molecules, complexes, as well as their rich conformational dynamics. Since different experimental methods produce information that is attributed to specific time and length scales, corresponding computational methods have to be tailored to these scales and experiments. These methods can be then combined and integrated in scales, hence producing a fuller picture of molecular structure and motion from the "puzzle pieces" offered by various experiments. Here, we describe a number of computational approaches to utilize experimental data to glance into structure of proteins and understand their dynamics. We will also discuss the limitations and the resolution of the constraints-based modeling approaches. SIGNIFICANCE: Experimentally-driven computational structure modeling and determination is a rapidly evolving alternative to traditional approaches for molecular structure determination. These new hybrid experimental-computational approaches are proving to be a powerful microscope to glance into the structural features of intrinsically or partially disordered proteins, dynamics of molecules and complexes. In this review, we describe various approaches in the field of experimentally-driven computational structure modeling.

Ligand-induced disorder-to-order transitions characterized by structural proteomics and molecular dynamics simulations

For disordered proteins, ligand binding can be a critical event that changes their structural dynamics. The ability to characterize such changes would facilitate the development of drugs designed to stabilize disordered proteins, whose mis-folding is important for a number of pathologies, including neurodegenerative diseases such as Parkinson's and Alzheimer's diseases. In this study, we used hydrogen/deuterium exchange, differential crosslinking, differential surface modification, and molecular dynamics (MD) simulations to characterize the structural changes in disordered proteins that result from ligand binding. We show here that both an ATP-independent protein chaperone, Spy L32P, and the FK506 binding domain of a prolyl isomerase, FKBP-25 F145A/I223P, are disordered, yet exhibit structures that are distinct from chemically denatured unfolded states in solution, and that they undergo transitions to a more structured state upon ligand binding. These systems may serve as models for the characterization of ligand-induced disorder-to-order transitions in proteins using structural proteomics approaches. SIGNIFICANCE: In this study, we used hydrogen/deuterium exchange, differential crosslinking, differential surface modification, and molecular-dynamics simulations to characterize the structural changes in disordered proteins that result from ligand binding. The protein-ligand systems studied here (the ATP-independent protein chaperone, Spy L32P, and the FK506 binding domain of a prolyl isomerase, FKBP-25 F145A/I223P) may serve as models for understanding ligand-induced disorder-to-order transitions in proteins. Additionally, the structural proteomic techniques demonstrated here are shown to be effective tools for the characterization of disorder-to-order transitions and can be used to facilitate study of other systems in which this class of structural transition can be used for modulating major pathological features of disease, such as the abnormal protein aggregation that occurs with Parkinson's disease and Alzheimer's disease.

Depicting Conformational Ensembles of α-Synuclein by Single Molecule Force Spectroscopy and Native Mass Spectroscopy

Description of heterogeneous molecular ensembles, such as intrinsically disordered proteins, represents a challenge in structural biology and an urgent question posed by biochemistry to interpret many physiologically important, regulatory mechanisms. Single-molecule techniques can provide a unique contribution to this field. This work applies single molecule force spectroscopy to probe conformational properties of α-synuclein in solution and its conformational changes induced by ligand binding. The goal is to compare data from such an approach with those obtained by native mass spectrometry. These two orthogonal, biophysical methods are found to deliver a complex picture, in which monomeric α-synuclein in solution spontaneously populates compact and partially compacted states, which are differently stabilized by binding to aggregation inhibitors, such as dopamine and epigallocatechin-3-gallate. Analyses by circular dichroism and Fourier-transform infrared spectroscopy show that these transitions do not involve formation of secondary structure. This comparative analysis provides support to structural interpretation of charge-state distributions obtained by native mass spectrometry and helps, in turn, defining the conformational components detected by single molecule force spectroscopy.

Insight into the Structure of the "Unstructured" Tau Protein

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 the conformational ensemble of tau protein in solution. The predicted structures were in agreement with surface modification and long-distance crosslinking data. Tau in solution was found as an ensemble of rather compact globular conformations with distinct topology, inter-residue contacts, and a number of transient secondary-structure elements. Regions important for pathological aggregation consistently were found to contain β strands. The determined structures are compatible with the tau protein in solution being a molten globule at near-ground state with persistent residual structural features which we were able to capture by CL-DMD. The predicted structure may facilitate an understanding of the misfolding and oligomerization pathways of the tau protein.

Mechanism of pharmacochaperoning in a mammalian KATP channel revealed by cryo-EM

ATP-sensitive potassium (K_ATP) channels composed of a pore-forming Kir6.2 potassium channel and a regulatory ABC transporter sulfonylurea receptor 1 (SUR1) regulate insulin secretion in pancreatic β-cells to maintain glucose homeostasis. Mutations that impair channel folding or assembly prevent cell surface expression and cause congenital hyperinsulinism. Structurally diverse K_ATP inhibitors are known to act as pharmacochaperones to correct mutant channel expression, but the mechanism is unknown. Here, we compare cryoEM structures of a mammalian K_ATP channel bound to pharmacochaperones glibenclamide, repaglinide, and carbamazepine. We found all three drugs bind within a common pocket in SUR1. Further, we found the N-terminus of Kir6.2 inserted within the central cavity of the SUR1 ABC core, adjacent the drug binding pocket. The findings reveal a common mechanism by which diverse compounds stabilize the Kir6.2 N-terminus within SUR1’s ABC core, allowing it to act as a firm ‘handle’ for the assembly of metastable mutant SUR1-Kir6.2 complexes.