Two-dimensional 1H nuclear magnetic resonance study of the (5-55) single-disulphide folding intermediate of bovine pancreatic trypsin inhibitor

An analogue of the bovine pancreatic trypsin inhibitor (BPTI) folding intermediate that contains only the disulphide bond between Cys5 and Cys55 has been prepared in Escherichia coli by protein engineering methods, with the other four Cys residues replaced by Ser. Two-dimensional 1H nuclear magnetic resonance studies of the analogue have resulted in essentially complete resonance assignments of the folded form of the protein. The folded protein has a compact conformation that is structurally very similar to that of native BPTI, although there are subtle differences and the folded conformation is not very stable. Approximately half of the protein molecules are unfolded at 3 degrees C, and this proportion increases at higher temperatures. The folded and unfolded conformations are in slow exchange. The conformational properties of the analogue can explain many aspects of the kinetic role that the normal (5-55) intermediate plays in the folding of BPTI.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Transient intermediary states with high and low folding probabilities in the apparent two-state folding equilibrium of ACBP at low pH

Measurements of the stability as a function of pH for the acyl-coenzyme A binding protein (ACBP) has shown a significant difference in the pH transition midpoint measured by NMR spectroscopy at pH 3.12 and the transition midpoint measured at pH 2.92 and 2.97 by circular dichroism and by fluorescence spectroscopy, respectively. A similar behavior has not been observed in other proteins. It is suggested that these differences arise because the population of the unfolded molecules still contains significant amounts of native like secondary and tertiary structure. NMR spectroscopy measures the concentration of the two components of the folding unfolding equilibrium individually, whereas circular dichroism and fluorescence measure the concentration of the conformations of the light-absorbing chromophores present in both the folded and the unfolded molecules. In the narrow pH range, nascent structure can be detected as the average amount of secondary structure per unfolded molecule and hydrophobic interactions in the population of unfolded molecules. These structures are not observable immediately by NMR spectroscopy; however, a chemical shift analysis of the peptide backbone (13)C chemical shift indicates strongly the existence of short-lived and transient helical structures at pH 2.3. Magnetization transfer studies have been applied to study the equilibrium between folded and unfolded ACBP near the pH transition point measured by NMR. This study has shown that there are two categories of subpopulations in the population of unfolded ACBP. One for which magnetization can be transferred to the folded form during the folding process, and one for which transfer is not observed. The molecules of the latter population of unfolded protein apparently, do not fold within the time-frame of the magnetization transfer experiment. This result suggests the existence of a subpopulation of the acid-unfolded protein molecules with a high propensity for folding. It is suggested that in this subpopulation, a particular set of native like interactions in the peptide backbone and between side-chains in the peptide chain have to be formed.     

Native protein sequences are designed to destabilize folding intermediates

Hydrophobic core mutants of sperm whale apomyoglobin were constructed to investigate the amino acid sequence features that determine the folding properties. Replacements of all of the Ile residues with Leu and of all of the Ile and Val residues with Leu decreased the thermodynamic stability of the folded states against the unfolded states but increased the stability of the folding intermediates against the unfolded states, indicating that the amino acid composition of the protein core is important for the protein stability and folding cooperativity. To examine the effect of the arrangement of these hydrophobic residues, mutant proteins were further constructed: 12 sites out of the 18 Leu, 9 Ile, and 8 Val residues of the wild-type myoglobin were randomly replaced with each other so that the amino acid compositions were similar to that of the wild-type protein. Four mutant proteins were obtained without selection of the protein properties. These residue replacements similarly resulted in the stabilization of both the intermediate and folded states against the unfolded states, as compared to the wild-type protein. Thus, the arrangements of the hydrophobic residues in the native amino acid sequence are selected to destabilize the folding intermediate rather than to stabilize the folded state. The present results suggest that the two-state transition of protein folding or the transient formation of the unstable intermediate, which seems to be required for effective production of the functional proteins, has been a major driving force in the molecular evolution of natural globular proteins.     

Structural characterization of partially folded intermediates of apomyoglobin H64F

We present a detailed investigation of unfolded and partially folded states of a mutant apomyoglobin (apoMb) where the distal histidine has been replaced by phenylalanine (H64F). Previous studies have shown that substitution of His64, located in the E helix of the native protein, stabilizes the equilibrium molten globule and native states and leads to an increase in folding rate and a change in the folding pathway. Analysis of changes in chemical shift and in backbone flexibility, detected via [1H]-15N heteronuclear nuclear Overhauser effect measurements, indicates that the phenylalanine substitution has only minor effects on the conformational ensemble in the acid- and urea-unfolded states, but has a substantial effect on the structure, dynamics, and stability of the equilibrium molten globule intermediate formed near pH 4. In H64F apomyoglobin, additional regions of the polypeptide chain are recruited into the compact core of the molten globule. Since the phenylalanine substitution has negligible effect on the unfolded ensemble, its influence on folding rate and stability comes entirely from interactions within the compact folded or partly folded states. Replacement of His64 with Phe leads to favorable hydrophobic packing between the helix E region and the molten globule core and leads to stabilization of helix E secondary structure and overall thermodynamic stabilization of the molten globule. The secondary structure of the equilibrium molten globule parallels that of the burst phase kinetic intermediate; both intermediates contain significant helical structure in regions of the polypeptide that comprise the A, B, E, G, and H helices of the fully folded protein.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

Fluorescence characterization of denatured proteins

Characterization of unfolded states, while critical to a complete understanding of protein folding, is inherently difficult due to structural heterogeneity and dynamic interchange between states. The growing body of work focusing on single molecule fluorescence techniques for the study of protein folding, also highlights their potential for studies of unfolded proteins. These methods can obtain conformational information about individual subpopulations of molecules in an ensemble, and measure dynamics without the need for synchronization. The studies highlighted here demonstrate the promise of these techniques for obtaining novel information about unfolded states in vitro and in more physiologically relevant milieu.     

8-anilino-1-naphthalene sulfonic acid (ANS) induces folding of acid unfolded cytochrome c to molten globule state as a result of electrostatic interactions.

Hydrophobic interaction of 8-anilino-1-naphthalene sulfonic acid (ANS) with proteins is one of the widely used methods for characterizing/detecting partially folded states of proteins. We have carried out a systematic investigation on the effect of ANS, a charged hydrophobic fluorescent dye, on structural properties of acid-unfolded horse heart cytochrome c at pH 2.0 by a combination of optical methods and electrospray ionization mass spectroscopy (ESI MS). ANS was found to induce, a secondary structure similar to native protein and quenching of fluorescence of tryptophan residue, in the acid-unfolded protein. However, the tertiary structure was found to be disrupted thus indicating that ANS stabilizes a molten globule state in acid-unfolded protein. To understand the mechanism of ANS-induced folding of acid-unfolded cytochrome c, comparative ESI MS, soret absorption, and tryptophan fluorescence studies using nile red, a neutral hydrophobic dye, and ANS were carried out. These studies suggested that, at low pH, electrostatic interactions between negatively charged ANS molecules and positively charged amino acid residues present in acid-unfolded cytochrome c are probably responsible for ANS-induced folding of acid-unfolded protein to partially folded compact state or molten globule state. This is the first experimental demonstration of ANS induced folding of unfolded protein and puts to question the usefulness of ANS for characterization/determination of partially folded intermediates of proteins observed under low pH conditions.     

Characterization of the hydrodynamic properties of the folding transition state of an SH3 domain by magnetization transfer NMR spectroscopy

Protein folding kinetic data have been obtained for the marginally stable N-terminal Src homology 3 domain of the Drosophila protein drk (drkN SH3) in an investigation of the hydrodynamic properties of its folding transition state. Due to the presence of NMR resonances of both folded and unfolded states at equilibrium, kinetic data can be derived from NMR magnetization transfer techniques under equilibrium conditions. Kinetic analysis as a function of urea (less than approximately 1 M) and glycerol enables determination of alpha values, measures of the energetic sensitivity of the transition state to the perturbation relative to the end states of the protein folding reaction (the folded and unfolded states). Both end states have previously been studied experimentally by NMR spectroscopic and other biophysical methods in great detail and under nondenaturing conditions. Combining these results with the kinetic folding data obtained here, we can characterize the folding transition state without requiring empirical models for the unfolded state structure. We are thus able to give a reliable measure of the solvent-accessible surface area of the transition state of the drkN SH3 domain (4730 +/- 360 A(2)) based on urea titration data. Glycerol titration data give similar results and additionally demonstrate that folding of this SH3 domain is dependent on solvent viscosity, which is indicative of at least partial hydration of the transition state. Because SH3 domains appear to fold by a common folding mechanism, the data presented here provide valuable insight into the transition states of the drkN and other SH3 domains.     

Molecular basis of co-operativity in protein folding.

The folding/unfolding transition of proteins is a highly co-operative process characterized by the presence of very few or no thermodynamically stable partially folded intermediate states. The purpose of this paper is to present a thermodynamic formalism aimed at describing quantitatively the co-operative folding behavior of proteins. In order to account for this behavior, a hierarchical algorithm aimed at evaluating the folding/unfolding partition function has been developed. This formalism defines the partition function in terms of multiple levels of interacting co-operative folding units. A co-operative folding unit is defined as a protein structural element that exhibits two-state folding/unfolding behavior. At the most fundamental level are those structural elements that behave co-operatively as a result of purely local interactions. Higher-order co-operative folding units are formed through interactions between different structural elements. The hierarchical formalism utilizes the crystallographic structure of the protein as a template to generate partially folded conformations defined in terms of co-operative folding units. The Gibbs free energy of those states and their corresponding statistical weights are then computed using experimental energetic parameters determined calorimetrically. This formalism has been applied to the case of myoglobin. It is shown that the hierarchical partition function correctly predicts the presence, energetics and co-operativity of the heat and cold denaturation transitions. The major contribution to the co-operative folding behavior arises from the solvent exposure of non-polar residues located in regions complementary to those that have undergone unfolding. This entropically uncompensated and energetically unfavorable solvent exposure characterizes all partially folded states but not the unfolded state, thus minimizing the population of partially folded intermediates throughout the folding/unfolding transition.     

Protein folding cooperativity: basic insights from minimalist models

Basic concepts about two-state, cooperative protein folding and its relation to first-order phase transitions are reviewed. Minimalist models capable of reproducing the required free energy barrier between folded and unfolded macroscopic states are described. A significantly more restrictive "calorimetric" criterion is also discussed, which is based on direct comparison between model and experimental heat capacities with additional assumptions about conformational enthalpy variation in the unfolded state.     

Replica exchange simulation of reversible folding/unfolding of the Trp-cage miniprotein in explicit solvent: on the structure and possible role of internal water

We simulate the folding/unfolding equilibrium of the 20-residue miniprotein Trp-cage. We use replica exchange molecular dynamics simulations of the AMBER94 atomic detail model of the protein explicitly solvated by water, starting from a completely unfolded configuration. We employ a total of 40 replicas, covering the temperature range between 280 and 538 K. Individual simulation lengths of 100 ns sum up to a total simulation time of about 4 micros. Without any bias, we observe the folding of the protein into the native state with an unfolding-transition temperature of about 440 K. The native state is characterized by a distribution of root mean square distances (RMSD) from the NMR data that peaks at 1.8A, and is as low as 0.4A. We show that equilibration times of about 40 ns are required to yield convergence. A folded configuration in the entire extended ensemble is found to have a lifetime of about 31 ns. In a clamp-like motion, the Trp-cage opens up during thermal denaturation. In line with fluorescence quenching experiments, the Trp-residue sidechain gets hydrated when the protein opens up, roughly doubling the number of water molecules in the first solvation shell. We find the helical propensity of the helical domain of Trp-cage rather well preserved even at very high temperatures. In the folded state, we can identify states with one and two buried internal water molecules interconnecting parts of the Trp-cage molecule by hydrogen bonds. The loss of hydrogen bonds of these buried water molecules in the folded state with increasing temperature is likely to destabilize the folded state at elevated temperatures.     

Water as a conformational editor in protein folding

As molecules approach one another in aqueous solution, desolvation free energy barriers to association are encountered. Experiments suggest these (de)solvation effects contribute to the free energy barriers separating the folded and unfolded states of protein molecules. To explore their influence on the energy landscapes of protein folding reactions, we have incorporated desolvation barriers into a semi-realistic, off-lattice protein model that uses a simplified physico-chemical force-field determined solely by the sequence of amino acids. Monte Carlo sampling techniques were used to study the effects on the thermodynamics and kinetics of folding of a number of systems, diverse in structure and sequence. In each case, desolvation barriers increase the stability of the native conformation and the cooperativity of the major folding/unfolding transition. The folding times of these systems are reduced significantly upon inclusion of desolvation barriers, demonstrating that the particulate nature of the solvent engenders a more defined route to the native fold.     

Reversible denaturation of the gene V protein of bacteriophage f1

The guanidine hydrochloride (GuHCl)-induced denaturation of the gene V protein of bacteriophage f1 has been studied, using the chemical reactivity of a cysteine residue that is buried in the folded protein and the circular dichroism (CD) at 211 and 229 nm as measures of the fraction of polypeptide chains in the folded form. It is found that this dimeric protein unfolds in a single cooperative transition from a folded dimer to two unfolded monomers. A folded, monomeric form of the gene V protein was not detected at equilibrium. The kinetics of unfolding of the gene V protein in 3 M GuHCl and the refolding in 2 M GuHCl are also consistent with a transition between a folded dimer and two unfolded monomers. The GuHCl concentration dependence of the rates of folding and unfolding suggests that the transition state for folding is near the folded conformation.     

Wild type, mutant protein unfolding and phase transition detected by single-nanopore recording

Understanding protein folding remains a challenge. A difficulty is to investigate experimentally all the conformations in the energy landscape. Only single molecule methods, fluorescence and force spectroscopy, allow observing individual molecules along their folding pathway. Here we observe that single-nanopore recording can be used as a new single molecule method to explore the unfolding transition and to examine the conformational space of native or variant proteins. We show that we can distinguish unfolded states from partially folded ones with the aerolysin pore. The unfolding transition curves of the destabilized variant are shifted toward the lower values of the denaturant agent compared to the wild type protein. The dynamics of the partially unfolded wild type protein follows a first-order transition. The denaturation curve obtained with the aerolysin pore is similar to that obtained with the α-hemolysin pore. The nanopore geometry or net charge does not influence the folding transition but changes the dynamics.     

Microscopic reversibility of protein folding in molecular dynamics simulations of the engrailed homeodomain

The principle of microscopic reversibility states that at equilibrium the number of molecules entering a state by a given path must equal those exiting the state via the same path under identical conditions or, in structural terms, that the conformations along the two pathways are the same. There has been some indirect evidence indicating that protein folding is such a process, but there have been few conclusive findings. In this study, we performed molecular dynamics simulations of an ultrafast unfolding and folding protein at its melting temperature to observe, on an atom-by-atom basis, the pathways the protein followed as it unfolded and folded within a continuous trajectory. In a total of 0.67 micros of simulation in water, we found six transient denaturing events near the melting temperature (323 and 330 K) and an additional refolding event following a previously identified unfolding event at a high temperature (373 K). In each case, unfolding and refolding transition state ensembles were identified, and they agreed well with experiment on the basis of a comparison of S and Phi values. On the basis of several structural properties, these 13 transition state ensembles agreed very well with each other and with four previously identified transition states from high-temperature denaturing simulations. Thus, not only were the unfolding and refolding transition states part of the same ensemble, but in five of the seven cases, the pathway the protein took as it unfolded was nearly identical to the subsequent refolding pathway. These events provide compelling evidence that protein folding is a microscopically reversible process. In the other two cases, the folding and unfolding transition states were remarkably similar to each other but the paths deviated.     

Conformation of peptide fragments of proteins in aqueous solution: implications for initiation of protein folding

Applications of sensitive new technologies, in particular, two-dimensional NMR spectroscopy, have allowed detection of folded structures in short peptide fragments of proteins in aqueous solution under conditions where native proteins fold. These structures are in rapid dynamic exchange with unfolded states. These observations provide evidence in support of models for protein folding which postulate localized regions of folded structure as initiation sites for the folding process. Since these initiation processes are expected to be rapid, such models are consistent with kinetic evidence that the rate-determining steps of protein folding occur late in the process and probably involve rearrangement of incorrectly folded intermediates.     

Local interactions and the optimization of protein folding

The role of local interactions in protein folding has recently been the subject of some controversy. Here we investigate an extension of Zwanzig's simple and general model of folding in which local and nonlocal interactions are represented by functions of single and multiple conformational degrees of freedom, respectively. The kinetics and thermodynamics of folding are studied for a series of energy functions in which the energy of the native structure is fixed, but the relative contributions of local and nonlocal interactions to this energy are varied over a broad range. For funnel shaped energy landscapes, we find that 1) the rate of folding increases, but the stability of the folded state decreases, as the contribution of local interactions to the energy of the native structure increases, and 2) the amount of native structure in the unfolded state and the transition state vary considerably with the local interaction strength. Simple exponential kinetics and a well-defined free energy barrier separating folded and unfolded states are observed when nonlocal interactions make an appreciable contribution to the energy of the native structure; in such cases a transition state theory type approximation yields reasonably accurate estimates of the folding rate. Bumps in the folding funnel near the native state, which could result from desolvation effects, side chain freezing, or the breaking of nonnative contacts, significantly alter the dependence of the folding rate on the local interaction strength: the rate of folding decreases when the local interaction strength is increased beyond a certain point. A survey of the distribution of strong contacts in the protein structure database suggests that evolutionary optimization has involved both kinetics and thermodynamics: strong contacts are enriched at both very short and very long sequence separations. Proteins 29:282–291, 1997. © 1997 Wiley-Liss, Inc.     

Long-range electrostatic effects on peptide folding.

The reversible folding/unfolding of a short peptide in solution is studied by molecular dynamics simulations. The effects of long-range electrostatic interactions are examined and found to be important both for the equilibrium between folded and unfolded states and the dynamics of the folding process. The neglect of long-range electrostatics leads to an increased population of unfolded states and increased structural fluctuations. When such interactions are taken into account, the peptide unfolds and folds to the experimentally determined structure several times during a 25 ns trajectory, with approximately equal populations of folded and unfolded states in the neighborhood of its proposed melting temperature. The effect of using spherical boundary conditions rather than periodic ones does not appear to have any major effect on the folding dynamics.