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.     

Assessing induced folding of an intrinsically disordered protein by site-directed spin-labeling electron paramagnetic resonance spectroscopy

We used site-directed spin-labeling electron paramagnetic resonance (EPR) spectroscopy to study the induced folding of the intrinsically disordered C-terminal domain of measles virus nucleoprotein (N(TAIL)). Four single-site N(TAIL) mutants (S407C, S488C, L496C, and V517C), located in three conserved regions, were prepared and labeled with a nitroxide paramagnetic probe. We could monitor the gain of rigidity that N(TAIL) undergoes in the presence of either the secondary structure stabilizer 2,2,2-trifluoroethanol (TFE) or one of its physiological partners, namely, the C-terminal domain (XD) of the viral phosphoprotein. The mobility of the spin label grafted at positions 488, 496, and 517 was significantly reduced upon addition of XD, contrary to that of the spin label bound to position 407, which was unaffected. Furthermore, the EPR spectra of spin-labeled S488C and L496C bound to XD in the presence of 30% sucrose are indicative of the formation of an alpha-helix in the proximity of the spin labels. Such an alpha-helix had been already identified by previous biochemical and structural studies. Using TFE we unveiled a previously undetected structural propensity within the N-terminal region of N(TAIL) and showed that its C-terminal region "resists" gaining structure even at high TFE concentrations. Finally, we for the first time showed the reversibility of the induced folding process that N(TAIL) undergoes in the presence of XD. These results highlight the suitability of site-directed spin-labeling EPR spectroscopy to identify protein regions involved in binding and folding events, while providing insights at the residue level.     

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.     

Long‐lived states in an intrinsically disordered protein domain

Long‐lived states (LLS) are relaxation‐favored spin population distributions of J‐coupled magnetic nuclei. LLS were measured, along with classical ¹H and ¹⁵N relaxation rate constants, in amino acids of the N‐terminal Unique domain of the c‐Src kinase, which is disordered in vitro under physiological conditions. The relaxation rates of LLS can probe motions and interactions in biomolecules. LLS of the aliphatic protons of glycines, with lifetimes approximately four times longer than their spin–lattice relaxation times, are reported for the first time in an intrinsically disordered protein domain. LLS relaxation experiments were integrated with 2D spectroscopy methods, further adapting them for studies on proteins. Copyright © 2013 John Wiley & Sons, Ltd.     

Rate constants and mechanisms of intrinsically disordered proteins binding to structured targets

The binding of intrinsically disordered proteins (IDPs) to structured targets is gaining increasing attention. Here we review experimental and computational studies on the binding kinetics of IDPs. Experiments have yielded both the binding rate constants and the binding mechanisms, the latter via mutation and deletion studies and NMR techniques. Most computational studies have aimed at qualitative understanding of the binding rate constants or at mapping the free energy surfaces after the IDPs are engaged with their targets. The experiments and computation show that IDPs generally gain structures after they are engaged with their targets; that is, interactions with the targets facilitate the IDPs' folding. It also seems clear that the initial contact of an IDP with the target is formed by just a segment, not the entire IDP. The docking of one segment to its sub-site followed by coalescing of other segments around the corresponding sub-sites emerges as a recurring feature in the binding of IDPs. Such a dock-and-coalesce model forms the basis for quantitative calculation of binding rate constants. For both disordered and ordered proteins, strong electrostatic attraction with their targets can enhance the binding rate constants by several orders of magnitude. There are now tremendous opportunities in narrowing the gap in our understanding of IDPs relative to ordered proteins with regard to binding kinetics.     

Differential dynamical effects of macromolecular crowding on an intrinsically disordered protein and a globular protein: implications for in-cell NMR spectroscopy

In-cell NMR provides a valuable means to assess how macromolecules, with concentrations up to 300 g/L in the cytoplasm, affect the structure and dynamics of proteins at atomic resolution. Here an intrinsically disordered protein, alpha-synuclein (alphaSN), and a globular protein, chymotrypsin inhibitor 2 (CI2) were examined by using in-cell NMR. High-resolution in-cell spectra of alphaSN can be obtained, but CI2 leaks from the cell and the remaining intracellular CI2 is not detectable. Even after stabilizing the cells from leakage by using alginate encapsulation, no CI2 signal is detected. From in vitro studies we conclude that this difference in detectability is the result of the differential dynamical response of disordered and ordered proteins to the changes of motion caused by the increased viscosity in cells.     

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).     

Thermal stability and folding kinetics analysis of disordered protein,securin

Lacking a stable tertiary structure, intrinsically disordered proteins (IDPs) possess particular functions in cell regulation, signaling, and controlling pathways. The study of their unique structural features, thermal stabilities, and folding kinetics is intriguing. In this study, an identified IDP, securin, was used as a model protein. By using a quasi-static five-step (on-path) folding process, the function of securin was restored and analyzed by isothermal titration calorimetry. Fluorescence spectroscopy and particle size analysis indicated that securin possessed a compact hydrophobic core and particle size. The glass transition of securin was characterized using differential scanning microcalorimetry. Furthermore, the folding/unfolding rates (k_obs) of securin were undetectable, implying that the folding/unfolding rate is very fast and that the conformation of securin is sensitive to solvent environmental change. Therefore, securin may fold properly under specific physiological conditions. In summary, the thermal glass transition behavior and undetectable k_obs of folding/unfolding reactions may be two of the indices of IDP.     

Proteus: a random forest classifier to predict disorder-to-order transitioning binding regions in intrinsically disordered proteins

The focus of the computational structural biology community has taken a dramatic shift over the past one-and-a-half decades from the classical protein structure prediction problem to the possible understanding of intrinsically disordered proteins (IDP) or proteins containing regions of disorder (IDPR). The current interest lies in the unraveling of a disorder-to-order transitioning code embedded in the amino acid sequences of IDPs/IDPRs. Disordered proteins are characterized by an enormous amount of structural plasticity which makes them promiscuous in binding to different partners, multi-functional in cellular activity and atypical in folding energy landscapes resembling partially folded molten globules. Also, their involvement in several deadly human diseases (e.g. cancer, cardiovascular and neurodegenerative diseases) makes them attractive drug targets, and important for a biochemical understanding of the disease(s). The study of the structural ensemble of IDPs is rather difficult, in particular for transient interactions. When bound to a structured partner, an IDPR adapts an ordered conformation in the complex. The residues that undergo this disorder-to-order transition are called protean residues, generally found in short contiguous stretches and the first step in understanding the modus operandi of an IDP/IDPR would be to predict these residues. There are a few available methods which predict these protean segments from their amino acid sequences; however, their performance reported in the literature leaves clear room for improvement. With this background, the current study presents ‘Proteus’, a random forest classifier that predicts the likelihood of a residue undergoing a disorder-to-order transition upon binding to a potential partner protein. The prediction is based on features that can be calculated using the amino acid sequence alone. Proteus compares favorably with existing methods predicting twice as many true positives as the second best method (55 vs. 27%) with a much higher precision on an independent data set. The current study also sheds some light on a possible ‘disorder-to-order’ transitioning consensus, untangled, yet embedded in the amino acid sequence of IDPs. Some guidelines have also been suggested for proceeding with a real-life structural modeling involving an IDPR using Proteus.     

The catalytic mechanism of an aspartic proteinase explored with neutron and X-ray diffraction

Hydrogen atoms play key roles in enzyme mechanism, but as this study shows, even high-quality X-ray data to a resolution of 1 A cannot directly visualize them. Neutron diffraction, however, can locate deuterium atoms even at resolutions around 2 A. Both neutron and X-ray diffraction data have been used to investigate the transition state of the aspartic proteinase endothiapepsin. The different techniques reveal a different part of the story, revealing the clearest picture yet of the catalytic mechanism by which the enzyme operates. Room temperature neutron and X-ray diffraction data were used in a newly developed joint refinement software package to visualize deuterium atoms within the active site of the enzyme when a gem-diol transition state analogue inhibitor is bound at the active site. These data were also used to estimate their individual occupancy, while analysis of the differences between the bond lengths of the catalytic aspartates was performed using atomic resolution X-ray data. The two methods are in agreement on the protonation state of the active site with a transition state analogue inhibitor bound confirming the catalytic mechanism at which the enzyme operates.     

The importance of the compact disordered state in the fuzzy interactions between intrinsically disordered proteins

The intrinsically disordered C-terminal domain (CTD) of protein 4.1G is able to specifically bind a 26-residue intrinsically disordered region of NuMA, forming a dynamic fuzzy complex. As one of a few cases of extremely fuzzy interactions between two intrinsically disordered proteins/regions (IDPs/IDRs) without induced folding, the principle of the binding is unknown. Here, we combined experimental and computational methods to explore the detailed mechanism of the interaction between 4.1G-CTD and NuMA. MD simulations suggest that the kinetic hub states in the structure ensemble of 4.1G-CTD are favorable in the fuzzy complex. The feature of these hub states is that the binding ‘hot spot’ motifs βA and βB exhibit β strand propensities and are well packed to each other. The binding between 4.1G-CTD and NuMA is disrupted at low pH, which changes the intramolecular packing of 4.1G-CTD and weakens the packing between βA and βB motifs. Low pH conditions also lead to increased hydrodynamic radius and acceleration of backbone dynamics of 4.1G-CTD. All these results underscore the importance of tertiary structural arrangements and overall compactness of 4.1G-CTD in its binding to NuMA, i.e. the compact disordered state of 4.1G-CTD is crucial for binding. Different from the short linear motifs (SLiMs) that are often found to mediate IDP interactions, 4.1G-CTD functions as an intrinsically disordered domain (IDD), which is a functional and structural unit similar to conventional protein domains. This work sheds light on the molecular recognition mechanism of IDPs/IDRs and expands the conventional structure-function paradigm in protein biochemistry.

Tertiary structural arrangements and overall compactness is important for interactions between IDPs.     

The acidic domains of the Toc159 chloroplast preprotein receptor family are intrinsically disordered protein domains

Background  

The Toc159 family of proteins serve as receptors for chloroplast-destined preproteins. They directly bind to transit peptides, and exhibit preprotein substrate selectivity conferred by an unknown mechanism. The Toc159 receptors each include three domains: C-terminal membrane, central GTPase, and N-terminal acidic (A-) domains. Although the function(s) of the A-domain remains largely unknown, the amino acid sequences are most variable within these domains, suggesting they may contribute to the functional specificity of the receptors.     

Continuous thermal collapse of the intrinsically disordered protein tau is driven by its entropic flexible domain

The tau protein belongs to the category of Intrinsically Disordered Proteins (IDP), which in their native state lack a folded structure and fluctuate between many conformations. In its physiological state, tau helps nucleating and stabilizing the microtubules' (MTs) surfaces in the axons of the neurons. Tau is mainly composed by two domains: (i) the binding domain that tightly bounds the MT surfaces and (ii) the projection domain that exerts a long-range entropic repulsive force and thus provides the proper spacing between adjacent MTs. Tau is also involved in the genesis and in the development of the Alzheimer disease when it detaches from MT surfaces and aggregates in paired helical filaments. Unfortunately, the molecular mechanisms behind these phenomena are still unclear. Temperature variation, rarely considered in biological studies, is here used to provide structural information on tau correlated to its role as an entropic spacer between adjacent MTs surfaces. In this paper, by means of small-angle X-ray scattering and molecular dynamics simulation, we demonstrate that tau undergoes a counterintuitive collapse phenomenon with increasing temperature. A detailed analysis of our results, performed by the Ensemble Optimization Method, shows that the thermal collapse is coupled to the occurrence of a transient long-range contact between a region encompassing the end of the proline-rich domain P2 and the first part of the repeats domain, and the region of the N-terminal domain entailing residues 80-150. Interestingly these two regions involved in the tau temperature collapse belong to the flexible projection domain that acts as an entropic bristle and regulates the MTs' architecture. Our results show that temperature is an important parameter that influences the dynamics of the tau projection domain, and hence its entropic behavior.     

Kinetics of protein binding in solid-phase immunoassays: theory

In a solid-phase immunoassay, binding between an antigen and its specific antibody takes place at the boundary of a liquid and a solid phase. One of the reactants (receptor) is immobilized on a surface. The other reactant (ligand) is initially free in solution. We present a theory describing the kinetics of immunochemical reaction in such a system. A single essential restriction of the theory is the assumption that the reaction conditions are uniform along the binding surface. In general, the reaction rate as a function of time can be obtained numerically as a solution of a nonlinear integral equation. For some special cases, analytical solutions are available. Various immunoassay geometries are considered, in particular, the case when the reaction is carried out on a microspot.     

Coupling ligand recognition to protein folding in an engineered variant of rabbit ileal lipid binding protein

We have engineered a variant of the beta-clam shell protein ILBP which lacks the alpha-helical motif that caps the central binding cavity; the mutant protein is sufficiently destabilised that it is unfolded under physiological conditions, however, it unexpectedly binds its natural bile acid substrates with high affinity forming a native-like beta-sheet rich structure and demonstrating strong thermodynamic coupling between ligand binding and protein folding.     

Disorder-to-order transition of an intrinsically disordered region of sortase revealed by multiscale enhanced sampling

Molecular functions of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs), such as molecular recognition and cellular signaling, are ascribed to dynamic changes in the conformational space in response to binding of target molecules. Sortase, a transpeptitase in Gram-positive bacteria, has an IDR in a loop which undergoes a disordered-to-ordered transition (called "disordered loop"), accompanying a tilt of another loop ("dynamic loop"), upon binding of a signal peptide and a calcium ion. In this study, all-atom conformational ensembles of sortase were calculated for the four different binding states (with/without the peptide and with/without a calcium ion) by the multiscale enhanced sampling (MSES) simulation to examine how the binding of the peptide and/or calcium influences the conformational ensemble. The MSES is a multiscale and multicopy simulation method that allows an enhanced sampling of the all-atom model of large proteins including explicit solvent. A 100 ns MSES simulation of the ligand-free sortase using 20 replicas (in total 2 μs) demonstrated large flexibility in both the disordered and dynamic loops; however, their distributions were not random but had a clear preference which populates the N-terminal part of the disordered loop near the bound form. The MSES simulations of the three binding states clarified the allosteric mechanism of sortase: the N- and C-terminal parts of the disordered loop undergo a disorder-to-order transition independently of each other upon binding of the peptide and a calcium ion, respectively; however, upon binding of both ligands, the two parts work cooperatively to stabilize the bound peptide.     

Computational strategy for intrinsically disordered protein ligand design leads to the discovery of p53 transactivation domain I binding compounds that activate the p53 pathway

Intrinsically disordered proteins or intrinsically disordered regions (IDPs) have gained much attention in recent years due to their vital roles in biology and prevalence in various human diseases. Although IDPs are perceived as attractive therapeutic targets, rational drug design targeting IDPs remains challenging because of their conformational heterogeneity. Here, we propose a hierarchical computational strategy for IDP drug virtual screening (IDPDVS) and applied it in the discovery of p53 transactivation domain I (TAD1) binding compounds. IDPDVS starts from conformation sampling of the IDP target, then it combines stepwise conformational clustering with druggability evaluation to identify potential ligand binding pockets, followed by multiple docking screening runs and selection of compounds that can bind multi-conformations. p53 is an important tumor suppressor and restoration of its function provides an opportunity to inhibit cancer cell growth. TAD1 locates at the N-terminus of p53 and plays key roles in regulating p53 function. No compounds that directly bind to TAD1 have been reported due to its highly disordered structure. We successfully used IDPDVS to identify two compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells. Our study demonstrates that IDPDVS is an efficient strategy for IDP drug discovery and p53 TAD1 can be directly targeted by small molecules.

A hierarchical computational strategy for IDP drug virtual screening (IDPDVS) was proposed and successfully applied to identify compounds that bind p53 TAD1 and restore wild-type p53 function in cancer cells.     

Kinetics of cooperative protein folding involving two separate conformational families

The master equation that describes the kinetics of protein folding is solved numerically for a portion of Staphylococcal Protein A by a Laplace transformation. The calculations are carried out with 50 local-minimum conformations belonging to two conformational families. The master equation allows for transitions among all the 50 conformations in the evolution toward the final folded equilibrium distribution of conformations. It is concluded that the native protein folds in a fast cooperative process. The global energy minimum of a native protein can be reached after a sufficiently long folding time regardless of the initial state and the existence of a large number of local energy minima. Conformations representing non-native states of the protein can transform to the native state even if they do not belong to the native conformational family. Given a starting conformation, the protein molecule can fold to its final conformation through different paths. Finally, when the folding reaches the equilibrium distribution, the protein molecule adopts a set of conformations in which the global minimum has the largest average probability.