Folding of a salivary intrinsically disordered protein upon binding to tannins

We used ion mobility spectrometry to explore conformational adaptability of intrinsically disordered proteins bound to their targets in complex mixtures. We investigated the interactions between a human salivary proline-rich protein IB5 and a model of wine and tea tannin: epigallocatechin gallate (EgCG). Collisional cross sections of naked IB5 and IB5 complexed with N = 1-15 tannins were recorded. The data demonstrate that IB5 undergoes an unfolded to folded structural transition upon binding with EgCG.     

A free-energy landscape for coupled folding and binding of an intrinsically disordered protein in explicit solvent from detailed all-atom computations

The N-terminal repressor domain of neural restrictive silencer factor (NRSF) is an intrinsically disordered protein (IDP) that binds to the paired amphipathic helix (PAH) domain of mSin3. An NMR experiment revealed that the minimal binding unit of NRSF is a 15-residue segment that adopts a helical structure upon binding to a cleft of mSin3. We computed a free-energy landscape of this system by an enhanced conformational sampling method, all-atom multicanonical molecular dynamics. The simulation started from a configuration where the NRSF segment was fully disordered and distant from mSin3 in explicit solvent. In the absence of mSin3, the disordered NRSF segment thermally fluctuated between hairpins, helices, and bent structures. In the presence of mSin3, the segment bound to mSin3 by adopting the structures involved in the isolated state, and non-native and native complexes were formed. The free-energy landscape comprised three superclusters, and free-energy barriers separated the superclusters. The native complex was located at the center of the lowest free-energy cluster. When NRSF landed in the largest supercluster, the generated non-native complex moved on the landscape to fold into the native complex, by increasing the interfacial hydrophobic contacts and the helix content. When NRSF landed in other superclusters, the non-native complex overcame the free-energy barriers between the various segment orientations in the binding cleft of mSin3. Both population-shift and induced-fit (or induced-folding) mechanisms work cooperatively in the coupled folding and binding. The diverse structural adaptability of NRSF may be related to the hub properties of the IDP.     

Long-range modulation of chain motions within the intrinsically disordered transactivation domain of tumor suppressor p53

The tumor suppressor p53 is a hub protein with a multitude of binding partners, many of which target its intrinsically disordered N-terminal domain, p53-TAD. Partners, such as the N-terminal domain of MDM2, induce formation of local structure and leave the remainder of the domain apparently disordered. We investigated segmental chain motions in p53-TAD using fluorescence quenching of an extrinsic label by tryptophan in combination with fluorescence correlation spectroscopy (PET-FCS). We studied the loop closure kinetics of four consecutive segments within p53-TAD and their response to protein binding and phosphorylation. The kinetics was multiexponential, showing that the conformational ensemble of the domain deviates from random coil, in agreement with previous findings from NMR spectroscopy. Phosphorylations or binding of MDM2 changed the pattern of intrachain kinetics. Unexpectedly, we found that upon binding and phosphorylation chain motions were altered not only within the targeted segments but also in remote regions. Long-range interactions can be induced in an intrinsically disordered domain by partner proteins that induce apparently only local structure or by post-translational modification.     

Scalable free energy calculation of proteins via multiscale essential sampling

A multiscale simulation method, "multiscale essential sampling (MSES)," is proposed for calculating free energy surface of proteins in a sizable dimensional space with good scalability. In MSES, the configurational sampling of a full-dimensional model is enhanced by coupling with the accelerated dynamics of the essential degrees of freedom. Applying the Hamiltonian exchange method to MSES can remove the biasing potential from the coupling term, deriving the free energy surface of the essential degrees of freedom. The form of the coupling term ensures good scalability in the Hamiltonian exchange. As a test application, the free energy surface of the folding process of a miniprotein, chignolin, was calculated in the continuum solvent model. Results agreed with the free energy surface derived from the multicanonical simulation. Significantly improved scalability with the MSES method was clearly shown in the free energy calculation of chignolin in explicit solvent, which was achieved without increasing the number of replicas in the Hamiltonian exchange.     

Direct observations of conformational distributions of intrinsically disordered p53 peptides using UV Raman and explicit solvent simulations

We report the first experimental measurements of Ramachandran Ψ-angle distributions for intrinsically disordered peptides: the N-terminal peptide fragment of tumor suppressor p53 and its P27S mutant form. To provide atomically detailed views of the conformational distributions, we performed classical, explicit-solvent molecular dynamics simulations on the microsecond time scale. Upon binding its partner protein, MDM2, wild-type p53 peptide adopts an α-helical conformation. Mutation of Pro27 to serine results in the highest affinity yet observed for MDM2-binding of the p53 peptide. Both UV resonance Raman spectroscopy (UVRR) and simulations reveal that the P27S mutation decreases the extent of PPII helical content and increases the probability for conformations that are similar to the α-helical MDM2-bound conformation. In addition, UVRR measurements were performed on peptides that were isotopically labeled at the Leu26 residue preceding the Pro27 in order to determine the conformational distributions of Leu26 in the wild-type and mutant peptides. The UVRR and simulation results are in quantitative agreement in terms of the change in the population of non-PPII conformations involving Leu26 upon mutation of Pro27 to serine. Finally, our simulations reveal that the MDM2-bound conformation of the peptide is significantly populated in both the wild-type and mutant isolated peptide ensembles in their unbound states, suggesting that MDM2 binding of the p53 peptides may involve conformational selection.     

Multitude of binding modes attainable by intrinsically disordered proteins: a portrait gallery of disorder-based complexes

Proteins are constantly involved in the multitude of various interactions creating sophisticated networks which define and control all (or almost all) the biological processes taking place in any living organism. Intrinsically disordered proteins or regions play a number of crucial roles in mediating protein interactions. The lack of fixed structure protruding to the high level of intrinsic dynamics and almost unrestricted flexibility at various structure levels, being the major characteristics of intrinsically disordered proteins, provides them with unprecedented advantages over the ordered proteins. The binding modes attainable by disordered proteins are highly diverse, creating a multitude of unusual complexes. Although the majority of studied to date intrinsic disorder-based complexes are ordered or static entities originating due to the global or local disorder-to-order transitions, a new development is the discovery of dynamic complexes in which intrinsically disordered proteins continue to sample an ensemble of rapidly interconverting conformations mostly devoid of structure even in their bound state. The goal of this critical review is to illustrate binding plasticity of intrinsically disordered proteins by representing a portrait gallery of the disorder-based complexes (119 references).     

Inherent flexibility and protein function: The open/closed conformational transition in the N-terminal domain of calmodulin

The key to understand a protein's function often lies in its conformational dynamics. We develop a coarse-grained variational model to investigate the interplay between structural transitions, conformational flexibility, and function of the N-terminal calmodulin domain (nCaM). In this model, two energy basins corresponding to the "closed" apo conformation and "open" holo conformation of nCaM are coupled by a uniform interpolation parameter. The resulting detailed transition route from our model is largely consistent with the recently proposed EFbeta-scaffold mechanism in EF-hand family proteins. We find that the N-terminal parts of the calcium binding loops shows higher flexibility than the C-terminal parts which form this EFbeta-scaffold structure. The structural transition of binding loops I and II are compared in detail. Our model predicts that binding loop II, with higher flexibility and earlier structural change than binding loop I, dominates the open/closed conformational transition in nCaM.     

Prediction of Protein Loop Conformations using the AGBNP Implicit Solvent Model and Torsion Angle Sampling

The OPLS-AA all-atom force field and the Analytical Generalized Born plus Non-Polar (AGBNP) implicit solvent model, in conjunction with torsion angle conformational search protocols based on the Protein Local Optimization Program (PLOP), are shown to be effective in predicting the native conformations of 57 9-residue and 35 13-residue loops of a diverse series of proteins with low sequence identity. The novel nonpolar solvation free energy estimator implemented in AGBNP augmented by correction terms aimed at reducing the occurrence of ion pairing are important to achieve the best prediction accuracy. Extended versions of the previously developed PLOP-based conformational search schemes based on calculations in the crystal environment are reported that are suitable for application to loop homology modeling without the crystal environment. Our results suggest that in general the loop backbone conformation is not strongly influenced by crystal packing. The application of the temperature Replica Exchange Molecular Dynamics (T-REMD) sampling method for a few examples where PLOP sampling is insufficient are also reported. The results reported indicate that the OPLS-AA/AGBNP effective potential is suitable for high-resolution modeling of proteins in the final stages of homology modeling and/or protein crystallographic refinement.     

Alteration of the Conformational Dynamics of a DNA Hairpin by α-Synuclein in the Presence of Aqueous Two-Phase Systems

The effect of an amyloidogenic intrinsically disordered protein, α-synuclein, which is associated with Parkinson's disease (PD), on the conformational dynamics of a DNA hairpin (DNA-HP) was studied by employing the single-molecule Förster resonance energy transfer method. The open-to-closed conformational equilibrium of the DNA-HP is drastically affected by binding of monomeric α-synuclein to the loop region of the DNA-HP. Formation of a protein-bound intermediate conformation is fostered in the presence of an aqueous two-phase system mimicking intracellular liquid-liquid phase separation. Using pressure modulation, additional mechanistic information about the binding complex could be retrieved. Hence, in addition to toxic amyloid formation, α-synuclein may alter expression profiles of disease-modifying genes in PD. Furthermore, these findings might also have significant bearings on the understanding of the physiology of organisms thriving at high pressures in the deep sea.     

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. © 2014 Wiley Periodicals, Inc.     

Helical Propensity in an Intrinsically Disordered Protein Accelerates Ligand Binding

Many intrinsically disordered proteins fold upon binding to other macromolecules. The secondary structure present in the well‐ordered complex is often formed transiently in the unbound state. The consequence of such transient structure for the binding process is, however, not clear. The activation domain of the activator for thyroid hormone and retinoid receptors (ACTR) is intrinsically disordered and folds upon binding to the nuclear coactivator binding domain (NCBD) of the CREB binding protein. A number of mutants was designed that selectively perturbs the amount of secondary structure in unbound ACTR without interfering with the intermolecular interactions between ACTR and NCBD. Using NMR spectroscopy and fluorescence‐monitored stopped‐flow kinetic measurements we show that the secondary structure content in helix 1 of ACTR indeed influences the binding kinetics. The results thus support the notion of preformed secondary structure as an important determinant for molecular recognition in intrinsically disordered proteins.     

Mapping the potential energy landscape of intrinsically disordered proteins at amino Acid resolution

Intrinsically disordered regions are predicted to exist in a significant fraction of proteins encoded in eukaryotic genomes. The high levels of conformational plasticity of this class of proteins endows them with unique capacities to act in functional modes not achievable by folded proteins, but also places their molecular characterization beyond the reach of classical structural biology. New techniques are therefore required to understand the relationship between primary sequence and biological function in this class of proteins. Although dependences of some NMR parameters such as chemical shifts (CSs) or residual dipolar couplings (RDCs) on structural propensity are known, so that sampling regimes are often inferred from experimental observation, there is currently no framework that allows for a statistical mapping of the available Ramachandran space of each amino acid in terms of conformational propensity. In this study we develop such an approach, combining highly efficient conformational sampling with ensemble selection to map the backbone conformational sampling of IDPs on a residue specific level. By systematically analyzing the ability of NMR data to map the conformational landscape of disordered proteins, we identify combinations of RDCs and CSs that can be used to raise conformational degeneracies inherent to different data types, and apply these approaches to characterize the conformational behavior of two intrinsically disordered proteins, the K18 domain from Tau protein and N(TAIL) from measles virus nucleoprotein. In both cases, we identify the enhanced populations of turn and helical regions in key regions of the proteins, as well as contiguous strands that show clear and enhanced polyproline II sampling.     

Two‐Dimensional Crowding Uncovers a Hidden Conformation of α‐Synuclein

The intrinsically disordered protein (IDP), α‐synuclein (αS), is well‐known for phospholipid membrane binding‐coupled folding into tunable helical conformers. Here, using single‐molecule experiments in conjunction with ensemble assays and a theoretical model, we present a unique case demonstrating that the interaction–folding landscape of αS can be tuned by two‐dimensional (2D) crowding through simultaneous binding of a second protein on the bilayer surface. Unexpectedly, the experimental data show a clear deviation from a simple competitive inhibition model, but are consistent with a bimodal inhibition mechanism wherein membrane binding of a second protein (a membrane interacting chaperone, Hsp27, in this case) differentially inhibits two distinct modules of αS–membrane interaction. As a consequence, αS molecules are forced to access a hidden conformational state on the phospholipid bilayer in which only the higher‐affinity module remains membrane‐bound. Our results demonstrate that macromolecular crowding in two dimensions can play a significant role in shaping the conformational landscape of membrane‐binding IDPs with multiple binding modes.     

Variation of free‐energy landscape of the p53 C‐terminal domain induced by acetylation: Enhanced conformational sampling

The C‐terminal domain (CTD) of tumor suppressor protein p53 is an intrinsically disordered region that binds to various partner proteins, where lysine of CTD is acetylated/nonacetylated and histidine neutralized/non‐neutralized. Because of the flexibility of the unbound CTD, a free‐energy landscape (FEL) is a useful quantity for determining its statistical properties. We conducted enhanced conformational sampling of CTD in the unbound state via virtual system coupled multicanonical molecular dynamics, in which the lysine was acetylated or nonacetylated and histidine was charged or neutralized. The fragments were expressed by an all‐atom model and were immersed in an explicit solvent. The acetylation and charge‐neutralization varied FEL greatly, which might be convenient to exert a hub property. The acetylation slightly enhanced alpha‐helix structures that are more compact than sheet/loop conformations. The charge‐neutralization produced hairpins. Additionally, circular dichroism experiments confirmed the computational results. We propose possible binding mechanisms of CTD to partners by investigating FEL. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Exhaustive search of the configurational space of heat-shock protein 90 with its inhibitor by multicanonical molecular dynamics based dynamic docking

Multicanonical molecular dynamics based dynamic docking was used to exhaustively search the configurational space of an inhibitor binding to the N-terminal domain of heat-shock protein 90 (Hsp90). The obtained structures at 300 K cover a wide structural ensemble, with the top two clusters ranked by their free energy coinciding with the native binding site. The representative structure of the most stable cluster reproduced the experimental binding configuration, but an interesting conformational change in Hsp90 could be observed. The combined effects of solvation and ligand binding shift the equilibrium from a preferred loop-in conformation in the unbound state to an α-helical one in the bound state for the flexible lid region of Hsp90. Thus, our dynamic docking method is effective at predicting the native binding site while exhaustively sampling a wide configurational space, modulating the protein structure upon binding.     

Mechanism‐Dependent Modulation of Ultrafast Interfacial Water Dynamics in Intrinsically Disordered Protein Complexes

The recognition of intrinsically disordered proteins (IDPs) is highly dependent on dynamics owing to the lack of structure. Here we studied the interplay between dynamics and molecular recognition in IDPs with a combination of time‐resolving tools on timescales ranging from femtoseconds to nanoseconds. We interrogated conformational dynamics and surface water dynamics and its attenuation upon partner binding using two IDPs, IBB and Nup153FG, both of central relevance to the nucleocytoplasmic transport machinery. These proteins bind the same nuclear transport receptor (Importinβ) with drastically different binding mechanisms, coupled folding–binding and fuzzy complex formation, respectively. Solvent fluctuations in the dynamic interface of the Nup153FG‐Importinβ fuzzy complex were largely unperturbed and slightly accelerated relative to the unbound state. In the IBB‐Importinβ complex, on the other hand, substantial relative slowdown of water dynamics was seen in a more rigid interface. These results show a correlation between interfacial water dynamics and the plasticity of IDP complexes, implicating functional relevance for such differential modulation in cellular processes, including nuclear transport.     

Adsorption of homopolypeptides on gold investigated using atomistic molecular dynamics

We investigate the role of dynamics on adsorption of peptides to gold surfaces using all-atom molecular dynamics simulations in explicit solvent. We choose six homopolypeptides [Ala(10), Ser(10), Thr(10), Arg(10), Lys(10), and Gln(10)], for which experimental surface coverages are not correlated with amino acid level affinities for gold, with the idea that dynamic properties may also play a role. To assess dynamics we determine both conformational movement and flexibility of the peptide within a given conformation. Low conformational movement indicates stability of a given conformation and leads to less adsorption than homopolypeptides with faster conformational movement. Likewise, low flexibility within a given conformation also leads to less adsorption. Neither amino acid affinities nor dynamic considerations alone predict surface coverage; rather both quantities must be considered in peptide adsorption to gold surfaces.     

How well does molecular simulation reproduce environment-specific conformations of the intrinsically disordered peptides PLP,TP2 and ONEG?

Understanding the conformational ensembles of intrinsically disordered proteins and peptides (IDPs) in their various biological environments is essential for understanding their mechanisms and functional roles in the proteome, leading to a greater knowledge of, and potential treatments for, a broad range of diseases. To determine whether molecular simulation is able to generate accurate conformational ensembles of IDPs, we explore the structural landscape of the PLP peptide (an intrinsically disordered region of the proteolipid membrane protein) in aqueous and membrane-mimicking solvents, using replica exchange with solute scaling (REST2), and examine the ability of four force fields (ff14SB, ff14IDPSFF, CHARMM36 and CHARMM36m) to reproduce literature circular dichroism (CD) data. Results from variable temperature (VT) ¹H and Rotating frame Overhauser Effect SpectroscopY (ROESY) nuclear magnetic resonance (NMR) experiments are also presented and are consistent with the structural observations obtained from the simulations and CD. We also apply the optimum simulation protocol to TP2 and ONEG (a cell-penetrating peptide (CPP) and a negative control peptide, respectively) to gain insight into the structural differences that may account for the observed difference in their membrane-penetrating abilities. Of the tested force fields, we find that CHARMM36 and CHARMM36m are best suited to the study of IDPs, and accurately predict a disordered to helical conformational transition of the PLP peptide accompanying the change from aqueous to membrane-mimicking solvents. We also identify an α-helical structure of TP2 in the membrane-mimicking solvents and provide a discussion of the mechanistic implications of this observation with reference to the previous literature on the peptide. From these results, we recommend the use of CHARMM36m with the REST2 protocol for the study of environment-specific IDP conformations. We believe that the simulation protocol will allow the study of a broad range of IDPs that undergo conformational transitions in different biological environments.

A protocol for simulating intrinsically disordered peptides in aqueous and hydrophobic solvents is proposed. Results from four force fields are compared with experiment. CHARMM36m performs the best for the simulated IDPs in all environments.