Investigating metal-binding in proteins by nuclear magnetic resonance.

Metal ions play a key role for the function of many proteins. The interaction of the metal ion with the protein and its involvement in the function of the protein vary widely. In some proteins, the metal ion is bound tightly to the ligand residues and may be the key player in the function of the protein, as in the case of blue copper proteins. In other proteins, the metal ion is bound only temporarily and loosely to the protein, as in the case of some metalloenzymes and other proteins where the metal ion acts as a cofactor necessary for the function of the protein. Such proteins are often known as metal ion-activated proteins. The review focuses on recent nuclear magnetic resonance (NMR) studies of a series of metal-dependent proteins and the characterization of the metal-binding sites. In particular, we focus on NMR techniques for studying metal binding to proteins such as chemical shift mapping, paramagnetic NMR and changes in backbone dynamics upon metal binding.

The (ProB27, ThrB28) human insulin analogue is more potent and more rapidly absorbed from subcutis than human insulin.

OBJECTIVE: To test the physiological properties of human insulin in which the amino acids Thr (B27) and Pro (B28) are interchanged (PT insulin). This was hypothesised to prevent dimerisation and accelerate the absorption from s.c. tissue without altering the affinity for the insulin receptor. DESIGN: PT insulin was expressed in Pichia pastoris and processed in vitro. The purified compound was used for physiological investigations. METHODS: Receptor binding activity to insulin and IGF receptors was evaluated in a competition assay using iodinated PT insulin and recombinant receptors while growth induction properties were evaluated by thymidine incorporation. Absorption kinetics from pig subcutis was investigated by measuring the disappearance of iodinated PT insulin. The potency was evaluated by measuring the blood glucose lowering activity in mice. RESULTS: The absorption of PT insulin was accelerated compared with human insulin, although still slower than Asp (B28) insulin. Human and PT insulin had similar affinities for the human insulin receptor (K(d)=3.6 x 10(-12) vs 5.2 x 10(-12) mol/l) while the affinity for the IGF receptor was four times higher for PT insulin than for human insulin (K(d)=3.4 x 10(-8) vs 1.3 x 10(-7) mol/l). This resulted in a slightly higher DNA synthesis when assayed in intermediary insulin concentrations. The blood glucose lowering effect in mice exceeded the effect of human insulin (integral 0-60 min: 61.4+/-7 vs 30+/-4, n=6, P=0.046). CONCLUSIONS: PT insulin is absorbed faster and is more potent than human insulin. Although PT insulin stimulates growth more than human insulin, this will not prevent its use in the clinic, but the main interest will probably focus on investigations to clarify the paradox of full biological activity in connection with the recently described lack of structure in the B-chain.

Measuring the longitudinal NMR relaxation rates of fast relaxing nuclei using a signal eliminating relaxation filter.

A new experiment for selective determination of the relaxation rates of fast relaxing NMR signals is presented. The experiment is derived from the conventional inversion recovery experiment by substituting the 180 degrees inversion pulse of this experiment with a signal eliminating relaxation filter (SERF) consisting of three 180 degrees pulses separated by two variable delays, Delta1 and Delta2. The SERF experiment allows a selective suppression of signals with relaxation rates below a given limit while monitoring the relaxation of faster relaxing signals. The experiment was tested on a sample of 20% oxidized plastocyanin from Anabaena variabilis, where the fast exchange of an electron between the reduced (diamagnetic) and the oxidized (paramagnetic) form results in a series of average signals with widely different relaxation rates. To ensure an optimum extraction of information from the experimental data, the relaxation rates were obtained from the SERF experiment by a simultaneous analysis of all the FIDs of the experiment using a fast linear prediction model method developed previously. The reliability of the relaxation rates obtained from the SERF experiment was confirmed by a comparison of the rates with the corresponding rates obtained from a conventional inversion recovery experiment.

Flexibility and bioactivity of insulin: an NMR investigation of the solution structure and folding of an unusually flexible human insulin mutant with increased biological activity.

The structure and folding of a novel human insulin mutant, [Thr(B27) --> Pro, Pro(B28) --> Thr]insulin (PT insulin), in aqueous solution and in mixtures of water and 2,2,2-trifluoroethanol (TFE) have been studied by NMR spectroscopy. It was found that PT insulin has a highly flexible structure in pure water and is present in at least two different conformations, although with an overall tertiary structure similar to that of native insulin. Furthermore, the native helical structures are poorly defined. Surprisingly, the mutant has a biological activity about 50% higher than native insulin. In contrast, in TFE/water solution the mutant reveals a propensity of forming a well-defined structure at the secondary structure level, similar to monomeric native insulin. Thus, as shown by a detailed determination of the structure from 208 distance restraints and 52 torsion angle restraints by distance geometry, simulated annealing, and restrained energy minimization, the native insulin helices (A2-A7, A13-A19, and B10-B19) as well as the beta-turn (B20-B23) are formed in 35% TFE. However, the amount of tertiary structure is decreased significantly in TFE/water solution. The obtained results suggest that only an overall tertiary fold, as observed for PT insulin in pure water, is necessary for expressing the biological activity of insulin, as long as the molecule is flexible and retains the propensity to form the secondary structure required for its receptor binding. In contrast, a compact secondary structure, as found for native insulin in solution, is unnecessary for the biological activity. A model for the receptor binding of insulin is suggested that relates the increased bioactivity to the enhanced flexibility of the mutant.

Determination of the electron self-exchange rates of blue copper proteins by super-WEFT NMR spectroscopy.

An NMR approach for determining the electron self-exchange (ESE) rate constants in blue copper proteins is presented. The approach uses the paramagnetic relaxation enhancement of resonances in 1D 1H super-WEFT spectra of partly oxidized (paramagnetic) proteins. These spectra allow a more precise determination of the relevant paramagnetic linebroadenings than conventional 1D 1H spectra and, thus, permit a more detailed investigation of the applicability of the linebroadenings for determining the electron exchange rates. The approach was used to estimate the ESE rate constant of plastocyanin from Anabaena variabilis. It was found that, although the rate constant can be determined accurately from a series of resonances, precise but erroneous constants are obtained from the resonances of the copper-bound residues, unless a narrow splitting of these resonances caused by the presence of two conformations is taken into account. As demonstrated here, this complication can be overcome by a correct analysis of the paramagnetic broadening of the combined double signals. Because of the high resolution and specific sensitivity of the approach it should be generally applicable to estimate electron transfer rates, k, if the paramagnetic relaxation enhancement R2p of the resonances can be determined, and the conditions k << R2p or delta omega(p) >> k >> R2p are fulfilled, delta omega(p) being the frequency separation between corresponding diamagnetic and paramagnetic sites.

Unraveling the symmetry ambiguity in a hexamer: calculation of the R6 human insulin structure.

Crystallographic and NMR studies of insulin have revealed a highly flexible molecule with a range of different aggregation and structural states; the importance of these states for the function of the hormone is still unclear. To address this question, we have studied the solution structure of the insulin R6 symmetric hexamer using NMR spectroscopy. Structure determination of symmetric oligomers by NMR is complicated due to 'symmetry ambiguity' between intra- and intermonomer NOEs, and between different classes of intermonomer NOEs. Hence, to date, only two symmetric tetramers and one symmetric pentamer (VTB, B subunit of verotoxin) have been solved by NMR: there has been no other symmetric hexamer or higher-order oligomer. Recently, we reported a solution structure for R6 insulin hexamer. However, in that study, a crystal structure was used as a reference to resolve ambiguities caused by the threefold symmetry; the same method was used in solving VTB. Here, we have successfully recalculated R6 insulin using the symmetry-ADR method, a computational strategy in which ambiguities are resolved using the NMR data alone. Thus the obtained structure is a refinement of the previous R6 solution structure. Correlated motions in the final structural ensemble were analysed using a recently developed principal component method; this suggests the presence of two major conformational substates. The study demonstrates that the solution structure of higher-order symmetric oligomers can be determined unambiguously from NMR data alone, using the symmetry-ADR method. This success bodes well for future NMR studies of higher-order symmetric oligomers. The correlated motions observed in the structural ensemble suggest a new insight into the mechanism of phenol exchange and the T6 <--> R6 transition of insulin in solution.

Assignment of side-chain conformation using adiabatic energy mapping, free energy perturbation, and molecular dynamic simulations.

NMR spectroscopic analysis of the C-terminal Kunitz domain fragment (alpha3(VI)) from the human alpha3-chain of type VI collagen has revealed that the side chain of Trp21 exists in two unequally populated conformations. The major conformation (M) is identical to the conformation observed in the X-ray crystallographic structure, while the minor conformation (m) cannot structurally be resolved in detail by NMR due to insufficient NOE data. In the present study, we have applied: (1) rigid and adiabatic mapping, (2) free energy simulations, and (3) molecular dynamic simulations to elucidate the structure of the m conformer and to provide a possible pathway of the Trp21 side chain between the two conformers. Adiabatic energy mapping of conformations of the Trp21 side chain obtained by energy minimization identified two energy minima: One corresponding to the conformation of Trp21 observed in the X-ray crystallographic structure and solution structure of alpha3(VI) (the M conformation) and the second corresponding to the m conformation predicted by NMR spectroscopy. A transition pathway between the M and m conformation is suggested. The free-energy difference between the two conformers obtained by the thermodynamic integration method is calculated to 1.77+/-0.7 kcal/mol in favor of the M form, which is in good agreement with NMR results. Structural and dynamic properties of the major and minor conformers of the alpha3(VI) molecule were investigated by molecular dynamic. Essential dynamics analysis of the two resulting 800 ps trajectories reveals that when going from the M to the m conformation only small, localized changes in the protein structure are induced. However, notable differences are observed in the mobility of the binding loop (residues Thr13-Ile18), which is more flexible in the m conformation than in the M conformation. This suggests that the reorientation of Trp2 might influence the inhibitory activity against trypsin, despite the relative large distance between the binding loop and Trp21.

Solution structures of the R6 human insulin hexamer,.

The three-dimensional solution structure of the phenol-stabilized 36 kDa R6 insulin hexamer was determined by NMR spectroscopy and restrained molecular dynamics. The hexamer structures were derived using a stepwise procedure. Initially, 60 monomers were obtained by distance geometry from 665 NOE-derived distance restraints and three disulfide bridges. Subsequently, the hexamer structures were calculated by simulated annealing, using 30 hexamers constructed from the best 36 monomer structures as the starting models. The NMR data show that the aromatic ring of residue Phe(B25) can take two different orientations in the solution hexamer: one in which it points inward (molecule 1, about 90%) and one in which it points outward from the surface of the monomer (molecule 2, about 10%). Therefore, two hexamer structures were calculated: a symmetric hexamer consisting of six molecule 1 monomers and a nonsymmetric hexamer consisting of five molecule 1 monomers and one molecule 2 monomer. For each of the six monomers, the restraints used in the calculations of the hexamer structures include, in addition to the intramonomeric restraints, 25 NOEs between insulin and phenol, 23 NOEs and two hydrogen bonds across the dimer interface, nine NOEs across the trimer interface, and five intramonomeric or two intermonomeric NOEs, respectively, specifying the different orientations of the Phe(B25) ring. The coordination of the two Zn atoms was defined by eight distance restraints. Thus, a total of 4394 and 4391 distance restraints, respectively, were used in the two hexamer calculations. The NOE restraints were classified in an iterative process as intra- or intermonomeric on the basis of their consistency or inconsistency with the structure of the monomer. The assignment of the dimer- and trimer-specific NOEs was made using the crystal structure of the R6 hexamer as the starting model. For both solution hexamers, the average backbone rms deviation is 0.81 A, if the less well-defined N- and C-terminal residues are excluded. The corresponding rms deviations for all heavy atoms are 1.17 and 1.19 A for the nonsymmetric and symmetric hexamer, respectively. The overall solution structure of the R6 insulin hexamer is compact, rigid, and symmetric and resembles the corresponding crystal structure. However, the extension of the B-chain alpha-helix, which characterizes the R state, is shorter in the solution structure than in the crystal structure. Also, the study shows that the orientation of the Phe(B25) ring has no effect on the structure of the rest of the molecule, within the uncertainty of the structure determination. The importance of these findings for the current model for the insulin-receptor interaction is discussed.

Solution structure and backbone dynamics of the human alpha3-chain type VI collagen C-terminal Kunitz domain,.

The solution structure and backbone dynamics of the 58-residue C-terminal Kunitz domain fragment [alpha3(VI)] of human alpha3-chain type VI collagen has been studied by two-dimensional 1H-1H and 1H-15N nuclear magnetic resonance spectroscopy at 303 K. The solution structure is represented by an ensemble of 20 structures calculated with X-PLOR using 612 distance and 47 dihedral angle restraints. The distance restraints were obtained by a complete relaxation matrix analysis using MARDIGRAS. The root mean squared (rms) deviation is 0.91 A for the backbone atoms of the residues Thr2(8)-Gly12(18), Arg15(21)-Tyr35(41), and Gly40(46)-Pro57(63). The central beta-sheet [residues Ile18(24)-Tyr35(41)] and the C-terminal alpha-helix [residues Gln48(54)-Cys55(61)] are better defined with a backbone rms deviation of 0.46 A. The solution structure of alpha3(VI) is virtually identical to the crystal structure of alpha3(VI) and to the solution structure of bovine pancreatic trypsin inhibitor (BPTI). The 15N spin-lattice and spin-spin relaxation rates and the 1H-15N heteronuclear nuclear Overhauser enhancement (NOE) were analyzed using both the "model-free" formalism [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4559 and 4559-4570] and the reduced spectral density mapping procedure [Farrow, N. A., Szabo, A., Torchia, D. A., & Kay, L. E. (1995) J. Biomol.NMR 6, 153-162]. The results obtained from the "model-free" analysis include an overall correlation time tauc of 3. 00 ns and backbone order parameters S2 in the range from 0.28 to 0. 93. The necessity of including an exchange term in the analysis of the relaxation data from 14 residues indicated that these residues are involved in motions on the micro- to millisecond time scale. The majority of the 14 residues are located in the vicinity of the Cys14(20)-Cys38(44) disulfide bond, suggesting the presence of a disulfide bond isomerization similar to the one observed in BPTI [Otting, G., Liepinsh, E., & Wüthrich, K. (1993) Biochemistry 32, 3571-3582]. It is suggested that this disulfide bond isomerization is the main reason for the surprisingly small effect on trypsin inhibition observed when Thr13(19) of alpha3(VI) is substituted with Pro.

Elucidation of the origin of multiple conformations of the human alpha 3-chain type VI collagen C-terminal Kunitz domain: the reorientation of the Trp21 ring.

The human alpha 3-chain type VI collagen C-terminal Kunitz domain fragment (alpha 3(VI)) has been studied by two dimensional 1H-1H and 1H-13C NMR spectroscopy at 303 K. It is shown that the secondary structure of the protein is strikingly similar to that of BPTI, and a number of unusual H alpha chemical shifts, which are highly conserved in Kunitz-domain proteins, are also observed for a alpha 3(VI). Furthermore, a series of exchange cross peaks observed in 1H-1H spectra shows that a large number of protons in the central beta-sheet exist in two different chemical environments, corresponding to two unequally populated conformations that are slowly exchanging on the NMR time scale. Several protons, including Ser47(53) H alpha, Arg32(28) H(gamma 1) and H(gamma 2), and GLN48(54) H(beta 2), all located in the vicinity of the Trp21(27) ring in the crystal structure of alpha 3(VI) [Arnoux, B. et al. (1995) J. Mol. Biol., 246, 609-617], have very different chemical shifts in the two conformations, the most affected being Gln48(54) H(beta 2) (delta sigma = 3 ppm), which is placed directly above the Trp21(27) ring in the crystal structure of alpha 3(VI). It should be concluded that the origin of the multiple conformations of the central beta-sheet is a reorientation of the Trp21(27) ring. From the intensities of corresponding signals in the two conformations, the populations, the population of the minor conformation was found to be 6.4 +/- 0.2% of that of the major conformation, while a rate constant kM = 1.01 +/- 0.05 s-1 for the major to minor interconversion was obtained from a series of NOESY spectra with different mixing times. In addition, it is shown that Cys14(20)-Cys38(44) disulfide bond isomerization, previously observed in BPTI [Otting, G. et al. (1993) Biochemistry, 32, 3571-3582], is also likely to occur in alpha 3(VI).

Solution structure of reduced plastocyanin from the blue-green alga Anabaena variabilis.

The three-dimensional solution structure of plastocyanin from Anabaena variabilis (A.v.PCu) has been determined by nuclear magnetic resonance spectroscopy. Sixty structures were calculated by distance geometry from 1141 distance restraints and 46 dihedral angle restraints. The distance geometry structures were optimized by simulated annealing and restrained energy minimization. The average rms deviation from the mean structure for the 20 structures with the lowest total energy is 1.25 A for the backbone atoms and 1.75 A for all heavy atoms. Overall, the global tertiary fold of A.v.PCu resembles those of other plastocyanins which have been structurally characterized by X-ray diffraction and NMR methods. This holds even though A.v.PCu is longer than any other known plastocyanins, contains far less invariant amino acid residues, and has an overall charge that differs considerably from those of other plastocyanins (+1 vs -9 +/- 1 at pH > or = 7). The most striking feature of the A.v. PCu structure is the absence of the beta-turn, formed at the remote site by residues (58)-(61) in most higher plant plastocyanins. The displacement caused by the absence of this turn is compensated for by an extension of the small helix [from Ala53(51) to Ser60(58) in A.v.PCu] found in other plastocyanins. Moreover, the extra residues of A.v.PCu from Pro77 to Asp79 form an appended loop. These two features allow A.v.PCu to retain almost the same global fold as observed in other plastocyanins. From a comparison with the structures of other plastocyanins it is concluded that the lack of negatively charged residues at the remote site, rather than the specific structure of A.v.PCu, is the main reason for the failure of the remote site of this plastocyanin to function as a significant electron transfer site.

Solution structure of the superactive monomeric des-[Phe(B25)] human insulin mutant: elucidation of the structural basis for the monomerization of des-[Phe(B25)] insulin and the dimerization of native insulin.

The three-dimensional solution structure of des-[Phe(B25)] human insulin has been determined by nuclear magnetic resonance spectroscopy and restrained molecular dynamics calculations. Thirty-five structures were calculated by distance geometry from 581 nuclear Overhauser enhancement-derived distance constraints, ten phi torsional angle restraints, the restraints from 16 helical hydrogen bonds, and three disulfide bridges. The distance geometry structures were optimized using simulated annealing and restrained energy minimization. The average root-mean-square (r.m.s.) deviation for the best 20 refined structures is 1.07 angstroms for the backbone and 1.92 angstroms for all atoms if the less well-defined N and C-terminal residues are excluded. The helical regions are more well defined, with r.m.s. deviations of 0.64 angstroms for the backbone and 1.51 angstroms for all atoms. It is found that the des-[Phe(B25)] insulin is a monomer under the applied conditions (4.6 to 4.7 mM, pH 3.0, 310 K), that the overall secondary and tertiary structures of the monomers in the 2Zn crystal hexamer of native insulin are preserved, and that the conformation-averaged NMR solution structure is close to the structure of molecule 1 in the hexamer. The structure reveals that the lost ability of des-[Phe(B25)] insulin to self-associate is caused by a conformational change of the C-terminal region of the B-chain, which results in an intra-molecular hydrophobic interaction between Pro(B28) and the hydrophobic region Leu(B11)-Leu(B15) of the B-chain alpha-helix. This interaction interferes with the inter-molecular hydrophobic interactions responsible for the dimerization of native insulin, depriving the mutant of the ability to dimerize. Further, the structure displays a series of features that may explain the high potency of the mutant on the basis of the current model for the insulin-receptor interaction. These features are: a change in conformation of the C-terminal region of the B-chain, the absence of strong hydrogen bonds between this region and the rest of the molecule, and a relatively easy accessibility to the Val(A3) residue.

Assignment of the backbone carbonyl resonances in (15)N-labelled proteins with (13)C at natural abundance by a 2D triple-resonance correlation technique.

A 2D NMR experiment for assignment of backbone carbon resonances in small and medium-sized (15)N-labelled proteins with (13)C at natural abundance is presented. The experiment is a two-dimensional variant of the HNCO triple-resonance experiment and is demonstrated by application to a 6 kDa protein at relatively low concentration (2 mM) and temperature (30°C). The experiment is particularly suitable for assignment of carbonyl resonances.

Structural details of Asp(B9) human insulin at low pH from two-dimensional NMR titration studies.

A titration study of the dimeric Asp(B9) mutant of human insulin was performed using two-dimensional NMR spectroscopy. Based on 10 NOESY spectra recorded in the pH range 1.73-3.93, the pKa values of the seven carboxyl groups in the mutant were determined, and the titration shifts of 46 pH-dependent protons in non-ionizable groups were investigated. Further, the pKa values of the two histidine imidazole rings were determined from a series of 1D spectra recorded in the pH range 6.65-10.0. The titration shifts of all pH-dependent protons were analyzed by a nonlinear least-squares fitting procedure, using an equation that describes a one-step titration. Also the pH dependence of the exchange rate of the amide proton of Phe(B24) was determined in the applied pH range. On the basis of the experimental results, it is concluded that the Asp(B9) residue forms an N-cap of the B-chain alpha-helix through an interaction between the side-chain carboxyl group of the residue and the dipole of the helix. Further, the titration data show that salt bridges are established between Glu(B13) and His(B10) and between Asn(A21) and Arg(B22) at pH values, where the interacting groups are ionized, and that a hydrogen bond exists between the amide proton of Val(A3) and the C-terminal carboxyl group of Thr(B30). Most surprisingly, the data analysis shows that the Asp(B9) insulin exists as a dimer throughout the investigated pH range, that is, also at pH values where there is a substantial negative charge repulsion in the monomer-monomer interface of the dimer.

A new 2D NMR method for measurement of JHH coupling constants.

A new 2D NMR pulse sequence for E.COSY-type measurement of J(HH) coupling constants is introduced. It exploits a heteronuclear spin, e.g., 13C, for displacement in the omega(1) frequency dimension via a large heteronuclear J coupling. The experiment is demonstrated by application to a heptapeptide at the natural abundance 13C level. It is suitable, for example, for measurement of 3J(HH) and 4J(HH) coupling constants in peptides and proteins.

Specific nitrogen-15 labelling of leucine residues in human growth hormone.

Biosynthetic human growth hormone (hGH) specifically 15N labelled in the leucine residues has been obtained by recombinant DNA technology, using 15N-labelled leucine and an E. coli strain that requires leucine. It is shown that, despite the possibility of minor transaminase activity, the labelling on the whole is specific, and that the two-dimensional 1H-15N correlation NMR spectra of hGH can be greatly simplified by this methodology.

Three-dimensional solution structure of an insulin dimer. A study of the B9(Asp) mutant of human insulin using nuclear magnetic resonance, distance geometry and restrained molecular dynamics.

The solution structure of the B9(Asp) mutant of human insulin has been determined by two-dimensional 1H nuclear magnetic resonance spectroscopy. Thirty structures were calculated by distance geometry from 451 interproton distance restraints based on intra-residue, sequential and long-range nuclear Overhauser enhancement data, 17 restraints on phi torsional angles obtained from 3JH alpha HN coupling constants, and the restraints from 17 hydrogen bonds, and the three disulphide bridges. The distance geometry structures were optimized using restrained molecular dynamics (RMD) and energy minimization. The average root-mean-square deviation for the best 20 RMD refined structures is 2.26 A for the backbone and 3.14 A for all atoms if the less well-defined N and C-terminal residues are excluded. The helical regions are better defined, with root-mean-square deviation values of 1.11 A for the backbone and 2.03 A for all atoms. The data analysis and the calculations show that B9(Asp) insulin, in water solution at the applied pH (1.8 to 1.9), is a well-defined dimer with no detectable difference between the two monomers. The association of the two monomers in the solution dimer is relatively loose as compared with the crystal dimer. The overall secondary and tertiary structures of the monomers in the 2Zn crystal hexamer is found to be preserved. The conformation-averaged NMR structures obtained for the monomer is close to the structure of molecule 1 in the hexamer of the 2Zn insulin crystal. However, minor, but significant deviations from this structure, as well as from the structure of monomeric insulin in solution, exist and are ascribed to the absence of the hexamer and crystal packing forces, and to the presence of monomer-monomer interactions, respectively. Thus, the monomer in the solution dimer shows a conformation similar to that of the crystal monomer in molecular regions close to the monomer-monomer interface, whereas it assumes a conformation similar to that of the solution structure of monomeric insulin in other regions, suggesting that B9(Asp) insulin adopts a monomer-like conformation when this is not inconsistent with the monomer-monomer arrangement in the dimer.

Characterization of tertiary interactions in a folded protein by NMR methods: studies of pH-induced structural changes in human growth hormone.

The pH-induced conformational changes in human growth hormone (hGH) have been studied, using a new quantitative NMR approach that combines 13C labeling of specific backbone carbonyl carbons with a complete spectral analysis of the corresponding 13C resonances. Thus, a complete analysis of the carbonyl resonances of the 26 Leu residues of hGH and their variation with pH provided detailed information about the equilibrium folding processes of the protein, including information about the kinetics of the folding. By combining this information with the pH dependence of readily identifiable 1H resonances, the pH-induced changes observed in the carbonyl carbon spectra can be associated with specific regions in the protein and can be ascribed to a series of localized adjustments in the tertiary structure, brought about by changes in the hydrogen bond interactions or electrostatic interactions between different residues in the globular folded protein. The preexchange lifetimes of these adjustments range from a fraction of a millisecond to a few milliseconds.