Early intermediates in the folding of dihydrofolate reductase from Escherichia coli detected by hydrogen exchange and NMR

The importance of the location of the N-terminus in successful protein folding in vivo and in vitro

Protein folding in the cell often begins during translation. Many proteins fold more efficiently cotranslationally than when refolding from a denatured state. Changing the vectorial synthesis of the polypeptide chain through circular permutation could impact functional, soluble protein expression and interactions with cellular proteostasis factors. Here, we measure the solubility and function of every possible circular permutant (CP) of HaloTag in Escherichia coli cell lysate using a gel-based assay, and in living E. coli cells via FACS-seq. We find that 78% of HaloTag CPs retain protein function, though a subset of these proteins are also highly aggregation-prone. We examine the function of each CP in E. coli cells lacking the cotranslational chaperone trigger factor and the intracellular protease Lon and find no significant changes in function as a result of modifying the cellular proteostasis network. Finally, we biophysically characterize two topologically interesting CPs in vitro via circular dichroism and hydrogen-deuterium exchange coupled with mass spectrometry to reveal changes in global stability and folding kinetics with circular permutation. For CP33, we identify a change in the refolding intermediate as compared to wild-type (WT) HaloTag. Finally, we show that the strongest predictor of aggregation-prone expression in cells is the introduction of termini within the refolding intermediate. These results, in addition to our finding that termini insertion within the conformationally restrained core is most disruptive to protein function, indicate that successful folding of circular permutants may depend more on changes in folding pathway and termini insertion in flexible regions than on the availability of proteostasis factors.

Resolving chaperone-assisted protein folding on the ribosome at the peptide level

Protein folding in vivo begins during synthesis on the ribosome and is modulated by molecular chaperones that engage the nascent polypeptide. How these features of protein biogenesis influence the maturation pathway of nascent proteins is incompletely understood. Here, we use hydrogen-deuterium exchange mass spectrometry to define, at peptide resolution, the cotranslational chaperone-assisted folding pathway of Escherichia coli dihydrofolate reductase. The nascent polypeptide folds along an unanticipated pathway through structured intermediates not populated during refolding from denaturant. Association with the ribosome allows these intermediates to form, as otherwise destabilizing carboxy-terminal sequences remain confined in the ribosome exit tunnel. Trigger factor binds partially folded states without disrupting their structure, and the nascent chain is poised to complete folding immediately upon emergence of the C terminus from the exit tunnel. By mapping interactions between the nascent chain and ribosomal proteins, we trace the path of the emerging polypeptide during synthesis. Our work reveals new mechanisms by which cellular factors shape the conformational search for the native state.

Differentiating between the sequence of structural events on alternative pathways of folding of a heterodimeric protein

Distinguishing between competing pathways of folding of a protein, on the basis of how they differ in their progress of structure acquisition, remains an important challenge in protein folding studies. A previous study had shown that the heterodimeric protein, double chain monellin (dcMN) switches between alternative folding pathways upon a change in guanidine hydrochloride (GdnHCl) concentration. In the current study, the folding of dcMN has been characterized by the pulsed hydrogen exchange (HX) labeling methodology used in conjunction with mass spectrometry. Quantification of the extent to which folding intermediates accumulate and then disappear with time of folding at both low and high GdnHCl concentrations, where the folding pathways are known to be different, shows that the folding mechanism is describable by a triangular three-state mechanism. Structural characterization of the productive folding intermediates populated on the alternative pathways has enabled the pathways to be differentiated on the basis of the progress of structure acquisition that occurs on them. The intermediates on the two pathways differ in the extent to which the α-helix and the rest of the β-sheet have acquired structure that is protective against HX. The major difference is, however, that β2 has not acquired any protective structure in the intermediate formed on one pathway, but it has acquired significant protective structure in the intermediate formed on the alternative pathway. Hence, the sequence of structural events is different on the two alternative pathways.

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.

DMSO-Quenched H/D-Exchange 2D NMR Spectroscopy and Its Applications in Protein Science

Hydrogen/deuterium (H/D) exchange combined with two-dimensional (2D) NMR spectroscopy has been widely used for studying the structure, stability, and dynamics of proteins. When we apply the H/D-exchange method to investigate non-native states of proteins such as equilibrium and kinetic folding intermediates, H/D-exchange quenching techniques are indispensable, because the exchange reaction is usually too fast to follow by 2D NMR. In this article, we will describe the dimethylsulfoxide (DMSO)-quenched H/D-exchange method and its applications in protein science. In this method, the H/D-exchange buffer is replaced by an aprotic DMSO solution, which quenches the exchange reaction. We have improved the DMSO-quenched method by using spin desalting columns, which are used for medium exchange from the H/D-exchange buffer to the DMSO solution. This improvement has allowed us to monitor the H/D exchange of proteins at a high concentration of salts or denaturants. We describe methodological details of the improved DMSO-quenched method and present a case study using the improved method on the H/D-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride.

Structure elements can be predicted using the contact volume among protein residues

Previously, the structure elements of dihydrofolate reductase (DHFR) were determined using comprehen­sive Ala-insertion mutation analysis, which is assumed to be a kind of protein “building blocks.” It is hypo­thesized that our comprehension of the structure elements could lead to understanding how an amino acid sequence dictates its tertiary structure. However, the comprehensive Ala-insertion mutation analysis is a time- and cost-consuming process and only a set of the DHFR structure elements have been reported so far. Therefore, developing a computational method to predict structure elements is an urgent necessity. We focused on intramolecular residue–residue contacts to predict the structure elements. We introduced a simple and effective parameter: the overlapped contact volume (CV) among the residues and calculated the CV along the DHFR sequence using the crystal structure. Our results indicate that the CV profile can recapitulate its precipitate ratio profile, which was used to define the structure elements in the Ala-insertion mutation analysis. The CV profile allowed us to predict structure elements like the experimentally determined structure elements. The strong correlation between the CV and precipitate ratio profiles indicates the importance of the intramolecular residue–residue contact in maintaining the tertiary structure. Additionally, the CVs between the structure elements are considerably more than those between a structure element and a linker or two linkers, indicating that the structure elements play a funda­mental role in increasing the intramolecular adhesion. Thus, we propose that the structure elements can be considered a type of “building blocks” that maintain and dictate the tertiary structures of proteins.

Early Folding Events, Local Interactions, and Conservation of Protein Backbone Rigidity

Protein folding is in its early stages largely determined by the protein sequence and complex local interactions between amino acids, resulting in lower energy conformations that provide the context for further folding into the native state. We compiled a comprehensive data set of early folding residues based on pulsed labeling hydrogen deuterium exchange experiments. These early folding residues have corresponding higher backbone rigidity as predicted by DynaMine from sequence, an effect also present when accounting for the secondary structures in the folded protein. We then show that the amino acids involved in early folding events are not more conserved than others, but rather, early folding fragments and the secondary structure elements they are part of show a clear trend toward conserving a rigid backbone. We therefore propose that backbone rigidity is a fundamental physical feature conserved by proteins that can provide important insights into their folding mechanisms and stability.

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.

Advances in turbulent mixing techniques to study microsecond protein folding reactions

Recent experimental and computational advances in the protein folding arena have shown that the readout of the one-dimensional sequence information into three-dimensional structure begins within the first few microseconds of folding. The initiation of refolding reactions has been achieved by several means, including temperature jumps, flash photolysis, pressure jumps, and rapid mixing methods. One of the most commonly used means of initiating refolding of chemically denatured proteins is by turbulent flow mixing with refolding dilution buffer, where greater than 99% mixing efficiency has been achieved within 10's of microseconds. Successful interfacing of turbulent flow mixers with complementary detection methods, including time-resolved Fluorescence Spectroscopy (trFL), Förster Resonance Energy Transfer, Circular Dichroism, Small-Angle X-ray Scattering, Hydrogen Exchange followed by Mass Spectrometry and Nuclear Magnetic Resonance Spectroscopy, Infrared Spectroscopy (IR), and Fourier Transform IR Spectroscopy, has made this technique very attractive for monitoring various aspects of structure formation during folding. Although continuous-flow (CF) mixing devices interfaced with trFL detection have a dead time of only 30 µs, burst phases have been detected in this time scale during folding of peptides and of large proteins (e.g., CheY and TIM barrels). Furthermore, a major limitation of the CF mixing technique has been the requirement of large quantities of sample. In this brief communication, we will discuss the recent flurry of activity in micromachining and microfluidics, guided by computational simulations, which are likely to lead to dramatic improvements in time resolution and sample consumption for CF mixers over the next few years.

Interplay between drying and stability of a TIM barrel protein: a combined simulation-experimental study

Recent molecular dynamics simulations have suggested important roles for nanoscale dewetting in the stability, function, and folding dynamics of proteins. Using a synergistic simulation-experimental approach on the αTS TIM barrel protein, we validated this hypothesis by revealing the occurrence of drying inside hydrophobic amino acid clusters and its manifestation in experimental measures of protein stability and structure. Cavities created within three clusters of branched aliphatic amino acids [isoleucine, leucine, and valine (ILV) clusters] were found to experience strong water density fluctuations or intermittent dewetting transitions in simulations. Individually substituting 10 residues in the large ILV cluster at the N-terminus with less hydrophobic alanines showed a weakening or diminishing effect on dewetting that depended on the site of the mutation. Our simulations also demonstrated that replacement of buried leucines with isosteric, polar asparagines enhanced the wetting of the N- and C-terminal clusters. The experimental results on the stability, secondary structure, and compactness of the native and intermediate states for the asparagine variants are consistent with the preferential drying of the large N-terminal cluster in the intermediate. By contrast, the region encompassing the small C-terminal cluster experiences only partial drying in the intermediate, and its structure and stability are unaffected by the asparagine substitution. Surprisingly, the structural distortions required to accommodate the replacement of leucine by asparagine in the N-terminal cluster revealed the existence of alternative stable folds in the native basin. This combined simulation-experimental study demonstrates the critical role of drying within hydrophobic ILV clusters in the folding and stability of the αTS TIM barrel.

Deciphering the preference and predicting the viability of circular permutations in proteins

Circular permutation (CP) refers to situations in which the termini of a protein are relocated to other positions in the structure. CP occurs naturally and has been artificially created to study protein function, stability and folding. Recently CP is increasingly applied to engineer enzyme structure and function, and to create bifunctional fusion proteins unachievable by tandem fusion. CP is a complicated and expensive technique. An intrinsic difficulty in its application lies in the fact that not every position in a protein is amenable for creating a viable permutant. To examine the preferences of CP and develop CP viability prediction methods, we carried out comprehensive analyses of the sequence, structural, and dynamical properties of known CP sites using a variety of statistics and simulation methods, such as the bootstrap aggregating, permutation test and molecular dynamics simulations. CP particularly favors Gly, Pro, Asp and Asn. Positions preferred by CP lie within coils, loops, turns, and at residues that are exposed to solvent, weakly hydrogen-bonded, environmentally unpacked, or flexible. Disfavored positions include Cys, bulky hydrophobic residues, and residues located within helices or near the protein's core. These results fostered the development of an effective viable CP site prediction system, which combined four machine learning methods, e.g., artificial neural networks, the support vector machine, a random forest, and a hierarchical feature integration procedure developed in this work. As assessed by using the hydrofolate reductase dataset as the independent evaluation dataset, this prediction system achieved an AUC of 0.9. Large-scale predictions have been performed for nine thousand representative protein structures; several new potential applications of CP were thus identified. Many unreported preferences of CP are revealed in this study. The developed system is the best CP viability prediction method currently available. This work will facilitate the application of CP in research and biotechnology.

β-Bulge triggers route-switching on the functional landscape of interleukin-1β

Proteins fold into three-dimensional structures in a funneled energy landscape. This landscape is also used for functional activity. Frustration in this landscape can arise from the competing evolutionary pressures of biological function and reliable folding. Thus, the ensemble of partially folded states can populate multiple routes on this journey to the native state. Although protein folding kinetics experiments have shown the presence of such routes for several proteins, there has been sparse information about the structural diversity of these routes. In addition, why a given protein populates a particular route more often than another protein of similar structure and sequence is not clear. Whereas multiple routes are observed in theoretical studies on the folding of interleukin-1β (IL-1β), experimental results indicate one dominant route where the central portion of the protein folds first, and is then followed by closure of the barrel in this β-trefoil fold. Here we show, using a combination of computation and experiment, that the presence of functionally important regions like the β-bulge in the signaling protein IL-1β strongly influences the choice of folding routes. By deleting the β-bulge, we directly observe the presence of route-switching. This route-switching provides a direct link between route selection and the folding and functional landscapes of a protein.

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.

Folding of an all-helical Greek-key protein monitored by quenched-flow hydrogen-deuterium exchange and NMR spectroscopy

To advance our understanding of the protein folding process, we use stopped-flow far-ultraviolet (far-UV) circular dichroism and quenched-flow hydrogen-deuterium exchange coupled with nuclear magnetic resonance (NMR) spectroscopy to monitor the formation of hydrogen-bonded secondary structure in the C-terminal domain of the Fas-associated death domain (Fadd-DD). The death domain superfamily fold consists of six α-helices arranged in a Greek-key topology, which is shared by the all-β-sheet immunoglobulin and mixed α/β-plait superfamilies. Fadd-DD is selected as our model death domain protein system because the structure of this protein has been solved by NMR spectroscopy, and both thermodynamic and kinetic analysis indicate it to be a stable, monomeric protein with a rapidly formed hydrophobic core. Stopped-flow far-UV circular dichroism spectroscopy revealed that the folding process was monophasic and the rate is 23.4 s(-1). Twenty-two amide hydrogens in the backbone of the helices and two in the backbone of the loops were monitored, and the folding of all six helices was determined to be monophasic with rate constants between 19 and 22 s(-1). These results indicate that the formation of secondary structure is largely cooperative and concomitant with the hydrophobic collapse. This study also provides unprecedented insight into the formation of secondary structure within the highly populated Greek-key fold more generally. Additional insights are gained by calculating the exchange rates of 23 residues from equilibrium hydrogen-deuterium exchange experiments. The majority of protected amide protons are found on helices 2, 4, and 5, which make up core structural elements of the Greek-key topology.

Minireview: structural insights into early folding events using continuous-flow time-resolved small-angle X-ray scattering

Small-angle X-ray scattering (SAXS) is a powerful method for obtaining quantitative structural information on the size and shape of proteins, and it is increasingly used in kinetic studies of folding and association reactions. In this minireview, we discuss recent developments in using SAXS to obtain structural information on the unfolded ensemble and early folding intermediates of proteins using continuous-flow mixing devices. Interfacing of these micromachined devices to SAXS beamlines has allowed access to the microsecond time regime. The experimental constraints in implementation of turbulence and laminar flow-based mixers with SAXS detection and a comparison of the two approaches are presented. Current improvements and future prospects of microsecond time-resolved SAXS and the synergy with ab initio structure prediction and molecular dynamics simulations are discussed.

Systematic alanine insertion reveals the essential regions that encode structure formation and activity of dihydrofolate reductase

Decoding sequence information is equivalent to elucidating the design principles of proteins. For this purpose, we conducted systematic alanine insertion analysis to reveal the regions in the primary structure where the sequence continuity cannot be disrupted. We applied this method to dihydrofolate reductase (DHFR), and examined the effects of alanine insertion on structure and the enzymatic activity by solubility assay and trimethoprim resistance, respectively. We revealed that DHFR is composed of “Structure Elements”, “Function Elements” and linkers connecting these elements. The “Elements” are defined as regions where the alanine insertion caused DHFR to become unstructured or inactive. Some “Structure Elements” overlap with “Function Elements”, indicating that loss of structure leads to loss of function. However, other “Structure Elements” are not “Function Elements”, in that alanine insertion mutants of these regions exhibit substrate- or inhibitor-induced folding. There are also some “Function Elements” which are not “Structure Elements”; alanine insertion into these elements deforms the catalytic site topology without the loss of tertiary structure. We hypothesize that these elements are involved essential interactions for structure formation and functional expression. The “Elements” are closely related to the module structure of DHFR. An “Element” belongs to a single module, and a single module is composed of some number of “Elements.” We propose that properties of a module are determined by the “Elements” it contains. Systematic alanine insertion analysis is an effective and unique method for deriving the regions of a sequence that are essential for structure formation and functional expression.

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.

Structural analysis of kinetic folding intermediates for a TIM barrel protein, indole-3-glycerol phosphate synthase, by hydrogen exchange mass spectrometry and Gō model simulation

The structures of partially folded states appearing during the folding of a (betaalpha)(8) TIM barrel protein, the indole-3-glycerol phosphate synthase from Sulfolobus solfataricus (sIGPS), was assessed by hydrogen exchange mass spectrometry (HX-MS) and Gō model simulations. HX-MS analysis of the peptic peptides derived from the pulse-labeled product of the sub-millisecond folding reaction from the urea-denatured state revealed strong protection in the (betaalpha)(4) region, modest protection in the neighboring (betaalpha)(1-3) and (betaalpha)(5)beta(6) segments and no significant protection in the remaining N and C-terminal segments. These results demonstrate that this species is not a collapsed form of the unfolded state under native-favoring conditions nor is it the native state formed via fast-track folding. However, the striking contrast of these results with the strong protection observed in the (betaalpha)(2-5)beta(6) region after 5 s of folding demonstrates that these species represent kinetically distinct folding intermediates that are not identical as previously thought. A re-examination of the kinetic folding mechanism by chevron analysis of fluorescence data confirmed distinct roles for these two species: the burst-phase intermediate is predicted to be a misfolded, off-pathway intermediate, while the subsequent 5 s intermediate corresponds to an on-pathway equilibrium intermediate. Comparison with the predictions using a C(alpha) Gō model simulation of the kinetic folding reaction for sIGPS shows good agreement with the core of the structure offering protection against exchange in the on-pathway intermediate(s). Because the native-centric Gō model simulations do not explicitly include sequence-specific information, the simulation results support the hypothesis that the topology of TIM barrel proteins is a primary determinant of the folding free energy surface for the productive folding reaction. The early misfolding reaction must involve aspects of non-native structure not detected by the Gō model simulation.

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.

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.

Prediction of viable circular permutants using a graph theoretic approach

MOTIVATION

In recent years graph-theoretic descriptions have been applied to aid the analysis of a number of complex biological systems. However, such an approach has only just begun to be applied to examine protein structures and the network of interactions between residues with promising results. Here we examine whether a graph measure known as closeness is capable of predicting regions where a protein can be split to form a viable circular permutant. Circular permutants are a powerful experimental tool to probe folding mechanisms and more recently have been used to design split enzyme reporter proteins.

RESULTS

We test our method on an extensive set of experiments carried out on dihydrofolate reductase in which circular permutants were constructed for every amino acid position in the sequence, together with partial data from studies on other proteins. Results show that closeness is capable of correctly identifying significantly more residues which are suitable for circular permutation than solvent accessibility. This has potential implications for the design of successful split enzymes having particular importance for the development of protein-protein interaction screening methods and offers new perspectives on protein folding. More generally, the method illustrates the success with which graph-theoretic measures encapsulate the variety of long and short range interactions between residues during the folding process.

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.

Specific collapse followed by slow hydrogen-bond formation of beta-sheet in the folding of single-chain monellin

Characterization of the conformational landscapes for proteins with different secondary structures is important in elucidating the mechanism of protein folding. The folding trajectory of single-chain monellin composed of a five-stranded beta-sheet and a helix was investigated by using a pH-jump from the alkaline unfolded to native state. The kinetic changes in the secondary structures and in the overall size and shape were measured by circular dichroism spectroscopy and small-angle x-ray scattering, respectively. The formation of the tertiary structure was monitored by intrinsic and extrinsic fluorescence. A significant collapse was observed within 300 micros after the pH-jump, leading to the intermediate with a small amount of secondary and tertiary structures but with an overall oblate shape. Subsequently, the stepwise formation of secondary and tertiary structures was detected. The current observation was consistent with the theoretical prediction that a more significant collapse precedes the formation of secondary structures in the folding of beta-sheet proteins than that of helical proteins [Shea, J. E., Onuchic, J. N. & Brooks, C. L., III (2002) Proc. Natl. Acad. Sci. USA 99, 16064-16068]. Furthermore, it was implied that the initial collapse was promoted by the formation of some specific structural elements, such as tight turns, to form the oblate shape.

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.

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.

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.

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.

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.

An essential intermediate in the folding of dihydrofolate reductase

The folding of Escherichia coli dihydrofolate reductase was examined at pH 7.8 and 15 degrees C by using stopped-flow fluorescence and absorbance spectroscopies. The formation of a highly fluorescent intermediate occurs with relaxation times ranging between 142 and 343 msec, whereas stopped-flow absorbance spectroscopy using methotrexate binding assays shows a distinct lag phase during these time frames for the native state. The lag in absorbance kinetics and the lack of fast-track folding events indicate that the formation of this ensemble of intermediates is an obligatory step in the folding reaction.

How native-state topology affects the folding of dihydrofolate reductase and interleukin-1beta

The overall structure of the transition-state and intermediate ensembles observed experimentally for dihydrofolate reductase and interleukin-1beta can be obtained by using simplified models that have almost no energetic frustration. The predictive power of these models suggests that, even for these very large proteins with completely different folding mechanisms and functions, real protein sequences are sufficiently well designed, and much of the structural heterogeneity observed in the intermediates and the transition-state ensembles is determined by topological effects.

Localized, stereochemically sensitive hydrophobic packing in an early folding intermediate of dihydrofolate reductase from Escherichia coli

Mutational analysis was performed to probe the development of hydrophobic clusters during the early events in the folding of dihydrofolate reductase. Replacements were made in several hydrophobic subdomains to examine the roles of hydrophobicity and stereochemistry in the formation of kinetic intermediates. Amide protons in two of these clusters, including residues I91, I94, and I155, have been shown to be protected against solvent exchange within 13 ms of folding. Additional hydrophobic clusters were probed by substitutions at residues I2, I61, and L112; these residues are not protected from exchange until later in the folding reaction. Valine and leucine replacements at positions I91, I94, and I155 significantly diminish the formation of the burst phase kinetic intermediate, relative to the wild-type protein. In contrast, I2 and I61 are insensitive to these substitutions in the first 5 ms of the folding reaction, as is the replacement of L112 with either isoleucine or valine. These results demonstrate that the tightly packed core of dihydrofolate reductase is acquired in a non-uniform fashion, beginning in the submillisecond time frame. The progressive development of specific side-chain packing in localized hydrophobic clusters may be a common theme for complex protein folding reactions.

Identifying the structural boundaries of independent folding domains in the alpha subunit of tryptophan synthase, a beta/alpha barrel protein

Two equilibrium intermediates have previously been observed in the urea denaturation of the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli, an eight-stranded beta/alpha barrel protein. In the current study, a series of amino-terminal fragments were characterized to probe the elementary folding units that may be in part responsible for this complex behavior. Stop-codon mutagenesis was used to produce eight fragments ranging in size from 105-214 residues and containing incremental elements of secondary structure. Equilibrium studies by circular dichroism indicate that all of these fragments are capable of adopting secondary structure. All except for the shortest fragment fold cooperatively. The addition of the fourth, sixth, and eighth beta-strands leads to distinct increases in structure, cooperativity, and/or stability, suggesting that folding involves the modular assembly of betaalphabeta supersecondary structural elements. One-dimensional NMR titrations at high concentrations of urea, probing the environment around His92, were also performed to test for the presence of residual structure in the fragments. All fragments that contained the first four betaalpha units of structure exhibited a cooperative unfolding transition at high concentrations of urea with significant but reduced stability relative to the full-length protein. These results suggest that the residual structure in alphaTS requires the participation of hydrophobic residues in multiple beta-strands that span the entire sequence.

The chaperonin GroEL binds to late-folding non-native conformations present in native Escherichia coli and murine dihydrofolate reductases

Dihydrofolate reductases from mouse (MuDHFR) or Escherichia coli (EcDHFR) are shown to refold via several intermediate forms, each of which can bind to the chaperonin GroEL. When stable complexes with GroEL are formed, they consist of late-folding intermediates. In addition, we find that late-folding intermediates that are present in the native enzyme bind to GroEL. For the E. coli and murine proteins, the extent of protein bound increases as the temperature is increased from 8 degreesC to 42 degreesC, at which temperature either protein is completely bound as the last (EcDHFR) or the last two (MuDHFR) folding intermediate(s). Thus for EcDHFR, the binding is transient at low temperature (<30 degreesC) and stable at high temperature (>35 degreesC). For MuDHFR, complex formation appears less temperature dependent. In general, the data demonstrate that the overall binding free energy for the interaction of GroEL with native DHFR is the sum of the free energy for the first step in DHFR unfolding, which is unfavorable, and the free energy of binding the non-native conformation, which is favorable. For EcDHFR, this results in an overall binding free energy that is unfavorable below 30 degreesC. Over the temperature range of 8 degreesC to 42 degreesC, GroEL binds MuDHFR more tightly than EcDHFR, due partially to a small free energy difference between two pre-existing non-native conformations of MuDHFR, resulting in binding to more than one folding intermediate.

Native Escherichia coli and murine dihydrofolate reductases contain late-folding non-native structures

We have examined the equilibrium and kinetic folding properties of two structurally homologous dihydrofolate reductases, Escherichia coli DHFR (EcDHFR) and murine DHFR (MuDHFR), as a function of temperature and ligand concentration. Conformational heterogeneity in native DHFR is well documented, and the results demonstrate that the non-native form(s) represents late intermediate(s) in the folding process. We have measured the concentrations of native and non-native forms and the rate constants for their interconversion over a temperature range of 3 degreesC to 49 degreesC, allowing characterization of the thermodynamic as well as the kinetic properties of the final folding step(s) relative to the overall folding reaction. Differences in ligand binding suggest that the intermediate structures for these two proteins may be different during refolding.

Amide protection in an early folding intermediate of cytochrome c

BACKGROUND

For many proteins, compact states appear long before the rate-limiting step in the formation of the native structure. A key issue is whether the initial collapse of the chain is driven by random or more specific hydrophobic interactions.

RESULTS

Hydrogen-exchange labeling coupled with NMR was used to monitor the formation of stable hydrogen-bonded and solvent-excluded structure in horse cytochrome c (cyt c). Protection was measured using a hydrogen exchange/folding competition protocol at variable pH and short competition time (2 ms). Protection factors of threefold to eightfold were observed in all three alpha helices of cyt c, whereas other regions showed no significant protection. This suggests that the compact states that are present contain segments of marginally stable hydrogen-bonded structure. When the intermediate(s) are destabilized, only amide protons from Cys14, Ala15 and His18 show significant protection, indicating a region of persistent residual structure near the covalently bound heme group in the unfolded protein. Fluorescence-detected stopped-flow studies showed that the maximum protection factor in the early intermediate is consistent with its unfolding equilibrium constant.

CONCLUSIONS

Together with previous fluorescence and CD results, the observed pattern of amide protection is consistent with the early formation of an alpha-helical core domain in an ensemble of compact states, indicating that efficient folding is facilitated by stepwise acquisition of native structural elements. These specific early interactions are established on the sub-millisecond time scale, prior to the rate-limiting step for folding.

Molecular dynamics simulations of hydrophobic collapse of ubiquitin

Nine nonnative conformations of ubiquitin, generated during two different thermal denaturation trajectories, were simulated under nearly native conditions (62 degrees C). The simulations included all protein and solvent atoms explicitly, and simulation times ranged from 1-2.4 ns. The starting structures had alpha-carbon root-mean-square deviations (RMSDs) from the crystal structure of 4-12 A and radii of gyration as high as 1.3 times that of the native state. In all but one case, the protein collapsed when the temperature was lowered and sampled conformations as compact as those reached in a control simulation beginning from the crystal structure. In contrast, the protein did not collapse when simulated in a 60% methanol:water mixture. The behavior of the protein depended on the starting structure: during simulation of the most native-like starting structures (<5 A RMSD to the crystal structure) the RMSD decreased, the number of native hydrogen bonds increased, and the secondary and tertiary structure increased. Intermediate starting structures (5-10 A RMSD) collapsed to the radius of gyration of the control simulation, hydrophobic residues were preferentially buried, and the protein acquired some native contacts. However, the protein did not refold. The least native starting structures (10-12 A RMSD) did not collapse as completely as the more native-like structures; instead, they experienced large fluctuations in radius of gyration and went through cycles of expansion and collapse, with improved burial of hydrophobic residues in successive collapsed states.

Prolyl isomerases do not catalyze isomerization of non-prolyl peptide bonds

Prolyl isomerases accelerate the cis <--> trans isomerization of prolyl peptide bonds during protein folding and probably also in folded proteins. We asked whether this catalytic function is in fact restricted to prolyl bonds or whether the isomerizations of 'normal' non-prolyl peptide bonds are catalyzed as well. By using the P39A variant of ribonuclease T1 as a substrate we find that the trans --> cis isomerization of the Tyr38-Ala39 bond in the refolding of this protein is not catalyzed by prolyl isomerases of the cyclophilin, FKBP and parvulin families. These enzymes are neither able to catalyze amide bond isomerizations in the proline-free model peptide Ala-Ala-Tyr-Ala-Ala.

Hydrogen exchange studies of protein structure

Hydrogen exchange techniques, with their residue-level specificity, exquisite sensitivity, and adaptability to many solution conditions, are becoming essential to the study of protein stability, folding and dynamics. Recent studies have elucidated the structures of intermediates formed transiently during protein folding and rare partially folded ensembles present at equilibrium. Analysis of hydrogen exchange mechanisms has revealed protein stability and kinetics at the level of individual residues.

Probing minimal independent folding units in dihydrofolate reductase by molecular dissection

Molecular dissection was employed to identify minimal independent folding units in dihydrofolate reductase (DHFR) from Escherichia coli. Eight overlapping fragments of DHFR, spanning the entire sequence and ranging in size from 36 to 123 amino acids, were constructed by chemical cleavage. These fragments were designed to examine the effect of tethering multiple elements of secondary structure on folding and to test if the secondary structural domains represent autonomous folding units. CD and fluorescence spectroscopy demonstrated that six fragments containing up to a total of seven alpha-helices or beta-strands and, in three cases, the adenine binding domain (residues 37-86), are largely disordered. A stoichiometric mixture of the two fragments comprising the large discontinuous domain, 1-36 and 87-159, also showed no evidence for folding beyond that observed for the isolated fragments. A fragment containing residues 1-107 appears to have secondary and tertiary structure; however, spontaneous self-association made it impossible to determine if this structure solely reflects the behavior of the monomeric form. In contrast, a monomeric fragment spanning residues 37-159 possesses significant secondary and tertiary structure. The urea-induced unfolding of fragment 37-159 in the presence of 0.5 M ammonium sulfate was found to be a well-defined, two-state process. The observation that fragment 37-159 can adopt a stable native fold with unique, aromatic side-chain packing is quite striking because residues 1-36 form an integral part of the structural core of the full-length protein.

Evidence for an obligatory intermediate in the folding of interleukin-1 beta

The folding of the beta-sheet protein, interleukin-1 beta, was examined at pH 5.0 and 25 degrees C using pulse-labelling hydrogen exchange and electrospray ionization mass spectrometric analysis, as well as stopped-flow circular dichroism and fluorescence spectroscopies. The first detectable event is the formation of a partially folded intermediate in a kinetic step with a relaxation time of 126 +/- 26 ms. There is a lag in native protein production of at least 400 ms. Optical studies indicate that the intermediate is converted to the native species in a reaction with a relaxation time of 43 +/- 5 s. The kinetic rates determined from stopped-flow fluorescence, circular dichroism and pulse-labelling experiments are similar and consistent with a simple sequential model for the folding pathway of interleukin-1 beta at pH 5.0 and 25 degrees C. Taken together, our data provide kinetic evidence that formation of the native state of interleukin-1 beta proceeds through an obligatory intermediate. We explain our results in terms of the classical and new views of protein folding.

Protein folding and intermediates

With the exception of the discovery of the rate of formation of the earliest intermediates, there have been no major conceptual leaps in our understanding of protein folding reactions over the past two years. Rather, this period has seen an extension of two established techniques: first, mutational analysis combined with a kinetic definition of the energy landscape of the reaction; and second, the use of hydrogen/deuterium exchange of backbone amide groups combined with NMR. Owing to the application of these methods to a wider range of proteins, it is now possible to draw some general conclusions about the physical processes that direct a protein to its native fold.

Comparison of the amide proton exchange behavior of the rapidly formed folding intermediate and the native state of an antibody scFv fragment

We have investigated the stability of backbone amide protons of the intermediate and the native state of the scFv fragment of an antibody. Stopped flow experiments analyzed by MS and NMR detected the formation of an exchange protected intermediate within the deadtime of the stopped flow apparatus (17 ms). H/D exchange rates of the native protein identified a number of very stable backbone amide protons in the V(L) and the V(H) domains. In the V(L) domain, this slowly exchanging core of the scFv fragment is similar to the folding core of the intermediate, while the V(H) domain possesses a great number of very stable amide protons which are not stabilized to a significant degree in the folding intermediate of the scFv fragment.