A reexamination of the folding mechanism of dihydrofolate reductase from Escherichia coli: verification and refinement of a four-channel model

How synonymous mutations alter enzyme structure and function over long timescales

The specific activity of enzymes can be altered over long timescales in cells by synonymous mutations that alter a messenger RNA molecule's sequence but not the encoded protein's primary structure. How this happens at the molecular level is unknown. Here, we use multiscale modelling of three Escherichia coli enzymes (type III chloramphenicol acetyltransferase, D-alanine-D-alanine ligase B and dihydrofolate reductase) to understand experimentally measured changes in specific activity due to synonymous mutations. The modelling involves coarse-grained simulations of protein synthesis and post-translational behaviour, all-atom simulations to test robustness and quantum mechanics/molecular mechanics calculations to characterize enzymatic function. We show that changes in codon translation rates induced by synonymous mutations cause shifts in co-translational and post-translational folding pathways that kinetically partition molecules into subpopulations that very slowly interconvert to the native, functional state. Structurally, these states resemble the native state, with localized misfolding near the active sites of the enzymes. These long-lived states exhibit reduced catalytic activity, as shown by their increased activation energies for the reactions they catalyse.

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Fast Protein Translation Can Promote Co- and Posttranslational Folding of Misfolding-Prone Proteins

Chemical kinetic modeling has previously been used to predict that fast-translating codons can enhance cotranslational protein folding by helping to avoid misfolded intermediates. Consistent with this prediction, protein aggregation in yeast and worms was observed to increase when translation was globally slowed down, possibly due to increased cotranslational misfolding. Observation of similar behavior in molecular simulations would confirm predictions from the simpler chemical kinetic model and provide a molecular perspective on cotranslational folding, misfolding, and the impact of translation speed on these processes. All-atom simulations cannot reach the timescales relevant to protein synthesis, and most conventional structure-based coarse-grained models do not allow for nonnative structure formation. Here, we introduce a protocol to incorporate misfolding using the functional forms of publicly available force fields. With this model we create two artificial proteins that are capable of undergoing structural transitions between a native and a misfolded conformation and simulate their synthesis by the ribosome. Consistent with the chemical kinetic predictions, we find that rapid synthesis of misfolding-prone nascent-chain segments increases the fraction of folded proteins by kinetically partitioning more molecules through on-pathway intermediates, decreasing the likelihood of sampling misfolded conformations. Novel to this study, to our knowledge, we observe that differences in protein dynamics, arising from different translation-elongation schedules, can persist long after the nascent protein has been released from the ribosome, and that a sufficient level of energetic frustration is needed for fast-translating codons to be beneficial for folding. These results provide further evidence that fast-translating codons can be as biologically important as pause sites in coordinating cotranslational folding.

Halophilic mechanism of the enzymatic function of a moderately halophilic dihydrofolate reductase from Haloarcula japonica strain TR-1

Dihydrofolate (DHF) reductase coded by a plasmid of the extremely halophilic archaeon Haloarcula japonica strain TR-1 (HjDHFR P1) shows moderate halophilicity on enzymatic activity at pH 6.0, although there is no significant effect of NaCl on its secondary structure. To elucidate the salt-activation and -inactivation mechanisms of this enzyme, we investigated the effects of pH and salt concentration, deuterium isotope effect, steady-state kinetics, and rapid-phase ligand-binding kinetics. Enzyme activity was increased eightfold by the addition of 500 mM NaCl at pH 6.0, fourfold by 250 mM at pH 8.0, and became independent of salt concentration at pH 10.0. Full isotope effects observed at pH 10.0 under 0-1000 mM NaCl indicated that the rate of hydride transfer, which was the rate-determining step at the basic pH region, was independent of salt concentration. Conversely, rapid-phase ligand-binding experiments showed that the amplitude of the DHF-binding reaction increased and the tetrahydrofolate (THF)-releasing rate decreased with increasing NaCl concentration. These results suggested that the salt-activation mechanism of HjDHFR P1 is via the population change of the anion-unbound and anion-bound conformers, which are binding-incompetent and -competent conformations for DHF, respectively, while that of salt inactivation is via deceleration of the THF-releasing rate, which is the rate-determining step at the neutral pH region.

Homodimeric Escherichia coli Toxin CcdB (Controller of Cell Division or Death B Protein) Folds via Parallel Pathways

The existence of parallel pathways in the folding of proteins seems intuitive, yet remains controversial. We explore the folding kinetics of the homodimeric Escherichia coli toxin CcdB (Controller of Cell Division or Death B protein) using multiple optical probes and approaches. Kinetic studies performed as a function of protein and denaturant concentrations demonstrate that the folding of CcdB is a four-state process. The two intermediates populated during folding are present on parallel pathways. Both form by rapid association of the monomers in a diffusion limited manner and appear to be largely unstructured, as they are silent to the optical probes employed in the current study. The existence of parallel pathways is supported by the insensitivity of the amplitudes of the refolding kinetic phases to the different probes used in the study. More importantly, interrupted refolding studies and ligand binding studies clearly demonstrate that the native state forms in a biexponential manner, implying the presence of at least two pathways. Our studies indicate that the CcdA antitoxin binds only to the folded CcdB dimer and not to any earlier folding intermediates. Thus, despite being part of the same operon, the antitoxin does not appear to modulate the folding pathway of the toxin encoded by the downstream cistron. This study highlights the utility of ligand binding in distinguishing between sequential and parallel pathways in protein folding studies, while also providing insights into molecular interactions during folding in Type II toxin-antitoxin systems.

Direct observation of parallel folding pathways revealed using a symmetric repeat protein system

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.

Structure of a partially unfolded form of Escherichia coli dihydrofolate reductase provides insight into its folding pathway

Proteins frequently fold via folding intermediates that correspond to local minima on the conformational energy landscape. Probing the structure of the partially unfolded forms in equilibrium under native conditions can provide insight into the properties of folding intermediates. To elucidate the structures of folding intermediates of Escherichia coli dihydrofolate reductase (DHFR), we investigated transient partial unfolding of DHFR under native conditions. We probed the structure of a high-energy conformation susceptible to proteolysis (cleavable form) using native-state proteolysis. The free energy for unfolding to the cleavable form is clearly less than that for global unfolding. The dependence of the free energy on urea concentration (m-value) also confirmed that the cleavable form is a partially unfolded form. By assessing the effect of mutations on the stability of the partially unfolded form, we found that native contacts in a hydrophobic cluster formed by the F-G and Met-20 loops on one face of the central β-sheet are mostly lost in the partially unfolded form. Also, the folded region of the partially unfolded form is likely to have some degree of structural heterogeneity. The structure of the partially unfolded form is fully consistent with spectroscopic properties of the near-native kinetic intermediate observed in previous folding studies of DHFR. The findings suggest that the last step of the folding of DHFR involves organization in the structure of two large loops, the F-G and Met-20 loops, which is coupled with compaction of the rest of the protein.

Folding pathway of a multidomain protein depends on its topology of domain connectivity

How do the folding mechanisms of multidomain proteins depend on protein topology? We addressed this question by developing an Ising-like structure-based model and applying it for the analysis of free-energy landscapes and folding kinetics of an example protein, Escherichia coli dihydrofolate reductase (DHFR). DHFR has two domains, one comprising discontinuous N- and C-terminal parts and the other comprising a continuous middle part of the chain. The simulated folding pathway of DHFR is a sequential process during which the continuous domain folds first, followed by the discontinuous domain, thereby avoiding the rapid decrease in conformation entropy caused by the association of the N- and C-terminal parts during the early phase of folding. Our simulated results consistently explain the observed experimental data on folding kinetics and predict an off-pathway structural fluctuation at equilibrium. For a circular permutant for which the topological complexity of wild-type DHFR is resolved, the balance between energy and entropy is modulated, resulting in the coexistence of the two folding pathways. This coexistence of pathways should account for the experimentally observed complex folding behavior of the circular permutant.

Autotransporters: The Cellular Environment Reshapes a Folding Mechanism to Promote Protein Transport

We know very little about how the cellular environment affects protein folding mechanisms. Here, we focus on one unique aspect of that environment that is difficult to recapitulate in the test tube: the effect of a folding vector. When protein folding is initiated at one end of the polypeptide chain, folding starts from a much smaller ensemble of conformations than during refolding of a full-length polypeptide chain. But to what extent can vectorial folding affect protein folding kinetics and the conformations of folding intermediates? We focus on recent studies of autotransporter proteins, the largest class of virulence proteins from pathogenic Gram-negative bacteria. Autotransporter proteins are secreted across the bacterial inner membrane from N→C-terminus, which, like refolding in vitro, retards folding. But in contrast, upon C→N-terminal secretion across the outer membrane autotransporter folding proceeds orders of magnitude faster. The potential impact of vectorial folding on the folding mechanisms of other proteins is also discussed.

ATP-independent control of autotransporter virulence protein transport via the folding properties of the secreted protein

Autotransporter (AT) proteins are the largest class of extracellular virulence proteins secreted from Gram-negative bacteria. The mechanism by which AT proteins cross the bacterial outer membrane (OM), in the absence of ATP or another external energy source, is unknown. Here we demonstrate a linear correlation between localized regions of stability (ΔG(folding)) in the mature virulence protein (the AT "passenger") and OM secretion efficiency. Destabilizing the C-terminal β-helical domain of a passenger reduced secretion efficiency. In contrast, destabilizing the globular N-terminal domain of a passenger produced a linearly correlated increase in secretion efficiency. Thus, C-terminal passenger stability facilitates OM secretion, whereas N-terminal stability hinders it. The contributions of regional passenger stability to OM secretion demonstrate a crucial role for the passenger itself in directing its secretion, suggesting a novel type of ATP-independent, folding-driven transporter.

Geofold: topology-based protein unfolding pathways capture the effects of engineered disulfides on kinetic stability

Protein unfolding is modeled as an ensemble of pathways, where each step in each pathway is the addition of one topologically possible conformational degree of freedom. Starting with a known protein structure, GeoFold hierarchically partitions (cuts) the native structure into substructures using revolute joints and translations. The energy of each cut and its activation barrier are calculated using buried solvent accessible surface area, side chain entropy, hydrogen bonding, buried cavities, and backbone degrees of freedom. A directed acyclic graph is constructed from the cuts, representing a network of simultaneous equilibria. Finite difference simulations on this graph simulate native unfolding pathways. Experimentally observed changes in the unfolding rates for disulfide mutants of barnase, T4 lysozyme, dihydrofolate reductase, and factor for inversion stimulation were qualitatively reproduced in these simulations. Detailed unfolding pathways for each case explain the effects of changes in the chain topology on the folding energy landscape. GeoFold is a useful tool for the inference of the effects of disulfide engineering on the energy landscape of protein unfolding.

The shape-shifting quasispecies of RNA: one sequence, many functional folds

E Unus pluribum, or "Of One, Many", may be at the root of decoding the RNA sequence-structure-function relationship. RNAs embody the large majority of genes in higher eukaryotes and fold in a sequence-directed fashion into three-dimensional structures that perform functions conserved across all cellular life forms, ranging from regulating to executing gene expression. While it is the most important determinant of the RNA structure, the nucleotide sequence is generally not sufficient to specify a unique set of secondary and tertiary interactions due to the highly frustrated nature of RNA folding. This frustration results in folding heterogeneity, a common phenomenon wherein a chemically homogeneous population of RNA molecules folds into multiple stable structures. Often, these alternative conformations constitute misfolds, lacking the biological activity of the natively folded RNA. Intriguingly, a number of RNAs have recently been described as capable of adopting multiple distinct conformations that all perform, or contribute to, the same function. Characteristically, these conformations interconvert slowly on the experimental timescale, suggesting that they should be regarded as distinct native states. We discuss how rugged folding free energy landscapes give rise to multiple native states in the Tetrahymena Group I intron ribozyme, hairpin ribozyme, sarcin-ricin loop, ribosome, and an in vitro selected aptamer. We further describe the varying degrees to which folding heterogeneity impacts function in these RNAs, and compare and contrast this impact with that of heterogeneities found in protein folding. Embracing that one sequence can give rise to multiple native folds, we hypothesize that this phenomenon imparts adaptive advantages on any functionally evolving RNA quasispecies.

Multiple native states reveal persistent ruggedness of an RNA folding landscape

According to the “thermodynamic hypothesis”, the sequence of a biological macromolecule defines its folded, active structure as a global energy minimum on the folding landscape.1,2 But the enormous complexity of folding landscapes of large macromolecules raises a question: Is there indeed a unique global energy minimum corresponding to a unique native conformation, or are there deep local minima corresponding to alternative active conformations?3 Folding of many proteins is well described by two-state models, leading to highly simplified representations of protein folding landscapes with a single native conformation.4,5 Nevertheless, accumulating experimental evidence suggests a more complex topology of folding landscapes with multiple active conformations that can take seconds or longer to interconvert.6,7,8 Here we employ single molecule experiments to demonstrate that an RNA enzyme folds into multiple distinct native states that interconvert much slower than the time scale of catalysis. These data demonstrate that the severe ruggedness of RNA folding landscapes extends into conformational space occupied by native conformations.

Multiple routes and structural heterogeneity in protein folding

Experimental studies show that many proteins fold along sequential pathways defined by folding intermediates. An intermediate may not always be a single population of molecules but may consist of subpopulations that differ in their average structure. These subpopulations are likely to fold via independent pathways. Parallel folding and unfolding pathways appear to arise because of structural heterogeneity. For some proteins, the folding pathways can effectively switch either because different subpopulations of an intermediate get populated under different folding conditions, or because intermediates on otherwise hidden pathways get stabilized, leading to their utilization becoming discernible, or because mutations stabilize different substructures. Therefore, the same protein may fold via different pathways in different folding conditions. Multiple folding pathways make folding robust, and evolution is likely to have selected for this robustness to ensure that a protein will fold under the varying conditions prevalent in different cellular contexts.

Kinetic folding of Haloferax volcanii and Escherichia coli dihydrofolate reductases: haloadaptation by unfolded state destabilization at high ionic strength

Salts affect protein stability by multiple mechanisms (e.g., the Hofmeister effect, preferential hydration, electrostatic effects and weak ion binding). These mechanisms can affect the stability of both the native state and the unfolded state. Previous equilibrium stability studies demonstrated that KCl stabilizes dihydrofolate reductases (DHFRs) from Escherichia coli (ecDHFR, E. coli DHFR) and Haloferax volcanii (hvDHFR1, H. volcanii DHFR encoded by the hdrA gene) with similar efficacies, despite adaptation to disparate physiological ionic strengths (0.2 M versus 2 M). Kinetic studies can provide insights on whether equilibrium effects reflect native state stabilization or unfolded state destabilization. Similar kinetic mechanisms describe the folding of urea-denatured ecDHFR and hvDHFR1: a 5-ms stopped-flow burst-phase species that folds to the native state through two sequential intermediates with relaxation times of 0.1-3 s and 25-100 s. The latter kinetic step is very similar to that observed for the refolding of hvDHFR1 from low ionic strength. The unfolding of hvDHFR1 at low ionic strength is relatively slow, suggesting kinetic stabilization as observed for some thermophilic enzymes. Increased KCl concentrations slow the urea-induced unfolding of ecDHFR and hvDHFR1, but much less than expected from equilibrium studies. Unfolding rates extrapolated to 0 M denaturant, k(unf)(H(2)O), are relatively independent of ionic strength, demonstrating that the KCl-induced stabilization of ecDHFR and hvDHFR1 results predominantly from destabilization of the unfolded state. This supports the hypothesis from previous equilibrium studies that haloadaptation harnesses the effects of elevated salt concentrations on the properties of the aqueous solvent to enhance protein stability.

Microsecond Hydrophobic Collapse in the Folding of Escherichia coli Dihydrofolate Reductase, an α/β-Type Protein

Using small-angle X-ray scattering combined with a continuous-flow mixing device, we monitored the microsecond compaction dynamics in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein. A significant collapse of the radius of gyration from 30 A to 23.2 A occurs within 300 micros after the initiation of refolding by a urea dilution jump. The subsequent folding after the major chain collapse occurs on a considerably longer time-scale. The protein folding trajectories constructed by comparing the development of the compactness and the secondary structure suggest that the specific hydrophobic collapse model rather than the framework model better explains the experimental observations. The folding trajectory of this alpha/beta-type protein is located between those of alpha-helical and beta-sheet proteins, suggesting that native structure determines the folding landscape.

Stabilization of Hyperactive Dihydrofolate Reductase by Cyanocysteine-mediated Backbone Cyclization

Stabilization of an enzyme while maintaining its activity has been a major challenge in protein chemistry. Although it is difficult to simultaneously improve stability and activity of a protein by amino acid substitutions due to the activity-stability trade-off, backbone cyclization by connecting the N and C termini with a linker is promising as a general method of stabilizing a protein without affecting its activity. Recently, we created a hyperactive, methionine- and cysteine-free mutant of dihydrofolate reductase from Escherichia coli, called ANLYF, by introducing seven amino acid substitutions, which, however, destabilized the protein. Here we show that ANLYF is stabilized without a loss of its high activity by a novel backbone cyclization method for unprotected proteins. The method is based on the in vitro cyanocysteine-mediated intramolecular ligation reaction, which can be conducted with relatively high efficiency by a simple procedure and under mild conditions. We also show that the reversibility of thermal denaturation is highly improved by the cyclization. Thus, activity and stability of the protein can be separately improved by amino acid substitutions and backbone cyclization, respectively. We suggest that the cyanocysteine-mediated cyclization method is complementary to the intein-mediated cyclization method in stabilizing a protein without affecting its activity.

Criteria for downhill protein folding: Calorimetry, chevron plot, kinetic relaxation, and single-molecule radius of gyration in chain models with subdued degrees of cooperativity

Recent investigations of possible downhill folding of small proteins such as BBL have focused on the thermodynamics of non-two-state, "barrierless" folding/denaturation transitions. Downhill folding is noncooperative and thermodynamically "one-state," a phenomenon underpinned by a unimodal conformational distribution over chain properties such as enthalpy, hydrophobic exposure, and conformational dimension. In contrast, corresponding distributions for cooperative two-state folding are bimodal with well-separated population peaks. Using simplified atomic modeling of a three-helix bundle-in a scheme that accounts for hydrophobic interactions and hydrogen bonding-and coarse-grained C(alpha) models of four real proteins with various degrees of cooperativity, we evaluate the effectiveness of several observables at defining the underlying distribution. Bimodal distributions generally lead to sharper transitions, with a higher heat capacity peak at the transition midpoint, compared with unimodal distributions. However, the observation of a sigmoidal transition is not a reliable criterion for two-state behavior, and the heat capacity baselines, used to determine the van't Hoff and calorimetric enthalpies of the transition, can introduce ambiguity. Interestingly we find that, if the distribution of the single-molecule radius of gyration were available, it would permit discrimination between unimodal and bimodal underlying distributions. We investigate kinetic implications of thermodynamic noncooperativity using Langevin dynamics. Despite substantial chevron rollovers, the relaxation of the models considered is essentially single-exponential over an extended range of native stabilities. Consistent with experiments, significant deviations from single-exponential behavior occur only under strongly folding conditions.

Population of On-pathway Intermediates in the Folding of Ubiquitin

The role that intermediate states play in protein folding is the subject of intense investigation and in the case of ubiquitin has been controversial. We present fluorescence-detected kinetic data derived from single and double mixing stopped-flow experiments to show that the F45W mutant of ubiquitin (WT*), a well-studied single-domain protein and most recently regarded as a simple two-state system, folds via on-pathway intermediates. To account for the discrepancy we observe between equilibrium and kinetic stabilities and m-values, we show that the polypeptide chain undergoes rapid collapse to an intermediate whose presence we infer from a fast lag phase in interrupted refolding experiments. Double-jump kinetic experiments identify two direct folding phases that are not associated with slow isomerisation reactions in the unfolded state. These two phases are explained by kinetic partitioning which allows molecules to reach the native state from the collapsed state via two possible competing routes, which we further examine using two destabilised ubiquitin mutants. Interrupted refolding experiments allow us to observe the formation and decay of an intermediate along one of these pathways. A plausible model for the folding pathway of ubiquitin is presented that demonstrates that obligatory intermediates and/or chain collapse are important events in restricting the conformational search for the native state of ubiquitin.

Evolutional Design of a Hyperactive Cysteine- and Methionine-free Mutant of Escherichia coli Dihydrofolate Reductase

We developed a strategy for finding out the adapted variants of enzymes, and we applied it to an enzyme, dihydrofolate reductase (DHFR), in terms of its catalytic activity so that we successfully obtained several hyperactive cysteine- and methionine-free variants of DHFR in which all five methionyl and two cysteinyl residues were replaced by other amino acid residues. Among them, a variant (M1A/M16N/M20L/M42Y/C85A/M92F/C152S), named as ANLYF, has an approximately seven times higher k(cat) value than wild type DHFR. Enzyme kinetics and crystal structures of the variant were investigated for elucidating the mechanism of the hyperactivity. Steady-state and transient binding kinetics of the variant indicated that the kinetic scheme of the catalytic cycle of ANLYF was essentially the same as that of wild type, showing that the hyperactivity was brought about by an increase of the dissociation rate constants of tetrahydrofolate from the enzyme-NADPH-tetrahydrofolate ternary complex. The crystal structure of the variant, solved and refined to an R factor of 0.205 at 1.9-angstroms resolution, indicated that an increased structural flexibility of the variant and an increased size of the N-(p-aminobenzoyl)-L-glutamate binding cleft induced the increase of the dissociation constant. This was consistent with a large compressibility (volume fluctuation) of the variant. A comparison of folding kinetics between wild type and the variant showed that the folding of these two enzymes was similar to each other, suggesting that the activity enhancement of the enzyme can be attained without drastic changes of the folding mechanism.

Folding rate prediction using n-order contact distance for proteins with two- and three-state folding kinetics

It is a challenging task to understand the relationship between sequences and folding rates of proteins. Previous studies are found that one of contact order (CO), long-range order (LRO), total contact distance (TCD), chain topology parameter (CTP), and effective length (Leff) has a significant correlation with folding rate of proteins. In this paper, we introduce a new parameter called n-order contact distance (nOCD) and use it to predict folding rate of proteins with two- and three-state folding kinetics. A good linear correlation between the folding rate logarithm lnkf and nOCD with n=1.2, alpha=0.6 is found for two-state folders (correlation coefficient is -0.809, P-value<0.0001) and n=2.8, alpha=1.5 for three-state folders (correlation coefficient is -0.816, P-value<0.0001). However, this correlation is completely absent for three-state folders with n=1.2, alpha=0.6 (correlation coefficient is 0.0943, P-value=0.661) and for two-state folders with n=2.8, alpha=1.5 (correlation coefficient is -0.235, P-value=0.2116). We also find that the average number of contacts per residue Pm in the interval of m for two-state folders is smaller than that for three-state folders. The probability distribution P(gamma) of residue having gamma pairs of contacts fits a Gaussian distribution for both two- and three-state folders. We observe that the correlations between square radius of gyration S2 and number of residues for two- and three-state folders are both good, and the correlation coefficient is 0.908 and 0.901, and the slope of the fitting line is 1.202 and 0.795, respectively. Maybe three-state folders are more compact than two-state folders. Comparisons with nTCD and nCTP are also made, and it is found that nOCD is the best one in folding rate prediction.

The relationship between chain connectivity and domain stability in the equilibrium and kinetic folding mechanisms of dihydrofolate reductase from E.coli

The role of domains in defining the equilibrium and kinetic folding properties of dihydrofolate reductase (DHFR) from Escherichia coli was probed by examining the thermodynamic and kinetic properties of a set of variants in which the chain connectivity in the discontinuous loop domain (DLD) and the adenosine-binding domain (ABD) was altered by permutation. To test the concept that chain cleavage can selectively destabilize the domain in which the N- and C-termini are resident, permutations were introduced at one position within the ABD, one within the DLD and one at a boundary between the domains. The results demonstrated that a continuous ABD is required for a stable thermal intermediate and a continuous DLD is required for a stable urea intermediate. The permutation at the domain interface had both a thermal and urea intermediate. Strikingly, the observable kinetic folding responses of all three permuted proteins were very similar to the wild-type protein. These results demonstrate a crucial role for stable domains in defining the energy surface for the equilibrium folding reaction of DHFR. If domain connectivity affects the kinetic mechanism, the effects must occur in the sub-millisecond time range.

Peptide fragment studies on the folding elements of dihydrofolate reductase from Escherichia coli

One of the necessary conditions for a protein to be foldable is the presence of a complete set of "folding elements" (FEs) that are short, contiguous peptide segments distributed over an amino acid sequence. The FE-assembly model of protein folding has been proposed, in which the FEs play a role in guiding structure formation through FE-FE interactions early in folding. However, two major issues remain to be clarified regarding the roles of the FEs in determining protein foldability. Are the FEs AFUs that can form nativelike structures in isolation? Is the presence of only the FEs without mutual connections a sufficient condition for a protein to be foldable? Here, we address these questions using peptide fragments corresponding to the FEs of DHFR from Escherichia coli. We show by CD measurement that the FE peptides are unfolded under the native conditions, and some of them have the propensities toward non-native helices. MD simulations also show the non-native helical propensities of the peptides, and the helix contents estimated from the simulations are well correlated with those estimated from the CD in TFE. Thus, the FEs of DHFR are not AFUs, suggesting the importance of the FEs in nonlocal interactions. We also show that equimolar mixtures of the FE peptides do not induce any structural formation. Therefore, mutual connections between the FEs, which should strengthen the nonlocal FE-FE interactions, are also one of the necessary conditions for a protein to be foldable.

Kinetic partitioning during folding of the p53 DNA binding domain

The DNA-binding domain (DBD) of wild-type p53 loses DNA binding activity spontaneously at 37 degrees C in vitro, despite being thermodynamically stable at this temperature. We test the hypothesis that this property is due to kinetic misfolding of DBD. Interrupted folding experiments and chevron analysis show that native molecules are formed via four tracks (a-d) under strongly native conditions. Folding half-lives of tracks a-d are 7.8 seconds, 50 seconds, 5.3 minutes and more than five hours, respectively, in 0.3M urea (10 degrees C). Approximately equal fractions of molecules fold through each track in zero denaturant, but above 2.0M urea approximately 90% fold via track c. A kinetic mechanism consisting of two parallel folding channels (fast and slow) is proposed. Each channel populates an on-pathway intermediate that can misfold to form an aggregation-prone, dead-end species. Track a represents direct folding through the fast channel. Track b proceeds through the fast channel but via the off-pathway state. Track c corresponds to folding via the slow channel, primarily through the off-pathway state. Track d proceeds by way of an even slower, uncharacterized route. We postulate that activity loss is caused by partitioning to the slower tracks, and that structural unfolding limits this process. In support of this view, tumorigenic hot-spot mutants G245S, R249S and R282Q accelerate unfolding rates but have no affect on folding kinetics. We suggest that these and other destabilizing mutants facilitate loss of p53 function by causing DBD to cycle unusually rapidly between folded and unfolded states. A significant fraction of DBD molecules become effectively trapped in a non-functional state with each unfolding-folding cycle.

Probing the interactions between the folding elements early in the folding of Escherichia coli dihydrofolate reductase by systematic sequence perturbation analysis

One of the necessary conditions for a protein to be foldable is the presence of a complete set of folding elements (FEs) that are short contiguous peptide segments distributed over an amino acid sequence. Previous studies indicated the FE assembly model of protein folding, in which the FEs interact with each other and coalesce to form an intermediate(s) early in the folding reaction. This suggests that a clue to the understanding of the determinants of protein foldability can be found by investigating how the FEs interact with each other early in the folding and thereby elucidating roles of the FEs in protein folding. To reveal the formation process of FE-FE interactions, we studied the early folding events of Escherichia coli dihydrofolate reductase (DHFR) utilizing systematic sequence perturbation analysis. Here, systematic single amino acid substitutions were introduced inside of the FEs (W30X in FE2, V40X in FE3, N59X in FE4, and I155X in FE10; X refers to various amino acid residues), and their kinetic refolding reactions were measured by stopped-flow circular dichroism and fluorescence. We show that the interactions around Trp30 and Ile155 are formed in the burst phase intermediate, while those around Val40 and Asn59 are formed in the transition state of the subsequent folding phase (tau5-phase) and in much later processes, respectively. These and previous results suggest that FE2 and FE10, and also FE1 and FE7, involved in the loop subdomain of DHFR, interact with each other within a millisecond time range, while the stable FE3-FE4 interactions are formed in the later processes. This may highlight the important roles of the FEs mainly inside of the loop subdomain in formation of the burst phase intermediate having a hydrophobic cluster and native-like overall topology and in acquisition of the foldability of DHFR.

SecA folding kinetics: a large dimeric protein rapidly forms multiple native states

SecA, a 202 kDa dimeric protein, is the ATPase for the Sec-dependent translocase of precursor proteins in vivo. SecA must undergo conformational changes, which may involve dissociation into a monomer, as it translocates the precursor protein across the inner membrane. To better understand the dynamics of SecA in vivo, protein folding studies to probe the native, intermediate, and unfolded species of SecA in vitro have been done. SecA folds through a stable dimeric intermediate and dimerizes in the dead-time of a manual-mixing kinetic experiment ( approximately 5-7 seconds). Here, stopped-flow fluorescence and CD, as well as ultra-rapid continuous flow fluorescence techniques, were used to further probe the rapid folding kinetics of SecA. In the absence of urea, rapid, near diffusion-limited ( approximately 10(9)M(-1)s(-1)) SecA dimerization occurs following a rate-limiting unimolecular rearrangement of a rapidly formed intermediate. Multiple kinetic folding and unfolding phases were observed and SecA was shown to have multiple native and unfolded states. Using sequential-mixing stopped-flow experiments, SecA was determined to fold via parallel channels with sequential intermediates. These results confirm that SecA is a highly dynamic protein, consistent with the rapid, major conformational changes it must undergo in vivo.

A concerted mechanism for the suppression of a folding defect through interactions with chaperones

Specific amino acid substitutions confer a temperature-sensitive-folding (tsf) phenotype to bacteriophage P22 coat protein. Additional amino acid substitutions, called suppressor substitutions (su), relieve the tsf phenotype. These su substitutions are proposed to increase the efficiency of procapsid assembly, favoring correct folding over improper aggregation. Our recent studies indicate that the molecular chaperones GroEL/ES are more effectively recruited in vivo for the folding of tsf:su coat proteins than their tsf parents. Here, the tsf:su coat proteins are studied with in vitro equilibrium and kinetic techniques to establish a molecular basis for suppression. The tsf:su coat proteins were monomeric, as determined by velocity sedimentation analytical ultracentrifugation. The stability of the tsf:su coat proteins was ascertained by equilibrium urea titrations, which were best described by a three-state folding model, N <--> I <--> U. The tsf:su coat proteins either had stabilized native or intermediate states as compared with their tsf coat protein parents. The kinetics of the I <--> U transition showed a decrease in the rate of unfolding and a small increase in the rate of refolding, thereby increasing the population of the intermediate state. The increased intermediate population may be the reason the tsf:su coat proteins are aggregation-prone and likely enhances GroEL-ES interactions. The N --> I unfolding rate was slower for the tsf:su proteins than their tsf coat parents, resulting in an increase in the native state population, which may allow more competent interactions with scaffolding protein, an assembly chaperone. Thus, the suppressor substitution likely improves folding in vivo through increased efficiency of coat protein-chaperone interactions.

Multiple parallel-pathway folding of proline-free Staphylococcal nuclease

When a protein exhibits complex kinetics of refolding, we often ascribe the complexity to slow isomerization events in the denatured protein, such as cis/trans isomerization of peptidyl prolyl bonds. Does the complex folding kinetics arise only from this well-known reason? Here, we have investigated the refolding of a proline-free variant of staphylococcal nuclease by stopped-flow, double-jump techniques, to examine the folding reactions without the slow prolyl isomerizations. As a result, the protein folds into the native state along at least two accessible parallel pathways, starting from a macroscopically single denatured-state ensemble. The presence of intermediates on the individual folding pathways has revealed the existence of multiple parallel pathways, and is characterized by multi-exponential folding kinetics with a lag phase. Therefore, a "single" amino acid sequence can fold along the multiple parallel pathways. This observation in staphylococcal nuclease suggests that the multiple folding may be more general than we have expected, because the multiple parallel-pathway folding cannot be excluded from proteins that show simpler kinetics.

Contact order revisited: influence of protein size on the folding rate

Guided by the recent success of empirical model predicting the folding rates of small two-state folding proteins from the relative contact order (CO) of their native structures, by a theoretical model of protein folding that predicts that logarithm of the folding rate decreases with the protein chain length L as L(2/3), and by the finding that the folding rates of multistate folding proteins strongly correlate with their sizes and have very bad correlation with CO, we reexamined the dependence of folding rate on CO and L in attempt to find a structural parameter that determines folding rates for the totality of proteins. We show that the Abs_CO = CO x L, is able to predict rather accurately folding rates for both two-state and multistate folding proteins, as well as short peptides, and that this Abs_CO scales with the protein chain length as L(0.70 +/- 0.07) for the totality of studied single-domain proteins and peptides.

Folding of phage P22 coat protein monomers: kinetic and thermodynamic properties

To assemble into a virus with icosahedral symmetry, capsid proteins must be able to attain multiple conformations. Whether this conformational diversity is achieved during folding of the subunit, or subsequently during assembly, is not clear. Phage P22 coat protein offers an ideal model to investigate the folding of a monomeric capsid subunit since its folding is independent of assembly. Our early studies indicated that P22 coat protein monomers could be folded into an assembly-competent state in vitro, with evidence of a kinetic intermediate. Using urea denaturation, coat protein monomers are shown to be marginally stable. The reversible folding of coat protein follows a three-state model, N if I if U, with an intermediate exhibiting most of the tryptophan fluorescence of the folded state, but little secondary structure. Folding and unfolding kinetics monitored by circular dichroism, tryptophan fluorescence, and bisANS fluorescence indicate that several kinetic intermediates are populated sequentially through parallel channels en route to the native state. Additionally, two native states were identified, suggesting that the several conformers required to assemble an icosahedral capsid may be found in solution before assembly ensues.

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.

Chain length is the main determinant of the folding rate for proteins with three-state folding kinetics

We demonstrate that chain length is the main determinant of the folding rate for proteins with the three-state folding kinetics. The logarithm of their folding rate in water (k(f)) strongly anticorrelates with their chain length L (the correlation coefficient being -0.80). At the same time, the chain length has no correlation with the folding rate for two-state folding proteins (the correlation coefficient is -0.07). Another significant difference of these two groups of proteins is a strong anticorrelation between the folding rate and Baker's "relative contact order" for the two-state folders and the complete absence of such correlation for the three-state folders.

Testing the relationship between foldability and the early folding events of dihydrofolate reductase from Escherichia coli

A "folding element" is a contiguous peptide segment crucial for a protein to be foldable and is a new concept that could assist in our understanding of the protein-folding problem. It is known that the presence of the complete set of folding elements of dihydrofolate reductase (DHFR) from Escherichia coli is essential for the protein to be foldable. Since almost all of the amino acid residues known to be involved in the early folding events of DHFR are located within the folding elements, a close relationship between the folding elements and early folding events is hypothesized. In order to test this hypothesis, we have investigated whether or not the early folding events are preserved in circular permutants and topological mutants of DHFR, in which the order of the folding elements is changed but the complete set of folding elements is present. The stopped-flow circular dichroism (CD) measurements show that the CD spectra at the early stages of folding are similar among the mutants and the wild-type DHFR, indicating that the presence of the complete set of folding elements is sufficient to preserve the early folding events. We have further examined whether or not sequence perturbation on the folding elements by a single amino acid substitution affects the early folding events of DHFR. The results show that the amino acid substitutions inside of the folding elements can affect the burst-phase CD spectra, whereas the substitutions outside do not. Taken together, these results indicate that the above hypothesis is true, suggesting a close relationship between the foldability of a protein and the early folding events. We propose that the folding elements interact with each other and coalesce to form a productive intermediate(s) early in the folding, and these early folding events are important for a protein to be foldable.

The coordination of the isomerization of a conserved non-prolyl cis peptide bond with the rate-limiting steps in the folding of dihydrofolate reductase

The propensity for peptide bonds to adopt the trans configuration in native and unfolded proteins, and the relatively slow rates of cis-trans isomerization reactions, imply that the formation of cis peptide bonds in native conformations are likely to limit folding reactions. The role of the conserved cis Gly95-Gly96 peptide bond in dihydrofolate reductase (DHFR) from Escherichia coli was examined by replacing Gly95 with alanine. The introduction of a beta carbon at position 95 is expected to increase the propensity for the trans isomer and perturb the isomerization reaction required to reach the native conformation. Although G95A DHFR is 1.30 kcal mol(-1) less stable than the wild-type protein, it adopts a well-folded structure that can be chemically denatured in a cooperative fashion. The mutant protein also retains the complex refolding kinetic pattern attributed to a parallel-channel mechanism in wild-type DHFR. The spectroscopic response upon refolding monitored by Trp fluorescence and the absence of a Trp/Trp exciton coupling apparent in the far-UV CD spectrum of the wild-type protein, however, indicated that the tertiary structure of the folded state for G95A DHFR is altered. The addition of methotrexate (MTX), a tight-binding inhibitor, to folded G95A DHFR restored the exciton coupling and the fluorescence properties through five slow kinetic events whose relaxation times are independent of the ligand and the denaturant concentrations. The results were interpreted to mean that MTX-binding drives the formation of the cis isomer of the peptide bond between Ala95 and Gly96 in five compact and stable but not wild-type-like conformations that contain the trans isomer. Folding studies in the presence of MTX for both wild-type and G95A DHFR support the notion that the cis peptide bond between Gly95 and Gly96 in the wild-type protein forms during four parallel rate-limiting steps, which are primarily controlled by folding reactions, and lead directly to a set of native, or native-like, conformers. The isomerization of the cis peptide bond is not a source of the parallel channels that characterize the complex folding mechanism for DHFR.

Sequential vs. parallel protein-folding mechanisms: experimental tests for complex folding reactions

The recent emphasis on rough energy landscapes for protein folding reactions by theoreticians, and the many observations of complex folding kinetics by experimentalists provide a rationale for a brief literature survey of various empirical approaches for validating the underlying mechanisms. The determination of the folding mechanism is a key step in defining the energy surface on which the folding reactions occurs and in interpreting the effects of amino acid replacements on this reaction. Case studies that illustrate methods for differentiating between sequential and parallel channel folding mechanisms are presented. The ultimate goal of such efforts is to understand how the one-dimensional information contained in the amino acid sequence is rapidly and efficiently translated into three-dimensional structure.

The effect of salts on the activity and stability of Escherichia coli and Haloferax volcanii dihydrofolate reductases

The extremely halophilic Archae require near-saturating concentrations of salt in the external environment and in their cytoplasm, potassium being the predominant intracellular cation. The proteins of these organisms have evolved to function in concentrations of salt that inactivate or precipitate homologous proteins from non-halophilic species. It has been proposed that haloadaptation is primarily due to clustering of acidic residues on the surface of the protein, and that these clusters bind networks of hydrated ions. The dihydrofolate reductases from Escherichia coli (ecDHFR) and two DHFR isozymes from Haloferax volcanii (hvDHFR1 and hvDHFR2) have been used as a model system to compare the effect of salts on a mesophilic and halophilic enzyme. The KCl-dependence of the activity and substrate affinity was investigated. ecDHFR is largely inactivated above 1M KCl, with no major effect on substrate affinity. hvDHFR1 and hvDHFR2 unfold at KCl concentrations below approximately 0.5M. Above approximately 1M, the KCl dependence of the hvDHFR activities can be attributed to the effect of salt on substrate affinity. The abilities of NaCl, KCl, and CsCl to enhance the stability to urea denaturation were determined, and similar efficacies of stabilization were observed for all three DHFR variants. The DeltaG degrees (H(2)O) values increased linearly with increasing KCl and CsCl concentrations. The increase of DeltaG degrees (H(2)O) as a function of the smallest cation, NaCl, is slightly curved, suggesting a minor stabilization from cation binding or screening of electrostatic repulsion. At their respective physiological ionic strengths, the DHFR variants exhibit similar stabilities. Salts stabilize ecDHFR and the hvDHFRs by a common mechanism, not a halophile-specific mechanism, such as the binding of hydrated salt networks. The primary mode of salt stabilization of the mesophilic and halophilic DHFRs appears to be through preferential hydration and the Hofmeister effect of salt on the activity and entropy of the aqueous solvent. In support of this conclusion, all three DHFRs are similarly stabilized by the non-ionic cosolute, sucrose.

Analysis of kinetics using a hybrid maximum-entropy/nonlinear-least-squares method: application to protein folding

A hybrid analysis that combines the maximum entropy method (MEM) with nonlinear least squares (NLS) fitting has been developed to interpret a general time-dependent signal. Data that include processes of opposite sign and a slow baseline drift can be inverted to obtain both a continuous distribution of lifetimes and a sum of discrete exponentials. Fits by discrete exponentials are performed with initial parameters determined from the distribution of lifetimes obtained with the MEM. The regularization of the parameter space achieved by the MEM stabilizes the introduction of each successive exponential in the NLS fits. This hybrid approach is particularly useful when fitting by a large number of exponentials. Revision of the MEM "prior" based on features in the data can improve the lifetime distribution obtained. Standard errors in the mean are estimated automatically for raw data. The results presented for simulated data and for fluorescence measurements of protein folding illustrate the utility and accuracy of the hybrid algorithm. Analysis of the folding of dihydrofolate reductase reveals six kinetic processes, one more than previously reported.

Thermal unfolding molecular dynamics simulation of Escherichia coli dihydrofolate reductase: thermal stability of protein domains and unfolding pathway

Temperature induced unfolding of Escherichia coli dihydrofolate reductase was carried out by using molecular dynamic simulations. The simulations show that the unfolding generally involves an initial end-to-end collapse of the adenine binding domain into partially extended loops, followed by a gradual breakdown of the remaining beta sheet core structure. The core, which consists of beta strands 5-7, was observed to be the most resistant to thermal unfolding. This region, which is made up of part of the N terminus domain and part of the large domain of the E. coli dihydrofolate reductase, may constitute the nucleation site for protein folding and may be important for the eventual formation of both domains. The unfolding of different domains at different stages of the unfolding process suggests that protein domains vary in stability and that the rate at which they unfold can affect the overall outcome of the unfolding pathway. This observation is compared with the recently proposed hierarchical folding model. Finally, the results of the simulation were found to be consistent with a previous experimental study (Frieden, Proc Natl Acad Sci USA 1990;87:4413-4416) which showed that the folding process of E. coli dihydrofolate reductase involves sequential formation of the substrate binding sites.

Highly divergent dihydrofolate reductases conserve complex folding mechanisms

To test the hypothesis that protein folding mechanisms are better conserved than amino acid sequences, the mechanisms for dihydrofolate reductases (DHFR) from human (hs), Escherichia coli (ec) and Lactobacillus casei (lc) were elucidated and compared using intrinsic Trp fluorescence and fluorescence-detected 8-anilino-1-naphthalenesulfonate (ANS) binding. The development of the native state was monitored using either methotrexate (absorbance at 380 nm) or NADPH (extrinsic fluorescence) binding. All three homologs displayed complex unfolding and refolding kinetic mechanisms that involved partially folded states and multiple energy barriers. Although the pairwise sequence identities are less than 30 %, folding to the native state occurs via parallel folding channels and involves two types of on-pathway kinetic intermediates for all three homologs. The first ensemble of kinetic intermediates, detected within a few milliseconds, has significant secondary structure and exposed hydrophobic cores. The second ensemble is obligatory and has native-like side-chain packing in a hydrophobic core; however, these intermediates are unable to bind active-site ligands. The formation of the ensemble of native states occurs via three channels for hsDHFR, and four channels for lcDHFR and ecDHFR. The binding of active-site ligands (methotrexate and NADPH) accompanies the rate-limiting formation of the native ensemble. The conservation of the fast, intermediate and slow-folding events for this complex alpha/beta motif provides convincing evidence for the hypothesis that evolutionarily related proteins achieve the same fold via similar pathways.

Exploring the folding landscape of a structured RNA

Structured RNAs achieve their active states by traversing complex, multidimensional energetic landscapes. Here we probe the folding landscape of the Tetrahymena ribozyme by using a powerful approach: the folding of single ribozyme molecules is followed beginning from distinct regions of the folding landscape. The experiments, combined with small-angle x-ray scattering results, show that the landscape contains discrete folding pathways. These pathways are separated by large free-energy barriers that prevent interconversion between them, indicating that the pathways lie in deep channels in the folding landscape. Chemical protection and mutagenesis experiments are then used to elucidate the structural features that determine which folding pathway is followed. Strikingly, a specific long-range tertiary contact can either help folding or hinder folding, depending on when it is formed during the process. Together these results provide an unprecedented view of the topology of an RNA folding landscape and the RNA structural features that underlie this multidimensional landscape.

Rough energy landscapes in protein folding: dimeric E. coli Trp repressor folds through three parallel channels

The folding mechanism of the dimeric Escherichia coli Trp repressor (TR) is a kinetically complex process that involves three distinguishable stages of development. Following the formation of a partially folded, monomeric ensemble of species, within 5 ms, folding to the native dimer is controlled by three kinetic phases. The rate-limiting step in each phase is either a non-proline isomerization reaction or a dimerization reaction, depending on the final denaturant concentration. Two approaches have been employed to test the previously proposed folding mechanism of TR through three parallel channels: (1) unfolding double-jump experiments demonstrate that all three folding channels lead directly to native dimer; and (2) the differential stabilization of the transition state for the final step in folding and the native dimer, by the addition of salt, shows that all three channels involve isomerization of a dimeric species. A refined model for the folding of Trp repressor is presented, in which all three channels involve a rapid dimerization reaction between partially folded monomers followed by the isomerization of the dimeric intermediates to yield native dimer. The ensemble of partially folded monomers can be captured at equilibrium by low pH; one-dimensional proton NMR spectra at pH 2.5 demonstrate that monomers exist in two distinct, slowly interconverting conformations. These data provide a potential structural explanation for the three-channel folding mechanism of TR: random association of two different monomeric forms, which are distinguished by alternative packing modes of the core dimerization domain and the DNA-binding, helix-turn-helix, domain. One, perhaps both, of these packing modes contains non-native contacts.