The early folding kinetics of apomyoglobin.

The folding pathway of apomyoglobin has been experimentally shown to have early kinetic intermediates involving the A, B, G, and H helices. The earliest detected kinetic events occur on a ns to micros time scale. We show that the early folding kinetics of apomyoglobin may be understood as the association of nascent helices through a network of diffusion-collision-coalescence steps G + H <--> GH + A <--> AGH + B <--> ABGH obtained by solving the diffusion-collision model in a chemical kinetics approximation. Our reproduction of the experimental results indicates that the model is a useful way to analyze folding data. One prediction from our fit is that the nascent A and H helices should be relatively more helix-like before coalescence than the other apomyoglobin helices.     

Resolution of the fluorescence equilibrium unfolding profile of trp aporepressor using single tryptophan mutants.

Single tryptophan mutants of the trp aporepressor, tryptophan 19-->phenylalanine (W19F) and tryptophan 99-->phenylalanine (W99F), were used in this study to resolve the individual steady-state and time-resolved fluorescence urea unfolding profiles of the two tryptophan residues in this highly intertwined, dimeric protein. The wild-type protein exhibits a large increase in fluorescence intensity and lifetime, as well as a large red shift in the steady-state fluorescence emission spectrum, upon unfolding by urea (Lane, A.N. & Jardetsky, O., 1987, Eur. J. Biochem. 164, 389-396; Gittelman, M.S. & Matthews, C.R., 1990, Biochemistry 29, 7011-7020; Fernando, T. & Royer, C.A., 1992, Biochemistry 31, 6683-6691). Unfolding of the W19F mutant demonstrated that Trp 99 undergoes a large increase in intensity and a red shift upon exposure to solvent. Lifetime studies revealed that the contribution of the dominant 0.5-ns component of this tryptophan tends toward zero with increasing urea, whereas the longer lifetime components increase in importance. This lifting of the quenching of Trp 99 may be due to disruption of the interaction between the two subunits upon denaturation, which abolishes the interaction of Trp 99 on one subunit with the amide quenching group of Asn 32 on the other subunit (Royer, C.A., 1992, Biophys. J. 63, 741-750). On the other hand, Trp 19 is quenched in response to unfolding in the W99F mutant. Exposure to solvent of Trp 19, which is buried at the hydrophobic dimer interface in the native protein, results in a large red shift of the average steady-state emission.(ABSTRACT TRUNCATED AT 250 WORDS)     

Effects of pressure on the structure of metmyoglobin: molecular dynamics predictions for pressure unfolding through a molten globule intermediate.

We investigated the pathway for pressure unfolding of metmyoglobin using molecular dynamics (MD) for a range of pressures (0.1 MPa to 1.2 GPa) and a temperature of 300 K. We find that the unfolding of metmyoglobin proceeds via a two-step mechanism native --> molten globule intermediate --> unfolded, where the molten globule forms at 700 MPa. The simulation describes qualitatively the experimental behavior of metmyoglobin under pressure. We find that unfolding of the alpha-helices follows the sequence of migrating hydrogen bonds (i,i + 4) --> (i,i + 2).     

Quantitative analysis of the kinetics of denaturation and renaturation of barstar in the folding transition zone.

The fluorescence-monitored kinetics of folding and unfolding of barstar by guanidine hydrochloride (GdnHCl) in the folding transition zone, at pH 7, 25 degrees C, have been quantitatively analyzed using a 3-state mechanism: U(S)<-->UF<-->N. U(S) and UF are slow-refolding and fast-refolding unfolded forms of barstar, and N is the native protein. U(S) and UF probably differ in possessing trans and cis conformations, respectively, of the Tyr 47-Pro 48 bond. The 3-state model could be used because the kinetics of folding and unfolding of barstar show 2 phases, a fast phase and a slow phase, and because the relative amplitudes of the 2 phases depend only on the final refolding conditions and not on the initial conditions. Analysis of the observed kinetics according to the 3-state model yields the values of the 4 microscopic rate constants that describe the transitions between the 3 states at different concentrations of GdnHCl. The value of the equilibrium unfolded ratio U(S):UF (K21) and the values of the rate constants of the U(S)-->UF and UF-->U(S) reactions, k12 and k21, respectively, are shown to be independent of the concentration of GdnHCl. K21 has a value of 2.1 +/- 0.1, and k12 and k21 have values of 5.3 x 10(-3) s-1 and 11.2 x 10(-3) s-1, respectively. Double-jump experiments that monitor reactions that are silent to fluorescence monitoring were used to confirm the values of K21, k12, and k21 obtained from the 3-state analysis and thereby the validity of the 3-state model.(ABSTRACT TRUNCATED AT 250 WORDS)     

Tryptophan replacements in the trp aporepressor from Escherichia coli: probing the equilibrium and kinetic folding models.

Mutants of the dimeric Escherichia coli trp aporepressor are constructed by replacement of the two tryptophan residues in each subunit in order to assess the effects on equilibrium and kinetic fluorescence properties of the folding reaction. The three kinetic phases detected by intrinsic tryptophan fluorescence in refolding of the wild-type aporepressor are also observed in folding of both Trp 19 to Phe and Trp 99 to Phe single mutants, demonstrating that these phases correspond to global rather than local conformational changes. Comparison of equilibrium fluorescence (Royer, C.A., Mann, C.J., & Matthews, C.R., 1993, Protein Sci. 2, 1844-1852) and circular dichroism transition curves induced by urea shows that replacement of either Trp 19 or Trp 99 results in noncoincident behavior. Unlike the wild-type protein (Gittelman, M.S. & Matthews, C.R., 1990, Biochemistry 29, 7011-7020), tertiary and/or quaternary structures are disrupted at lower denaturant concentration than is secondary structure. The equilibrium results can be interpreted in terms of enhancement in the population of a monomeric folding intermediate in which the lone tryptophan residue is highly exposed to solvent, but in which substantial secondary structure is retained. The location of both mutations at the interface between the two subunits (Zhang, R.G., et al., 1987, Nature 327, 591-597) provides a simple explanation for this phenomenon.     

Helical peptides with three pairs of Asp-Arg and Glu-Arg residues in different orientations and spacings.

The helix-stabilizing effects of repeating pairs of Asp-Arg and Glu-Arg residues have been characterized using a peptide system of the same design used earlier to study Glu-Lys (Marqusee, S. & Baldwin, R.L., 1987, Proc. Natl. Acad. Sci. USA 84, 8898-8902) and Asp-Lys ion pairs (Marqusee, S. & Baldwin, R.L., 1990, In Protein Folding [Gierasch, L.M. & King, J., Eds.], pp. 85-94, AAAS, Washington, D.C.). The consequences of breaking ion pair and charge-helix dipole interactions by titration to pH 2 have been compared with the results of screening these interactions with NaCl at pH 7.0 and pH 2.5. The four peptides in each set contain three pairs of acidic (A) and basic (B) residues spaced either i, i + 4 or i, i + 3 apart. In one peptide of each kind the pairwise order of residues is AB, with the charges oriented favorably to the helix macrodipole, and in the other peptide the order is BA. The results are as follows: (1) Remarkably, both Asp-Arg and Glu-Arg peptides show the same pattern of helix stabilization at pH 7.0 found earlier for Glu-Lys and Asp-Lys peptides: i + 4 AB > i + 4 BA approximately i + 3 AB > i + 3 BA. (2) The ion pairs and charge-helix dipole interactions cannot be cleanly separated, but the results suggest that both interactions make important contributions to helix stability.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular basis of cooperativity in protein folding. V. Thermodynamic and structural conditions for the stabilization of compact denatured states

The heat-denatured state of proteins has been usually assumed to be a fully hydrated random coil. It is now evident that under certain solvent conditions or after chemical or genetic modifications, the protein molecule may exhibit a hydrophobic core and residual secondary structure after thermal denaturation. This state of the protein has been called the “compact denatured” or “molten globule” state. Recently is has been shown that α-lactalbumin at pH < 5 denatures into a molten globule state upon increasing the temperature (Griko, Y., Freire, E., Privalov, P. L. Biochemistry 33:1889–1899, 1994). This state has a lower heat capacity and a higher enthalpy at low temperatures than the unfolded state. At those temperatures the stabilization of the molten globule state is of an entropic origin since the enthalpy contributes unfavorably to the Gibbs free energy. Since the molten globule is more structured than the unfolded state and, therefore, is expected to have a lower configurational entropy, the net entropic gain must originate primarily from solvent related entropy arising from the hydrophobic effect, and to a lesser extent from protonation or electrostatic effects. In this work, we have examined a large ensemble of partly folded states derived from the native structure of α-lactalbumin in order to identify those states that satisfy the energetic criteria of the molten globule. It was found that only few states satisfied the experimental constraints and that, furthermore, those states were part of the same structural family. In particular, the regions corresponding to the A, B, and C helices were found to be folded, while the β sheet and the D helix were found to be unfolded. At temperatures below 45°C the states exhibiting those structural characteristics are enthalpically higher than the unfolded state in agreement with the experimental data. Interestingly, those states have a heat capacity close to that observed for the acid pH compact denatured state of α-lactalbumin [980 cal (mol.K)^?l]. In addition, the folded regions of these states include those residues found to be highly protected by NMR hydrogen exchange experiments. This work represents an initial attempt to model the structural origin of the thermodynamic properties of partly folded states. The results suggest a number of structural features that are consistent with experimental data. © 1994 Wiley-Liss, Inc.     

Cysteine scanning mutagenesis of putative transmembrane helices IX and X in the lactose permease of Escherichia coli.

Using a functional lactose permease mutant devoid of Cys residues (C-less permease), each amino-acid residue in putative transmembrane helices IX and X and the short intervening loop was systematically replaced with Cys (from Asn-290 to Lys-335). Thirty-four of 46 mutants accumulate lactose to high levels (70-100% or more of C-less), and an additional 7 mutants exhibit lower but highly significant lactose accumulation. As expected (see Kaback, H.R., 1992, Int. Rev. Cytol. 137A, 97-125), Cys substitution for Arg-302, His-322, or Glu-325 results in inactive permease molecules. Although Cys replacement for Lys-319 or Phe-334 also inactivates lactose accumulation, Lys-319 is not essential for active lactose transport (Sahin-Tóth, M., Dunten, R.L., Gonzalez, A., & Kaback, H.R., 1992, Proc. Natl. Acad. Sci. USA 89, 10547-10551), and replacement of Phe-334 with leucine yields permease with considerable activity. All single-Cys mutants except Gly-296 --> Cys are present in the membrane in amounts comparable to C-less permease, as judged by immunological techniques. In contrast, mutant Gly-296 --> Cys is hardly detectable when expressed at a relatively low rate from the lac promoter/operator but present in the membrane in stable form when expressed at a high rate from T7 promoter. Finally, studies with N-ethylmaleimide (NEM) show that only a few mutants are inactivated significantly. Remarkably, the rate of inactivation of Val-315 --> Cys permease is enhanced at least 10-fold in the presence of beta-galactopyranosyl 1-thio-beta-D-galactopyranoside (TDG) or an H+ electrochemical gradient (delta mu-H+). The results demonstrate that only three residues in this region of the permease -Arg-302, His-322, and Glu-325-are essential for active lactose transport. Furthermore, the enhanced reactivity of the Val-315 --> Cys mutant toward NEM in the presence of TDG or delta mu-H+ probably reflects a conformational alteration induced by either substrate binding or delta mu-H+.     

Mechanism of polyethylene glycol interaction with the molten globule folding intermediate of bovine carbonic anhydrase B.

Polyethylene glycol has been shown to bind to the molten globule intermediate on the bovine carbonic anhydrase B folding pathway. The mechanism of this interaction has been extensively probed. Polyethylene glycol (PEG) binds weakly to the molten globule first intermediate as measured by hydrophobic interaction chromatography, but PEG does not bind to either the native state or the second intermediate. The binding of PEG to the molten globule has been confirmed with both intrinsic fluorescence and fluorescence quenching experiments which indicate a single PEG-binding site on the molten globule. Electron paramagnetic resonance spectroscopic studies with nitroxide-labeled PEG also indicate a single binding site. Additional electron paramagnetic resonance studies with spin-labeled carbonic anhydrase B suggest that a conformational change occurs in the molten globule intermediate after PEG binds to the surface. The formation of a PEG-molten globule complex results in a reduction in self-association of this compact hydrophobic structure. PEG-molten globule complex formation is analogous to the observed interaction between chaperonins and a molten globule intermediate (Martin, J., Langer, T., Boteva, R., Schramel, A., Horwich, A.L., and Hartl, F.U. (1991) Nature 352, 36-42).     

Fourier transform infrared spectroscopy characterization of the lamellar and nonlamellar structures of free lipid A and Re lipopolysaccharides from Salmonella minnesota and Escherichia coli.

The structural polymorphism of free lipid A and deep rough mutant lipopolysaccharide (LPS Re) from Salmonella minnesota strain R595 and Escherichia coli strain F515 was characterized by Fourier transform infrared (IR) spectroscopy. For this, the beta <--> alpha phase states and the three-dimensional supramolecular structures, the latter deduced from small-angle synchrotron radiation x-ray diffraction, were investigated at different water contents, Mg2+ concentrations, and temperatures. The analysis of the IR data for vibrations originating from the hydrophobic moiety shows that the beta <--> alpha acyl chain melting is strongly expressed only for the stretching and scissoring modes of the methylene groups. Vibrational groups originating from the interface region sense the acyl chain melting well (ester carbonyl bands) or only weakly (amide bands), and those resulting from the pure polar moiety not at all. From the x-ray data, the existence of lamellar (L), different cubic, and, for lipid A and LPS R595, also inverted hexagonal (HII) structures could be proven in the temperature range 20-80 degrees C with cubic <--> cubic and cubic <--> HII transitions for the Mg(2+)-free and L <--> HII transitions for the Mg(2+)-containing samples. These structural transitions can be characterized most readily by specific changes of the vibrational bands resulting from the interface region: the ester carbonyl and the amide bands. The magnitude of the changes corresponds to that of the structural rearrangement, i.e., is highest for the L <--> HII, lower for the cubic <--> HII, and lowest for the cubic <--> cubic transitions. The structural transitions are only marginally expressed for vibrational bands of the hydrophobic moiety. Similarly, the band contours of vibrations from the hydrophilic region are no indicators of the structural reorientations except for the carboxylate bands of LPS Re. Particularly the stretching vibrations of the phosphate groups are nearly completely invariant; the absolute values of their half bandwidths, however, differ significantly for lipid A and LPS Re, which seems to be of biological relevance. The ability of IR spectroscopy to detect supramolecular changes also beyond the measurability by x-ray diffraction, i.e., at water contents > 95 to 99.5%, is demonstrated.     

Peptide-protein interaction markedly alters the functional properties of the catalytic subunit of aspartate transcarbamoylase.

Interaction of a 70-amino acid zinc-binding polypeptide from the regulatory chain of aspartate transcarbamoylase (ATCase) with the catalytic (C) subunit leads to dramatic changes in enzyme activity and affinity for ligand binding at the active sites. The complex between the polypeptide (zinc domain) and wild-type C trimer exhibits hyperbolic kinetics in contrast to the sigmoidal kinetics observed with the intact holoenzyme. Moreover, the Scatchard plot for binding N-(phosphonacetyl)-L-aspartate (PALA) to the complex is linear with a Kd corresponding to that evaluated for the holoenzyme converted to the relaxed (R) state. Additional evidence that the binding of the zinc domain to the C trimer converts it to the R state was attained with a mutant form of ATCase in which Lys 164 in the catalytic chain is replaced by Glu. As shown previously (Newell, J.O. & Schachman, H.K., 1990, Biophys. Chem. 37, 183-196), this mutant holoenzyme, which exists in the R conformation even in the absence of active site ligands, has a 50-fold greater affinity for PALA than the free C subunit. Adding the zinc domain to the C trimer containing the Lys 164-->Glu substitution leads to a 50-fold enhancement in the affinity for the bisubstrate analog yielding a value of Kd equal to that for the holoenzyme. A different mutant ATCase containing the Gln 231 to Ile replacement was shown (Peterson, C.B., Burman, D.L., & Schachman, H.K., 1992, Biochemistry 31, 8508-8515) to be much less active as a holoenzyme than as the free C trimer. For this mutant holoenzyme, the addition of substrates does not cause its conversion to the R state. However, the addition of the zinc domain to the Gln 231-->Ile C trimer leads to a marked increase in enzyme activity, and PALA binding data indicate that the complex resembles the R state of the holoenzyme. This interaction leading to a more active conformation serves as a model of intergenic complementation in which peptide binding to a protein causes a conformational correction at a site remote from the interacting surfaces resulting in activation of the protein. This linkage was also demonstrated by difference spectroscopy using a chromophore covalently bound at the active site, which served as a spectral probe for a local conformational change. The binding of ligands at the active sites was shown also to lead to a strengthening of the interaction between the zinc domain and the C trimer.     

Deriving the intermediate spectra and photocycle kinetics from time-resolved difference spectra of bacteriorhodopsin. The simpler case of the recombinant D96N protein.

The bacteriorhodopsin photocycle contains more than five spectrally distinct intermediates, and the complexity of their interconversions has precluded a rigorous solution of the kinetics. A representation of the photocycle of mutated D96N bacteriorhodopsin near neutral pH was given earlier (Váró, G., and J. K. Lanyi. 1991. Biochemistry. 30:5008-5015) as BRhv-->K<==>L<==>M1-->M2--> BR. Here we have reduced a set of time-resolved difference spectra for this simpler system to three base spectra, each assumed to consist of an unknown mixture of the pure K, L, and M difference spectra represented by a 3 x 3 matrix of concentration values between 0 and 1. After generating all allowed sets of spectra for K, L, and M (i.e., M1 + M2) at a 1:50 resolution of the matrix elements, invalid solutions were eliminated progressively in a search based on what is expected, empirically and from the theory of polyene excited states, for rhodopsin spectra. Significantly, the average matrix values changed little after the first and simplest of the search criteria that disallowed negative absorptions and more than one maximum for the M intermediate. We conclude from the statistics that during the search the solutions strongly converged into a narrow region of the multidimensional space of the concentration matrix. The data at three temperatures between 5 and 25 degrees C yielded a single set of spectra for K, L, and M; their fits are consistent with the earlier derived photocycle model for the D96N protein.     

The allosteric activator ATP induces a substrate-dependent alteration of the quaternary structure of a mutant aspartate transcarbamoylase impaired in active site closure.

Aspartate transcarbamoylase from Escherichia coli shows homotropic cooperativity for aspartate as well as heterotropic regulation by nucleotides. Structurally, it consists of two trimeric catalytic subunits and three dimeric regulatory subunits, each chain being comprised of two domains. Glu-50 and Ser-171 are involved in stabilizing the closed conformation of the catalytic chain. Replacement of Glu-50 or Ser-171 by Ala in the holoenzyme has been shown previously to result in marked decreases in the maximal observed specific activity, homotropic cooperativity, and affinity for aspartate (Dembowski NJ, Newton CJ, Kantrowitz ER, 1990, Biochemistry 29:3716-3723; Newton CJ, Kantrowitz ER, 1990, Biochemistry 29:1444-1451). We have constructed a double mutant enzyme combining both mutations. The resulting Glu-50/ser-171-->Ala enzyme is 9-fold less active than the Ser-171-->Ala enzyme, 69-fold less active than the Glu-50-->Ala enzyme, and shows 1.3-fold and 1.6-fold increases in the [S]0.5Asp as compared to the Ser-171-->Ala and Glu-50-->Ala enzymes, respectively. However, the double mutant enzyme exhibits some enhancement of homotropic cooperativity with respect to aspartate, relative to the single mutant enzymes. At subsaturating concentrations of aspartate, the Glu-50/Ser-171 -->Ala enzyme is activated less by ATP than either the Glu-50-->Ala or Ser-171-->Ala enzyme, whereas CTP inhibition is intermediate between that of the two single mutants. As opposed to the wild-type enzyme, the Glu-50/Ser-171 -->Ala enzyme is activated by ATP and inhibited by CTP at saturating concentrations of aspartate. Structural analysis of the Ser-171-->Ala and Glu-50/Ser-171-->Ala enzymes by solution X-ray scattering indicates that both mutants exist in the same T quaternary structure as the wild-type enzyme in the absence of ligands, and in the same R quaternary structure in the presence of saturating N-(phosphonoacetyl)-L-aspartate. However, saturating concentrations of carbamoyl phosphate and succinate are unable to convert a significant fraction of either mutant enzyme population to the R quaternary structure, as has been observed previously for the Glu-50-->Ala enzyme. The curves for both the Ser-171-->Ala and Glu-50/Ser-171-->Ala enzymes obtained in the presence of substoichiometric amounts of PALA are linear combinations of the two extreme T and R states. The structural consequences of nucleotide binding to these two enzymes were also investigated. Most surprisingly, the direction and amplitude of the effect of ATP upon the double mutant enzyme were shown to vary depending upon the substrate analogue used.     

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.     

The apomyoglobin folding pathway revisited: structural heterogeneity in the kinetic burst phase intermediate

Extensive analysis of accurate quench-flow hydrogen exchange results indicates that the burst phase kinetic intermediate in the folding of apomyoglobin (apoMb) from urea is structurally heterogeneous. The structural variability is associated with the partial folding of the E helix during the burst phase (<6.4ms) of the folding process. Analysis of the effects of exchange-out of amide proton labels during the labeling pulse ( approximately pH 10) of the quench-flow process indicates that three of the amide protons in the E helix are in fact largely protected in the burst phase of folding, while the remainder of the E helix has a substantial complement of amide protons that show biphasic kinetics, i.e. are protected partly during the burst phase and partly during the slow phase of folding. The locations of these amide protons can be used to map the sites of structural heterogeneity in the kinetic molten globule. These sites include, besides the E helix, the ends of the A and B helices and part of the C helix. Our results give significant support to the hypothesis that the kinetic molten globule intermediate of apoMb is native-like.     

Transient association of the first intermediate during the refolding of bovine carbonic anhydrase B.

Many proteins which aggregate during refolding may form transiently populated aggregated states which do not reduce the final recovery of active species. However, the transient association of a folding intermediate will result in reduced refolding rates if the dissociation process occurs slowly. Previous studies on the refolding and aggregation of bovine carbonic anhydrase B (CAB) have shown that the molten globule first intermediate on the CAB folding pathway will form dimers and trimers prior to the formation of large aggregates (Cleland, J. L.; Wang, D. I. C. Biochemistry 1990, 29, 11072-11078; Cleland, J. L.; Wang, D. I. C. In Protein Refolding; Georgiou, G., De-Bernardez-Clark, E., Eds.; ACS Symposium Series 470; American Chemical Society: Washington, DC, 1991; pp 169-179). Refolding of CAB from 5 M guanidine hydrochloride (GuHCl) was achieved at conditions ([CAB]f = 10-33 microM, [GuHCl]f = 1.0 M) which allowed complete recovery of active protein as well as the formation of a transiently populated dimer of the molten globule intermediate on the refolding pathway. A kinetic analysis of CAB refolding provided insight into the mechanism of the association phenomenon. Using the kinetic results, a model of the refolding with transient association was constructed. By adjusting a single variable, the dimer dissociation rate constant, the model prediction fit both the experimentally determined active protein and dimer concentrations. The model developed in this analysis should also be applicable to the refolding of proteins which have been observed to form aggregates during refolding. In particular, the transient association of hydrophobic folding intermediates may also occur during the refolding of other proteins.(ABSTRACT TRUNCATED AT 250 WORDS)     

Illuminating the Off-Pathway Nature of the Molten Globule Folding Intermediate of an α-β Parallel Protein

Partially folded protein species transiently form during folding of most proteins. Often, these species are molten globules, which may be on- or off-pathway to the native state. Molten globules are ensembles of interconverting protein conformers that have a substantial amount of secondary structure, but lack virtually all tertiary side-chain packing characteristics of natively folded proteins. Due to solvent-exposed hydrophobic groups, molten globules are prone to aggregation, which can have detrimental effects on organisms. The molten globule observed during folding of the 179-residue apoflavodoxin from Azotobacter vinelandii is off-pathway, as it has to unfold before native protein can form. Here, we study folding of apoflavodoxin and characterize its molten globule using fluorescence spectroscopy and Förster Resonance Energy Transfer (FRET). Apoflavodoxin is site-specifically labeled with fluorescent donor and acceptor dyes, utilizing dye-inaccessibility of Cys69 in cofactor-bound protein. Donor (i.e., Alexa Fluor 488) is covalently attached to Cys69 in all apoflavodoxin variants used. Acceptor (i.e., Alexa Fluor 568) is coupled to Cys1, Cys131 and Cys178, respectively. Our FRET data show that apoflavodoxin’s molten globule forms in a non-cooperative manner and that its N-terminal 69 residues fold last. In addition, striking conformational differences between molten globule and native protein are revealed, because the inter-label distances sampled in the 111-residue C-terminal segment of the molten globule are shorter than observed for native apoflavodoxin. Thus, FRET sheds light on the off-pathway nature of the molten globule during folding of an α-β parallel protein.     

Kinetics of cytochrome C folding: atomically detailed simulations

The vast range of time scales (from nanoseconds to seconds) during protein folding is a challenge for experiments and computations. To make concrete predictions on folding mechanisms, atomically detailed simulations of protein folding, using potentials derived from chemical physics principles, are desired. However, due to their computational complexity, straightforward molecular dynamics simulations of protein folding are impossible today. An alternative algorithm is used that makes it possible to compute approximate atomically detailed long time trajectories (the Stochastic Difference Equation in Length). This algorithm is used to compute 26 atomically detailed folding trajectories of cytochrome c (a millisecond process). The early collapse of the protein chain (with marginal formation of secondary structure), and the earlier formation of the N and C helices (compare to the 60's helix) are consistent with the experiment. The existence of an energy barrier upon entry to the molten globule is examined as well. In addition to (favorable) comparison to experiments, we show that non-native contacts drive the formation of the molten globule. In contrast to popular folding models, the non-native contacts do not form off-pathway kinetic traps in cytochrome c.     

Fast compaction of alpha-lactalbumin during folding studied by stopped-flow X-ray scattering

To monitor the fast compaction process during protein folding, we have used a stopped-flow small-angle X-ray scattering technique combined with a two-dimensional charge-coupled device-based X-ray detector that makes it possible to improve the signal-to-noise ratio of data dramatically, and measured the kinetic refolding reaction of alpha-lactalbumin. The results clearly show that the radius of gyration and the overall shape of the kinetic folding intermediate of alpha-lactalbumin are the same as those of the molten globule state observed at equilibrium. Thus, the identity between the kinetic folding intermediate and the equilibrium molten globule state is firmly established. The present results also suggest that the folding intermediate is more hydrated than the native state and that the hydrated water molecules are dehydrated when specific side-chain packing is formed during the change from the molten globule to the native state.     

Folding of bovine growth hormone is consistent with the molten globule hypothesis

Previous results from equilibrium and kinetic studies of the folding of bovine growth hormone (bGH) have demonstrated that bGH does not follow a simple two-step folding mechanism. These results are summarized and interpreted according to the "molten globule" model. The molten globule state of bGH is characterized as a folding intermediate which is largely alpha-helical, retains a compact hydrodynamic radius, has packing of the aromatic side chains that is similar to the unfolded state, and possesses a solvent-exposed hydrophobic surface along helix 106-127 that readily leads to association.