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.     

Homogeneous and heterogeneous tertiary structure ensembles of amyloid-β peptides

The interplay of modern molecular simulation and high-quality nuclear magnetic resonance (NMR) experiments has reached a fruitful stage for quantitative characterization of structural ensembles of disordered peptides. Amyloid-β 1-42 (Aβ42), the primary peptide associated with Alzheimer's disease, and fragments such as Aβ21-30 are both classified as intrinsically disordered peptides (IDPs). We use a variety of NMR observables to validate de novo molecular dynamics simulations in explicit water to characterize the tertiary structure ensemble of Aβ42 and Aβ21-30 from the perspective of their classification as IDPs. Unlike the Aβ21-30 fragment that conforms to expectations of an IDP that is primarily extended, we find that Aβ42 samples conformations reflecting all possible secondary structure categories and spans the range of IDP classifications from collapsed structured states to highly extended conformations, making it an IDP with a far more heterogeneous tertiary ensemble.     

Compaction properties of an intrinsically disordered protein: Sic1 and its kinase-inhibitor domain

IDPs in their unbound state can transiently acquire secondary and tertiary structure. Describing such intrinsic structure is important to understand the transition between free and bound state, leading to supramolecular complexes with physiological interactors. IDP structure is highly dynamic and, therefore, difficult to study by conventional techniques. This work focuses on conformational analysis of the KID fragment of the Sic1 protein, an IDP with a key regulatory role in the cell-cycle of Saccharomyces cerevisiae. FT-IR spectroscopy, ESI-MS, and IM measurements are used to capture dynamic and short-lived conformational states, probing both secondary and tertiary protein structure. The results indicate that the isolated Sic1 KID retains dynamic helical structure and populates collapsed states of different compactness. A metastable, highly compact species is detected. Comparison between the fragment and the full-length protein suggests that chain length is crucial to the stabilization of compact states of this IDP. The two proteins are compared by a length-independent compaction index.     

Conformational properties of intrinsically disordered proteins bound to the surface of silica nanoparticles

Background

Protein-nanoparticle (NP) interactions dictate properties of nanoconjugates relevant to bionanotechnology. Non-covalent adsorption generates a protein corona (PC) formed by an inner and an outer layer, the hard and soft corona (HC, SC). Intrinsically disordered proteins (IDPs) exist in solution as conformational ensembles, whose response to the presence of NPs is not known.

Not Only Expansion: Proline Content and Density Also Induce Disordered Protein Conformation Compaction

Intrinsically disordered proteins (IDPs) adopt a wide array of different conformations that can be constrained by the presence of proline residues, which are frequently found in IDPs. To assess the effects of proline, we designed a series of peptides that differ with respect to the number of prolines in the sequence and their organization. Using high-resolution atomistic molecular dynamics simulations, we found that accounting for whether the proline residues are clustered or isolated contributed significantly to explaining deviations in the experimentally-determined gyration radii of IDPs from the values expected based on the Flory scaling-law. By contrast, total proline content makes smaller contribution to explaining the effect of prolines on IDP conformation. Proline residues exhibit opposing effects depending on their organizational pattern in the IDP sequence. Clustered prolines (i.e., prolines with ≤2 intervening non-proline residues) result in expanded peptide conformations whereas isolated prolines (i.e., prolines with >2 intervening non-proline residues) impose compacted conformations. Clustered prolines were estimated to induce an expansion of ∼20% in IDP dimension (via formation of PPII structural elements) whereas isolated prolines were estimated to induce a compaction of ∼10% in IDP dimension (via the formation of backbone turns). This dual role of prolines provides a mechanism for conformational switching that does not rely on the kinetically much slower isomerization of cis proline to the trans form. Bioinformatic analysis demonstrates high populations of both isolated and clustered prolines and implementing them in coarse-grained molecular dynamics models illustrates that they improve the characterization of the conformational ensembles of IDPs.     

Quantitative characterization of intrinsic disorder in polyglutamine: insights from analysis based on polymer theories

Intrinsically disordered proteins (IDPs) are unfolded under physiological conditions. Here we ask if archetypal IDPs in aqueous milieus are best described as swollen disordered coils in a good solvent or collapsed disordered globules in a poor solvent. To answer this question, we analyzed data from molecular simulations for a 20-residue polyglutamine peptide and concluded, in accord with experimental results, that water is a poor solvent for this system. The relevance of monomeric polyglutamine is twofold: It is an archetypal IDP sequence and its aggregation is associated with nine neurodegenerative diseases. The main advance in this work lies in our ability to make accurate assessments of solvent quality from analysis of simulations for a single, rather than multiple chain lengths. We achieved this through the proper design of simulations and analysis of order parameters that are used to describe conformational equilibria in polymer physics theories. Despite the preference for collapsed structures, we find that polyglutamine is disordered because a heterogeneous ensemble of conformations of equivalent compactness is populated at equilibrium. It is surprising that water is a poor solvent for polar polyglutamine and the question is: why? Our preliminary analysis suggests that intrabackbone interactions provide at least part of the driving force for the collapse of polyglutamine in water. We also show that dynamics for conversion between distinct conformations resemble structural relaxation in disordered, glassy systems, i.e., the energy landscape for monomeric polyglutamine is rugged. We end by discussing generalizations of our methods to quantitative studies of conformational equilibria of other low-complexity IDP sequences.     

DNP-Enhanced MAS NMR: A Tool to Snapshot Conformational Ensembles of α-Synuclein in Different States

Intrinsically disordered proteins dynamically sample a wide conformational space and therefore do not adopt a stable and defined three-dimensional conformation. The structural heterogeneity is related to their proper functioning in physiological processes. Knowledge of the conformational ensemble is crucial for a complete comprehension of this kind of proteins. We here present an approach that utilizes dynamic nuclear polarization-enhanced solid-state NMR spectroscopy of sparsely isotope-labeled proteins in frozen solution to take snapshots of the complete structural ensembles by exploiting the inhomogeneously broadened line-shapes. We investigated the intrinsically disordered protein α-synuclein (α-syn), which plays a key role in the etiology of Parkinson’s disease, in three different physiologically relevant states. For the free monomer in frozen solution we could see that the so-called “random coil conformation” consists of α-helical and β-sheet-like conformations, and that secondary chemical shifts of neighboring amino acids tend to be correlated, indicative of frequent formation of secondary structure elements. Based on these results, we could estimate the number of disordered regions in fibrillar α-syn as well as in α-syn bound to membranes in different protein-to-lipid ratios. Our approach thus provides quantitative information on the propensity to sample transient secondary structures in different functional states. Molecular dynamics simulations rationalize the results.     

Conformational response to charge clustering in synthetic intrinsically disordered proteins

Background

Recent theoretical and computational studies have shown that the charge content and, most importantly, the linear distribution of opposite charges are major determinants of conformational properties of intrinsically disordered proteins (IDPs). Charge segregation in a sequence can be measured through κ, which represents a normalized measure of charge asymmetry. A strong inverse correlation between κ and radius of gyration has been previously demonstrated for two independent sets of permutated IDP sequences.

Contribution of proline to the pre-structuring tendency of transient helical secondary structure elements in intrinsically disordered proteins

Background

IDPs function without relying on three-dimensional structures. No clear rationale for such a behavior is available yet. PreSMos are transient secondary structures observed in the target-free IDPs and serve as the target-binding “active” motifs in IDPs. Prolines are frequently found in the flanking regions of PreSMos. Contribution of prolines to the conformational stability of the helical PreSMos in IDPs is investigated.

Methods

MD simulations are performed for several IDP segments containing a helical PreSMo and the flanking prolines. To measure the influence of flanking-prolines on the structural content of a helical PreSMo calculations were done for wild type as well as for mutant segments with Pro → Asp, His, Lys, or Ala. The change in the helicity due to removal of a proline was measured both for the PreSMo region and for the flanking regions.

Results

The α-helical content in ~ 70% of the helical PreSMos at the early stage of simulation decreases due to replacement of an N-terminal flanking proline by other residues whereas the helix content in nearly all PreSMos increases when the same replacements occur at the C-terminal flanking region. The helix destabilizing/terminating role of the C-terminal flanking prolines is more pronounced than the helix promoting effect of the N-terminal flanking prolines.

General significance

This work represents a novel example demonstrating that a proline is encoded in an IDP with a defined purpose. The helical PreSMos presage their target-bound conformations. As they most likely mediate IDP-target binding via conformational selection their helical content can be an important feature for IDP function.     

CIDER: Resources to Analyze Sequence-Ensemble Relationships of Intrinsically Disordered Proteins

Intrinsically disordered proteins and regions (IDPs) represent a large class of proteins that are defined by conformational heterogeneity and lack of persistent tertiary/secondary structure. IDPs play important roles in a range of biological functions, and their dysregulation is central to numerous diseases, including neurodegeneration and cancer. The conformational ensembles of IDPs are encoded by their amino acid sequences. Here, we present two computational tools that are designed to enable rapid and high-throughput analyses of a wide range of physicochemical properties encoded by IDP sequences. The first, CIDER, is a user-friendly webserver that enables rapid analysis of IDP sequences. The second, localCIDER, is a high-performance software package that enables a wide range of analyses relevant to IDP sequences. In addition to introducing the two packages, we demonstrate the utility of these resources using examples where sequence analysis offers biophysical insights.     

Chain collapse of an amyloidogenic intrinsically disordered protein

Natively unfolded or intrinsically disordered proteins (IDPs) are under intense scrutiny due to their involvement in both normal biological functions and abnormal protein misfolding disorders. Polypeptide chain collapse of amyloidogenic IDPs is believed to play a key role in protein misfolding, oligomerization, and aggregation leading to amyloid fibril formation, which is implicated in a number of human diseases. In this work, we used bovine κ-casein, which serves as an archetypal model protein for amyloidogenic IDPs. Using a variety of biophysical tools involving both prediction and spectroscopic techniques, we first established that monomeric κ-casein adopts a collapsed premolten-globule-like conformational ensemble under physiological conditions. Our time-resolved fluorescence and light-scattering data indicate a change in the mean hydrodynamic radius from ∼4.6 nm to ∼1.9 nm upon chain collapse. We then took the advantage of two cysteines separated by 77 amino-acid residues and covalently labeled them using thiol-reactive pyrene maleimide. This dual-labeled protein demonstrated a strong excimer formation upon renaturation from urea- and acid-denatured states under both equilibrium and kinetic conditions, providing compelling evidence of polypeptide chain collapse under physiological conditions. The implication of the IDP chain collapse in protein aggregation and amyloid formation is also discussed.     

Effects of Molecular Crowding on the Dynamics of Intrinsically Disordered Proteins

Inside cells, the concentration of macromolecules can reach up to 400 g/L. In such crowded environments, proteins are expected to behave differently than in vitro. It has been shown that the stability and the folding rate of a globular protein can be altered by the excluded volume effect produced by a high density of macromolecules. However, macromolecular crowding effects on intrinsically disordered proteins (IDPs) are less explored. These proteins can be extremely dynamic and potentially sample a wide ensemble of conformations under non-denaturing conditions. The dynamic properties of IDPs are intimately related to the timescale of conformational exchange within the ensemble, which govern target recognition and how these proteins function. In this work, we investigated the macromolecular crowding effects on the dynamics of several IDPs by measuring the NMR spin relaxation parameters of three disordered proteins (ProTα, TC1, and α-synuclein) with different extents of residual structures. To aid the interpretation of experimental results, we also performed an MD simulation of ProTα. Based on the MD analysis, a simple model to correlate the observed changes in relaxation rates to the alteration in protein motions under crowding conditions was proposed. Our results show that 1) IDPs remain at least partially disordered despite the presence of high concentration of other macromolecules, 2) the crowded environment has differential effects on the conformational propensity of distinct regions of an IDP, which may lead to selective stabilization of certain target-binding motifs, and 3) the segmental motions of IDPs on the nanosecond timescale are retained under crowded conditions. These findings strongly suggest that IDPs function as dynamic structural ensembles in cellular environments.     

How multisite phosphorylation impacts the conformations of intrinsically disordered proteins

Phosphorylation of intrinsically disordered proteins (IDPs) can produce changes in structural and dynamical properties and thereby mediate critical biological functions. How phosphorylation effects intrinsically disordered proteins has been studied for an increasing number of IDPs, but a systematic understanding is still lacking. Here, we compare the collapse propensity of four disordered proteins, Ash1, the C-terminal domain of RNA polymerase (CTD2’), the cytosolic domain of E-Cadherin, and a fragment of the p130Cas, in unphosphorylated and phosphorylated forms using extensive all-atom molecular dynamics (MD) simulations. We find all proteins to show V-shape changes in their collapse propensity upon multi-site phosphorylation according to their initial net charge: phosphorylation expands neutral or overall negatively charged IDPs and shrinks positively charged IDPs. However, force fields including those tailored towards and commonly used for IDPs overestimate these changes. We find quantitative agreement of MD results with SAXS and NMR data for Ash1 and CTD2’ only when attenuating protein electrostatic interactions by using a higher salt concentration (e.g. 350 mM), highlighting the overstabilization of salt bridges in current force fields. We show that phosphorylation of IDPs also has a strong impact on the solvation of the protein, a factor that in addition to the actual collapse or expansion of the IDP should be considered when analyzing SAXS data. Compared to the overall mild change in global IDP dimension, the exposure of active sites can change significantly upon phosphorylation, underlining the large susceptibility of IDP ensembles to regulation through post-translational modifications.     

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.     

Towards the physical basis of how intrinsic disorder mediates protein function

Intrinsically disordered proteins (IDPs) are an important class of functional proteins that is highly prevalent in biology and has broad association with human diseases. In contrast to structured proteins, free IDPs exist as heterogeneous and dynamical conformational ensembles under physiological conditions. Many concepts have been discussed on how such intrinsic disorder may provide crucial functional advantages, particularly in cellular signaling and regulation. Establishing the physical basis of these proposed phenomena requires not only detailed characterization of the disordered conformational ensembles, but also mechanistic understanding of the roles of various ensemble properties in IDP interaction and regulation. Here, we review the experimental and computational approaches that may be integrated to address many important challenges of establishing a "structural" basis of IDP function, and discuss some of the key emerging ideas on how the conformational ensembles of IDPs may mediate function, especially in coupled binding and folding interactions.     

Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins

Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 μs. Four combinations of protein and water force fields were tested: ff19SB/OPC , ff19SB/TIP4P-D , ff03CMAP/TIP4P-D , and a99SB-disp/TIP4P-disp , with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 μs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and ³J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.     

Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.     

A J-modulated protonless NMR experiment characterizes the conformational ensemble of the intrinsically disordered protein WIP

Intrinsically disordered proteins (IDPs) are multi-conformational polypeptides that lack a single stable three-dimensional structure. It has become increasingly clear that the versatile IDPs play key roles in a multitude of biological processes, and, given their flexible nature, NMR is a leading method to investigate IDP behavior on the molecular level. Here we present an IDP-tailored J-modulated experiment designed to monitor changes in the conformational ensemble characteristic of IDPs by accurately measuring backbone one- and two-bond J(¹⁵N,¹³Cα) couplings. This concept was realized using a unidirectional (H)NCO ¹³C-detected experiment suitable for poor spectral dispersion and optimized for maximum coverage of amino acid types. To demonstrate the utility of this approach we applied it to the disordered actin-binding N-terminal domain of WASp interacting protein (WIP), a ubiquitous key modulator of cytoskeletal changes in a range of biological systems. One- and two-bond J(¹⁵N,¹³Cα) couplings were acquired for WIP residues 2–65 at various temperatures, and in denaturing and crowding environments. Under native conditions fitted J-couplings identified in the WIP conformational ensemble a propensity for extended conformation at residues 16–23 and 45–60, and a helical tendency at residues 28–42. These findings are consistent with a previous study of the based upon chemical shift and RDC data and confirm that the WIP^2–65 conformational ensemble is biased towards the structure assumed by this fragment in its actin-bound form. The effects of environmental changes upon this ensemble were readily apparent in the J-coupling data, which reflected a significant decrease in structural propensity at higher temperatures, in the presence of 8 M urea, and under the influence of a bacterial cell lysate. The latter suggests that crowding can cause protein unfolding through protein–protein interactions that stabilize the unfolded state. We conclude that J-couplings are a useful measureable in characterizing structural ensembles in IDPs, and that the proposed experiment provides a practical method for accurately performing such measurements, once again emphasizing the power of NMR in studying IDP behavior.