Theoretical study of the conformation of the lipoamide arm in a mutant H protein

The lipoamide arm of the H protein plays a pivotal role in the catalytic cycle of the glycine decarboxylase complex (GDC) by being successively methylamine loaded (Hmet), reduced (Hred), and oxidized (Hox). In a previous study, we calculated free-energy surfaces as a function of the lipoamide arm position of the three forms of the wild-type protein and found close agreement with the available experimental data. Our simulations, together with crystallographic and NMR data, showed that the methylamine-loaded arm is locked in a cavity by interaction with Ser12, Glu14, and Asp67. In this work, we investigate the behavior of the methylamine-loaded form of a mutant H protein (HEA) where Glu14 has been replaced by Ala. We find that the arm can still be held in the cavity but that the energy barrier to release of the arm is halved from approximately 40 kcal mol(-1) for Hmet to approximately 12 kcal mol(-1) for HEA. To compensate for the loss of Glu14, the methylamine group shifts toward Ser66 in the mutant form. These results provide a structural basis for the equilibrium between the loaded and the unloaded forms of the arm observed by Gueguen et al. (Gueguen et al., J Biol Chem 1999;274:26344-26352) in HEA.     

Combined structural and biochemical analysis of the H-T complex in the glycine decarboxylase cycle: evidence for a destabilization mechanism of the H-protein

The lipoate containing H-protein plays a pivotal role in the catalytic cycle of the glycine decarboxylase complex (GDC), undergoing reducing methylamination, methylene transfer, and oxidation. The transfer of the CH(2) group is catalyzed by the T-protein, which forms a 1:1 complex with the methylamine-loaded H-protein (Hmet). The methylamine group is then deaminated and transferred to the tetrahydrofolate-polyglutamate (H(4)FGlu(n)) cofactor of T-protein, forming methylenetetrahydrofolate-polyglutamate. The methylamine group is buried inside the protein structure and highly stable. Experimental data show that the H(4)FGlu(n) alone does not induce transfer of the methylene group, and molecular modeling also indicates that the reaction cannot take place without significant structural perturbations of the H-protein. We have, therefore, investigated the effect of the presence of the T-protein on the stability of Hmet. Addition of T-protein without H(4)FGlu(n) greatly increases the rate of the unloading reaction of Hmet, reducing the activation energy by about 20 kcal mol(-1). Differences of the (1)H and (15)N chemical shifts of the H-protein in its isolated form and in the complex with the T-protein show that the interaction surface for the H-protein is localized on one side of the cleft where the lipoate arm is positioned. This suggests that the role of the T-protein is not only to locate the tetrahydrofolate cofactor in a position favorable for a nucleophilic attack on the methylene carbon but also to destabilize the H-protein in order to facilitate the unlocking of the arm and initiate the reaction.     

Theoretical study of the conformation of the H-protein lipoamide arm as a function of its terminal group.

The glycine decarboxylase complex consists of four different proteins (the L-, P-, H-, and T-proteins). The H-protein plays a central role in communication among the other enzymes, as its lipoamide arm interacts successively with each of the components of the complex. The crystal structures of two states of the H-protein have been resolved: the oxidized form, Hox at 2 A and the methylamine-loaded form, Hmet at 2.2 A. However, the position of the arm for the reduced form, Hred, is still unknown. We have performed numerical free-energy calculations in order to better understand the differences in the structures and to elucidate the conformation of the arm in Hred. The results of the simulations are in agreement with the crystallographic results, as the minima of the free energy surface for Hox and Hmet correspond to the crystal structures. For Hred, we observe a single minimum in which the arm is on the surface of the H-protein, close to its position in the Hox structure. In all of our simulations, the lower, lysine portion of the arm remains bound to the protein, which substantially reduces the number of accessible arm configurations. An analysis of the stability of Hmet in the cavity shows that electrostatic interactions are crucial for locking the arm in the bottom of the cavity, especially near Glu14. In addition, the analysis shows that there is a water molecule, also observed in the crystallographic structure, that binds to the arm's terminal NH3+ group and helps to fix it in the cavity. In conclusion, because of the close agreement of the results of our calculations with the available experimental evidence, we are able to suggest a structural basis for the observed behavior. Proteins 1999;36: 228-237.     

NMR structural and dynamic characterization of the acid-unfolded state of apomyoglobin provides insights into the early events in protein folding

Apomyoglobin forms a denatured state under low-salt conditions at pH 2.3. The conformational propensities and polypeptide backbone dynamics of this state have been characterized by NMR. Nearly complete backbone and some side chain resonance assignments have been obtained, using a triple-resonance assignment strategy tailored to low protein concentration (0.2 mM) and poor chemical shift dispersion. An estimate of the population and location of residual secondary structure has been made by examining deviations of (13)C(alpha), (13)CO, and (1)H(alpha) chemical shifts from random coil values, scalar (3)J(HN,H)(alpha) coupling constants and (1)H-(1)H NOEs. Chemical shifts constitute a highly reliable indicator of secondary structural preferences, provided the appropriate random coil chemical shift references are used, but in the case of acid-unfolded apomyoglobin, (3)J(HN,H)(alpha) coupling constants are poor diagnostics of secondary structure formation. Substantial populations of helical structure, in dynamic equilibrium with unfolded states, are formed in regions corresponding to the A and H helices of the folded protein. In addition, the deviation of the chemical shifts from random coil values indicates the presence of helical structure encompassing the D helix and extending into the first turn of the E helix. The polypeptide backbone dynamics of acid-unfolded apomyoglobin have been investigated using reduced spectral density function analysis of (15)N relaxation data. The spectral density J(omega(N)) is particularly sensitive to variations in backbone fluctuations on the picosecond to nanosecond time scale. The central region of the polypeptide spanning the C-terminal half of the E helix, the EF turn, and the F helix behaves as a free-flight random coil chain, but there is evidence from J(omega(N)) of restricted motions on the picosecond to nanosecond time scale in the A and H helix regions where there is a propensity to populate helical secondary structure in the acid-unfolded state. Backbone fluctuations are also restricted in parts of the B and G helices due to formation of local hydrophobic clusters. Regions of restricted backbone flexibility are generally associated with large buried surface area. A significant increase in J(0) is observed for the NH resonances of some residues located in the A and G helices of the folded protein and is associated with fluctuations on a microsecond to millisecond time scale that probably arise from transient contacts between these distant regions of the polypeptide chain. Our results indicate that the equilibrium unfolded state of apomyoglobin formed at pH 2.3 is an excellent model for the events that are expected to occur in the earliest stages of protein folding, providing insights into the regions of the polypeptide that spontaneously undergo local hydrophobic collapse and sample nativelike secondary structure.     

Native and non-native secondary structure and dynamics in the pH 4 intermediate of apomyoglobin

The partly folded state of apomyoglobin at pH 4 represents an excellent model for an obligatory kinetic folding intermediate. The structure and dynamics of this intermediate state have been extensively examined using NMR spectroscopy. Secondary chemical shifts, (1)H-(1)H NOEs, and amide proton temperature coefficients have been used to probe residual structure in the intermediate state, and NMR relaxation parameters T(1) and T(2) and ?(1)H?-(15)N NOE have been analyzed using spectral densities to correlate motion of the polypeptide chain with these structural observations. A significant amount of helical structure remains in the pH 4 state, indicated by the secondary chemical shifts of the (13)C(alpha), (13)CO, (1)H(alpha), and (13)C(beta) nuclei, and the boundaries of this helical structure are confirmed by the locations of (1)H-(1)H NOEs. Hydrogen bonding in the structured regions is predominantly native-like according to the amide proton chemical shifts and their temperature dependence. The locations of the A, G, and H helix segments and the C-terminal part of the B helix are similar to those in native apomyoglobin, consistent with the early, complete protection of the amides of residues in these helices in quench-flow experiments. These results confirm the similarity of the equilibrium form of apoMb at pH 4 and the kinetic intermediate observed at short times in the quench-flow experiment. Flexibility in this structured core is severely curtailed compared with the remainder of the protein, as indicated by the analysis of the NMR relaxation parameters. Regions with relatively high values of J(0) and low values of J(750) correspond well with the A, B, G, and H helices, an indication that nanosecond time scale backbone fluctuations in these regions of the sequence are restricted. Other parts of the protein show much greater flexibility and much reduced secondary chemical shifts. Nevertheless, several regions show evidence of the beginnings of helical structure, including stretches encompassing the C helix-CD loop, the boundary of the D and E helices, and the C-terminal half of the E helix. These regions are clearly not well-structured in the pH 4 state, unlike the A, B, G, and H helices, which form a native-like structured core. However, the proximity of this structured core most likely influences the region between the B and F helices, inducing at least transient helical structure.     

Characterization of the binding interface between the copper chaperone Atx1 and the first cytosolic domain of Ccc2 ATPase

The interaction of the copper chaperone Atx1 and the first cytosolic domain of Ccc2 ATPase, Ccc2a, was investigated by NMR in solution. In particular, a solution of Cu(I)-15NAtx1 was titrated with apo-Ccc2a, and, vice versa, a solution of Cu(I)-15NCcc2a was titrated with apo-Atx1. By following the 15N and 1H chemical shifts, a new species is detected in both experiments. This species is the same in both titrations and is in fast exchange with the parent species on the NMR time scale. Nuclear relaxation data are consistent with the formation of an adduct. Judging from the nuclear Overhauser effect spectroscopy patterns, the structure of Cu(I)-15NCcc2a in the presence of apo-Atx1 is not significantly altered, whereas Cu(I)-15NAtx1 in the presence of apo-Ccc2a experiences some changes with respect to both the apoproteins and the Cu(I)-loaded proteins. The structure of the Cu(I)-15NAtx1 moiety in the adduct was obtained from 1137 nuclear Overhauser effects to a final root mean square deviation to the mean structure of 0.76 +/- 0.13 A for the backbone and 1.11 +/- 0.11 A for the heavy atoms. 15N and 1H chemical shifts suggest the regions of interaction that, together with independent information, allow a structural model of the adduct to be proposed. The apo form of Atx1 displays significant mobility in loops 1 and 5, the N-terminal part of helix alpha1, and the C-terminal part of helix alpha2 on the ms-micros time scale. These regions correspond to the metal binding site. Such mobility is largely reduced in the free Cu(I)-Atx1 and in the adduct with apo-Ccc2a. The analogous mobility of Ccc2a in both Cu(I) and apo forms is reduced with respect to Atx1. Such an adduct is relevant as a structural and kinetic model for copper transfer from Atx1 to Ccc2a in physiological conditions.     

Structural investigations of chromatin core protein by nuclear magnetic resonance

A complex derived from chromatin containing one molecule of each of histones H2A, H2B, H3, and H4, termed core protein, was studied by 13C and 1H nuclear magnetic resonance. 13C line widths, when analyzed and compared with those of native and thermally unfolded representative globular proteins, showed that regions of the core protein possess considerable mobility. Studies of Calpha and Cbeta line widths, and Calpha spin-spin relaxation times, show that this mobility arises from sections of random-coil polypeptide. It is argued that these regions are N-terminal "tails", attached to C-terminal globular polypeptides. The 270-MHz 1H nuclear magnetic resonance spectrum shows numerous ring current shifted resonances, indicating that the C-terminal globular domain has a precise tertiary structure. The globular domain most likely forms the histone "core" of the chromatin monomer particle, whilst the basic tails probably wind around the grooves of the double helix, enabling the basic side chains to interact with the DNA phosphate groups. Some biological implications of this model are considered.     

Internal dynamics of the tryptophan repressor (TrpR) and two functionally distinct TrpR variants, L75F-TrpR and A77V-TrpR, in their l-Trp-bound forms

Backbone amide dynamics of the Escherichia coli tryptophan repressor protein (WT-TrpR) and two functionally distinct variants, L75F-TrpR and A77V-TrpR, in their holo (l-tryptophan corepressor-bound) form have been characterized using (15)N nuclear magnetic resonance (NMR) relaxation. The three proteins possess very similar structures, ruling out major conformational differences as the source of their functional differences, and suggest that changes in protein flexibility are at the origin of their distinct functional properties. Comparison of site specific (15)N-T(1), (15)N-T(2), (15)N-{(1)H} nuclear Overhauser effect, reduced spectral density, and generalized order (S(2)) parameters indicates that backbone dynamics in the three holo-repressors are overall very similar with a few notable and significant exceptions for backbone atoms residing within the proteins' DNA-binding domain. We find that flexibility is highly restricted for amides in core α-helices (i.e., helices A-C and F), and a comparable "stiffening" is observed for residues in the DNA recognition helix (helix E) of the helix D-turn-helix E (HTH) DNA-binding domain of the three holo-repressors. Unexpectedly, amides located in helix D and in adjacent turn regions remain flexible. These data support the concept that residual flexibility in TrpR is essential for repressor function, DNA binding, and molecular recognition of target operators. Comparison of the (15)N NMR relaxation parameters of the holo-TrpRs with those of the apo-TrpRs indicates that the single-point amino acid substitutions, L75F and A77V, perturb the flexibility of backbone amides of TrpR in very different ways and are most pronounced in the apo forms of the three repressors. Finally, we present these findings in the context of other DNA-binding proteins and the role of protein flexibility in molecular recognition.     

Backbone dynamics of calmodulin studied by 15N relaxation using inverse detected two-dimensional NMR spectroscopy: the central helix is flexible.

The backbone dynamics of Ca(2+)-saturated recombinant Drosophila calmodulin has been studied by 15N longitudinal and transverse relaxation experiments, combined with 15N(1H) NOE measurements. Results indicate a high degree of mobility near the middle of the central helix of calmodulin, from residue K77 through S81, with order parameters (S2) in the 0.5-0.6 range. The anisotropy observed in the motion of the two globular calmodulin domains is much smaller than expected on the basis of hydrodynamic calculations for a rigid dumbbell type structure. This indicates that, for the purposes of 15N relaxation, the tumbling of the N-terminal (L4-K77) and C-terminal (E82-S147) lobes of calmodulin is effectively independent. A slightly shorter motional correlation time (tau c approximately 6.3 ns) is obtained for the C-terminal domain compared to the N-terminal domain (tau c approximately 7.1 ns), in agreement with the smaller size of the C-terminal domain. A high degree of mobility, with order parameters of approximately 0.5, is also observed in the loop that connects the first with the second EF-hand type calcium binding domain and in the loop connecting the third and fourth calcium binding domain.     

Conformational dynamics of yeast calmodulin in the Ca(2+)-bound state probed using NMR relaxation dispersion

Most calmodulin (CaM) in apo and Ca(2+)-bound states show a dumb-bell-like structure, involving the N- and C-terminal domains, connected with a flexible linker. However, Ca(2+)-bound yeast calmodulin (yCaM) takes on a unique globular structure; the target-binding site of this protein is autoinhibited. We applied NMR relaxation dispersion experiments to yCaM in the Ca(2+)-bound state. The amide (15)N and (1)H(N) relaxation dispersion profiles indicated the presence of conformational dynamics for specific residues at the interface between the N- and C-terminal domains. We conclude that these conformational dynamics were derived from the mobility of the C-terminal domain.     

Assignment of amide 1H and 15N NMR resonances in detergent-solubilized M13 coat protein: a model for the coat protein dimer.

The major coat protein of the filamentous coliphage M13 is a 50-residue integral membrane protein. Detergent-solubilized M13 coat protein is a promising candidate for structure determination by nuclear magnetic resonance methods as the protein can be prepared in large quantities and the protein-containing micelle is reasonably small. Under the conditions of our experiments, SDS-bound coat protein exists as a dimer with an apparent molecular weight of 27,000. Broad lines and poor resolution in the 1H spectrum have led us to adopt an 15N-directed approach, in which the coat protein was labeled both uniformly with 15N and selectively with [alpha-15N]alanine, -glycine, -valine, -leucine, -isoleucine, phenylalanine, -lysine, -tyrosine, and -methionine. Nitrogen resonances were assigned as far as possible using carboxypeptidase digestion, double-labeling, and an independent knowledge of the amide proton exchange rates determined from neighboring assigned 13C-labeled carbonyl carbons. 1H/15N heteronuclear multiple quantum coherence (HMQC) spectroscopy of both uniform and site-selectively-labeled proteins subsequently correlated amide nitrogen with amide proton chemical shifts, and the assignments were completed sequentially from homonuclear NOESY and HMQC-NOESY spectra. The most slowly exchanging amide protons were shown to occur in a continuous stretch extending from methionine-28 to phenylalanine-42. This sequence includes most of the resonances of the hydrophobic core, although it is shifted toward the C-terminal end of the protein. Strong NH to NH (i,i+1) nuclear Overhauser enhancements are a feature of the coat protein, which appears to be largely helical. Between 20 and 25 residues give rise to 2 juxtaposed resonances which can be seen clearly in the HMQC spectrum of uniform 15N-labeled coat protein. These residues are concentrated in a region extending from the beginning of the membrane-spanning sequence through to the disordered region near the C-terminus. We propose that dodecyl sulfate-bound M13 coat protein consists of two independent domains, an N-terminal helix which is in a state of moderately fast dynamic flux and a long, stable, C-terminal membrane-spanning helix, which undergoes extensive interactions with a second monomer. Amide 1H chemical shifts are consistent with this picture; in addition, a marked periodicity is observed at the C-terminal end of the molecule.     

Hsc70-interacting HPD loop of the J domain of polyomavirus T antigens fluctuates in ps to ns and micros to ms

The backbone dynamics of the J domain from polyomavirus T antigens have been investigated using 15N NMR relaxation and molecular dynamics simulation. Model-free relaxation analysis revealed picosecond to nanosecond motions in the N terminus, the I-II loop, the C-terminal end of helix II through the HPD loop to the beginning of helix III, and the C-terminal end of helix III to the C terminus. The backbone dynamics of the HPD loop and termini are dominated by motions with moderately large amplitudes and correlation times of the order of a nanosecond or longer. Conformational exchange on the microsecond to millisecond timescale was identified in the HPD loop, the N and C termini, and the I-II loop. A 9.7ns MD trajectory manifested concerted swings of the HPD loop. Transitions between major and minor conformations of the HPD loop featured distinct patterns of change in backbone dihedral angles and hydrogen bonds. Fraying of the C-terminal end of helix II and the N-terminal end of helix III correlated with displacements of the HPD loop. Correlation of crankshaft motions of Gly46 and Gly47 with the collective motions of the HPD loop suggested an important role of the two glycine residues in the mobility of the loop. Fluctuations of the HPD loop correlated with relative reorientation of side-chains of Lys35 and Asp44 that interact with Hsc70.     

Direct demonstration of the flexibility of the glycosylated proline-threonine linker in the Cellulomonas fimi Xylanase Cex through NMR spectroscopic analysis

The modular xylanase Cex (or CfXyn10A) from Cellulomonas fimi consists of an N-terminal catalytic domain and a C-terminal cellulose-binding domain, joined by a glycosylated proline-threonine (PT) linker. To characterize the conformation and dynamics of the Cex linker and the consequences of its modification, we have used NMR spectroscopy to study full-length Cex in its nonglycosylated ( approximately 47 kDa) and glycosylated ( approximately 51 kDa) forms. The PT linker lacks any predominant structure in either form as indicated by random coil amide chemical shifts. Furthermore, heteronuclear (1)H-(15)N nuclear Overhauser effect relaxation measurements demonstrate that the linker is flexible on the ns-to-ps time scale and that glycosylation partially dampens this flexibility. The catalytic and cellulose-binding domains also exhibit identical amide chemical shifts whether in isolation or in the context of either unmodified or glycosylated full-length Cex. Therefore, there are no noncovalent interactions between the two domains of Cex or between either domain and the linker. This conclusion is supported by the distinct (15)N relaxation properties of the two domains, as well as their differential alignment within a magnetic field by Pf1 phage particles. These data demonstrate that the PT linker is a flexible tether, joining the structurally independent catalytic and cellulose-binding domains of Cex in an ensemble of conformations; however, more extended forms may predominate because of restrictions imparted by the alternating proline residues. This supports the postulate that the binding-domain anchors Cex to the surface of cellulose, whereas the linker provides flexibility for the catalytic domain to hydrolyze nearby hemicellulose (xylan) chains.     

The distal cytoplasmic tail of the influenza A M2 protein dynamically extends from the membrane

The influenza A M2 protein is a multifunctional membrane-associated homotetramer that orchestrates several essential events in the viral infection cycle. The monomeric subunits of the M2 homotetramer consist of an N-terminal ectodomain, a transmembrane domain, and a C-terminal cytoplasmic domain. The transmembrane domain forms a four-helix proton channel that promotes uncoating of virions upon host cell entry. The membrane-proximal region of the C-terminal domain forms a surface-associated amphipathic helix necessary for viral budding. The structure of the remaining ~34 residues of the distal cytoplasmic tail has yet to be fully characterized despite the functional significance of this region for influenza infectivity. Here, we extend structural and dynamic studies of the poorly characterized M2 cytoplasmic tail. We used SDSL-EPR to collect site-specific information on the mobility, solvent accessibility, and conformational properties of residues 61–70 of the full-length, cell-expressed M2 protein reconstituted into liposomes. Our analysis is consistent with the predominant population of the C-terminal tail dynamically extending away from the membranes surface into the aqueous medium. These findings provide insight into the hypothesis that the C-terminal domain serves as a sensor that regulates how M2 protein participates in critical events in the viral infection cycle.     

A further clue to understanding the mobility of mitochondrial yeast cytochrome c: a (15)N T1rho investigation of the oxidized and reduced species.

A new approach was developed to overproduce 15N-enriched yeast iso-1-cytochrome c in the periplasm of Escherichia coli in order to perform a study of the motions in the ms-micros time scale on the oxidized and reduced forms through rotating frame 15N relaxation rates and proton/deuterium exchange studies. It is confirmed that the reduced protein is rather rigid whereas the oxidized species is more flexible. The regions of the protein that display increased internal mobility upon oxidation are easily identified by the number of residues experiencing conformational equilibria and by their exchange rates. These data complement the information already available in the literature and provide a comprehensive picture of the mobility in the protein. In particular, oxidation mobilizes the loop containing Met80 and, through specific contacts, affects the mobility of helix 3 and possibly of helix 5, and of a section of protein connecting the heme propionates to helix 2. The relevance of internal motions to molecular recognition and to the early steps of the unfolding process of the oxidized species is also discussed. In agreement with the reported data, subnanosecond mobility is found to be less informative than the ms-micros with respect to redox dependent properties.     

The C terminus of apocytochrome b562 undergoes fast motions and slow exchange among ordered conformations resembling the folded state

The present work describes the dynamics of the apo form of cytochrome b(562), a small soluble protein consisting of 106 amino acid residues [Itagaki, E., and Hager, L. P. (1966) J. Biol. Chem. 241, 3687-3695]. The presence of exchange in the millisecond time scale is demonstrated for the last part of helix IV (residues 95-105 in the holo form). The chemical shift index analysis [Wishart, D. S., and Sykes, B. D. (1994) J. Biomol. NMR 4, 171-180] based on H(alpha), C(alpha), C(beta), and C' chemical shifts suggests a larger helical content than shown in the NMR structure based on NOEs. These results indicate the presence of helical-like conformations participating in the exchange process. This hypothesis is consistent with amide deuterium exchange rates and the presence of some hydrogen bonds identified from amide chemical shift temperature coefficients [Baxter, N. J., and Williamson, M. P. (1997) J. Biomol. NMR 9, 359-369]. (15)N relaxation indicates limited mobility for the amide protons of this part of the helix in the picosecond time scale. A 30 ns stochastic dynamics simulation shows small fluctuations around the helical conformation on this time scale. These fluctuations, however, do not result in a significant decrease of the calculated order parameters which are consistent with the experimental (15)N relaxation data. These results resolve an apparent discrepancy in the NMR structures between the disorder observed in helix IV due to a lack of NOEs and the secondary structure predictions based on H(alpha) chemical shifts [Feng, Y., Wand, A. J., and Sligar, S. G. (1994) Struct. Biol. 1, 30-35].     

Effects of varying the local propensity to form secondary structure on the stability and folding kinetics of a rapid folding mixed alpha/beta protein: characterization of a truncation mutant of the N-terminal domain of the ribosomal protein L9.

The N-terminal domain of the ribosomal protein L9 forms a split betaalphabeta structure with a long C-terminal helix. The folding transitions of a 56 residue version of this protein have previously been characterized, here we report the results of a study of a truncation mutant corresponding to residues 1-51. The 51 residue protein adopts the same fold as the 56 residue protein as judged by CD and two-dimensional NMR, but it is less stable as judged by chemical and thermal denaturation experiments. Studies with synthetic peptides demonstrate that the C-terminal helix of the 51 residue version has very little propensity to fold in isolation in contrast to the C-terminal helix of the 56 residue variant. The folding rates of the two proteins, as measured by stopped-flow fluorescence, are essentially identical, indicating that formation of local structure in the C-terminal helix is not involved in the rate-limiting step of folding.     

Conformational and dynamic characterization of the molten globule state of an apomyoglobin mutant with an altered folding pathway.

Kinetic and equilibrium studies of apomyoglobin folding pathways and intermediates have provided important insights into the mechanism of protein folding. To investigate the role of intrinsic helical propensities in the apomyoglobin folding process, a mutant has been prepared in which Asn132 and Glu136 have been substituted with glycine to destabilize the H helix. The structure and dynamics of the equilibrium molten globule state formed at pH 4.1 have been examined using NMR spectroscopy. Deviations of backbone (13)C(alpha) and (13)CO chemical shifts from random coil values reveal high populations of helical structure in the A and G helix regions and in part of the B helix. However, the H helix is significantly destabilized compared to the wild-type molten globule. Heteronuclear [(1)H]-(15)N NOEs show that, although the polypeptide backbone in the H helix region is more flexible than in the wild-type protein, its motions are restricted by transient hydrophobic interactions with the molten globule core. Quench flow hydrogen exchange measurements reveal stable helical structure in the A and G helices and part of the B helix in the burst phase kinetic intermediate and confirm that the H helix is largely unstructured. Stabilization of structure in the H helix occurs during the slow folding phases, in synchrony with the C and E helices and the CD region. The kinetic and equilibrium molten globule intermediates formed by N132G/E136G are similar in structure. Although both the wild-type apomyoglobin and the mutant fold via compact helical intermediates, the structures of the intermediates and consequently the detailed folding pathways differ. Apomyoglobin is therefore capable of compensating for mutations by using alternative folding pathways within a common basic framework. Tertiary hydrophobic interactions appear to play an important role in the formation and stabilization of secondary structure in the H helix of the N132G/E136G mutant. These studies provide important insights into the interplay between secondary and tertiary structure formation in protein folding.