Residual Structure Accelerates Binding of Intrinsically Disordered ACTR by Promoting Efficient Folding upon Encounter

Intrinsically disordered proteins (IDPs) often fold into stable structures upon specific binding. The roles of residual structure of unbound IDPs in coupling binding and folding have been under much debate. While many studies emphasize the importance of conformational flexibility for IDP recognition, it was recently demonstrated that stabilization the N-terminal helix of intrinsically disordered ACTR accelerated its binding to another IDP, NCBD of the CREB-binding protein. To understand how enhancing ACTR helicity accelerates binding, we derived a series of topology-based coarse-grained models that mimicked various ACTR mutants with increasing helical contents and reproduced their NCBD binding affinities. Molecular dynamics simulations were then performed to sample hundreds of reversible coupled binding and folding transitions. The results show that increasing ACTR helicity does not alter the baseline mechanism of synergistic folding, which continues to follow “extended conformational selection” with multiple stages of selection and induced folding. Importantly, these coarse-grained models, while only calibrated based on binding thermodynamics, recapitulate the observed kinetic acceleration with increasing ACTR helicity. However, the residual helices do not enhance the association kinetics via more efficient seeding of productive collisions. Instead, they allow the nonspecific collision complexes to evolve more efficiently into the final bound and folded state, which is the primary source of accelerated association kinetics. Meanwhile, reduced dissociation kinetics with increasing ACTR helicity can be directly attributed to smaller entropic cost of forming the bound state. Altogether, this study provides important mechanistic insights into how residual structure may modulate thermodynamics and kinetics of IDP interactions.     

Dynamical coupling of intrinsically disordered proteins and their hydration water: comparison with folded soluble and membrane proteins

Hydration water is vital for various macromolecular biological activities, such as specific ligand recognition, enzyme activity, response to receptor binding, and energy transduction. Without hydration water, proteins would not fold correctly and would lack the conformational flexibility that animates their three-dimensional structures. Motions in globular, soluble proteins are thought to be governed to a certain extent by hydration-water dynamics, yet it is not known whether this relationship holds true for other protein classes in general and whether, in turn, the structural nature of a protein also influences water motions. Here, we provide insight into the coupling between hydration-water dynamics and atomic motions in intrinsically disordered proteins (IDP), a largely unexplored class of proteins that, in contrast to folded proteins, lack a well-defined three-dimensional structure. We investigated the human IDP tau, which is involved in the pathogenic processes accompanying Alzheimer disease. Combining neutron scattering and protein perdeuteration, we found similar atomic mean-square displacements over a large temperature range for the tau protein and its hydration water, indicating intimate coupling between them. This is in contrast to the behavior of folded proteins of similar molecular weight, such as the globular, soluble maltose-binding protein and the membrane protein bacteriorhodopsin, which display moderate to weak coupling, respectively. The extracted mean square displacements also reveal a greater motional flexibility of IDP compared with globular, folded proteins and more restricted water motions on the IDP surface. The results provide evidence that protein and hydration-water motions mutually affect and shape each other, and that there is a gradient of coupling across different protein classes that may play a functional role in macromolecular activity in a cellular context.     

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.     

Multi-scaled explorations of binding-induced folding of intrinsically disordered protein inhibitor IA3 to its target enzyme

Biomolecular function is realized by recognition, and increasing evidence shows that recognition is determined not only by structure but also by flexibility and dynamics. We explored a biomolecular recognition process that involves a major conformational change - protein folding. In particular, we explore the binding-induced folding of IA3, an intrinsically disordered protein that blocks the active site cleft of the yeast aspartic proteinase saccharopepsin (YPrA) by folding its own N-terminal residues into an amphipathic alpha helix. We developed a multi-scaled approach that explores the underlying mechanism by combining structure-based molecular dynamics simulations at the residue level with a stochastic path method at the atomic level. Both the free energy profile and the associated kinetic paths reveal a common scheme whereby IA3 binds to its target enzyme prior to folding itself into a helix. This theoretical result is consistent with recent time-resolved experiments. Furthermore, exploration of the detailed trajectories reveals the important roles of non-native interactions in the initial binding that occurs prior to IA3 folding. In contrast to the common view that non-native interactions contribute only to the roughness of landscapes and impede binding, the non-native interactions here facilitate binding by reducing significantly the entropic search space in the landscape. The information gained from multi-scaled simulations of the folding of this intrinsically disordered protein in the presence of its binding target may prove useful in the design of novel inhibitors of aspartic proteinases.     

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.     

Protein domain definition should allow for conditional disorder

Abstract: Proteins are often classified in a binary fashion as either structured or disordered. However this approach has several deficits. Firstly, protein folding is always conditional on the physiochemical environment. A protein which is structured in some circumstances will be disordered in others. Secondly, it hides a fundamental asymmetry in behavior. While all structured proteins can be unfolded through a change in environment, not all disordered proteins have the capacity for folding. Failure to accommodate these complexities confuses the definition of both protein structural domains and intrinsically disordered regions. We illustrate these points with an experimental study of a family of small binding domains, drawn from the RNA polymerase of mumps virus and its closest relatives. Assessed at face value the domains fall on a structural continuum, with folded, partially folded, and near unstructured members. Yet the disorder present in the family is conditional, and these closely related polypeptides can access the same folded state under appropriate conditions. Any heuristic definition of the protein domain emphasizing conformational stability divides this domain family in two, in a way that makes no biological sense. Structural domains would be better defined by their ability to adopt a specific tertiary structure: a structure that may or may not be realized, dependent on the circumstances. This explicitly allows for the conditional nature of protein folding, and more clearly demarcates structural domains from intrinsically disordered regions that may function without folding.     

Atomic-level characterization of disordered protein ensembles

The roles of unfolded states of proteins in normal folding and in diseases involving aggregation, as well as the prevalence and regulatory functions of intrinsically disordered proteins, have become increasingly recognized. The structural representation of these disordered states as ensembles of interconverting conformers can therefore provide critical insights. Experimental methods can be used to probe ensemble-averaged structural properties of disordered states and computational approaches generate representative ensembles of conformers using experimental restraints. In particular, NMR and small-angle X-ray scattering provide quantitative data that can readily be incorporated into calculations. These techniques have gleaned structural information about denatured, unfolded and intrinsically disordered proteins. The use of experimental data in different computational approaches, including ensemble molecular dynamics simulations and algorithms that assign populations to pregenerated conformers, has highlighted the presence of both local and long-range structure, and the occurrence of native-like and non-native interactions in unfolded and denatured states. Analysis of the resulting ensembles has suggested important implications of this fluctuating structure for folding, aggregation and binding.     

Theory of protein folding

Protein folding should be complex. Proteins organize themselves into specific three-dimensional structures, through a myriad of conformational changes. The classical view of protein folding describes this process as a nearly sequential series of discrete intermediates. In contrast, the energy landscape theory of folding considers folding as the progressive organization of an ensemble of partially folded structures through which the protein passes on its way to the natively folded structure. As a result of evolution, proteins have a rugged funnel-like landscape biased toward the native structure. Connecting theory and simulations of minimalist models with experiments has completely revolutionized our understanding of the underlying mechanisms that control protein folding.     

Leucine-rich hydrophobic clusters promote folding of the N-terminus of the intrinsically disordered transactivation domain of p53

Molecular dynamics simulations have been performed on the intrinsically disordered 39-residue N-terminal transactivation domain of p53 (p53_1-39). Simulations not only revealed that p53_1-39 is natively compact, but also possesses a folded structure. Furthermore, leucine-rich hydrophobic clusters were found to play a crucial role in the formation and stabilization of the folded structure of p53_1-39. Collapsing in the sub-microsecond timescale might allow for rapid conformational turnovers of p53_1-39, necessary for its efficient transactivation activity and modulation. Fast collapsing might be the result of unique conformational landscapes, featuring several energy minima separated by small energy barriers. It is suggested that IDPs with highly specialized functions in the cell, such as transactivation, possibly display more ordered patterns than their less specialized counterparts.     

Topology-based modeling of intrinsically disordered proteins: balancing intrinsic folding and intermolecular interactions

Coupled binding and folding is frequently involved in specific recognition of so-called intrinsically disordered proteins (IDPs), a newly recognized class of proteins that rely on a lack of stable tertiary fold for function. Here, we exploit topology-based Gō-like modeling as an effective tool for the mechanism of IDP recognition within the theoretical framework of minimally frustrated energy landscape. Importantly, substantial differences exist between IDPs and globular proteins in both amino acid sequence and binding interface characteristics. We demonstrate that established Gō-like models designed for folded proteins tend to over-estimate the level of residual structures in unbound IDPs, whereas under-estimating the strength of intermolecular interactions. Such systematic biases have important consequences in the predicted mechanism of interaction. A strategy is proposed to recalibrate topology-derived models to balance intrinsic folding propensities and intermolecular interactions, based on experimental knowledge of the overall residual structure level and binding affinity. Applied to pKID/KIX, the calibrated Gō-like model predicts a dominant multistep sequential pathway for binding-induced folding of pKID that is initiated by KIX binding via the C-terminus in disordered conformations, followed by binding and folding of the rest of C-terminal helix and finally the N-terminal helix. This novel mechanism is consistent with key observations derived from a recent NMR titration and relaxation dispersion study and provides a molecular-level interpretation of kinetic rates derived from dispersion curve analysis. These case studies provide important insight into the applicability and potential pitfalls of topology-based modeling for studying IDP folding and interaction in general.     

Conformational characteristics of unstructured peptides: alpha-synuclein

We have performed replica-exchange molecular dynamics simulations on 41 residue peptides containing NAC region of alpha-synuclein in various force fields and solvent conditions. Alpha-synuclein is known to be the major cause of Parkinson's disease by amyloid-like aggregation, and one of the natively unfolded proteins. To investigate conformational characteristics of intrinsically unstructured peptides, we carried out structural analysis by introducing 'representative structure' for ensemble of structures occurring during the overall trajectory. Representative structures may be defined by using either coordinate averaging or distance averaging. When applied to the natively folded proteins such as villin headpiece, structural analysis based on representative structure was found to yield consistent results with those obtained from conventional analysis. Individual conformations obtained from the simulations of NAC peptide for various conditions show flexible structures close to random coil. Secondary structure contents and free energy surfaces showed dependency on solvent conditions, which may be interpreted as another manifestation of structural diversity. It is found that representative structures can provide useful information about structural characteristics of intrinsically unstructured proteins.     

Electrostatically Accelerated Encounter and Folding for Facile Recognition of Intrinsically Disordered Proteins

Achieving facile specific recognition is essential for intrinsically disordered proteins (IDPs) that are involved in cellular signaling and regulation. Consideration of the physical time scales of protein folding and diffusion-limited protein-protein encounter has suggested that the frequent requirement of protein folding for specific IDP recognition could lead to kinetic bottlenecks. How IDPs overcome such potential kinetic bottlenecks to viably function in signaling and regulation in general is poorly understood. Our recent computational and experimental study of cell-cycle regulator p27 (Ganguly et al., J. Mol. Biol. (2012)) demonstrated that long-range electrostatic forces exerted on enriched charges of IDPs could accelerate protein-protein encounter via “electrostatic steering” and at the same time promote “folding-competent” encounter topologies to enhance the efficiency of IDP folding upon encounter. Here, we further investigated the coupled binding and folding mechanisms and the roles of electrostatic forces in the formation of three IDP complexes with more complex folded topologies. The surface electrostatic potentials of these complexes lack prominent features like those observed for the p27/Cdk2/cyclin A complex to directly suggest the ability of electrostatic forces to facilitate folding upon encounter. Nonetheless, similar electrostatically accelerated encounter and folding mechanisms were consistently predicted for all three complexes using topology-based coarse-grained simulations. Together with our previous analysis of charge distributions in known IDP complexes, our results support a prevalent role of electrostatic interactions in promoting efficient coupled binding and folding for facile specific recognition. These results also suggest that there is likely a co-evolution of IDP folded topology, charge characteristics, and coupled binding and folding mechanisms, driven at least partially by the need to achieve fast association kinetics for cellular signaling and regulation.     

AlphaFold and Implications for Intrinsically Disordered Proteins

Accurate predictions of the three-dimensional structures of proteins from their amino acid sequences have come of age. AlphaFold, a deep learning-based approach to protein structure prediction, shows remarkable success in independent assessments of prediction accuracy. A significant epoch in structural bioinformatics was the structural annotation of over 98% of protein sequences in the human proteome. Interestingly, many predictions feature regions of very low confidence, and these regions largely overlap with intrinsically disordered regions (IDRs). That over 30% of regions within the proteome are disordered is congruent with estimates that have been made over the past two decades, as intense efforts have been undertaken to generalize the structure–function paradigm to include the importance of conformational heterogeneity and dynamics. With structural annotations from AlphaFold in hand, there is the temptation to draw inferences regarding the “structures” of IDRs and their interactomes. Here, we offer a cautionary note regarding the misinterpretations that might ensue and highlight efforts that provide concrete understanding of sequence-ensemble-function relationships of IDRs. This perspective is intended to emphasize the importance of IDRs in sequence-function relationships (SERs) and to highlight how one might go about extracting quantitative SERs to make sense of how IDRs function.     

The structural context of posttranslational modifications at a proteome-wide scale

The recent revolution in computational protein structure prediction provides folding models for entire proteomes, which can now be integrated with large-scale experimental data. Mass spectrometry (MS)-based proteomics has identified and quantified tens of thousands of posttranslational modifications (PTMs), most of them of uncertain functional relevance. In this study, we determine the structural context of these PTMs and investigate how this information can be leveraged to pinpoint potential regulatory sites. Our analysis uncovers global patterns of PTM occurrence across folded and intrinsically disordered regions. We found that this information can help to distinguish regulatory PTMs from those marking improperly folded proteins. Interestingly, the human proteome contains thousands of proteins that have large folded domains linked by short, disordered regions that are strongly enriched in regulatory phosphosites. These include well-known kinase activation loops that induce protein conformational changes upon phosphorylation. This regulatory mechanism appears to be widespread in kinases but also occurs in other protein families such as solute carriers. It is not limited to phosphorylation but includes ubiquitination and acetylation sites as well. Furthermore, we performed three-dimensional proximity analysis, which revealed examples of spatial coregulation of different PTM types and potential PTM crosstalk. To enable the community to build upon these first analyses, we provide tools for 3D visualization of proteomics data and PTMs as well as python libraries for data accession and processing.

A combination of the comprehensive structural predictions of AlphaFold2 and large-scale proteomics data on post-translational modifications (PTMs) reveals novel insights into the functional importance of PTMs, based on their structural context.     

Studying backbone torsional dynamics of intrinsically disordered proteins using fluorescence depolarization kinetics

Intrinsically disordered proteins (IDPs) do not autonomously adopt a stable unique 3D structure and exist as an ensemble of rapidly interconverting structures. They are characterized by significant conformational plasticity and are associated with several biological functions and dysfunctions. The rapid conformational fluctuation is governed by the backbone segmental dynamics arising due to the dihedral angle fluctuation on the Ramachandran ?–ψ conformational space. We discovered that the intrinsic backbone torsional mobility can be monitored by a sensitive fluorescence readout, namely fluorescence depolarization kinetics, of tryptophan in an archetypal IDP such as α-synuclein. This methodology allows us to map the site-specific torsional mobility in the dihedral space within picosecond-nanosecond time range at a low protein concentration under the native condition. The characteristic timescale of ~?1.4 ns, independent of residue position, represents collective torsional dynamics of dihedral angles (? and ψ) of several residues from tryptophan and is independent of overall global tumbling of the protein. We believe that fluorescence depolarization kinetics methodology will find broad application to study both short-range and long-range correlated motions, internal friction, binding-induced folding, disorder-to-order transition, misfolding and aggregation of IDPs.     

Compromise in competition between free energy and binding effect of intrinsically disordered protein p53 C-terminal domain

The C-terminal domain (CTD) of tumour suppressor p53 is an intrinsically disordered protein which has been shown to be able to bind multiple partner proteins and exercise diverse physiological functions in the cell. In this study, we performed molecular dynamics simulations on the isolated p53 CTD, as well as three regulatory binding complexes to investigate the conformational ensemble of isolated p53 CTD and its dynamic structures when different binding partner present. The results demonstrate that the isolated p53 CTD resembles a molten globule rather than extended structure. It mainly adopts random coil conformations with some tendency to form helical structures, which is consistent with experimental observations. For isolated p53 CTD, the dynamics is exclusively dominated by the intrinsic free energy and the p53 CTD could not folded spontaneously to each binding competent state which is located in high free energy region. However, when the binding partners present, the dynamics of p53 CTD are dominated by two mechanisms, the p53 CTD tending to adopt the structure with minimum free energy as isolate existed and the binding energy from partner protein tending to minimum. Each of them has an extreme tendency and corresponds to a possible characteristic state, the random coil state and each binding competent state. The compromise in competition between these two mechanisms results in alternate realisation of different characteristic states, while the relative strength of each mechanism determines the sampling frequency of each characteristic state.