Cancer-Associated Mutations Perturb the Disordered Ensemble and Interactions of the Intrinsically Disordered p53 Transactivation Domain

Intrinsically disordered proteins (IDPs) are key components of regulatory networks that control crucial aspects of cell decision making. The intrinsically disordered transactivation domain (TAD) of tumor suppressor p53 mediates its interactions with multiple regulatory pathways to control the p53 homeostasis during the cellular response to genotoxic stress. Many cancer-associated mutations have been discovered in p53-TAD, but their structural and functional consequences are poorly understood. Here, by combining atomistic simulations, NMR spectroscopy, and binding assays, we demonstrate that cancer-associated mutations can significantly perturb the balance of p53 interactions with key activation and degradation regulators. Importantly, the four mutations studied in this work do not all directly disrupt the known interaction interfaces. Instead, at least three of these mutations likely modulate the disordered state of p53-TAD to perturb its interactions with regulators. Specifically, NMR and simulation analysis together suggest that these mutations can modulate the level of conformational expansion as well as rigidity of the disordered state. Our work suggests that the disordered conformational ensemble of p53-TAD can serve as a central conduit in regulating the response to various cellular stimuli at the protein–protein interaction level. Understanding how the disordered state of IDPs may be modulated by regulatory signals and/or disease associated perturbations will be essential in the studies on the role of IDPs in biology and diseases.     

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.     

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.     

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.     

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.     

Criticality in the conformational phase transition among self-similar groups in intrinsically disordered proteins: Probed by salt-bridge dynamics

Intrinsically disordered proteins (IDP) serve as one of the key components in the global proteome. In contrast to globular proteins, they harbor an enormous amount of physical flexibility enforcing them to be retained in conformational ensembles rather than stable folds. Previous studies in an aligned direction have revealed the importance of transient dynamical phenomena like that of salt-bridge formation in IDPs to support their physical flexibility and have further highlighted their functional relevance. For this characteristic flexibility, IDPs remain amenable and accessible to different ordered binding partners, supporting their potential multi-functionality. The current study further addresses this complex structure-functional interplay in IDPs using phase transition dynamics to conceptualize the underlying (avalanche type) mechanism of their being distributed across and hopping around degenerate structural states (conformational ensembles). For this purpose, extensive molecular dynamics simulations have been done and the data analyzed from a statistical physics perspective. Investigation of the plausible scope of 'self-organized criticality' (SOC) to fit into the complex dynamics of IDPs was found to be assertive, relating the conformational degeneracy of these proteins to their functional multiplicity. In accordance with the transient nature of 'salt-bridge dynamics', the study further uses it as a probe to explain the structural basis of the proposed criticality in the conformational phase transition among self-similar groups in IDPs. The analysis reveal scale-invariant self-similar fractal geometries in the structural conformations of different IDPs. The insights from the study has the potential to be extended further to benefit structural tinkering of IDPs in their functional characterization and drugging.     

Discrete Molecular Dynamics Can Predict Helical Prestructured Motifs in Disordered Proteins

Intrinsically disordered proteins (IDPs) lack a stable tertiary structure, but their short binding regions termed Pre-Structured Motifs (PreSMo) can form transient secondary structure elements in solution. Although disordered proteins are crucial in many biological processes and designing strategies to modulate their function is highly important, both experimental and computational tools to describe their conformational ensembles and the initial steps of folding are sparse. Here we report that discrete molecular dynamics (DMD) simulations combined with replica exchange (RX) method efficiently samples the conformational space and detects regions populating α-helical conformational states in disordered protein regions. While the available computational methods predict secondary structural propensities in IDPs based on the observation of protein-protein interactions, our ab initio method rests on physical principles of protein folding and dynamics. We show that RX-DMD predicts α-PreSMos with high confidence confirmed by comparison to experimental NMR data. Moreover, the method also can dissect α-PreSMos in close vicinity to each other and indicate helix stability. Importantly, simulations with disordered regions forming helices in X-ray structures of complexes indicate that a preformed helix is frequently the binding element itself, while in other cases it may have a role in initiating the binding process. Our results indicate that RX-DMD provides a breakthrough in the structural and dynamical characterization of disordered proteins by generating the structural ensembles of IDPs even when experimental data are not available.     

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.     

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.     

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.     

Intrinsically Disordered Protein Exhibits Both Compaction and Expansion under Macromolecular Crowding

Conformational malleability allows intrinsically disordered proteins (IDPs) to respond agilely to their environments, such as nonspecifically interacting with in vivo bystander macromolecules (or crowders). Previous studies have emphasized conformational compaction of IDPs due to steric repulsion by macromolecular crowders, but effects of soft attraction are largely unexplored. Here we studied the conformational ensembles of the IDP FlgM in both polymer and protein crowders by small-angle neutron scattering. As crowder concentrations increased, the mean radius of gyration of FlgM first decreased but then exhibited an uptick. Ensemble optimization modeling indicated that FlgM conformations under protein crowding segregated into two distinct populations, one compacted and one extended. Coarse-grained simulations showed that compacted conformers fit into an interstitial void and occasionally bind to a surrounding crowder, whereas extended conformers snake through interstitial crevices and bind multiple crowders simultaneously. Crowder-induced conformational segregation may facilitate various cellular functions of IDPs.     

Molecular Dynamics Ensemble Refinement of Intrinsically Disordered Peptides According to Deconvoluted Spectra from Circular Dichroism

We have developed a computational method of atomistically refining the structural ensemble of intrinsically disordered peptides (IDPs) facilitated by experimental measurements using circular dichroism spectroscopy (CD). A major challenge surrounding this approach stems from the deconvolution of experimental CD spectra into secondary structure features of the IDP ensemble. Currently available algorithms for CD deconvolution were designed to analyze the spectra of proteins with stable secondary structures. Herein, our work aims to minimize any bias from the peptide deconvolution analysis by implementing a non-negative linear least-squares fitting algorithm in conjunction with a CD reference data set that contains soluble and denatured proteins (SDP48). The non-negative linear least-squares method yields the best results for deconvolution of proteins with higher disordered content than currently available methods, according to a validation analysis of a set of protein spectra with Protein Data Bank entries. We subsequently used this analysis to deconvolute our experimental CD data to refine our computational model of the peptide secondary structure ensemble produced by all-atom molecular dynamics simulations with implicit solvent. We applied this approach to determine the ensemble structures of a set of short IDPs, that mimic the calmodulin binding domain of calcium/calmodulin-dependent protein kinase II and its 1-amino-acid and 3-amino-acid mutants. Our study offers a, to our knowledge, novel way to solve the ensemble secondary structures of IDPs in solution, which is important to advance the understanding of their roles in regulating signaling pathways through the formation of complexes with multiple partners.     

Modulation of Correlated Segment Fluctuations in IDPs upon Complex Formation as an Allosteric Regulatory Mechanism

Molecular recognition of and by intrinsically disordered proteins (IDPs) is an intriguing and still largely elusive phenomenon. Typically, protein recognition involving IDPs requires either folding upon binding or, alternatively, the formation of “fuzzy complexes.” Here we show via correlation analyses of paramagnetic relaxation enhancement data unprecedented and striking alterations of the concerted fluctuations within the conformational ensemble of IDPs upon ligand binding. We study the binding of α-synuclein to calmodulin, a ubiquitous calcium-binding protein, and the binding of the extracellular matrix IDP osteopontin to heparin, a mimic of the extracellular matrix ligand hyaluronic acid. In both cases, binding leads to reduction of correlated long-range motions in these two IDPs and thus indicates a loosening of structural compaction upon binding. Most importantly, however, the simultaneous presence of correlated and anti-correlated fluctuations in IDPs suggests the prevalence of “energetic frustration” and provides an explanation for the puzzling observation of disordered allostery in IDPs.     

NMR illuminates intrinsic disorder

Nuclear magnetic resonance (NMR) has long been instrumental in the characterization of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs). This method continues to offer rich insights into the nature of IDPs in solution, especially in combination with other biophysical methods such as small-angle scattering, single-molecule fluorescence, electron paramagnetic resonance (EPR), and mass spectrometry. Substantial advances have been made in recent years in studies of proteins containing both ordered and disordered domains and in the characterization of problematic sequences containing repeated tracts of a single or a few amino acids. These sequences are relevant to disease states such as Alzheimer's, Parkinson's, and Huntington's diseases, where disordered proteins misfold into harmful amyloid. Innovative applications of NMR are providing novel insights into mechanisms of protein aggregation and the complexity of IDP interactions with their targets. As a basis for understanding the solution structural ensembles, dynamic behavior, and functional mechanisms of IDPs and IDRs, NMR continues to prove invaluable.