Two-state Folding Research Articles

Kinetic folding studies of the P22 tailspike beta-helix domain reveal multiple unfolded states

The beta-helix is an important protein fold in many pathogens, and is a challenging system for folding pathway prediction because it primarily is stabilized by non-local interactions along the primary sequence. A useful experimental model of this fold is a monomeric truncation of P22 tailspike protein, the beta-helix domain (bhx). This report describes a systematic in vitro study of the chemical denaturation and refolding of bhx. Results from equilibrium chemical denaturation experiments were consistent with a two-state folding mechanism, but showed only partial reversibility. Stopped-flow fluorescence studies showed a single unfolding step, but two refolding steps. The slow refolding step could be partly attributed to proline isomerization, based on an increased rate during refolding in the presence of PPIase and an increased relative amplitude of this step with increasing delay time in double-jump refolding experiments observed over delays of 5–100 s. However, double-jump refolding experiments with delay times longer than 100 s along with size exclusion chromatography and dynamic light scattering of refolding samples showed that the overall refolding yield decreased as bhx was unfolded for longer periods of time. Furthermore, the losses resulted from aggregate formation during refolding. This suggests that a change occurs over time in the unfolded or denatured state ensemble that increases the propensity for aggregation upon the shift to more native-favoring conditions. Alternatively aggregate nuclei may be able to form even under high denaturant conditions, and these subsequently grow and consume monomer when placed under native-favoring conditions.

What have we learned from the studies of two-state folders, and what are the unanswered questions about two-state protein folding?

Small proteins with globular structures often fold by simple all-or-none mechanisms, both in an equilibrium and a kinetic sense, despite the very large number of partly folded conformations available. This type of ‘two-state’ folding will be discussed in terms of experimental tests, underlying molecular mechanisms, and limits to two-state behavior. Factors that appear to be important for two-state folding include topology (sequence distance of contacts in the native structure), molecular cooperativity and local energy distribution. Because their local stability distributions and cooperativities can be dissected and analyzed separately from topological features, recent studies of the folding of symmetric proteins will be discussed as a means to better understand the origins of two-state folding.

The perspectives of studying multi-domain protein folding.

Most of fundamental studies on protein folding have been performed with small globular proteins consisting of a single domain. In vitro many of these proteins are well characterized by a reversible two-state folding scheme. However, the majority of proteins in the cell belong to the class of larger multi-domain proteins that often unfold irreversibly under in vitro conditions. This makes folding studies difficult or even impossible. In spite of these problems for many multi-domain proteins, folding has been investigated by classical refolding. Co-translational folding of nascent polypeptide chains when synthesized by ribosomes has also been studied. Single molecule techniques represent a promising approach for future studies on the folding of multi-domain proteins, and tremendous advances have been made in these techniques in recent years. In particular, fluorescence-based methods can contribute significantly to an understanding of the fundamental principles of multi-domain protein folding.

Measurement of Single Molecule Folding/unfolding Trajectories

We have measured folding/unfolding trajectories of single protein G (B1 domain) molecules, a simple two-state folder, by simultaneously measuring the fluorescence intensity, lifetime, and spectrum at various concentrations of denaturant. Protein molecules were labeled by a fluorescence resonance energy transfer (FRET) pair, Alexa Fluor 488 and Alexa Fluor 594 and were immobilized on a glass surface coated with polyethyleneglycol via streptavidin-biotin linkage. The vast majority of molecules (∼ 85%) exhibits simple two-state trajectories, with either high or low values of the FRET efficiency, corresponding to the folded and unfolded states, respectively, with unresolvable jumps between them. About 10% of the trajectories show transitions in the unfolded state that can be attributed to a ∼20 nm spectral shift of the donor, as revealed by measurements of their emission spectra. The mean FRET efficiency of immobilized molecules matches the value measured in free diffusion experiments. There is a distribution of these values beyond the width expected from shot noise, which can, however, be quantitatively accounted for by the distribution of acceptor lifetimes. In spite of these complications from photophysics, rate coefficients obtained from the exponential distribution of residence times in either the folded or unfolded state yield relaxation times that agree within a factor of 2 with those measured on the dye-labeled protein by stopped flow kinetics. In addition, no correlation is observed between the donor and acceptor intensity in the unfolded state from microseconds to seconds suggesting that structural averaging between unfolded conformations occurs on the nanosecond timescale, as expected from previous measurements by B. Schuler and coworkers (PNAS:104,2655,2007). All these results indicate that we have successfully immobilized the protein without significantly altering its structure, kinetics, or dynamics, and represent a major step forward toward the goal of “watching” individual molecules fold.

Experimental Studies on Protein Folding in the Presence of the Hsp70 Chaperone System

The Hsp70 chaperone system (including DnaK, DnaJ and GrpE in bacteria) plays a vital role in preventing aggregation and assisting protein folding. While a lot of the chaperone biochemistry has already been worked out, very little is know about the interaction of DnaK and DnaJ with substrates. So far, such interaction had only been probed with small peptide or protein substrates incapable of independent/efficient refolding. This paucity of information has precluded studies on the competition between inter-molecular binding and intra-molecular folding. Such competition is important for cell viability because the accessibility of DnaK to substrate binding is modulated by the intrinsic rates of substrate folding/unfolding and by the substrate's thermodynamic stability.We developed a computational kinetic model to predict the interaction of protein substrates exhibiting two-state folding behavior with DnaK, DnaJ and GrpE. We found that, under physiological conditions, only proteins that fold slowly and/or have moderate thermodynamic stability are predicted to bind chaperones during their folding cycle.Experimental protein folding kinetics followed by stopped-flow in the presence and absence of the Hsp70 chaperone system shows good agreement with the predictions by the computational model. Furthermore, gel filtration and reverse phase chromatography data further support the stopped-flow results by providing evidence for DnaK-DnaJ-substrate interactions at equilibrium. This result is in agreement with the predictions of the computational model. In summary, the combination of experiments and computational predictions developed in this work is a powerful tool to help unveiling the relations between protein folding and chaperone binding.

Monte Carlo Simulations of Proteins in Cages: Influence of Confinement on the Stability of Intermediate States

We study the folding of small proteins inside confining potentials. Proteins are described using an effective potential model that contains the Ramachandran angles as degrees of freedom and does not need any a priori information about the native state. Hydrogen bonds, dipole-dipole-, and hydrophobic interactions are taken explicitly into account. An interesting feature displayed by this potential is the presence of metastable intermediates between the unfolded and native states. We consider different types of confining potentials to describe proteins folding inside cages with repulsive or attractive walls. Using the Wang-Landau algorithm, we determine the density of states and analyze in detail the thermodynamical properties of the confined proteins for different sizes of the cages. We show that confinement dramatically reduces the phase space available to the protein and that the presence of intermediate states can be controlled by varying the properties of the confining potential. Cages with strongly attractive walls destabilize the intermediate states and lead to a two-state folding into a configuration that is less stable than the native structure. However, cages with slightly attractive walls enhance the stability of native structure and induce a folding process, which occurs through intermediate configurations.

Folding rate prediction using complex network analysis for proteins with two- and three-state folding kinetics

It is a challenging task to investigate the different in- fluence of long-range and short-range interactions on two-state and three-state folding kinetics of protein. The networks of the 30 two-state proteins and 15 three-state proteins were constructed by complex networks analysis at three length scales: Protein Contact Networks, Long-range Interaction Networks and Short-range Interaction Networks. To uncover the relationship between structural properties and folding kinetics of the proteins, the correlations of protein network parameters with protein folding rate and topology parameters contact order were analyzed. The results show that Protein Contact Networks and Short-range Interaction Networks (for both two-state and three-state proteins) exhibit the “small-world” property and Long-range Interaction networks indicate “scale-free” behavior. Our results further indicate that all Protein Contact Networks and Short- range Interaction networks are assortative type. While some of Long-range Interaction Networks are of assortative type, the others are of disassortative type. For two-state proteins, the clustering coefficients of Short-range Interaction Networks show prominent correlation with folding rate and contact order. The assortativity coefficients of Short-range Interaction Networks also show remarkable correlation with folding rate and contact order. Similar correlations exist in Protein Contact Networks of three-state proteins. For two-state proteins, the correlation between contact order and folding rate is determined by the numbers of local contacts. Short- range interactions play a key role in determining the connecting trend among amino acids and they impact the folding rate of two-state proteins directly. For three-state proteins, the folding rate is determined by short-range and long-range interactions among residues together.

Structural ultrafast dynamics of macromolecules: diffraction of free DNA and effect of hydration

Of special interest in molecular biology is the study of structural and conformational changes which are free of the additional effects of the environment. In the present contribution, we report on the ultrafast unfolding dynamics of a large DNA macromolecular ensemble in vacuo for a number of temperature jumps, and make a comparison with the unfolding dynamics of the DNA in aqueous solution. A number of coarse-graining approaches, such as kinetic intermediate structure (KIS) model and ensemble-averaged radial distribution functions, are used to account for the transitional dynamics of the DNA without sacrificing the structural resolution. The studied ensembles of DNA macromolecules were generated using distributed molecular dynamics (MD) simulations, and the ensemble convergence was ensured by monitoring the ensemble-averaged radial distribution functions and KIS unfolding trajectories. Because the order-disorder transition in free DNA implies unzipping, coiling, and strand-separation processes which occur consecutively or competitively depending on the initial and final temperature of the ensemble, DNA order-disorder transition in vacuo cannot be described as a two-state (un)folding process.

Two-state protein folding kinetics through all-atom molecular dynamics based sampling

This review focuses on advanced computational techniques that employ all atom molecular dynamics to study the folding of small two state proteins. As protein folding is a rare event process, special sampling techniques are required to overcome high folding free energy barriers. Several biased sampling methods enable computation of the free energy landscape. Trajectory based sampling methods can assess the kinetics and the dynamical folding mechanisms. Proper sampling is only the first step, and further analysis is required to obtain the folding mechanisms reaction coordinate. Only a combination of several simulation techniques can solve the sampling problems connected with all-atom protein folding, and allow computation of experimental observables that can validate the force fields and simulation techniques. Several of the involved issues are illustrated with folding of small protein (fragments) such as beta hairpins and the Trp-cage mini protein.

Ruggedness in the folding landscape of protein L

By exploring the folding pathways of the B1 domain of protein L with a series of equilibrium and rapid kinetic experiments, we have found its unfolded state to be more complex than suggested by two-state folding models. Using an ultrarapid mixer to initiate protein folding within approximately 2-4 microseconds, we observe folding kinetics by intrinsic tryptophan fluorescence and fluorescence resonance energy transfer. We detect at least two processes faster than 100 mus that would be hidden within the burst phase of a stopped-flow instrument measuring tryptophan fluorescence. Previously reported measurements of slow intramolecular diffusion are commensurate with the slower of the two observed fast phases. These results suggest that a multidimensional energy landscape is necessary to describe the folding of protein L, and that the dynamics of the unfolded state is dominated by multiple small energy barriers.

Calculation of Protein Heat Capacity from Replica-Exchange Molecular Dynamics Simulations with Different Implicit Solvent Models

The heat capacity has played a major role in relating microscopic and macroscopic properties of proteins and their disorder-order phase transition of folding. Its calculation by atomistic simulation methods remains a significant challenge due to the complex and dynamic nature of protein structures, their solvent environment, and configurational averaging. To better understand these factors on calculating a protein heat capacity, we provide a comparative analysis of simulation models that differ in their implicit solvent description and force-field resolution. Our model protein system is the src Homology 3 (SH3) domain of alpha-spectrin, and we report a series of 10 ns replica-exchange molecular dynamics simulations performed at temperatures ranging from 298 to 550 K, starting from the SH3 native structure. We apply the all-atom CHARMM22 force field with different modified analytical generalized Born solvent models (GBSW and GBMV2) and compare these simulation models with the distance-dependent dielectric screening of charge-charge interactions. A further comparison is provided with the united-atom CHARMM19 plus a pairwise GB model. Unfolding-folding transition temperatures of SH3 were estimated from the temperature-dependent profiles of the heat capacity, root-mean-square distance from the native structure, and the fraction of native contacts, each calculated from the density of states by using the weighted histogram analysis method. We observed that, for CHARMM22, the unfolding transition and energy probability density were quite sensitive to the implicit solvent description, in particular, the treatment of the protein-solvent dielectric boundary in GB models and their surface-area-based hydrophobic term. Among the solvent models tested, the calculated melting temperature varied in the range 353-438 K and was higher than the experimental value near 340 K. A reformulated GBMV2 model of employing a smoother molecular-volume dielectric interface was the most accurate in reproducing the native conformation and a two-state folding landscape, although the melting transition temperature did not show the smallest deviation from experiment. For the lower-resolution CHARMM19/GB model, the simulations failed to yield a bimodal energy distribution, yet the melting temperature was observed to be a good estimate of higher-resolution simulation models. We also demonstrate that a careful analysis of a relatively long simulation is necessary to avoid trapping in local minima and to find a true thermodynamic transition temperature.

Determination of the structure of the energy landscape for coarse-grained off-lattice models of folding heteropolymers

The structure of energy landscapes for a minimalist coarse-grained off-lattice protein model is presented to investigate the folding behaviors of heteropolymers. The obtained energy landscapes serve as a useful tool for visualization of the funnel-like structure of a considered system in the configuration space. Despite the simplicity of the model, the knowledge of the free-energy landscape enables us to show different folding characteristics known from real proteins and synthetic peptides, such as two-state folding and metastability.

Using an Amino Acid Fluorescence Resonance Energy Transfer Pair To Probe Protein Unfolding: Application to the Villin Headpiece Subdomain and the LysM Domain

Previously, we have shown that p-cyanophenylalanine (Phe CN) and tryptophan (Trp) constitute an efficient fluorescence resonance energy transfer (FRET) pair that has several advantages over commonly used dye pairs. Here, we aim to examine the general applicability of this FRET pair in protein folding-unfolding studies by applying it to the urea-induced unfolding transitions of two small proteins, the villin headpiece subdomain (HP35) and the lysin motif (LysM) domain. Depending on whether Phe CN is exposed to solvent, we are able to extract either qualitative information about the folding pathway, as demonstrated by HP35, which has been suggested to unfold in a stepwise manner, or quantitative thermodynamic and structural information, as demonstrated by LysM, which has been shown to be an ideal two-state folder. Our results show that the unfolding transition of HP35 reported by FRET occurs at a denaturant concentration lower than that measured by circular dichroism (CD) and that the loop linking helix 2 and helix 3 remains compact in the denatured state, which are consistent with the notion that HP35 unfolds in discrete steps and that its unfolded state contains residual structures. On the other hand, our FRET results on the LysM domain allow us to develop a model for extracting structural and thermodynamic parameters about its unfolding, and we find that our results are in agreement with those obtained by other methods. Given the fact that Phe CN is a non-natural amino acid and, thus, amenable to incorporation into peptides and proteins via existing peptide synthesis and protein expression methods, we believe that the FRET method demonstrated here is widely applicable to protein conformational studies, especially to the study of relatively small proteins.

The Unfolded Ensemble and Folding Mechanism of the C-Terminal GB1 β-Hairpin

In this work, we shed new light on a much-studied case of beta-hairpin folding by means of advanced molecular dynamics simulations. A fully atomistic description of the protein and the solvent molecule is used, together with metadynamics, to accelerate the sampling and estimate free-energy landscapes. This is achieved using the path collective variables approach, which provides an adaptive description of the mechanism under study. We discover that the folding mechanism is a multiscale process where the turn region conformation leads to two different energy pathways that are connected by elongated structures. The former displays a stable 2:4 native-like structure in which an optimal hydrophobic packing and hydrogen bond pattern leads to 8 kcal/mol of stabilization. The latter shows a less-structured 3:5 beta-sheet, where hydrogen bonds and hydrophobic packing provide only 2.5 kcal/mol of stability. This perspective is fully consistent with experimental evidence that shows this to be a prototypical two-state folder, while it redefines the nature of the unfolded state.

Peptide Folding Using Multiscale Coarse-Grained Models

The multiscale coarse-graining (MS-CG) method has been previously used to describe the equilibrium properties of peptides. The present study reveals that MS-CG models of alpha-helical polyalanine and the beta-hairpin V 5PGV 5 possess the capacity to efficiently refold in simulations initiated from unfolded configurations. The MS-CG peptides exhibit free energy landscapes that are funneled toward folded configurations and two-state folding behavior, consistent with the known characteristics of small, rapidly folding peptides. Moreover, the models demonstrate enhanced sampling capabilities when compared to systems with full atomic detail. The significance of these observations with respect to the theoretical basis of the MS-CG approach is discussed. The MS-CG peptides were used to reconstruct atomically detailed configurations in order to evaluate the extent to which MS-CG ensembles embody all-atom peptide free energy landscapes. Ensembles obtained from these reconstructed configurations display good agreement with the all-atom simulation data used to generate the MS-CG models and also corroborate the presence of features observed in the MS-CG peptide free energy landscapes. These findings suggest that MS-CG models may be of significant utility in the study of peptide folding.

ProbeCheck – a central resource for evaluating oligonucleotide probe coverage and specificity

The web server probeCheck, freely accessible at http://www.microbial-ecology.net/probecheck, provides a pivotal forum for rapid specificity and coverage evaluations of probes and primers against selected databases of phylogenetic and functional marker genes. Currently, 24 widely used sequence collections including the Ribosomal Database Project (RDP) II, Greengenes, SILVA and the Functional Gene Pipeline/Repository can be queried. For this purpose, probeCheck integrates a new online version of the popular ARB probe match tool with free energy (ΔG) calculations for each perfectly matched and mismatched probe-target hybrid, allowing assessment of the theoretical binding stabilities of oligo-target and non-target hybrids. For each output sequence, the accession number, the GenBank taxonomy and a link to the respective entry at GenBank, EMBL and, if applicable, the query database are displayed. Filtering options allow customizing results on the output page. In addition, probeCheck is linked with probe match tools of RDP II and Greengenes, NCBI blast, the Oligonucleotide Properties Calculator, the two-state folding tool of the DINAMelt server and the rRNA-targeted probe database probeBase. Taken together, these features provide a multifunctional platform with maximal flexibility for the user in the choice of databases and options for the evaluation of published and newly developed probes and primers.

Thermodynamics, Kinetics, and Salt dependence of Folding of YopM, a Large Leucine-rich Repeat Protein

Small globular proteins have many contacts between residues that are distant in primary sequence. These contacts create a complex network between sequence-distant segments of secondary structure, which may be expected to promote the cooperative folding of globular proteins. Although repeat proteins, which are composed of tandem modular units, lack sequence-distant contacts, several of considerable length have been shown to undergo cooperative two-state folding. To explore the limits of cooperativity in repeat proteins, we have studied the unfolding of YopM, a leucine-rich repeat (LRR) protein of over 400 residues. Despite its large size and modular architecture (15 repeats), YopM equilibrium unfolding is highly cooperative, and shows a very strong dependence on the concentration of urea. In contrast, kinetic studies of YopM folding indicate a mechanism that includes one or more transient intermediates. The urea dependence of the folding and unfolding rates suggests a relatively small transition state ensemble. As with the urea dependence, we have found an extreme dependence of the free energy of unfolding on the concentration of salt. This salt dependence likely results from general screening of a large number of unfavorable columbic interactions in the folded state, rather than from specific cation binding.

Predicting protein folding rates from geometric contact and amino acid sequence

Protein folding speeds are known to vary over more than eight orders of magnitude. Plaxco, Simons, and Baker (see References) first showed a correlation of folding speed with the topology of the native protein. That and subsequent studies showed, if the native structure of a protein is known, its folding speed can be predicted reasonably well through a correlation with the "localness" of the contacts in the protein. In the present work, we develop a related measure, the geometric contact number, N (alpha), which is the number of nonlocal contacts that are well-packed, by a Voronoi criterion. We find, first, that in 80 proteins, the largest such database of proteins yet studied, N (alpha) is a consistently excellent predictor of folding speeds of both two-state fast folders and more complex multistate folders. Second, we show that folding rates can also be predicted from amino acid sequences directly, without the need to know the native topology or other structural properties.

On the Pathway of Forming Enzymatically Productive Ligand-Protein Complexes in Lactate Dehydrogenase

We have carried out a series of studies on the binding of a substrate mimic to the enzyme lactate dehydrogenase (LDH) using advanced kinetic approaches, which begin to provide a molecular picture of the dynamics of ligand binding for this protein. Binding proceeds via a binding-competent subpopulation of the nonligated form of the protein (the LDH/NADH binary complex) to form a protein-ligand encounter complex. The work here describes the collapse of the encounter complex to form the catalytically competent Michaelis complex. Isotope-edited static Fourier transform infrared studies on the bound oxamate protein complex reveal two kinds of oxamate environments: 1), a major populated structure wherein all significant hydrogen-bonding patterns are formed at the active site between protein and bound ligand necessary for the catalytically productive Michaelis complex and 2), a minor structure in a configuration of the active site that is unfavorable to carry out catalyzed chemistry. This latter structure likely simulates a dead-end complex in the reaction mixture. Temperature jump isotope-edited transient infrared studies on the binding of oxamate with LDH/NADH suggest that the evolution of the encounter complex between LDH/NADH and oxamate collapses via a branched reaction pathway to form the major and minor bound species. The production of the catalytically competent protein-substrate complex has strong similarities to kinetic pathways found in two-state protein folding processes. Once the encounter complex is formed between LDH/NADH and substrate, the ternary protein-ligand complex appears to “fold” to form a compact productive complex in an all or nothing like fashion with all the important molecular interactions coming together at the same time.

Folding kinetics and thermodynamics of Pseudomonas syringae effector protein AvrPto provide insight into translocation via the type III secretion system

In order to infect their hosts, many Gram-negative bacteria translocate agents of infection, called effector proteins, through the type III secretion system (TTSS) into the host cytoplasm. This process is thought to require at least partial unfolding of these agents, raising the question of how an effector protein might unfold to enable its translocation and then refold once it reaches the host cytoplasm. AvrPto is a well-studied effector protein of Pseudomonas syringae pv tomato. The presence of a readily observed unfolded population of AvrPto in aqueous solution and the lack of a known secretion chaperone make it ideal for studying the kinetic and thermodynamic characteristics that facilitate translocation. Application of Nzz exchange spectroscopy revealed a global, two-state folding equilibrium with 16% unfolded population, a folding rate of 1.8 s(-1), and an unfolding rate of 0.33 s(-1) at pH 6.1. TrAvrPto stability increases with increasing pH, with only 2% unfolded population observed at pH 7.0. The R(1) relaxation of TrAvrPto, which is sensitive to both the global anisotropy of folded TrAvrPto and slow exchange between folded and unfolded conformations, provided independent verification of the global kinetic rate constants. Given the acidic apoplast in which the pathogen resides and the more basic host cytoplasm, these results offer an intriguing mechanism by which the pH dependence of stability and slow folding kinetics of AvrPto would allow efficient translocation of the unfolded form through the TTSS and refolding into its functional folded form once inside the host.