Beyond the random coil: stochastic conformational switching in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) participate in critical cellular functions that exploit the flexibility and rapid conformational fluctuations of their native state. Limited information about the native state of IDPs can be gained by the averaging over many heterogeneous molecules that is unavoidable in ensemble approaches. We used single molecule fluorescence to characterize native state conformational dynamics in five synaptic proteins confirmed to be disordered by other techniques. For three of the proteins, SNAP-25, synaptobrevin and complexin, their conformational dynamics could be described with a simple semiflexible polymer model. Surprisingly, two proteins, neuroligin and the NMDAR-2B glutamate receptor, were observed to stochastically switch among distinct conformational states despite the fact that they appeared intrinsically disordered by other measures. The hop-like intramolecular diffusion found in these proteins is suggested to define a class of functionality previously unrecognized for IDPs.     

Structural and dynamic perturbations induced by heme binding in cytochrome b5

The water-soluble domain of rat hepatic cytochrome b(5) undergoes marked structural changes upon heme removal. The solution structure of apocytochrome b(5) shows that the protein is partially folded in the absence of the heme group, exhibiting a stable module and a disordered heme-binding loop. The quality of the apoprotein structure in solution was improved with the use of heteronuclear NMR data. Backbone amide hydrogen exchange was studied to characterize cooperative units in the protein. It was found that this criterion distinguished the folded module from the heme-binding loop in the apoprotein, in contrast to the holoprotein. The osmolyte trimethylamine N-oxide (TMAO) did not affect the structure of the apoprotein in the disordered region. TMAO imparted a small stabilization consistent with an unfolded state effect correlating with the extent of buried surface area in the folded region of the native apoprotein. The failure of the osmolyte to cause large conformational shifts in the disordered loop supported the view that the specificity of the local sequence for the holoprotein fold was best developed with the stabilization of the native state through heme binding. To dissect the role of the heme prosthetic group in forcing the disordered region into the holoprotein conformation, the axial histidine belonging to the flexible loop (His63) was replaced with an alanine, and the structural properties of the protein with carbon-monoxide-ligated reduced iron were studied. The His63Ala substitution resulted in a protein with lower heme affinity but nevertheless capable of complete refolding. This indicated that the coordination bond was not necessary to establish the structural features of the holoprotein. In addition, the weak binding of the heme in this protein resulted in conformational shifts at a location distant from the binding site. The data suggested an uneven distribution of cooperative elements in the structure of the cytochrome.     

Hierarchical Conformational Analysis of Native Lysozyme Based on Sub-Millisecond Molecular Dynamics Simulations

Hierarchical organization of free energy landscape (FEL) for native globular proteins has been widely accepted by the biophysics community. However, FEL of native proteins is usually projected onto one or a few dimensions. Here we generated collectively 0.2 milli-second molecular dynamics simulation trajectories in explicit solvent for hen egg white lysozyme (HEWL), and carried out detailed conformational analysis based on backbone torsional degrees of freedom (DOF). Our results demonstrated that at micro-second and coarser temporal resolutions, FEL of HEWL exhibits hub-like topology with crystal structures occupying the dominant structural ensemble that serves as the hub of conformational transitions. However, at 100ns and finer temporal resolutions, conformational substates of HEWL exhibit network-like topology, crystal structures are associated with kinetic traps that are important but not dominant ensembles. Backbone torsional state transitions on time scales ranging from nanoseconds to beyond microseconds were found to be associated with various types of molecular interactions. Even at nanoseconds temporal resolution, the number of conformational substates that are of statistical significance is quite limited. These observations suggest that detailed analysis of conformational substates at multiple temporal resolutions is both important and feasible. Transition state ensembles among various conformational substates at microsecond temporal resolution were observed to be considerably disordered. Life times of these transition state ensembles are found to be nearly independent of the time scales of the participating torsional DOFs.     

Structural characterization of the native NH2-terminal transactivation domain of the human androgen receptor: a collapsed disordered conformation underlies structural plasticity and protein-induced folding

The androgen receptor (AR) mediates the action of the steroid hormones testosterone and dihydrotestosterone. The protein contains two globular alpha-helical domains responsible for binding hormone and DNA. In contrast, the N-terminal domain is less well structurally defined and contains the main determinants for receptor-dependent transactivation, termed AF1. Previously, we have shown this region has the propensity to form alpha-helix structure. Significantly, the binding of specific protein targets or a natural osmolyte resulted in a more protease resistant conformation for the AF1 domain, consistent with an increase in conformational stability. Computational and experimental analyses were used to investigate the conformational properties of the native AF1 domain. This region of the receptor is predicted to contain significant regions of natural disordered structure, when analyzed by amino acid composition, PONDR (Predictor of Natural Disordered Regions), RONN (Regional Order Neural Network), and GlobPlot, but is grouped with ordered proteins on a charge-hydropathy plot. The binding of a hydrophobic fluorescence probe, 8-anilinonaphthalene-1-sulfonic acid (ANS), together with size-exclusion chromatography suggests that native AR-AF1 exists in a collapsed disordered conformation, distinct from extended disordered (random coil) and a stable globular fold. This state has also been described as premolten or molten globule-like. These findings are discussed in terms of the functional importance of the intrinsic plasticity of the AF1 domain.     

Improved structural characterizations of the drkN SH3 domain unfolded state suggest a compact ensemble with native-like and non-native structure

Due to their dynamic ensemble nature and a deficiency of experimental restraints, disordered states of proteins are difficult to characterize structurally. Here, we have expanded upon our previous work on the unfolded state of the Drosophila drk N-terminal (drkN) SH3 domain with our program ENSEMBLE, which assigns population weights to pregenerated conformers in order to calculate ensembles of structures whose properties are collectively consistent with experimental measurements. The experimental restraint set has been enlarged with newly measured paramagnetic relaxation enhancements from Cu(2+) bound to an amino terminal Cu(2+)-Ni(2+) binding (ATCUN) motif as well as nuclear Overhauser effect (NOE) and hydrogen exchange data from recent studies. In addition, two new pseudo-energy minimization algorithms have been implemented that have dramatically improved the speed of ENSEMBLE population weight assignment. Finally, we have greatly improved our conformational sampling by utilizing a variety of techniques to generate both random structures and structures that are biased to contain elements of native-like or non-native structure. Although it is not possible to uniquely define a representative structural ensemble, we have been able to assess various properties of the drkN SH3 domain unfolded state by performing ENSEMBLE minimizations of different conformer pools. Specifically, we have found that the experimental restraint set enforces a compact structural distribution that is not consistent with an overall native-like topology but shows preference for local non-native structure in the regions corresponding to the diverging turn and the beta5 strand of the folded state and for local native-like structure in the region corresponding to the beta6 and beta7 strands. We suggest that this approach could be generally useful for the structural characterization of disordered states.     

Osmolyte‐induced conformational changes in the Hsp90 molecular chaperone

Osmolytes are small molecules that play a central role in cellular homeostasis and the stress response by maintaining protein thermodynamic stability at controlled levels. The underlying physical chemistry that describes how different osmolytes impact folding free energy is well understood, however little is known about their influence on other crucial aspects of protein behavior, such as native‐state conformational changes. Here we investigate this issue with the Hsp90 molecular chaperone, a large dimeric protein that populates a complex conformational equilibrium. Using small angle X‐ray scattering we observe dramatic osmolyte‐dependent structural changes within the native ensemble. The degree to which different osmolytes affect the Hsp90 conformation strongly correlates with thermodynamic metrics of their influence on stability. This observation suggests that the well‐established osmolyte principles that govern stability also apply to large‐scale conformational changes, a proposition that is corroborated by structure‐based fitting of the scattering data, surface area comparisons and m‐value analysis. This approach shows how osmolytes affect a highly cooperative open/closed structural transition between two conformations that differ by a domain‐domain interaction. Hsp90 adopts an additional ligand‐specific conformation in the presence of ATP and we find that osmolytes do not significantly affect this conformational change. Together, these results extend the scope of osmolytes by suggesting that they can maintain protein conformational heterogeneity at controlled levels using similar underlying principles that allow them to maintain protein stability; however the relative impact of osmolytes on different structural states can vary significantly.     

Automated use of mutagenesis data in structure prediction

In the absence of experimental structural determination, numerous methods are available to indirectly predict or probe the structure of a target molecule. Genetic modification of a protein sequence is a powerful tool for identifying key residues involved in binding reactions or protein stability. Mutagenesis data is usually incorporated into the modeling process either through manual inspection of model compatibility with empirical data, or through the generation of geometric constraints linking sensitive residues to a binding interface. We present an approach derived from statistical studies of lattice models for introducing mutation information directly into the fitness score. The approach takes into account the phenotype of mutation (neutral or disruptive) and calculates the energy for a given structure over an ensemble of sequences. The structure prediction procedure searches for the optimal conformation where neutral sequences either have no impact or improve stability and disruptive sequences reduce stability relative to wild type. We examine three types of sequence ensembles: information from saturation mutagenesis, scanning mutagenesis, and homologous proteins. Incorporating multiple sequences into a statistical ensemble serves to energetically separate the native state and misfolded structures. As a result, the prediction of structure with a poor force field is sufficiently enhanced by mutational information to improve accuracy. Furthermore, by separating misfolded conformations from the target score, the ensemble energy serves to speed up conformational search algorithms such as Monte Carlo-based methods.     

Conformational heterogeneity of an equilibrium folding intermediate quantified and mapped by scanning mutagenesis

It is challenging to experimentally define an energy landscape for protein folding that comprises multiple partially unfolded states. Experimental results are often ambiguous as to whether a non-native state is conformationally homogeneous. Here, we tested an approach combining systematic mutagenesis and a Br?nsted-like analysis to reveal and quantify conformational heterogeneity of folding intermediate states. Using this method, we resolved an otherwise apparently homogeneous equilibrium folding intermediate of Borrelia burgdorferi OspA into two conformationally distinct species and determined their relative populations. Furthermore, we mapped the structural differences between these intermediate species, which are consistent with the non-native species that we previously proposed based on native-state hydrogen exchange studies. When treated as a single state, the intermediate ensemble exhibited fractional Phi-values for mutations and Hammond-type behaviors that are often observed for folding transition states. We found that a change in relative population of the two species within the intermediate ensemble explains these properties well, suggesting that fractional Phi-values and Hammond-type behaviors exhibited by folding intermediates and transition states may arise more often from conformational heterogeneity than from a single partial structure. Our results are consistent with the presence of multiple minima in a rugged energy landscape predicted from theoretical studies. The method described here provides a promising means to probe a complex folding energy landscape.     

Effects of cell volume regulating osmolytes on glycerol 3-phosphate binding to triosephosphate isomerase

During cell volume regulation, intracellular concentration changes occur in both inorganic and organic osmolytes in order to balance the extracellular osmotic stress and maintain cell volume homeostasis. Generally, salt and urea increase the Km's of enzymes and trimethylamine N-oxide (TMAO) counteracts these effects by decreasing Km's. The hypothesis to account for these effects is that urea and salt shift the native state ensemble of the enzyme toward conformers that are substrate-binding incompetent (BI), while TMAO shifts the ensemble toward binding competent (BC) species. Km's are often complex assemblies of rate constants involving several elementary steps in catalysis, so to better understand osmolyte effects we have focused on a single elementary event, substrate binding. We test the conformational shift hypothesis by evaluating the effects of salt, urea, and TMAO on the mechanism of binding glycerol 3-phosphate, a substrate analogue, to yeast triosephosphate isomerase. Temperature-jump kinetic measurements promote a mechanism consistent with osmolyte-induced shifts in the [BI]/[BC] ratio of enzyme conformers. Importantly, salt significantly affects the binding constant through its effect on the activity coefficients of substrate, enzyme, and enzyme-substrate complex, and it is likely that TMAO and urea affect activity coefficients as well. Results indicate that the conformational shift hypothesis alone does not account for the effects of osmolytes on Km's.     

Exploring the Minimally Frustrated Energy Landscape of Unfolded ACBP

The unfolded state of globular proteins is not well described by a simple statistical coil due to residual structural features, such as secondary structure or transiently formed long-range contacts. The principle of minimal frustration predicts that the unfolded ensemble is biased toward productive regions in the conformational space determined by the native structure. Transient long-range contacts, both native-like and non-native-like, have previously been shown to be present in the unfolded state of the four-helix-bundle protein acyl co-enzyme binding protein (ACBP) as seen from both perturbations in nuclear magnetic resonance (NMR) chemical shifts and structural ensembles generated from NMR paramagnetic relaxation data. To study the nature of the contacts in detail, we used paramagnetic NMR relaxation enhancements, in combination with single-point mutations, to obtain distance constraints for the acid-unfolded ensemble of ACBP. We show that, even in the acid-unfolded state, long-range contacts are specific in nature and single-point mutations affect the free-energy landscape of the unfolded protein. Using this approach, we were able to map out concerted, interconnected, and productive long-range contacts. The correlation between the native-state stability and compactness of the denatured state provides further evidence for native-like contact formation in the denatured state. Overall, these results imply that, even in the earliest stages of folding, ACBP dynamics are governed by native-like contacts on a minimally frustrated energy landscape.     

Structured disorder and conformational selection

Traditionally, molecular disorder has been viewed as local or global instability. Molecules or regions displaying disorder have been considered inherently unstructured. The term has been routinely applied to cases for which no atomic coordinates can be derived from crystallized molecules. Yet, even when it appears that the molecules are disordered, prevailing conformations exist, with population times higher than those of all alternate conformations. Disordered molecules are the outcome of rugged energy landscapes away from the native state around the bottom of the funnel. Ruggedness has a biological function, creating a distribution of structured conformers that bind via conformational selection, driving association and multimolecular complex formation, whether chain-linked in folding or unlinked in binding. We classify disordered molecules into two types. The first type possesses a hydrophobic core. Here, even if the native conformation is unstable, it still has a large enough population time, enabling its experimental detection. In the second type, no such hydrophobic core exists. Hence, the native conformations of molecules belonging to this category have shorter population times, hindering their experimental detection. Although there is a continuum of distribution of hydrophobic cores in proteins, an empirical, statistically based hydrophobicity function may be used as a guideline for distinguishing the two disordered molecule types. Furthermore, the two types relate to steps in the protein folding reaction. With respect to protein design, this leads us to propose that engineering-optimized specific electrostatic interactions to avoid electrostatic repulsion would reduce the type I disordered state, driving the molten globule (MG) --> native (N) state. In contrast, for overcoming the type II disordered state, in addition to specific interactions, a stronger hydrophobic core is also indicated, leading to the denatured --> MG --> N state.     

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.     

Compaction of ribosomal protein S6 by sucrose occurs only under native conditions

The effect of osmolyte sucrose on the stability and compaction of the folded and unfolded states of ribosomal protein S6 from Thermus thermophilus was analyzed. Confirming previous results obtained with sodium sulfate and trehalose, refolding stopped-flow measurements of S6 show that sucrose favors the conversion of the unfolded state ensemble to a highly compact structure (75% as compact as the folded state). This conversion occurs when the unfolded state is suddenly placed under native conditions and the compact state accumulates in a transient off-folding pathway. This effect of sucrose on the compaction of the unfolded state ensemble is counteracted by guanidinium hydrochloride. The compact state does not accumulate at higher guanidinium concentrations and the unfolded state ensemble does not display increased compaction in the presence of 6 M guanidinium as evaluated by collisional quenching of tryptophan fluorescence. In contrast, accessibility of the tryptophan residue of folded S6 above 1 M sucrose concentration decreased as a result of an increased compaction of the folded state. Unfolding stopped-flow measurements of S6 reflect this increased compaction of the folded state, but the unfolding pathway is not affected by sucrose. Compaction of folded and unfolded S6 induced by sucrose occurs under native conditions indicating that decreased protein conformational entropy significantly contributes to the mechanism of protein stabilization by osmolytes.     

Direct evidence for a two-state protein unfolding transition from hydrogen-deuterium exchange, mass spectrometry, and NMR.

We use mass spectrometry in conjunction with hydrogen-deuterium exchange and NMR to characterize the conformational dynamics of the 62-residue IgG binding domain of protein L under conditions in which the native state is marginally stable. Mass spectra of protein L after short incubations in D2O reveal the presence of two distinct populations containing different numbers of protected protons. NMR experiments indicate that protons in the hydrophobic core are protected in one population, whereas all protons are exchanged for deuterons in the other. As the exchange period is increased, molecules are transferred from the former population to the latter. The absence of molecules with a subset of the core protons protected suggests that exchange occurs in part via a highly concerted transition to an excited state in which all protons exchange rapidly with deuterons. A steady increase in the molecular weight of the population with protected protons, and variation in the exchange rates of the individual protected protons indicates the presence of an additional exchange mechanism. A simple model in which exchange results from rapid (> 10(5)/s) local fluctuations around the native state superimposed upon transitions to an unfolded excited state at approximately 0.06/s is supported by qualitative agreement between the observed mass spectra and the mass spectra simulated according to the model using NMR-derived estimates of the proton exchange rates.     

The transition state of coupled folding and binding for a flexible β-finger

Flexible and fully disordered protein regions that fold upon binding mediate numerous protein-protein interactions. However, little is known about their mechanism of interaction. One such coupled folding and binding occurs when a flexible region of neuronal nitric oxide synthase adopts a β-finger structure upon binding to its protein ligand, a PDZ [PSD-95 (postsynaptic density protein-95)/Discs large/ZO-1] domain from PSD-95. We have analyzed this binding reaction by protein engineering combined with kinetic experiments. Mutational destabilization of the β-finger changed mainly the dissociation rate constant of the proteins and, to a lesser extent, the association rate constant. Thus, mutation affected late events in the coupled folding and binding reaction. Our results therefore suggest that the native binding interactions of the β-finger are not present in the rate-limiting transition state for binding but form on the downhill side in a cooperative manner. However, by mutation, we could destabilize the β-finger further and change the rate-limiting step such that an initial conformational change becomes rate limiting. This switch in rate-limiting step shows that multistep binding mechanisms are likely to be found among flexible and intrinsically disordered regions of proteins.     

Observing the osmophobic effect in action at the single molecule level

Protecting osmolytes are widespread small organic molecules able to stabilize the folded state of most proteins against various denaturing stresses in vivo. The osmophobic model explains thermodynamically their action through a preferential exclusion of the osmolyte molecules from the protein surface, thus favoring the formation of intrapeptide hydrogen bonds. Few works addressed the influence of protecting osmolytes on the protein unfolding transition state and kinetics. Among those, previous single molecule force spectroscopy experiments evidenced a complexation of the protecting osmolyte molecules at the unfolding transition state of the protein, in apparent contradiction with the osmophobic nature of the protein backbone. We present single-molecule evidence that glycerol, which is a ubiquitous protecting osmolyte, stabilizes a globular protein against mechanical unfolding without binding into its unfolding transition state structure. We show experimentally that glycerol does not change the position of the unfolding transition state as projected onto the mechanical reaction coordinate. Moreover, we compute theoretically the projection of the unfolding transition state onto two other common reaction coordinates, that is, the number of native peptide bonds and the weighted number of native contacts. To that end, we augment an analytic Ising-like protein model with support for group-transfer free energies. Using this model, we find again that the position of the unfolding transition state does not change in the presence of glycerol, giving further support to the conclusions based on the single-molecule experiments.