Residual structures, conformational fluctuations, and electrostatic interactions in the synergistic folding of two intrinsically disordered proteins

Coupled binding and folding of disordered SPIN N-terminal region in myeloperoxidase inhibition

Gram-positive pathogenic bacteria Staphylococcus express and secret staphylococcal peroxidase inhibitor (SPIN) proteins to help evade neutrophil-mediated immunity by inhibiting the activity of the main oxidative-defense player myeloperoxidase (MPO) enzyme. SPIN contains a structured 3-helix bundle C-terminal domain, which can specifically bind to MPO with high affinity, and an intrinsically disordered N-terminal domain (NTD), which folds into a structured β-hairpin and inserts itself into the active site of MPO for inhibition. Mechanistic insights of the coupled folding and binding process are needed in order to better understand how residual structures and/or conformational flexibility of NTD contribute to the different strengths of inhibition of SPIN homologs. In this work, we applied atomistic molecular dynamics simulations on two SPIN homologs, from S. aureus and S. delphini, respectively, which share high sequence identity and similarity, to explore the possible mechanistic basis for their different inhibition efficacies on human MPO. Direct simulations of the unfolding and unbinding processes at 450 K reveal that these two SPIN/MPO complexes systems follow surprisingly different mechanisms of coupled binding and folding. While coupled binding and folding of SPIN-aureus NTD is highly cooperative, SPIN-delphini NTD appears to mainly utilize a conformational selection-like mechanism. These observations are in contrast to an overwhelming prevalence of induced folding-like mechanisms for intrinsically disordered proteins that fold into helical structures upon binding. Further simulations of unbound SPIN NTDs at room temperature reveal that SPIN-delphini NTD has a much stronger propensity of forming β-hairpin like structures, consistent with its preference to fold and then bind. These may help explain why the inhibition strength is not well correlated with binding affinity for different SPIN homologs. Altogether, our work establishes the relationship between the residual conformational stability of SPIN-NTD and their inhibitory function, which can help us develop new strategies towards treating Staphylococcal infections.

A modular approach to map out the conformational landscapes of unbound intrinsically disordered proteins

SignificanceIntrinsically disordered proteins have the unique ability to morph in response to multiple partners and thereby process sophisticated inputs and outputs. It is, however, a mystery whether their response is passive, that is, entirely determined by the partner, or controlled via an internal, yet unknown, folding mechanism. Here we introduce an approach to examine this key question and demonstrate its potential by dissecting the conformational properties of the partially disordered protein NCBD and obtaining important clues about how it performs its biological function.

NMR Provides Unique Insight into the Functional Dynamics and Interactions of Intrinsically Disordered Proteins

Intrinsically disordered proteins are ubiquitous throughout all known proteomes, playing essential roles in all aspects of cellular and extracellular biochemistry. To understand their function, it is necessary to determine their structural and dynamic behavior and to describe the physical chemistry of their interaction trajectories. Nuclear magnetic resonance is perfectly adapted to this task, providing ensemble averaged structural and dynamic parameters that report on each assigned resonance in the molecule, unveiling otherwise inaccessible insight into the reaction kinetics and thermodynamics that are essential for function. In this review, we describe recent applications of NMR-based approaches to understanding the conformational energy landscape, the nature and time scales of local and long-range dynamics and how they depend on the environment, even in the cell. Finally, we illustrate the ability of NMR to uncover the mechanistic basis of functional disordered molecular assemblies that are important for human health.

Topological frustration leading to backtracking in a coupled folding-binding process

Intrinsically disordered proteins (IDPs) are abundant in all species. Their discovery challenges the traditional "sequence-structure-function" paradigm of protein science because IDPs play important roles in various biological processes without preformed folded structures. Bioinformatic analysis reveals that the intrinsically conformational disorder of IDPs as well as their conformational transition upon binding to their targets is encoded by their amino acid sequences. The rRNase domain of colicin E3 and the immunity protein Im3 are a pair of proteins involved in bacterial survival. While the N-terminal segment and the central segment of E3 make comparable intermolecular contacts with Im3 in the bound state, binding of E3 with Im3 is dominantly triggered by the central segment of E3. In this work, to further investigate the binding mechanism of disordered E3 with Im3, we performed systematic free energy and transition path analysis through coarse-grained molecular dynamics simulations. We observed backtracking of the N-terminal segment of E3 in the binding process, whose occurrence depends on salt concentration. Conformational analysis revealed that initial binding of the N-terminal segment of E3 to Im3 usually leads to misorientation of a central hairpin of E3 on Im3, which generates topological frustration and results in backtracking of the N-terminal segment. Our results not only provide deeper mechanistic insights into the coupled folding-binding process of the E3/Im3 complex, but also suggest that topological frustration could be present in the coupled folding-binding process of IDPs and play an important role in regulating the binding transition pathways.

Combining enhanced sampling and deep learning dimensionality reduction for the study of the heat shock protein B8 and its pathological mutant K141E

The biological functions of proteins closely depend on their conformational dynamics. This aspect is especially relevant for intrinsically disordered proteins (IDP) for which structural ensembles often offer more useful representations than individual conformations. Here we employ extensive enhanced sampling temperature replica-exchange atomistic simulations (TREMD) and deep learning dimensionality reduction to study the conformational ensembles of the human heat shock protein B8 and its pathological mutant K141E, for which no experimental 3D structures are available. First, we combined homology modelling with TREMD to generate high-dimensional data sets of 3D structures. Then, we employed a recently developed machine learning based post-processing algorithm, EncoderMap, to project the large conformational data sets into meaningful two-dimensional maps that helped us interpret the data and extract the most significant conformations adopted by both proteins during TREMD. These studies provide the first 3D structural characterization of HSPB8 and reveal the effects of the pathogenic K141E mutation on its conformational ensembles. In particular, this missense mutation appears to increase the compactness of the protein and its structural variability, at the same time rearranging the hydrophobic patches exposed on the protein surface. These results offer the possibility of rationalizing the pathogenic effects of the K141E mutation in terms of conformational changes.

The study provides the first 3D structural characterization of HSPB8 and its K141E mutant: extensive TREMD are combined with a deep learning algorithm to rationalize the disordered ensemble of structures adopted by each variant.

Advanced Sampling Methods for Multiscale Simulation of Disordered Proteins and Dynamic Interactions

Intrinsically disordered proteins (IDPs) are highly prevalent and play important roles in biology and human diseases. It is now also recognized that many IDPs remain dynamic even in specific complexes and functional assemblies. Computer simulations are essential for deriving a molecular description of the disordered protein ensembles and dynamic interactions for a mechanistic understanding of IDPs in biology, diseases, and therapeutics. Here, we provide an in-depth review of recent advances in the multi-scale simulation of disordered protein states, with a particular emphasis on the development and application of advanced sampling techniques for studying IDPs. These techniques are critical for adequate sampling of the manifold functionally relevant conformational spaces of IDPs. Together with dramatically improved protein force fields, these advanced simulation approaches have achieved substantial success and demonstrated significant promise towards the quantitative and predictive modeling of IDPs and their dynamic interactions. We will also discuss important challenges remaining in the atomistic simulation of larger systems and how various coarse-grained approaches may help to bridge the remaining gaps in the accessible time- and length-scales of IDP simulations.

Progressive Phosphorylation Modulates the Self-Association of a Variably Modified Histone H3 Peptide

Protein phosphorylation is a key regulatory mechanism in eukaryotic cells. In the intrinsically disordered histone tails, phosphorylation is often a part of combinatorial post-translational modifications and an integral part of the “histone code” that regulates gene expression. Here, we study the association between two histone H3 tail peptides modified to different degrees, using fully atomistic molecular dynamics simulations. Assuming that the initial conformations are either α-helical or fully extended, we compare the propensity of the two peptides to associate with one another when both are unmodified, one modified and the other unmodified, or both modified. The simulations lead to the identification of distinct inter- and intramolecular interactions in the peptide dimer, highlighting a prominent role of a fine-tuned phosphorylation rheostat in peptide association. Progressive phosphorylation appears to modulate peptide charge, inducing strong and specific intermolecular interactions between the monomers, which do not result in the formation of amorphous or ordered aggregates, as documented by experimental evidence derived from Circular Dichroism and NMR spectroscopy. However, upon complete saturation of positive charges by phosphate groups, this effect is reversed: intramolecular interactions prevail and dimerization of zero-charge peptides is markedly reduced. These findings underscore the role of phosphorylation thresholds in the dynamics of intrinsically disordered proteins. Phosphorylation rheostats might account for the divergent effects of histone modifications on the modulation of chromatin structure.

Accelerating atomistic simulations of proteins using multiscale enhanced sampling with independent tempering

Efficient sampling of the conformational space is essential for quantitative simulations of proteins. The multiscale enhanced sampling (MSES) method accelerates atomistic sampling by coupling it to a coarse-grained (CG) simulation. Bias from coupling to the CG model is removed using Hamiltonian replica exchange, such that one could benefit simultaneously from the high accuracy of atomistic models and fast dynamics of CG ones. Here, we extend MSES to allow independent control of the effective temperatures of atomistic and CG simulations, by directly scaling the atomistic and CG Hamiltonians. The new algorithm, named MSES with independent tempering (MSES-IT), supports more sophisticated Hamiltonian and temperature replica exchange protocols to further improve the sampling efficiency. Using a small but nontrivial β-hairpin, we show that setting the effective temperature of CG model in all conditions to its melting temperature maximizes structural transition rates at the CG level and promotes more efficient replica exchange and diffusion in the condition space. As the result, MSES-IT drive faster reversible transitions at the atomic level and leads to significant improvement in generating converged conformational ensembles compared to the original MSES scheme.

Recent advances in atomic molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) play important roles in cellular functions. The inherent structural heterogeneity of IDPs makes the high-resolution experimental characterization of IDPs extremely difficult. Molecular dynamics (MD) simulation could provide the atomic-level description of the structural and dynamic properties of IDPs. This perspective reviews the recent progress in atomic MD simulation studies of IDPs, including the development of force fields and sampling methods, as well as applications in IDP-involved protein-protein interactions. The employment of large-scale simulations and advanced sampling techniques allows more accurate estimation of the thermodynamics and kinetics of IDP-mediated protein interactions, and the holistic landscape of the binding process of IDPs is emerging.

Molecular simulations of IDPs: From ensemble generation to IDP interactions leading to disorder-to-order transitions

Intrinsically disordered proteins (IDPs) lack a well-defined three-dimensional structure but do exhibit some dynamical and structural ordering. The structural plasticity of IDPs indicates that entropy-driven motions are crucial for their function. Many IDPs undergo function-related disorder-to-order transitions upon by their interaction with specific binding partners. Approaches that are based on both experimental and theoretical tools enable the biophysical characterization of IDPs. Molecular simulations provide insights into IDP structural ensembles and disorder-to-order transition mechanisms. However, such studies depend strongly on the chosen force field parameters and simulation techniques. In this chapter, we provide an overview of IDP characteristics, review all-atom force fields recently developed for IDPs, and present molecular dynamics-based simulation methods that allow IDP ensemble generation as well as the characterization of disorder-to-order transitions. In particular, we introduce metadynamics, replica exchange molecular dynamics simulations, and also kinetic models resulting from Markov State modeling, and provide various examples for the successful application of these simulation methods to IDPs.

A disordered encounter complex is central to the yeast Abp1p SH3 domain binding pathway

Protein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general.

Identification of an α-MoRF in the Intrinsically Disordered Region of the Escargot Transcription Factor

Molecular recognition features (MoRFs) are common in intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs). MoRFs are in constant order-disorder structural transitions and adopt well-defined structures once they are bound to their targets. Here, we study Escargot (Esg), a transcription factor in Drosophila melanogaster that regulates multiple cellular functions, and consists of a disordered N-terminal domain and a group of zinc fingers at its C-terminal domain. We analyzed the N-terminal domain of Esg with disorder predictors and identified a region of 45 amino acids with high probability to form ordered structures, which we named S2. Through 54 μs of molecular dynamics (MD) simulations using CHARMM36 and implicit solvent (generalized Born/surface area (GBSA)), we characterized the conformational landscape of S2 and found an α-MoRF of ∼16 amino acids stabilized by key contacts within the helix. To test the importance of these contacts in the stability of the α-MoRF, we evaluated the effect of point mutations that would impair these interactions, running 24 μs of MD for each mutation. The mutations had mild effects on the MoRF, and in some cases, led to gain of residual structure through long-range contacts of the α-MoRF and the rest of the S2 region. As this could be an effect of the force field and solvent model we used, we benchmarked our simulation protocol by carrying out 32 μs of MD for the (AAQAA)₃ peptide. The results of the benchmark indicate that the global amount of helix in shorter peptides like (AAQAA)₃ is reasonably predicted. Careful analysis of the runs of S2 and its mutants suggests that the mutation to hydrophobic residues may have nucleated long-range hydrophobic and aromatic interactions that stabilize the MoRF. Finally, we have identified a set of residues that stabilize an α-MoRF in a region still without functional annotations in Esg.

Targeting Intrinsically Disordered Proteins through Dynamic Interactions

Intrinsically disordered proteins (IDPs) are over-represented in major disease pathways and have attracted significant interest in understanding if and how they may be targeted using small molecules for therapeutic purposes. While most existing studies have focused on extending the traditional structure-centric drug design strategies and emphasized exploring pre-existing structure features of IDPs for specific binding, several examples have also emerged to suggest that small molecules could achieve specificity in binding IDPs and affect their function through dynamic and transient interactions. These dynamic interactions can modulate the disordered conformational ensemble and often lead to modest compaction to shield functionally important interaction sites. Much work remains to be done on further elucidation of the molecular basis of the dynamic small molecule-IDP interaction and determining how it can be exploited for targeting IDPs in practice. These efforts will rely critically on an integrated experimental and computational framework for disordered protein ensemble characterization. In particular, exciting advances have been made in recent years in enhanced sampling techniques, Graphic Processing Unit (GPU)-computing, and protein force field optimization, which have now allowed rigorous physics-based atomistic simulations to generate reliable structure ensembles for nontrivial IDPs of modest sizes. Such de novo atomistic simulations will play crucial roles in exploring the exciting opportunity of targeting IDPs through dynamic interactions.

Coupled Binding and Helix Formation Monitored by Synchrotron-Radiation Circular Dichroism

Intrinsically disordered proteins organize interaction networks in the cell in many regulation and signaling processes. These proteins often gain structure upon binding to their target proteins in multistep reactions involving the formation of both secondary and tertiary structure. To understand the interactions of disordered proteins, we need to understand the mechanisms of these coupled folding and binding reactions. We studied helix formation in the binding of the molten globule-like nuclear coactivator binding domain and the disordered interaction domain from activator of thyroid hormone and retinoid receptors. We demonstrate that helix formation in a rapid binding reaction can be followed by stopped-flow synchrotron-radiation circular dichroism (CD) spectroscopy and describe the design of such a beamline. Fluorescence-monitored binding experiments of activator of thyroid hormone and retinoid receptors and nuclear coactivator binding domain display several kinetic phases, including one concentration-independent phase, which is consistent with an intermediate stabilized at high ionic strength. Time-resolved CD experiments show that almost all helicity is formed upon initial association of the proteins or separated from the encounter complex by only a small energy barrier. Through simulation of mechanistic models, we show that the intermediate observed at high ionic strength likely involves a structural rearrangement with minor overall changes in helicity. Our experiments provide a benchmark for simulations of coupled binding reactions and demonstrate the feasibility of using synchrotron-radiation CD for mechanistic studies of protein-protein interactions.

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.

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.

A structurally heterogeneous transition state underlies coupled binding and folding of disordered proteins

Many intrinsically disordered proteins (IDPs) attain a well-defined structure in a coupled folding and binding reaction with another protein. Such reactions may involve early to late formation of different native structural regions along the reaction pathway. To obtain insights into the transition state for a coupled binding and folding reaction, we performed restrained molecular dynamics simulations using previously determined experimental binding Φ_b values of the interaction between two IDP domains: the activation domain from the p160 transcriptional co-activator for thyroid hormone and retinoid receptors (ACTR) and the nuclear co-activator binding domain (NCBD) of CREB-binding protein, each forming three well-defined α-helices upon binding. These simulations revealed that both proteins are largely disordered in the transition state for complex formation, except for two helices, one from each domain, that display a native-like structure. The overall transition state structure was extended and largely dynamic with many weakly populated contacts. To test the transition state model, we combined site-directed mutagenesis with kinetic experiments, yielding results consistent with overall diffuse interactions and formation of native intramolecular interactions in the third NCBD helix during the binding reaction. Our findings support the view that the transition state and, by inference, any encounter complex in coupled binding and folding reactions are structurally heterogeneous and largely independent of specific interactions. Furthermore, experimental Φ_b values and Brønsted plots suggested that the transition state is globally robust with respect to most mutations but can display more native-like features for some highly destabilizing mutations, possibly because of Hammond behavior or ground-state effects.

Atomistic Peptide Folding Simulations Reveal Interplay of Entropy and Long-Range Interactions in Folding Cooperativity

Understanding how proteins fold has remained a problem of great interest in biophysical research. Atomistic computer simulations using physics-based force fields can provide important insights on the interplay of different interactions and energetics and their roles in governing the folding thermodynamics and mechanism. In particular, generalized Born (GB)-based implicit solvent force fields can be optimized to provide an appropriate balance between solvation and intramolecular interactions and successfully recapitulate experimental conformational equilibria for a set of helical and β-hairpin peptides. Here, we further demonstrate that key thermodynamic properties and their temperature dependence obtained from replica exchange molecular dynamics simulations of these peptides are in quantitative agreement with experimental results. Useful lessons can be learned on how the interplay of entropy and sequentially long-range interactions governs the mechanism and cooperativity of folding. These results highlight the great potential of high-quality implicit solvent force fields for studying protein folding and large-scale conformational transitions.

HyRes: a coarse-grained model for multi-scale enhanced sampling of disordered protein conformations

Efficient coarse-grained (CG) models can be coupled with atomistic force fields to accelerate the sampling of atomistic energy landscapes in the multi-scale enhanced sampling (MSES) framework. This approach may be particularly suitable for generating atomistic conformational ensembles of intrinsically disordered proteins (IDPs). While MSES is relatively robust to inherent CG artifacts, achieving optimal sampling efficiency requires CG modeling to generate the local and long-range fluctuations that are largely consistent with those at the atomistic level. Here, we describe a new hybrid resolution CG model (HyRes) for MSES simulations of disordered protein states, which is specifically designed to provide semi-quantitative secondary structure propensities together with a qualitative description of long-range nonspecific interactions. The HyRes model contains an atomistic description of the backbone with intermediate resolution side chains. The secondary structure propensities are tuned by adjusting the backbone hydrogen-bonding strength and the ϕ/ψ torsion profile. The sizes and covalent geometries of the side chains are parameterized to reproduce distributions derived from atomistic simulations. Lennard-Jones parameters for sidechain beads are assigned to reproduce statistical potentials derived from the protein structural database, and then globally parameterized with nonspecific electrostatic interactions to reproduce the free energy profiles of pair wise interactions and the key conformational properties of model peptides. Application of HyRes to MSES simulations of small IDPs suggests that it is capable of driving faster structural transitions at the atomistic level and increasing the convergence rate compared to the Cα-only Gō-like models previously utilized. With further optimization, we believe that the new CG model could greatly improve the efficiency of MSES simulations of the larger and more complex IDPs frequently involved in cellular signalling and regulation.

Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental approaches for investigating the conformational behaviour of intrinsically disordered proteins (IDPs). IDPs represent a significant fraction of all proteomes, and, despite their importance for understanding fundamental biological processes, the molecular basis of their activity still remains largely unknown. The functional mechanisms exploited by IDPs in their interactions with other biomolecules are defined by their intrinsic dynamic modes and associated timescales, justifying the considerable interest over recent years in the development of technologies adapted to measure and describe this behaviour. NMR spin relaxation delivers information-rich, site-specific data reporting on conformational fluctuations occurring throughout the molecule. Here we review recent progress in the use of ¹⁵N relaxation to identify local backbone dynamics and long-range chain-like motions in unfolded proteins.

Deciphering the promiscuous interactions between intrinsically disordered transactivation domains and the KIX domain

The kinase-inducible domain interacting (KIX) domain of the transcriptional coactivator CBP protein carries 2 isolated binding sites (designated as the c-Myb site and the MLL site) and is capable of binding numerous intrinsically disordered transactivation domains (TADs), including c-Myb and pKID via the c-Myb site, and MLL, E2A and c-Jun via the MLL site. In this study we compared the kinetics for binding of various disordered TADs to the KIX domain via computational biophysical analyses. We found that the binding rates are heavily affected by long-range electrostatic interactions. The basal rate constants for forming the encounter complexes are similar for different KIX binding peptides, favorable electrostatic interactions between the MLL site and the peptides result in greater association rates when peptides bind to the MLL site. FOXO3a and p53 TAD each contains 2 copies of KIX binding motif and each motif interacts with both the MLL site and the c-Myb site. Our kinetics studies suggest that binding of FOXO3a or p53 TAD to the KIX domain is via a sequential mechanism, where one KIX binding motif binds to the MLL site first and then the other KIX binding motif binds to the c-Myb site. Considering the promiscuous interactions between FOXO3a and KIX, and p53 TAD and KIX, electrostatic steering simplifies the binding mechanism. This study highlights the importance of long-range electrostatic interactions in molecular recognition process involving multi-motif intrinsically disordered proteins and promiscuous interactions.

Probing the dynamics of disorder

Intrinsically disordered proteins (IDPs) play an important role in many diseases. IDPs are a large and important class of proteins; estimated to represent a significant fraction of many genomes. In contrast to protein-protein interactions between well-folded proteins, IDPs typically bind to targets using short consecutive stretches of amino acids. Structures of IDPs complexed with a target have shown great diversity in binding modes. However, how this binding diversity is achieved at the molecular level is not well understood. Unfortunately, the prediction and detailed characterization of IDPs experimentally is still a very challenging task; however molecular mechanics based molecular dynamics simulation are well suited for studying the dynamic behavior of IDPs. We look into the current state for force fields for simulating IDPs and an example of how these methods have been applied to the p53 protein. p53 is one of the most extensively studied IDPs, with multiple intrinsically disordered regulatory domains that mediate its interactions with many other proteins engaged in multiple biological pathways. We show how molecular dynamics simulations can be used to elucidate on the mechanisms involved in selection of the different binding partners.

Optimization of the GBMV2 implicit solvent force field for accurate simulation of protein conformational equilibria

Accurate treatment of solvent environment is critical for reliable simulations of protein conformational equilibria. Implicit treatment of solvation, such as using the generalized Born (GB) class of models arguably provides an optimal balance between computational efficiency and physical accuracy. Yet, GB models are frequently plagued by a tendency to generate overly compact structures. The physical origins of this drawback are relatively well understood, and the key to a balanced implicit solvent protein force field is careful optimization of physical parameters to achieve a sufficient level of cancellation of errors. The latter has been hampered by the difficulty of generating converged conformational ensembles of non-trivial model proteins using the popular replica exchange sampling technique. Here, we leverage improved sampling efficiency of a newly developed multi-scale enhanced sampling technique to re-optimize the generalized-Born with molecular volume (GBMV2) implicit solvent model with the CHARMM36 protein force field. Recursive optimization of key GBMV2 parameters (such as input radii) and protein torsion profiles (via the CMAP torsion cross terms) has led to a more balanced GBMV2 protein force field that recapitulates the structures and stabilities of both helical and β-hairpin model peptides. Importantly, this force field appears to be free of the over-compaction bias, and can generate structural ensembles of several intrinsically disordered proteins of various lengths that seem highly consistent with available experimental data. © 2017 Wiley Periodicals, Inc.

Sampling conformational space of intrinsically disordered proteins in explicit solvent: Comparison between well-tempered ensemble approach and solute tempering method

Intrinsically disordered proteins (IDPs) are a class of proteins that expected to be largely unstructured under physiological conditions. Due to their heterogeneous nature, experimental characterization of IDP is challenging. Temperature replica exchange molecular dynamics (T-REMD) is a widely used enhanced sampling method to probe structural characteristics of these proteins. However, its application has been hindered due to its tremendous computational cost, especially when simulating large systems in explicit solvent. Two methods, parallel tempering well-tempered ensemble (PT-WTE) and replica exchange with solute tempering (REST), have been proposed to alleviate the computational expense of T-REMD. In this work, we select three different IDP systems to compare the sampling characteristics and efficiencies of the two methods Both the two methods could efficiently sample the conformational space of IDP and yield highly consistent results for all the three IDPs. The efficiencies of the two methods: are compatible, with about 5-6 times better than the plain T-REMD. Besides, the advantages and disadvantages of each method are also discussed. Specially, the PT-WTE method could provide temperature dependent data of the system which could not be achieved by REST, while the REST method could readily be used to a part of the system, which is quite efficient to simulate some biological processes.

Analytical methods for structural ensembles and dynamics of intrinsically disordered proteins

Intrinsically disordered proteins, proteins that do not have a well-defined three-dimensional structure, make up a significant proportion of our proteome and are particularly prevalent in signaling and regulation. Although their importance has been realized for two decades, there is a lack of high-resolution experimental data. Molecular dynamics simulations have been crucial in reaching our current understanding of the dynamical structural ensemble sampled by intrinsically disordered proteins. In this review, we discuss enhanced sampling simulation methods that are particularly suitable to characterize the structural ensemble, along with examples of applications and limitations. The dynamics within the ensemble can be rigorously analyzed using Markov state models. We discuss recent developments that make Markov state modeling a viable approach for studying intrinsically disordered proteins. Finally, we briefly discuss challenges and future directions when applying molecular dynamics simulations to study intrinsically disordered proteins.

Enhancing backbone sampling in Monte Carlo simulations using internal coordinates normal mode analysis

Normal mode methods are becoming a popular alternative to sample the conformational landscape of proteins. In this study, we describe the implementation of an internal coordinate normal mode analysis method and its application in exploring protein flexibility by using the Monte Carlo method PELE. This new method alternates two different stages, a perturbation of the backbone through the application of torsional normal modes, and a resampling of the side chains. We have evaluated the new approach using two test systems, ubiquitin and c-Src kinase, and the differences to the original ANM method are assessed by comparing both results to reference molecular dynamics simulations. The results suggest that the sampled phase space in the internal coordinate approach is closer to the molecular dynamics phase space than the one coming from a Cartesian coordinate anisotropic network model. In addition, the new method shows a great speedup (∼5-7×), making it a good candidate for future normal mode implementations in Monte Carlo methods.

Investigating the Role of Large-Scale Domain Dynamics in Protein-Protein Interactions

Intrinsically disordered linkers provide multi-domain proteins with degrees of conformational freedom that are often essential for function. These highly dynamic assemblies represent a significant fraction of all proteomes, and deciphering the physical basis of their interactions represents a considerable challenge. Here we describe the difficulties associated with mapping the large-scale domain dynamics and describe two recent examples where solution state methods, in particular NMR spectroscopy, are used to investigate conformational exchange on very different timescales.

The Neurite Outgrowth Inhibitory Nogo-A-Δ20 Region Is an Intrinsically Disordered Segment Harbouring Three Stretches with Helical Propensity

Functional recovery from central neurotrauma, such as spinal cord injury, is limited by myelin-associated inhibitory proteins. The most prominent example, Nogo-A, imposes an inhibitory cue for nerve fibre growth via two independent domains: Nogo-A-Δ20 (residues 544-725 of the rat Nogo-A sequence) and Nogo-66 (residues 1026-1091). Inhibitory signalling from these domains causes a collapse of the neuronal growth cone via individual receptor complexes, centred around sphingosine 1-phosphate receptor 2 (S1PR2) for Nogo-A-Δ20 and Nogo receptor 1 (NgR1) for Nogo-66. Whereas the helical conformation of Nogo-66 has been studied extensively, only little structural information is available for the Nogo-A-Δ20 region. We used nuclear magnetic resonance (NMR) spectroscopy to assess potential residual structural propensities of the intrinsically disordered Nogo-A-Δ20. Using triple resonance experiments, we were able to assign 94% of the non-proline backbone residues. While secondary structure analysis and relaxation measurements highlighted the intrinsically disordered character of Nogo-A-Δ20, three stretches comprising residues 561EAIQESL567, 639EAMNVALKALGT650, and 693SNYSEIAK700 form transient α-helical structures. Interestingly, 561EAIQESL567 is situated directly adjacent to one of the most conserved regions of Nogo-A-Δ20 that contains a binding motif for β1-integrin. Likewise, 639EAMNVALKALGT650 partially overlaps with the epitope recognized by 11C7, a Nogo-A-neutralizing antibody that promotes functional recovery from spinal cord injury. Diffusion measurements by pulse-field gradient NMR spectroscopy suggest concentration- and oxidation state-dependent dimerisation of Nogo-A-Δ20. Surprisingly, NMR and isothermal titration calorimetry (ITC) data could not validate previously shown binding of extracellular loops of S1PR2 to Nogo-A-Δ20.

Disorder-to-Order Transition of an Active-Site Loop Mediates the Allosteric Activation of Sortase A

Intrinsically disordered proteins and intrinsically disordered regions are implicated in many biological functions and associated with many diseases, but their conformational characterizations are challenging. The disordered β6/β7 loop of Staphylococcus aureus sortase A is involved in the binding of both sorting signals and calcium. Calcium binding allosterically activates the enzyme, but the detailed mechanism has been unclear. Here we adapted the replica exchange with solute tempering method to sample the conformations of the β6/β7 loop, in apo form and in three liganded forms (bound with a sorting signal or calcium or both). The extensive molecular dynamics simulations yield atomic details of the disordered-to-order transition of the loop and suggest a mechanism for allosteric activation: calcium binding results in partial closure and ordering of the loop, thereby leading to preorganization of the binding pocket for the sorting signal. The approach has general applicability to the study of intrinsically disordered regions.

Simulation of coupled folding and binding of an intrinsically disordered protein in explicit solvent with metadynamics

The C-terminal domain of measles virus nucleoprotein is an intrinsically disordered protein that could bind to the X domain (XD) of phosphoprotein P to exert its physiological function. Experiments reveal that the minimal binding unit is a 21-residue α-helical molecular recognition element (α-MoRE-MeV), which adopts a fully helical conformation upon binding to XD. Due to currently limited computing power, direct simulation of this coupled folding and binding process with atomic force field in explicit solvent cannot be achieved. In this work, two advanced sampling methods, metadynamics and parallel tempering, are combined to characterize the free energy surface of this process and investigate the underlying mechanism. Starting from an unbound and partially folded state of α-MoRE-MeV, multiple folding and binding events are observed during the simulation and the energy landscape was well estimated. The results demonstrate that the isolated α-MoRE-MeV resembles a molten globule and rapidly interconverts between random coil and multiple partially helical states in solution. The coupled folding and binding process occurs through the induced fit mechanism, with the residual helical conformations providing the initial binding sites. Upon binding, α-MoRE-MeV can easily fold into helical conformation without obvious energy barriers. Two mechanisms, namely, the system tending to adopt the structure in which the free energy of isolated α-MoRE-MeV is the minimum, and the binding energy of α-MoRE-MeV to its partner protein XD tending to the minimum, jointly dominate the coupled folding and binding process. With the advanced sampling approach, more IDP systems could be simulated and common mechanisms concerning the coupled folding and binding process could be investigated in the future.

Free energy surface of an intrinsically disordered protein: comparison between temperature replica exchange molecular dynamics and bias-exchange metadynamics

Intrinsically disordered proteins (IDPs), which are expected to be largely unstructured under physiological conditions, make up a large fraction of eukaryotic proteins. Molecular dynamics simulations have been utilized to probe structural characteristics of these proteins, which are not always easily accessible to experiments. However, exploration of the conformational space by brute force molecular dynamics simulations is often limited by short time scales. Present literature provides a number of enhanced sampling methods to explore protein conformational space in molecular simulations more efficiently. In this work, we present a comparison of two enhanced sampling methods: temperature replica exchange molecular dynamics and bias exchange metadynamics. By investigating both the free energy landscape as a function of pertinent order parameters and the per-residue secondary structures of an IDP, namely, human islet amyloid polypeptide, we found that the two methods yield similar results as expected. We also highlight the practical difference between the two methods by describing the path that we followed to obtain both sets of data.

Identifying intrinsically disordered protein regions likely to undergo binding-induced helical transitions

Many proteins contain intrinsically disordered regions (IDRs) lacking stable secondary and ordered tertiary structure. IDRs are often implicated in macromolecular interactions, and may undergo structural transitions upon binding to interaction partners. However, as binding partners of many protein IDRs are unknown, these structural transitions are difficult to verify and often are poorly understood. In this study we describe a method to identify IDRs that are likely to undergo helical transitions upon binding. This method combines bioinformatics analyses followed by circular dichroism spectroscopy to monitor 2,2,2-trifluoroethanol (TFE)-induced changes in secondary structure content of these IDRs. Our results demonstrate that there is no significant change in the helicity of IDRs that are not predicted to fold upon binding. IDRs that are predicted to fold fall into two groups: one group does not become helical in the presence of TFE and includes examples of IDRs that form β-strands upon binding, while the other group becomes more helical and includes examples that are known to fold into helices upon binding. Therefore, we propose that bioinformatics analyses combined with experimental evaluation using TFE may provide a general method to identify IDRs that undergo binding-induced disorder-to-helix transitions.

Identification of Dynamic Modes in an Intrinsically Disordered Protein Using Temperature-Dependent NMR Relaxation

The dynamic modes and time scales sampled by intrinsically disordered proteins (IDPs) define their function. Nuclear magnetic resonance (NMR) spin relaxation is probably the most powerful tool for investigating these motions delivering site-specific descriptions of conformational fluctuations from throughout the molecule. Despite the abundance of experimental measurement of relaxation in IDPs, the physical origin of the measured relaxation rates remains poorly understood. Here we measure an extensive range of auto- and cross-correlated spin relaxation rates at multiple magnetic field strengths on the C-terminal domain of the nucleoprotein of Sendai virus, over a large range of temperatures (268-298 K), and combine these data to describe the dynamic behavior of this archetypal IDP. An Arrhenius-type relationship is used to simultaneously analyze up to 61 relaxation rates per amino acid over the entire temperature range, allowing the measurement of local activation energies along the chain, and the assignment of physically distinct dynamic modes. Fast (τ ≤ 50 ps) components report on librational motions, a dominant mode occurs on time scales around 1 ns, apparently reporting on backbone sampling within Ramachandran substates, while a slower component (5-25 ns) reports on segmental dynamics dominated by the chain-like nature of the protein. Extending the study to three protein constructs of different lengths (59, 81, and 124 amino acids) substantiates the assignment of these contributions. The analysis is shown to be remarkably robust, accurately predicting a broad range of relaxation data measured at different magnetic field strengths and temperatures. The ability to delineate intrinsic modes and time scales from NMR spin relaxation will improve our understanding of the behavior and function of IDPs, adding a new and essential dimension to the description of this biologically important and ubiquitous class of proteins.

Long range recognition and selection in IDPs: the interactions of the C-terminus of p53

The C-terminal domain of p53 is an extensively studied IDP, interacting with different partners through multiple distinct conformations. To explore the interplay between preformed structural elements and intrinsic fluctuations in its folding and binding we combine extensive atomistic equilibrium and non-equilibrium simulations. We find that the free peptide segment rapidly interconverts between ordered and disordered states with significant populations of the conformations that are seen in the complexed states. The underlying global folding-binding landscape points to a synergistic mechanism in which recognition is dictated via long range electrostatic recognition which results in the formation of reactive structures as far away as 10 Å, and binding proceeds with the steering of selected conformations followed by induced folding at the target surface or within a close range.

Role of Intrinsic Protein Disorder in the Function and Interactions of the Transcriptional Coactivators CREB-binding Protein (CBP) and p300

The transcriptional coactivators CREB-binding protein (CBP) and p300 undergo a particularly rich set of interactions with disordered and partly ordered partners, as a part of their ubiquitous role in facilitating transcription of genes. CBP and p300 contain a number of small structured domains that provide scaffolds for the interaction of disordered transactivation domains from a wide variety of partners, including p53, hypoxia-inducible factor 1α (HIF-1α), NF-κB, and STAT proteins, and are the targets for the interactions of disordered viral proteins that compete with cellular factors to disrupt signaling and subvert the cell cycle. The functional diversity of the CBP/p300 interactome provides an excellent example of the power of intrinsic disorder to facilitate the complexity of living systems.

Effect of O-Linked Glycosylation on the Equilibrium Structural Ensemble of Intrinsically Disordered Polypeptides

Glycosylation is one of the most common post-translational modifications (PTMs), which provides a large proteome diversity. Previous work on glycosylation of globular proteins has revealed remarkable effects of glycosylation on protein function, altering the folding stability and structure and/or altering the protein surface which affects their binding characteristics. Intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs) of large proteins are also frequently glycosylated, yet how glycosylation affects their function remains to be elucidated. An important open question is, does glycosylation affect IDP structure or binding characteristics or both? In this work, we particularly address the structural effects of O-linked glycosylation by investigating glycosylated and unglycosylated forms of two different IDPs, tau174-183 and human islet amyloid polypeptide (hIAPP), by all-atom explicit solvent simulations. We simulate these IDPs in aqueous solution for O-linked glycosylated and unglycosylated forms by employing two modern all-atom force fields for which glycan parameters are also available. We find that O-linked glycosylation only has a modest effect on equilibrium structural ensembles of IDPs, for the cases studied here, which suggests that the functional role of glycosylation may be primarily exerted by modulation of the protein binding characteristics rather than structure.

Sequence- and Temperature-Dependent Properties of Unfolded and Disordered Proteins from Atomistic Simulations

We use all-atom molecular simulation with explicit solvent to study the properties of selected intrinsically disordered proteins and unfolded states of foldable proteins, which include chain dimensions and shape, secondary structure propensity, solvent accessible surface area, and contact formation. We find that the qualitative scaling behavior of the chains matches expectations from theory under ambient conditions. In particular, unfolded globular proteins tend to be more collapsed under the same conditions than charged disordered sequences of the same length. However, inclusion of explicit solvent in addition naturally captures temperature-dependent solvation effects, which results in an initial collapse of the chains as temperature is increased, in qualitative agreement with experiment. There is a universal origin to the collapse, revealed in the change of hydration of individual residues as a function of temperature: namely, that the initial collapse is driven by unfavorable solvation free energy of individual residues, which in turn has a strong temperature dependence. We also observe that in unfolded globular proteins, increased temperature also initially favors formation of native-like (rather than non-native-like) structure. Our results help to establish how sequence encodes the degree of intrinsic disorder or order as well as its response to changes in environmental conditions.

Conformational propensities of intrinsically disordered proteins influence the mechanism of binding and folding

Intrinsically disordered proteins (IDPs) frequently function in protein interaction networks that regulate crucial cellular signaling pathways. Many IDPs undergo transitions from disordered conformational ensembles to folded structures upon binding to their cellular targets. Several possible binding mechanisms for coupled folding and binding have been identified: folding of the IDP after association with the target ("induced fit"), or binding of a prefolded state in the conformational ensemble of the IDP to the target protein ("conformational selection"), or some combination of these two extremes. The interaction of the intrinsically disordered phosphorylated kinase-inducible domain (pKID) of the cAMP-response element binding (CREB) protein with the KIX domain of a general transcriptional coactivator CREB-binding protein (CBP) provides an example of the induced-fit mechanism. Here we show by NMR relaxation dispersion experiments that a different intrinsically disordered ligand, the transactivation domain of the transcription factor c-Myb, interacts with KIX at the same site as pKID but via a different binding mechanism that involves elements of conformational selection and induced fit. In contrast to pKID, the c-Myb activation domain has a strong propensity for spontaneous helix formation in its N-terminal region, which binds to KIX in a predominantly folded conformation. The C-terminal region of c-Myb exhibits a much smaller helical propensity and likely folds via an induced-fit process after binding to KIX. We propose that the intrinsic secondary structure propensities of pKID and c-Myb determine their binding mechanisms, consistent with their functions as inducible and constitutive transcriptional activators.

Peptide Binding to a PDZ Domain by Electrostatic Steering via Nonnative Salt Bridges

We have captured the binding of a peptide to a PDZ domain by unbiased molecular dynamics simulations. Analysis of the trajectories reveals on-pathway encounter complex formation, which is driven by electrostatic interactions between negatively charged carboxylate groups in the peptide and positively charged side chains surrounding the binding site. In contrast, the final stereospecific complex, which matches the crystal structure, features completely different interactions, namely the burial of the hydrophobic side chain of the peptide C-terminal residue and backbone hydrogen bonds. The simulations show that nonnative salt bridges stabilize kinetically the encounter complex during binding. Unbinding follows the inverse sequence of events with the same nonnative salt bridges in the encounter complex. Thus, in contrast to protein folding, which is driven by native interactions, the binding of charged peptides can be steered by nonnative interactions, which might be a general mechanism, e.g., in the recognition of histone tails by bromodomains.

Modulation of the disordered conformational ensembles of the p53 transactivation domain by cancer-associated mutations

Intrinsically disordered proteins (IDPs) are frequently associated with human diseases such as cancers, and about one-fourth of disease-associated missense mutations have been mapped into predicted disordered regions. Understanding how these mutations affect the structure-function relationship of IDPs is a formidable task that requires detailed characterization of the disordered conformational ensembles. Implicit solvent coupled with enhanced sampling has been proposed to provide a balance between accuracy and efficiency necessary for systematic and comparative assessments of the effects of mutations as well as post-translational modifications on IDP structure and interaction. Here, we utilize a recently developed replica exchange with guided annealing enhanced sampling technique to calculate well-converged atomistic conformational ensembles of the intrinsically disordered transactivation domain (TAD) of tumor suppressor p53 and several cancer-associated mutants in implicit solvent. The simulations are critically assessed by quantitative comparisons with several types of experimental data that provide structural information on both secondary and tertiary levels. The results show that the calculated ensembles reproduce local structural features of wild-type p53-TAD and the effects of K24N mutation quantitatively. On the tertiary level, the simulated ensembles are overly compact, even though they appear to recapitulate the overall features of transient long-range contacts qualitatively. A key finding is that, while p53-TAD and its cancer mutants sample a similar set of conformational states, cancer mutants could introduce both local and long-range structural modulations to potentially perturb the balance of p53 binding to various regulatory proteins and further alter how this balance is regulated by multisite phosphorylation of p53-TAD. The current study clearly demonstrates the promise of atomistic simulations for detailed characterization of IDP conformations, and at the same time reveals important limitations in the current implicit solvent protein force field that must be sufficiently addressed for reliable description of long-range structural features of the disordered ensembles.

Tumor suppressor p53 is the most frequently mutated protein in human cancers. Clinical studies have suggested that the type of p53 mutation can be linked to cancer prognosis, response to drug treatment, and patient survival. It is thus crucial to understand the molecular basis of p53 inactivation by various types of mutations, so as to understand the biological outcomes and assess potential cancer intervention strategies. Here, we utilize a recently developed replica exchange with guided annealing enhanced sampling technique to calculate well-converged atomistic conformational ensembles of the intrinsically disordered transactivation domain (TAD) of tumor suppressor p53 and several cancer-associated mutants in an implicit solvent protein force field. The calculated ensembles are in quantitative agreement with several types of existing NMR data on the wild-type protein and the K24N mutant. The results suggest that, while all sequences sample a similar set of conformational substates, cancer mutants could introduce both local and long-range structural modulations and in turn perturb the balance of p53 binding to various regulatory proteins and further alter how this balance is regulated by multisite phosphorylation of p53-TAD. The study also reveals important limitations in implicit solvent for simulations of disordered proteins like p53-TAD.

Visualizing the molecular recognition trajectory of an intrinsically disordered protein using multinuclear relaxation dispersion NMR

Despite playing important roles throughout biology, molecular recognition mechanisms in intrinsically disordered proteins remain poorly understood. We present a combination of (1)H(N), (13)C', and (15)N relaxation dispersion NMR, measured at multiple titration points, to map the interaction between the disordered domain of Sendai virus nucleoprotein (NT) and the C-terminal domain of the phosphoprotein (PX). Interaction with PX funnels the free-state equilibrium of NT by stabilizing one of the previously identified helical substates present in the prerecognition ensemble in a nonspecific and dynamic encounter complex on the surface of PX. This helix then locates into the binding site at a rate coincident with intrinsic breathing motions of the helical groove on the surface of PX. The binding kinetics of complex formation are thus regulated by the intrinsic free-state conformational dynamics of both proteins. This approach, providing high-resolution structural and kinetic information about a complex folding and binding interaction trajectory, can be applied to a number of experimental systems to provide a general framework for understanding conformational disorder in biomolecular function.

Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding

The diverse biological functions of intrinsically disordered proteins (IDPs) have markedly raised our appreciation of protein conformational versatility, whereas the existence of energetically favorable yet functional detrimental nonnative interactions underscores the physical limitations of evolutionary optimization. Here we survey recent advances in using biophysical modeling to gain insight into experimentally observed nonnative behaviors and IDP properties. Simulations of IDP interactions to date focus mostly on coupled folding-binding, which follows essentially the same organizing principle as the local-nonlocal coupling mechanism in cooperative folding of monomeric globular proteins. By contrast, more innovative theories of electrostatic and aromatic interactions are needed for the conceptually novel but less-explored 'fuzzy' complexes in which the functionally bound IDPs remain largely disordered.

A relationship between the transient structure in the monomeric state and the aggregation propensities of α-synuclein and β-synuclein

α-Synuclein is an intrinsically disordered protein whose aggregation is implicated in Parkinson's disease. A second member of the synuclein family, β-synuclein, shares significant sequence similarity with α-synuclein but is much more resistant to aggregation. β-Synuclein is missing an 11-residue stretch in the central non-β-amyloid component region that forms the core of α-synuclein amyloid fibrils, yet insertion of these residues into β-synuclein to produce the βSHC construct does not markedly increase the aggregation propensity. To investigate the structural basis of these different behaviors, quantitative nuclear magnetic resonance data, in the form of paramagnetic relaxation enhancement-derived interatomic distances, are combined with molecular dynamics simulations to generate ensembles of structures representative of the solution states of α-synuclein, β-synuclein, and βSHC. Comparison of these ensembles reveals that the differing aggregation propensities of α-synuclein and β-synuclein are associated with differences in the degree of residual structure in the C-terminus coupled to the shorter separation between the N- and C-termini in β-synuclein and βSHC, making protective intramolecular contacts more likely.

Replica exchange with guided annealing for accelerated sampling of disordered protein conformations

We critically examine a recently proposed convective replica exchange (cRE) method for enhanced sampling of protein conformation based on theoretical and numerical analysis. The results demonstrate that cRE and related replica exchange with guided annealing (RE-GA) schemes lead to unbalanced exchange attempt probabilities and break detailed balance whenever the system undergoes slow conformational transitions (relative to the temperature diffusion timescale). Nonetheless, numerical simulations suggest that approximate canonical ensembles can be generated for systems with small conformational transition barriers. This suggests that RE-GA maybe suitable for simulating intrinsically disordered proteins, an important class of newly recognized functional proteins. The efficacy of RE-GA is demonstrated by calculating the conformational ensembles of intrinsically disordered kinase inducible domain protein. The results show that RE-GA helps the protein to escape nonspecific compact states more efficiently and provides several fold speedups in generating converged and largely correct ensembles compared to the standard temperature RE.

Quantifying the Sources of Kinetic Frustration in Folding Simulations of Small Proteins

Experiments and atomistic simulations of polypeptides have revealed structural intermediates that promote or inhibit conformational transitions to the native state during folding. We invoke a concept of “kinetic frustration” to quantify the prevalence and impact of these behaviors on folding rates within a large set of atomistic simulation data for 10 fast-folding proteins, where each protein’s conformational space is represented as a Markov state model of conformational transitions. Our graph theoretic approach addresses what conformational features correlate with folding inhibition and therefore permits comparison among features within a single protein network and also more generally between proteins. Nonnative contacts and nonnative secondary structure formation can thus be quantitatively implicated in inhibiting folding for several of the tested peptides.