High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

Universality in the timescales of internal loop formation in unfolded proteins and single-stranded oligonucleotides

Understanding the rate at which various parts of a molecular chain come together to facilitate the folding of a biopolymer (e.g., a protein or RNA) into its functional form remains an elusive goal. Here we use experiments, simulations, and theory to study the kinetics of internal loop closure in disordered biopolymers such as single-stranded oligonucleotides and unfolded proteins. We present theoretical arguments and computer simulation data to show that the relationship between the timescale of internal loop formation and the positions of the monomers enclosing the loop can be recast in a form of a universal master dependence. We also perform experimental measurements of the loop closure times of single-stranded oligonucleotides and show that both these and previously reported internal loop closure kinetics of unfolded proteins are well described by this theoretically predicted dependence. Finally, we propose that experimental deviations from the master dependence can then be used as a sensitive probe of dynamical and structural order in unfolded proteins and other biopolymers.     

Mechanism of induced folding: Both folding before binding and binding before folding can be realized in staphylococcal nuclease mutants

A considerable number of functional proteins are unstructured under physiological condition. These "intrinsically disordered" proteins exhibit induced folding when they bind their targets. The induced folding comprises two elementary processes: folding and binding. Two mechanisms are possible for the induced folding: either folding before binding or binding before folding. We found that these two mechanisms can be distinguished by the target-concentration dependence of folding kinetics. We also created two types of mutants of staphylococcal nuclease showing the different inhibitor-concentration dependence of induced folding kinetics. One mutant obeys the scheme of binding before folding, while the other the folding before binding. This is the first experimental evidence demonstrating that both mechanisms are realized for a single protein. Binding before folding is possible, when the protein lacks essential nonlocal interaction to stabilize the native conformation. The results cast light on the protein folding mechanism involved in the intrinsically disordered proteins.     

Folding and stability of sweet protein single-chain monellin. An insight to protein engineering

Engineered single-chain monellin (SCM) proteins were constructed by recombinant technology without disrupting the topology and sweet activity of native protein. Data from 8-anilinonaphthalene-1-sulfonic acid fluorescence, size-exclusion chromatography, and heteronuclear NMR strongly suggest the presence of a folding intermediate at 1.5 m GdnHCl for SCM protein. The structural feature of the folding intermediate from NMR data reveals that the secondary structures became mostly unstable, and protein experiences a dynamic equilibrium between native and unfolded state. All backbone amide protons exchange within 10 min, which imply that no stable hydrogen bonds exist in the secondary structural regions in the folding intermediate. From equilibrium unfolding and mutagenesis studies, the unfolding transition midpoints of mutant proteins gradually shifted toward lower denaturant concentration, indicating stability reductions of mutant proteins. Our results suggest that stability and folding pathways of SCM proteins could be regulated by a combined study of spectroscopy and mutagenesis, and these studies will provide useful information for understanding the folding kinetics of novel engineered proteins.     

Atomic-level characterization of disordered protein ensembles

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

Protein folding kinetics beyond the phi value: using multiple amino acid substitutions to investigate the structure of the SH3 domain folding transition state

The SH3 domain folding transition state structure contains two well-ordered turn regions, known as the diverging turn and the distal loop. In the Src SH3 domain transition state, these regions are stabilized by a hydrogen bond between Glu30 in the diverging turn and Ser47 in the distal loop. We have examined the effects on folding kinetics of amino acid substitutions at the homologous positions (Glu24 and Ser41) in the Fyn SH3 domain. In contrast to most other folding kinetics studies which have focused primarily on non-disruptive substitutions with Ala or Gly, here we have examined the effects of substitutions with diverse amino acid residues. Using this approach, we demonstrate that the transition state structure is generally tolerant to amino acid substitutions. We also uncover a unique role for Ser at position 41 in facilitating folding of the distal loop, which can only be replicated by Asp at the same position. Both these residues appear to accelerate folding through the formation of short-range side-chain to backbone hydrogen bonds. The folding of the diverging turn region is shown to be driven primarily by local interactions. The diverging turn and distal loop regions are found to interact in the transition state structure, but only in the context of particular mutant backgrounds. This work demonstrates that studying the effects of a variety of amino acid substitutions on protein folding kinetics can provide unique insights into folding mechanisms which cannot be obtained by standard Phi value analysis.     

Intrinsically disordered protein from a pathogenic mesophile Mycobacterium tuberculosis adopts structured conformation at high temperature

Compared to eukaryotes, the occurrence of "intrinsically disordered" or "natively unfolded" proteins in prokaryotes has not been explored extensively. Here, we report the occurrence of an intrinsically disordered protein from the mesophilic human pathogen Mycobacterium tuberculosis. The Histidine-tagged recombinant Rv3221c biotin-binding protein is intrinsically disordered at ambient and physiological growth temperatures as revealed by circular dichroism and Fourier transform infrared (FTIR) spectroscopic studies. However, an increase in temperature induces a transition from disordered to structured state with a folding temperature of approximately 53 degrees C. Addition of a structure inducing solvent trifluoroethanol (TFE) causes the protein to fold at lower temperatures suggesting that TFE fosters hydrophobic interactions, which drives protein folding. Differential Scanning Calorimetry studies revealed that folding is endothermic and the transition from a disordered to structured state is continuous (higher-order), implying existence of intermediates during folding process. Secondary structure analysis revealed that the protein has propensity to form beta-sheets. This is in conformity with FTIR spectrum that showed an absorption peak at wave number of 1636 cm(-1), indicative of disordered beta-sheet conformation in the native state. These data suggest that although Rv3221c may be disordered under ambient or optimal growth temperature conditions, it has the potential to fold into ordered structure at high temperature driven by increased hydrophobic interactions. In contrast to the generally known behavior of other intrinsically disordered proteins folding at high temperature, Rv3221c does not appear to oligomerize or aggregate as revealed through numerous experiments including Congo red binding, Thioflavin T-binding, turbidity measurements, and examining molar ellipticity as a function of protein concentration. The amino acid composition of Rv3221c reveals that it has 24% charged and 54.9% hydrophobic amino acid residues. In this respect, this protein, although belonging to the class of intrinsically disordered proteins, has distinct features. The intrinsically disordered state and the biotin-binding feature of this protein suggest that it may participate in many biochemical processes requiring biotin as a cofactor and adopt suitable conformations upon binding other folded targets.     

Hydrogen-Bond Driven Loop-Closure Kinetics in Unfolded Polypeptide Chains

Characterization of the length dependence of end-to-end loop-closure kinetics in unfolded polypeptide chains provides an understanding of early steps in protein folding. Here, loop-closure in poly-glycine-serine peptides is investigated by combining single-molecule fluorescence spectroscopy with molecular dynamics simulation. For chains containing more than 10 peptide bonds loop-closing rate constants on the 20–100 nanosecond time range exhibit a power-law length dependence. However, this scaling breaks down for shorter peptides, which exhibit slower kinetics arising from a perturbation induced by the dye reporter system used in the experimental setup. The loop-closure kinetics in the longer peptides is found to be determined by the formation of intra-peptide hydrogen bonds and transient β-sheet structure, that accelerate the search for contacts among residues distant in sequence relative to the case of a polypeptide chain in which hydrogen bonds cannot form. Hydrogen-bond-driven polypeptide-chain collapse in unfolded peptides under physiological conditions found here is not only consistent with hierarchical models of protein folding, that highlights the importance of secondary structure formation early in the folding process, but is also shown to speed up the search for productive folding events.     

Dimensions,energetics, and denaturant effects of the protein unstructured state

Determining the energetics of the unfolded state of a protein is essential for understanding the folding mechanics of ordered proteins and the structure–function relation of intrinsically disordered proteins. Here, we adopt a coil‐globule transition theory to develop a general scheme to extract interaction and free energy information from single‐molecule fluorescence resonance energy transfer spectroscopy. By combining protein stability data, we have determined the free energy difference between the native state and the maximally collapsed denatured state in a number of systems, providing insight on the specific/nonspecific interactions in protein folding. Both the transfer and binding models of the denaturant effects are demonstrated to account for the revealed linear dependence of inter‐residue interactions on the denaturant concentration, and are thus compatible under the coil‐globule transition theory to further determine the dimension and free energy of the conformational ensemble of the unfolded state. The scaling behaviors and the effective θ‐state are also discussed.     

Refolding simulations of an isolated fragment of barnase into a native-like β hairpin: Evidence for compactness and hydrogen bonding as concurrent stabilizing factors

Experimental evidence and theoretical models both suggest that protein folding is initiated within specific fragments intermittently adopting conformations close to that found in the protein native structure. These folding initiation sites encompassing short portions of the protein are ideally suited for study in isolation by computational methods aimed at peering into the very early events of folding. We have used Molecular Dynamics (MD) technique to investigate the behavior of an isolated protein fragment formed by residues 85 to 102 of barnase that folds into a β hairpin in the protein native structure. Three independent MD simulations of 1.3 to 1.8 ns starting from unfolded conformations of the peptide portrayed with an all-atom model in water were carried out at gradually decreasing temperature. A detailed analysis of the conformational preferences adopted by this peptide in the course of the simulations is presented. Two of the unfolded peptide conformations fold into a hairpin characterized by native and a larger bulk of nonnative interactions. Both refolding simulations substantiate the close relationship between interstrand compactness and hydrogen bonding network involving backbone atoms. Persistent compactness witnessed by side-chain interactions always occurs concomitantly with the formation of backbone hydrogen bonds. No highly populated conformations generated in a third simulation starting from the remotest unfolded conformer relative to the native structure are observed. However, nonnative long-range and medium-range contacts with the aromatic moiety of Trp94 are spotted, which are in fair agreement with a former nuclear magnetic resonance study of a denaturing solution of an isolated barnase fragment encompassing the β hairpin. All this lends reason to believe that the 85–102 barnase fragment is a strong initiation site for folding. Proteins 29:212–227, 1997. © 1997 Wiley-Liss, Inc.     

Delineation of a new structural motif involving NHN γ-turn

Macromolecules are characterized by distinctive arrangement of hydrogen bonds. Different patterns of hydrogen bonds give rise to distinct and stable structural motifs. An analysis of 4114 non-redundant protein chains reveals the existence of a three-residue, (i − 1) to (i + 1), structural motif, having two hydrogen-bonded five-membered pseudo rings (the first, an N H···OC involving the first residue, and the second being N H∙∙∙N involving the last two residues), separated by a peptide bond. There could be an additional hydrogen bond between the side-chain at (i-1) and the main-chain NH of (i + 1). The average backbone torsion angles of −76(±21)° and – 12(±17)° at i creates a tight turn in the polypeptide chain, akin to a γ-turn. Indeed, a search of three-residue fragments with restriction on the terminal C^α···C^α distance and the existence of the two pseudo rings on either side revealed the presence 14 846 cases of a variant, termed NHN γ-turn, distinct from the NHO γ-turn (2032 cases) that has traditionally been characterized by the presence of NHO hydrogen bond linking the terminal main-chain atoms. As in the latter, the newly identified γ-turns are also of two types—classical and inverse, occurring in the ratio of 1:6. The propensities of residues to occur in these turns and their secondary structural features have been enumerated. An understanding of these turns would be useful for structure prediction and loop modeling, and may serve as models to represent some of the unfolded state or disordered region in proteins.     

On the unyielding hydrophobic core of villin headpiece

Villin headpiece (HP67) is a small, autonomously-folding domain that has become a model system for understanding the fundamental tenets governing protein folding. In this communication, we explore the role that Leu61 plays in the structure and stability of the construct. Deletion of Leu61 results in a completely unfolded protein that cannot be expressed in Escherichia coli. Omission of only the aliphatic leucine side chain (HP67 L61G) perturbed neither the backbone conformation nor the orientation of local hydrophobic side chains. As a result, a large, solvent-exposed hydrophobic pocket, a negative replica of the leucine side-chain, was created on the surface. The loss of the hydrophobic interface between leucine 61 and the hydrophobic pocket destabilized the construct by ~3.3 kcal/mol. Insertion of a single glycine residue immediately before Leu61 (HP67 L61[GL]) was also highly destabilizing and had the effect of altering the backbone conformation (α-helix to π-helix) in order to precisely preserve the wild-type position and conformation of all hydrophobic residues, including Leu61. In addition to demonstrating that the hydrophobic side-chain of Leu61 is critically important for the stability of villin headpiece, our results are consistent with the notion that the precise interactions present within the hydrophobic core, rather than the hydrogen bonds that define the secondary structure, specify a protein's fold.     

The Role of Backbone Hydrogen Bonds in the Transition State for Protein Folding of a PDZ Domain

Backbone hydrogen bonds are important for the structure and stability of proteins. However, since conventional site-directed mutagenesis cannot be applied to perturb the backbone, the contribution of these hydrogen bonds in protein folding and stability has been assessed only for a very limited set of small proteins. We have here investigated effects of five amide-to-ester mutations in the backbone of a PDZ domain, a 90-residue globular protein domain, to probe the influence of hydrogen bonds in a β-sheet for folding and stability. The amide-to-ester mutation removes NH-mediated hydrogen bonds and destabilizes hydrogen bonds formed by the carbonyl oxygen. The overall stability of the PDZ domain generally decreased for all amide-to-ester mutants due to an increase in the unfolding rate constant. For this particular region of the PDZ domain, it is therefore clear that native hydrogen bonds are formed after crossing of the rate-limiting barrier for folding. Moreover, three of the five amide-to-ester mutants displayed an increase in the folding rate constant suggesting that the hydrogen bonds are involved in non-native interactions in the transition state for folding.     

Limited proteolysis of natively unfolded protein 4E‐BP1 in the presence of trifluoroethanol

Natively unfolded (intrinsically disordered (ID) proteins) have been attracting an increasing attention due to their involvement in many regulatory processes. Natively unfolded proteins can fold upon binding to their metabolic partners. Coupled folding and binding events usually involve only relatively short motifs (binding motifs). These binding motifs which are able to fold should have an increased propensity to form a secondary structure. The aim of the present work was to probe the conformation of the intrinsically disordered protein 4E‐BP1 in the native and partly folded states by limited proteolysis and to reveal regions with a high propensity to form an ordered structure. Trifuoroethanol (TFE) in low concentrations (up to 15 vol%) was applied to increase the helical population of protein regions with a high intrinsic propensity to fold. When forming helical structures, these regions lose mobility and become more protected from proteases than random/unfolded protein regions. Limited proteolysis followed by mass spectrometry analysis allows identification of the regions with decreased mobility in TFE solutions. Trypsin and V8 proteases were used to perform limited proteolysis of the 4E‐BP1 protein in buffer and in solutions with low TFE concentrations at 37°C and at elevated temperatures (42 and 50°C). Comparison of the results obtained with the previously established 4E‐BP1 structure and the binding motif illustrates the ability of limited proteolysis in the presence of a folding assistant (TFE) to map the regions with high and low propensities to form a secondary structure revealing potential binding motifs inside the intrinsically disordered protein. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 591–602, 2014.     

The contribution of electrostatics to hydrogen exchange in the unfolded protein state

Although electrostatics have long been recognized to play an important role in hydrogen exchange (HX) with solvent, the quantitative assessment of its magnitude in the unfolded state has hitherto been lacking. This limits the utility of HX as a quantitative method to study protein stability, folding, and dynamics. Using the intrinsically disordered human protein α-synuclein as a proxy for the unfolded state, we show that a hybrid mean-field approach can effectively compute the electrostatic potential at all backbone amide positions along the chain. From the electrochemical potential, a fourfold reduction in hydroxide concentration near the protein backbone is predicted for the C-terminal domain, a prognosis that is in direct agreement with experimentally derived protection factors from NMR spectroscopy. Thus, impeded HX for the C-terminal region of α-synuclein is not the result of intramolecular hydrogen bonding and/or structure formation.     

Protein domain definition should allow for conditional disorder

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

Equilibrium hydrogen exchange reveals extensive hydrogen bonded secondary structure in the on-pathway intermediate of Im7

The four-helical immunity protein Im7 folds through an on-pathway intermediate that has a specific, but partially misfolded, hydrophobic core. In order to gain further insight into the structure of this species, we have identified the backbone hydrogen bonds formed in the ensemble by measuring the amide exchange rates (under EX2 conditions) of the wild-type protein and a variant, I72V. In this mutant the intermediate is significantly destabilised relative to the unfolded state (deltadeltaG(ui) = 4.4 kJ/mol) but the native state is only slightly destabilised (deltadeltaG(nu) = 1.8 kJ/mol) at 10 degrees C in 2H2O, pH* 7.0 containing 0.4 M Na2SO4, consistent with the view that this residue forms significant non-native stabilising interactions in the intermediate state. Comparison of the hydrogen exchange rates of the two proteins, therefore, enables the state from which hydrogen exchange occurs to be identified. The data show that amides in helices I, II and IV in both proteins exchange slowly with a free energy similar to that associated with global unfolding, suggesting that these helices form highly protected hydrogen-bonded helical structure in the intermediate. By contrast, amides in helix III exchange rapidly in both proteins. Importantly, the rate of exchange of amides in helix III are slowed substantially in the Im7* variant, I72V, compared with the wild-type protein, whilst other amides exchange more rapidly in the mutant protein, in accord with the kinetics of folding/unfolding measured using chevron analysis. These data demonstrate, therefore, that local fluctuations do not dominate the exchange mechanism and confirm that helix III does not form stable secondary structure in the intermediate. By combining these results with previously obtained Phi-values, we show that the on-pathway folding intermediate of Im7 contains extensive, stable hydrogen-bonded structure in helices I, II and IV, and that this structure is stabilised by both native and non-native interactions involving amino acid side-chains in these helices.     

Forces stabilizing proteins

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site-directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a –CH₂− group on folding contributes 1.1 ± 0.5 kcal/mol to protein stability. (2) The burial of non-polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1 ± 0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.     

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.