HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Abstract Full Text (PDF) All Versions of this Article: ps.062398606v1 15/12/2697    most recent Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Butterwick, J. A. Articles by Palmer, A. G. Search for Related Content PubMed PubMed Citation Articles by Butterwick, J. A. Articles by Palmer, A. G., III Social Bookmarking                   What's this?

An inserted Gly residue fine tunes dynamics between mesophilic and thermophilic ribonucleases H

Department of Biochemistry and Molecular Biophysics, Columbia University, New York, New York 10032, USA

(RECEIVED June 13, 2006; FINAL REVISION September 12, 2006; ACCEPTED September 18, 2006)

	Abstract

TOP Abstract Introduction Results Discussion Materials and methods Acknowledgments References

 
Dynamic processes are inherent properties of proteins and are crucial for a wide range of biological functions. To address how changes in protein sequence and structure affect dynamic processes, a quantitative comparison of microsecond-to-microsecond time scale conformational changes, measured by solution NMR spectroscopy, within homologous mesophilic and thermophilic ribonuclease H (RNase H) enzymes is presented. Kinetic transitions between the observed major state (high population) and alternate (low population) conformational state(s) of the substrate-binding handle region in RNase H from the mesophile Escherichia coli (ecRNH) and thermophile Thermus thermophilus (ttRNH) occur with similar kinetic exchange rate constants, but the difference in stability between exchanging conformers is smaller in ttRNH compared to ecRNH. The altered thermodynamic equilibrium between kinetically exchanging conformers in the thermophile is recapitulated in ecRNH by the insertion of a Gly residue within a putative hinge between -helices B and C. This Gly insertion is conserved among thermophilic RNases H, and allows the formation of additional intrahelical hydrogen bonds. A Gly residue inserted between -helices B and C appears to relieve unfavorable interactions in the transition state and alternate conformer(s) and represents an important adaptation to adjust conformational changes within RNase H for activity at high temperatures.

Keywords: nuclear magnetic resonance; ribonuclease H; protein dynamics; thermal stability

	Introduction

TOP Abstract Introduction Results Discussion Materials and methods Acknowledgments References

 
Protein function relies on a combination of time-average static structure and time-dependent conformational fluctuations. Intrinsic dynamic properties influence many biological processes, including folding, ligand binding, catalysis, and allosterism (Wand 2001; Akke 2002; Benkovic and Hammes-Schiffer 2003; Kern and Zuiderweg 2003; Palmer 2004; Tousignant and Pelletier 2004; Kay 2005). Dynamic processes have been shown to be (partially) rate-limiting for catalysis, and differences in the time scale of motion are directly correlated with differences in the catalytic efficiency of homologous enzymes (Williams and McDermott 1995; Kohen and Klinman 2000; Codreanu et al. 2002; Cole and Loria 2002; Wolf-Watz et al. 2004). Several recent studies have highlighted the importance of individual amino acid residues in modulating dynamic processes within proteins (Eisenmesser et al. 2005; Kovrigin and Loria 2006). However, detailed mechanisms by which local interactions within the protein structure control the kinetic and thermodynamic properties of dynamic processes remain largely unknown.

Ribonucleases H (RNases H, EC 3.1.26.4 [EC] ) from the mesophile Escherichia coli (ecRNH) and thermophile Thermus thermophilus (ttRNH) share significant amino acid sequence and structural similarity, but display differences in thermal stability and catalytic activity. Typical of other mesophilic and thermophilic proteins (Jaenicke and Böhm 1998; Kumar et al. 2001; Vieille and Zeikus 2001), ttRNH is more resistant to thermal (and chemical) denaturation (Kanaya and Itaya 1992; Hollien and Marqusee 1999a,b) and is significantly less catalytically active than ecRNH at mesophilic temperatures (Kanaya and Itaya 1992). One source for these differences may stem from altered dynamic properties between these enzymes: Several regions display localized differences in dynamics on the picosecond-to-nanosecond (ps–ns) and microsecond-to-millisecond (µs–ms) time scales (Butterwick et al. 2004).

Chemical exchange line broadening is observed for nuclear spins in a highly basic loop between -helices C (_C) and D, termed the handle (see Fig. 1), in both ecRNH and ttRNH (Mandel et al. 1995, 1996; Kroenke et al. 1998; Yamasaki et al. 1998; Butterwick et al. 2004), suggesting that µs–ms time scale conformational changes in this region are an integral property of both enzymes. Residues in _C and the handle are important for RNA:DNA hybrid recognition (Kanaya et al. 1991), although they are not essential for catalytic activity (Keck and Marqusee 1995). Nuclear magnetic resonance (NMR) spectroscopy and circular dichroism studies suggest that the handle changes conformation with respect to the catalytic domain upon complex formation (Oda et al. 1993; Iwai et al. 1996). The structural details of these transitions remain elusive; however, based on chemical shift changes upon RNA:DNA hybrid titration and models of the substrate-bound complex, residues between _B and _C have been suggested as a possible hinge in ecRNH (Oda et al. 1993). This region contains a single amino acid insertion in ttRNH relative to ecRNH: Gly85 is in a left-handed helical conformation (backbone dihedral angles = 70° and = –9°), and allows an additional hydrogen bond between the Gly85 carboxyl and Gly89 amide, forming the N-terminal cap of _C (Ishikawa et al. 1993d).

View larger version (54K): [in this window] [in a new window]   Figure 1. Sequence and structural alignments of RNases H. (A) ClustalW (Thompson et al. 1994) multiple sequence alignment of RNase H from the mesophilic bacterium E. coli (37°C) and several thermophilic bacteria: Chlorobium tepidum (47°C), Moorella thermoacetica (58°C), T. thermophilus (70°C), and Thermoanaerobacter tengcongensis (75°C). The optimal growth temperature for each organism is indicated in parentheses. Conserved active site and substrate-binding residues are highlighted in blue; the inserted Gly is shown in red. -Helical (A–E) and -strand (1–5) secondary structural elements and the handle in WT ecRNH are indicated above the alignment. (B) Structural superposition of WT ecRNH (blue, PDB 2RN2), WT ttRNH (orange, PDB 1RIL), and iG80b ecRNH (gray, PDB 1GOA) using backbone atoms along core helices A and D. Structural elements from (A) and the N and C termini are indicated. Ribbon diagrams for this and subsequent figures were drawn using MOLMOL (Koradi et al. 1996).

1

1

	Results

TOP Abstract Introduction Results Discussion Materials and methods Acknowledgments References

 
Gly insertion is conserved among thermophilic RNases H
Sequence alignments for mesophilic and thermophilic RNases H highlight several amino acid substitutions that may confer thermal stability or alter catalytic activity. More than 100 bacterial RNase H sequences with a conserved handle are found in the Swiss-Prot database, including four from thermophilic organisms (Fig. 1A). In addition to several amino acid substitutions known to promote increased thermal stability in ecRNH, for example H62P and V74L (Kimura et al. 1992b; Ishikawa et al. 1993b), the Gly insertion between _B and _C is conserved in thermophilic RNases H. Furthermore, this insertion is largely absent from mesophilic sequences: Only RNase H enzymes from mesophiles Lactobacillus plantarum and Deinococcus radiodurans contain this insertion.

Structural and functional consequences of Gly insertion
Although the Gly insertion is conserved in thermophilic RNases H, insertion of Gly80b into ecRNH provides only a moderate increase in thermal stability of 0.1 kJ/mol (Kimura et al. 1992b). However, this insertion leads to a drastic decrease in catalytic efficiency, increasing the binding affinity by 50-fold and decreasing the catalytic rate constant (k _cat) by 10-fold at 30°C (Ishikawa et al. 1993c), such that the activity of iG80b ecRNH is similar to the activity of ttRNH (Kanaya and Itaya 1992). Structurally, Gly80b is nicely accommodated between _B and _C, and adopts a similar left-handed helical conformation ( = 52°, = 46°) as Gly85 in ttRNH (Fig. 1B). In addition to the expected N-terminal cap hydrogen bond between Gly80b and Asn84, another hydrogen bond is present between the Ile78 carboxyl and Gly80b amide forming the C-terminal cap of _B (Ishikawa et al. 1993c). The additional Gly also shifts the Gly77 carboxyl–Ile82 amide hydrogen bond to Gly77–Trp81. Given the minimal structural perturbations and marginal stabilization induced by the Gly80b insertion, the drastic decrease in activity may be due to modification of dynamic processes involving the substrate-binding handle.

Backbone amide ¹H and ¹⁵N resonances for the iG80b ecRNH and dG85 ttRNH mutant enzymes were assigned by comparison to WT spectra complemented with an HMQC-NOESY-HSQC spectrum recorded at 300 K. Chemical shift differences between WT and mutant spectra are small outside the site of mutation, suggesting that insertion or removal of the Gly residue does not significantly alter the protein structure. In addition, NOESY cross-peaks are observed between backbone amide protons separated by one or two residues throughout both _B and _C, confirming the helical nature of these regions (data not shown). The largest ¹⁵N chemical shift differences between WT and iG80b ecRNH are for residues Thr79 and Ile82: _iG80b–_WT = 4.94 ppm and –3.77 ppm, respectively. These shift changes are the correct sign and magnitude as expected for gain of a hydrogen bond involving the Ile78 carboxyl and switch of the hydrogen bond involving the Ile82 amide (Xu and Case 2002).

Interestingly, the equivalent residues also exhibit similar ¹⁵N chemical shift differences between WT and dG85 ttRNH, but with the opposite sign: _dG85–_WT = –4.78 ppm for Thr83 and 3.89 ppm for Leu87. Chemical shift changes for residues in _B, _C and the handle between WT and mutant ecRNH and ttRNH are negatively correlated overall (Fig. 2A). Not only are the changes in chemical shift similar in magnitude for these proteins, but the iG80b ecRNH chemical shifts resemble WT ttRNH, and dG85 ttRNH shifts resemble WT ecRNH (Fig. 2B). Thus, the hydrogen bonding pattern and backbone geometry is similar between WT ttRNH and iG80b ecRNH and between WT ecRNH and dG85 ttRNH. These results are in slight disagreement with the crystal structure of WT ttRNH (Fig. 2C) (Ishikawa et al. 1993d): Based on chemical shift changes, hydrogen bonds would be expected between carboxyls of Phe82 and Gly77 and amides of Gly85 and Leu87. These interactions appear to be present in solution and slight packing rearrangements may prevent these hydrogen bonds from being recognized in the crystalline state.

View larger version (36K): [in this window] [in a new window]   Figure 2. Structural consequences of the Gly insertion. (A) Chemical shift changes between WT and iG80b ecRNH (ec) and WT and dG85 ttRNH (tt) are compared as = mutant–WT for backbone amide (left) 1H and (right) 15N nuclei. Residues in B, C and the handle (70–100 in ecRNH) are indicated by open circles; the correlation coefficient for these residues are –0.63 for 1H and –0.84 for 15N. (B) Overlay of 1H–15N correlation spectra for iG80b ecRNH and dG85 ttRNH with (left) WT ttRNH and (right) WT ecRNH highlighting the positions of Thr79 (ecRNH) and Thr83 (ttRNH). (C) Main-chain hydrogen bonds in the B/C hinge and handle are shown in green for WT ecRNH (blue, PDB 2RN2), iG80b ecRNH (gray, PDB 1GOA), and WT ttRNH (orange, PDB 1RIL). Additional hydrogen bonds seen in the crystal structures of iG80b ecRNH and WT ttRNH are shown in magenta.

View larger version (28K): [in this window] [in a new window]   Figure 3. Chemical exchange line broadening in WT and mutant RNases H. Backbone amide 15N R ex at 300 K, 14.1 T for WT ecRNH, WT ttRNH, iG80b ecRNH, and dG85 ttRNH are shown. Residues with significant chemical exchange line broadening (R ex 2.5 sec–1) are indicated by open circles. Vertical dashed lines indicate the sites of mutation. Secondary structural elements are diagrammed at the top of the figure. Data for WT ttRNH are from Butterwick et al. (2004).

Temperature dependence of chemical exchange
The temperature dependence of R _ex provides information on the time scale and energetics of conformational transitions (Mandel et al. 1996; Evenäs et al. 1999; Butterwick et al. 2004). R _ex measurements were repeated at temperatures between 285 K and 310 K for WT ecRNH and mutant RNases H; data for WT ttRNH have been reported previously (Butterwick et al. 2004). Values of R _ex increase with decreasing temperature in both WT enzymes, indicating that chemical exchange is intermediate to fast on the chemical shift time scale (Fig. 4A; Mandel et al. 1996; Butterwick et al. 2004). In addition, although R _ex is not a kinetic rate constant per se, the temperature dependence of R _ex displays Arrhenius behavior over the temperature ranges studied, where the thermophilic enzyme is characterized by a larger "apparent" activation energy barrier (Ê _a) (Fig. 4B; Butterwick et al. 2004). At T 290 K, Ile82 and Asn100 have significant R _ex in WT ecRNH, but accurate Ê _a could not be determined for these residues (see Materials and Methods).

View larger version (25K): [in this window] [in a new window]   Figure 4. Temperature dependence of chemical exchange line broadening. (A) R ex measured at multiple temperatures and 14.1 T and (B) Arrhenius plots are shown for (left) Trp90 in WT and iG80b ecRNH and (right) Trp95 in WT and dG85 ttRNH. Open and closed symbols represent data for WT and mutant RNases H, respectively. Apparent activation energies, Ê a, derived from the slope of the lines in B are listed in Table 1. Square symbols represent data that deviate from a linear Arrhenius relationship, and are not included for determining Ê a. Data for WT ttRNH are from Butterwick et al. (2004).

Relaxation dispersion of backbone amide nuclei
To further characterize the exchange kinetics within WT and mutant RNases H, backbone amide ¹⁵N R ₂( _e) relaxation dispersion was measured at a temperature of 300 K for WT ecRNH, WT ttRNH, and iG80b ecRNH using R ₁ spin relaxation techniques. The magnitude of exchange broadening in the dG85 ttRNH mutant is too small for accurate relaxation dispersion measurements (see Materials and Methods). All dispersion data are well described by a two-state exchange model

characterized by a kinetic exchange rate constant k _ex = k _AB + k _BA, where k _AB and k _BA are the forward and reverse rate constants between exchanging conformers A and B. On- and off-resonance R ₁ experiments (Akke and Palmer 1996; Mulder et al. 1998; Massi et al. 2004) were combined to span effective fields of _e 10²–10⁴ sec^–1 at static magnetic fields of 11.7 and 18.8 T. For all exchanging residues in WT ecRNH and ttRNH, R _ex( _e) is barely suppressed at the highest _e obtained, such that k _ex >> 10⁴ sec^–1 for residues in both enzymes (Fig. 5; Table 2). In this regime, the individual forward and reverse rate constants and the thermodynamic and structural characteristics of the exchange process cannot be separated and only k _ex and _ex = p _A p _B ² can be determined, where p _A and p _B are the fractional populations of the conformers A and B, respectively (p _A p _B), and is the chemical shift difference between exchanging conformers (Palmer 2004).

View larger version (24K): [in this window] [in a new window]   Figure 5. R 1 relaxation dispersion for backbone amide 15N nuclei. R 2( e) dispersion results are shown for residues Trp90 and Lys60 in WT and iG80b ecRNH, and Trp95 and Val126 in WT ttRNH at 300 K. Closed and open symbols represent data collected at 11.7 and 18.8 T, respectively.

In a similar manner to the WT enzymes, chemical exchange for all residues in iG80b ecRNH are fast on the chemical shift time scale (Fig. 5; Table 2). Residues in _C and the handle region were generally individually fit with slower kinetics than WT ecRNH with k _ex 15,000–40,000 sec^–1; subsequent global analysis of Trp81, Ile82, Arg88, Trp90, Thr92, Lys95, Lys99, and Asn100 yielded k _ex = (27 ± 2) x 10³ sec^–1. The reduced k _ex does not completely account for the increased R _ex relative to WT ecRNH, and residues in iG80b ecRNH are fit with approximately fourfold larger _ex. Relaxation dispersion for Lys60, however, yields k _ex = (101 ± 30) x 10³ sec^–1, which is identical within experimental uncertainties to the value obtained for WT ecRNH. Thus, reduced R _ex for Lys60 in iG80b ecRNH compared to WT is due to an 0.8-fold reduction in _ex.

Interpreting chemical exchange effects
In the fast exchange regime, only k _ex and _ex can be determined from relaxation dispersion measurements; consequently, the underlying energetics and chemical shift changes of the conformational transition cannot be determined. However, ratios of k _ex and _ex provide information on the relative changes in the free energy difference between exchanging conformers (G) and between transition states (G ) for two proteins. Assuming that k _AB and k _BA are described by transition state theory, k _AB = k _oexp[–G /RT] and k _BA = k _oexp[–(G –G)/RT], and that site populations follow a Boltzmann distribution, p _B/p _A = exp[–G/RT], then

where the subscripts 1 and 2 represent the two proteins of interest, k _o is the preexponential factor for k _AB and k _BA, R is the universal gas constant, and G and G are the energy differences of the minor conformer and transition state relative to the major conformer, respectively, and G = G ₁ – G ₂ and G = G ₁ – G ₂ . The second equalities in Equations 2 and 3 arise if p _A – p _B 1 for either protein.

Given the small chemical shift changes between WT and iG80b ecRNH outside the site of mutation, the exchanging conformer structures probably are similar for both proteins, such that _{WT iG80b}. In this case, using R ₁ results for Trp90 (_N,iG80b – _N,WT = 0.41 ppm), G = G _iG80b – G _WT = –3.6 ± 0.1 kJ/mol and G = G _iG80b – G _WT = –2.5 ± 0.5 kJ/mol at 300 K. Thus, insertion of Gly80b reduces the stability difference between exchanging conformers, increasing p _B relative to WT ecRNH (Fig. 6). The decreased stability difference between conformers is slightly larger than the decrease in activation energy, leading to slower kinetics in the mutant enzyme. Virtually identical results are obtained with data from Asn100.

View larger version (12K): [in this window] [in a new window]   Figure 6. Summary of conformational changes in RNases H. Energy diagram comparing the relative stabilities and activation energies for exchanging conformers in WT ecRNH (blue), WT ttRNH (orange), and iG80b ecRNH (gray) at 300 K. Diagram is drawn to scale assuming p B = 1% for WT ecRNH and k o = k B T/h, where k B and h are Boltzmann's and Planck's constants, respectively.

1

	Discussion

TOP Abstract Introduction Results Discussion Materials and methods Acknowledgments References

 
Conformational changes are ubiquitous in biology. While structural studies have identified numerous cases in which structural changes occur during a biological process, the roles of specific residues and local interactions in controlling thermodynamics and kinetics of structural transitions remain largely uncharacterized. In the present work, NMR spectroscopy was used to characterize µs–ms time scale motional processes in homologous mesophilic and thermophilic RNase H enzymes. Comparison of the conformational dynamic properties of WT and Gly insertion/deletion mutants was used to probe the role of local interactions surrounding a putative hinge in modulating dynamics within the substrate-binding handle.

A Gly insertion between _B and _C is conserved in thermophilic RNases H that utilize a handle region for RNA:DNA hybrid recognition, suggesting that this residue is important for biological function at high temperatures. This insertion is a focal point for additional intrahelical hydrogen bonds that marginally stabilize the protein. Chemical shifts surrounding this region are similar between WT ttRNH and iG80b ecRNH and between WT ecRNH and dG85 ttRNH, suggesting that interactions within the iG80b ecRNH mutant are more "thermophilic-like" while those in dG85 ttRNH are more "mesophilic-like." The similarity between WT ttRNH and iG80b ecRNH is further reflected in their reduced catalytic efficiencies relative to WT ecRNH (Kanaya and Itaya 1992; Ishikawa et al. 1993c). dG85 ttRNH, however, does not show significantly increased catalytic activity at low temperatures (to resemble WT ecRNH), possibly due to the uncoupling of conformational changes between the handle and ₅ near the active site (J.A. Butterwick and A.G. Palmer III, unpubl.). Thus, while the crystallographic structures of WT ttRNH and iG80b ecRNH display slight differences, both proteins likely incorporate the same additional hydrogen bonds in solution, compared to WT ecRNH and dG85 ttRNH.

Currently, the structural details of alternative conformations sampled in solution and leading to chemical exchange line broadening are unknown. Differences in R _ex between WT and mutant RNases H are spread throughout residues in _C and the handle, suggesting that conformational changes within these regions are highly coupled. Large R _ex are not observed for most residues, indicating that are small and the local structure in _C and the handle is maintained. Because chemical exchange line broadening is not observed for residues Thr79, Gly80b and Asn84 in iG80b ecRNH, the additional hydrogen bonds most likely are maintained during conformational changes about this region. Residues Trp90 and Asn100 may serve as hinges for motions of the handle: Trp90 is the first and Asn100 is the last residue in the handle, both Trp90 and Asn100 show significant increase in R _ex in the iG80b ecRNH mutant, and both Trp90 and Asn100 have left-handed helical backbone dihedral angles in the observed major conformation.

The hypothesis that chemical exchange reflects movement of the entire handle as a whole is supported by the dynamic properties of a thermostabilized ecRNH quintuple G23A/H62P/V74L/K95G/D134H mutant (Yamasaki et al. 1998). Mutation of the left-handed Lys95 to Gly is well accommodated within the handle of ecRNH (Ishikawa et al. 1993a) and significantly stabilizes the protein (Kimura et al. 1992a). While no global changes in dynamics were observed in the quintuple mutant compared to the WT protein, several residues surrounding the sites of mutation experienced a change in dynamics across multiple time scales. Interestingly, an increase in R _ex was observed for residues Ile82, Trp90, Lys91, and Asn100 in the quintuple mutant (Yamasaki et al. 1998), although to a lesser extent than seen in iG80b ecRNH. Lys95 is 5–14 Å away from these residues, yet removal of its bulky side chain likely alters the thermodynamics or kinetics of structural transitions across the handle.

The remarkable similarity between the free energy profiles of iG80b ecRNH and WT ttRNH (Fig. 6) suggests that the Gly insertion plays a major role in modulating conformational dynamics of the handle region of RNase H. In iG80b ecRNH, both the transition and minor states are stabilized by a similar extent relative to WT ecRNH. Thus, the Gly insertion does not drastically alter the kinetic rate constants, but rather adjusts the thermodynamic equilibrium between conformers. Lacking a C atom, Gly occupies approximately threefold more Ramachandran space compared to other amino acids. This structural plasticity appears to alleviate unfavorable interactions in the transition state and minor conformer(s) in RNase H.

For RNases H, the kinetics of the chemical exchange process are likely too fast to be rate limiting for enzymatic activity. Previous characterization of the WT proteins using a nonanucleotide hybrid substrate at 30°C yielded k _cat 1.5 and 0.3 sec^–1 for WT ecRNH and WT ttRNH, respectively (Kanaya et al. 1990; Kanaya and Itaya 1992). Given the current results, the lower limit for WT ecRNH is p _B 1%, corresponding to k _AB 400 sec^–1 and k _BA 41,600 sec^–1, which are much faster than k _cat. Thus, the altered thermodynamic distribution rather than changes in the kinetic rate constants afforded by the Gly insertion likely represents an important adaptation to fine tune conformational changes within RNase H for activity at high temperatures. The Gly80b insertion into ecRNH does not significantly stabilize the protein; thus, the mutant ecRNH enzyme is unfolded and inactive at the optimal temperature for activity in the wild-type thermophilic ttRNH enzyme (343 K) (Kanaya and Itaya 1992). Additional amino acid substitutions, such as H62P and V74L (Kimura et al. 1992b; Ishikawa et al. 1993b), may be required to confer increased thermal stability in order for the effects of the Gly insertion to be seen at higher temperatures.

Insertion of Gly80b into ecRNH has few effects on the ground-state structure outside of the _B/_C hinge region and insignificantly affects overall thermodynamic stability, but renders the structure of the _B/_C hinge region and enzymatic activity more similar to the hinge structure and enzymatic activity observed for WT ttRNH (Ishikawa et al. 1993c). The crystallographically observed structures of apo WT ecRNH, WT ttRNH, and iG80b ecRNH are consistent with the dominant conformations observed in solution by NMR spectroscopy, but these structures appear incapable of productive interactions between the handle and substrate (Iwai et al. 1996). The present results demonstrate that the potential energy surface controlling µs–ms time scale transitions to alternate higher energy conformations also is more similar in iG80b ecRNH to WT ttRNH than WT ecRNH. These observations imply that enzymatic activity of RNases H depends upon the structures and energetics of the alternative conformations that are sampled at equilibrium in solution.

	Materials and methods

TOP Abstract Introduction Results Discussion Materials and methods Acknowledgments References

 
Sample preparation
Plasmids containing the WT RNase H coding sequences have been described previously (Mandel et al. 1995; Hollien and Marqusee 1999b). For improved solubility, ttRNH contains four Cys mutations (Hollien and Marqusee 1999b): C17A, C45S, C67A, and C154S. Gly insertion (GGT codon) and deletion mutations were created using the Quik-Change protocol (Stratagene). WT and mutant RNases H were expressed in E. coli BL21(DE3) (Stratagene) using M9 minimal media (Sambrook et al. 1989) prepared with 99% (v/v) deuterium oxide and [98% ¹⁵N]-ammonium chloride (Cambridge Isotope Laboratories, Inc.) and purified as previously described (Kroenke et al. 1998; Hollien and Marqusee 1999b). ²H incorporation at nonexchangeable proton sites was estimated to be >85% for all samples by mass spectrometry. Final protein concentrations were 1.0 mM for WT and iG80b ecRNH, and 0.3 mM for WT and dG85 ttRNH in 100 mM d ₃-sodium acetate at pH 5.5, 10% (v/v) deuterium oxide, 0.02% (w/v) sodium azide, and 3 mM sodium 2,2-dimethyl-2-silapentane-sulfonate (DSS); ecRNH samples also contained 1 mM d ₁₀-dithiothreitol.

NMR spectroscopy
All NMR data were collected at Columbia University (11.7 and 14.1 T) or the New York Structural Biology Center (16.4 and 18.8 T). Data at static magnetic fields of 11.7 and 18.8 T were acquired on Bruker DRX500 and AV800, respectively, equipped with triple-resonance probes and three-axis gradients. Data at static magnetic fields of 14.1 and 16.4 T were acquired on Bruker DRX600 and AV700, respectively, equipped with triple-resonance cryoprobes with Z-axis gradients. For all experiments, the temperature was calibrated using a 100% (v/v) methanol standard (Cavanagh et al. 1996). Data processing and analysis were performed with NMRPipe (Delaglio et al. 1995) and Sparky (T.D. Goddard and D.G. Kneller, Sparky 3, University of California, San Francisco) software, along with in-house written programs.

Resonance assignments for mutant RNases H
Backbone ¹H and ¹⁵N resonance assignments for iG80b ecRNH and dG85 ttRNH were determined by comparison with WT spectra supplemented by HMQC-NOESY-HSQC spectra (Cavanagh et al. 1996) acquired at a static magnetic field of 16.4 T and a temperature of 300 K. The spectra were recorded using with (t ₁ x t ₂ x t ₃) 48 x 32 x 768 complex points and 2.857 x 2.857 x 10 kHz spectral widths. Assignments, referenced to internal DSS, have been deposited to the BioMagResBank (http://www.bmrb.wisc.edu) under accession codes 7277 (dG85 ttRNH) and 7278 (iG80b ecRNH).

Analysis of chemical exchange line broadening
For residues subject to chemical exchange line broadening, R ₂ = R ₂° + R _ex, where R ₂° is the intrinsic transverse relaxation rate constant which is sensitive to motion on the ps–ns time scale and the protein size. A TROSY-based Hahn-echo sequence (Wang et al. 2003) was used to measure R ₂° and R _ex at a static magnetic field of 14.1 T and temperatures between 285 K and 310 K. This method uses the transverse ¹⁵N chemical shift anisotropy/¹H–¹⁵N dipolar interference rate constant ( _xy) to determine R ₂° = _xy, where is the trimmed mean R ₂/ _xy ratio for nonexchanging residues. Values of were between 1.31 and 1.47 at all temperatures for all proteins. Analysis of the data with a weighted mean = 1.40 does not significantly change the results. Relaxation data were collected in an interleaved manner with (t ₁ x t ₂) 150 x 1024 complex points and 2.5 x 10 kHz spectral widths. A minimum of four experiments were acquired at each temperature for each protein. The average and standard deviation were taken as the relaxation rate constant and uncertainty, respectively.

Apparent activation energies, Ê _a, were determined from the temperature dependence of chemical exchange assuming an Arrhenius relationship for R _ex:

where H is the transition state enthalpy for the forward reaction and H is the difference in enthalpy between the two exchanging states (Mandel et al. 1996; Evenäs et al. 1999; Butterwick et al. 2004). The second equality assumes that exchange is fast (k _ex >> ). For each residue, Ê _a was determined from a linear fit of ln(R _ex) versus 1/T using R _ex data acquired at a minimum of three temperatures. At low temperatures, deviations from linearity may become apparent if the exchange kinetics slow to the intermediate regime (k _ex ); these data are not included in the determination of Ê _a.

R₁ relaxation dispersion analysis
Chemical exchange line broadening is suppressed in the presence of an effective field, _e. In an R ₁ experiment, magnetization is spin-locked parallel to _e by the application of a continuous radiofrequency (rf) field with frequency _rf and amplitude ₁. The magnitude of the effective field is _e ² = ₁ ² + ², where = _o – _rf is the observed resonance offset from the applied field. The relaxation rate constant, R ₁, for magnetization aligned parallel to _e is given by

where R ₁ is the longitudinal relaxation rate constant, = arctan( ₁/) is the tilt angle of the effective field, and

is the effective transverse relaxation constant where R _ex( _e) is the residual R _ex.

R ₁ and on-resonance R ₁ relaxation rate constants were measured at a temperature of 300 K and static magnetic fields of 11.7 and 18.8 T using relaxation time delays between 0 and 1.2 sec and 0 and 100 msec, respectively (Farrow et al. 1994; Massi et al. 2004). Seven to nine unique time points were sampled with three duplicates for estimation of the error in signal intensity. The relaxation rate constants were determined from the best single exponential fit to the data with the error estimated using a jackknife algorithm. Off-resonance R ₁ relaxation rate constants were measured from pairs of interleaved spectra with relaxation delays of 2 msec and 60–200 msec assuming a single exponential decay (Akke and Palmer 1996; Mulder et al. 1998); errors were estimated based on repeat experiments. For on-resonance experiments, magnetization was aligned using a free precession period of 1/ ₁ (Yamazaki et al. 1994; Massi et al. 2004), while off-resonance experiments utilized 10 msec tanh/tan adiabatic pulses each with a frequency sweep of 15 kHz (Mulder et al. 1998). For all R ₁ experiments, ₁ was calibrated using the residual scalar coupling in the presence and absence of the applied field (Palmer et al. 2001).

Analysis of relaxation dispersion data followed a similar protocol previously described (Grey et al. 2006). Initially, R ₁ and R ₁ relaxation data at both static magnetic fields along with R ₂° and free precession R _ex estimated from TROSY Hahn-echo experiments acquired at 11.7 and 18.8 T were analyzed on a per-residue basis using Equations 5 and 6. F-statistical testing (at = 0.05 confidence level) was used to select the appropriate functional form for R _ex( _e): either chemical exchange is not present, with R _ex( _e) = 0; or exchange is fast on the chemical shift time scale (k _ex >> _e), with

Additional R _ex( _e) models that include exchange outside the fast time scale regime (Trott and Palmer 2002), or a three-site exchange mechanism (Grey et al. 2003) did not significantly improve the fitted parameters. Optimization of R ₁ and R ₂° rate constants and _ex and k _ex exchange parameters were accomplished with routines written in Mathematica (v4.1, Wolfram Research, Inc.) utilizing the Levenberg-Marquart algorithm. Off-resonance R ₁ data were included only if 20° 65° while on-resonance R ₁ data required |0.4 ₁/2|. Residues characterized by similar k _ex values by inspection were subsequently analyzed globally with a single k _ex, maintaining individual _ex.

	Footnotes

 
¹ Present address: Laboratory of Molecular Neurobiology and Biophysics, The Rockefeller University, 1230 York Avenue, New York, NY 10021, USA.

Reprint requests to: Arthur G. Palmer III, Department of Biochemistry and Molecular Biophysics, Columbia University, 650 West 168 Street, Black Building 507, New York, NY 10032, USA; e-mail: agp6{at}columbia.edu; fax: (212) 305-6949.

Article published online ahead of print. Article and publication date are at http://www.proteinscience.org/cgi/doi/10.1110/ps.062398606.

	Acknowledgments

TOP Abstract Introduction Results Discussion Materials and methods Acknowledgments References

 
We thank members of the Palmer laboratory for helpful discussions, and M. Grey and F. Massi additionally for assistance with R ₁ measurements and critical reviews of the manuscript. This work was supported by National Institutes of Health (NIH) Grant GM50291 (A.G.P.). A.G.P. is a member of the New York Structural Biology Center supported by NIH Grant GM66354.

	References

TOP Abstract Introduction Results Discussion Materials and methods Acknowledgments References

 
Akke, M. 2002. NMR methods for characterizing microsecond to millisecond dynamics in recognition and catalysis. Curr. Opin. Struct. Biol. 12: 642–647.[CrossRef][Medline]

Akke, M. and Palmer III, A.G. 1996. Monitoring macromolecular motions on microsecond to millisecond time scales by R ₁ – R ₁ constant relaxation time NMR spectroscopy. J. Am. Chem. Soc. 118: 911–912.[CrossRef]

Benkovic, S.J. and Hammes-Schiffer, S. 2003. A perspective on enzyme catalysis. Science 301: 1196–1202.[Abstract/Free Full Text]

Butterwick, J.A., Loria, J.P., Astrof, N.S., Kroenke, C.D., Cole, R., Rance, M., and Palmer III, A.G. 2004. Multiple time scale backbone dynamics of homologous thermophilic and mesophilic ribonuclease HI enzymes. J. Mol. Biol. 339: 855–871.[CrossRef][Medline]

Cavanagh, J., Fairbrother, W.J., Palmer III, A.G., and Skelton, N.J. 1996. Protein NMR spectroscopy: Principles and practice. Academic Press, San Diego, CA.

Codreanu, S.G., Ladner, J.E., Xiao, G., Stourman, N.V., Hachey, D.L., Gilliland, G.L., and Armstrong, R.N. 2002. Local protein dynamics and catalysis: Detection of segmental motion associated with rate-limiting product release by a glutathione transferase. Biochemistry 41: 15161–15172.[CrossRef][Medline]

Cole, R. and Loria, J.P. 2002. Evidence for flexibility in the function of ribonuclease A. Biochemistry 41: 6072–6081.[CrossRef][Medline]

Delaglio, F., Grzesiek, S., Vuister, G.W., Zhu, G., Pfeifer, J., and Bax, A. 1995. NMRPipe: A multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 6: 277–293.[Medline]

Eisenmesser, E.Z., Millet, O., Labeikovsky, W., Korzhnev, D.M., Wolf-Watz, M., Bosco, D.A., Skalicky, J.J., Kay, L.E., and Kern, D. 2005. Intrinsic dynamics of an enzyme underlies catalysis. Nature 438: 117–121.[CrossRef][Medline]

Evenäs, J., Forsén, S., Malmendal, A., and Akke, M. 1999. Backbone dynamics and energetics of a calmodulin domain mutant exchanging between closed and open conformations. J. Mol. Biol. 289: 603–607.[CrossRef][Medline]

Farrow, N.A., Muhandiram, R., Singer, A.U., Pascal, S.M., Kay, C.M., Gish, G., Shoelson, S.E., Pawson, T., Forman-Kay, J.D., and Kay, L.E. 1994. Backbone dynamics of a free and phosphopeptide-complexed Src homology 2 domain studied by ¹⁵N NMR relaxation. Biochemistry 33: 5984–6003.[CrossRef][Medline]

Grey, M.J., Wang, C.Y., and Palmer III, A.G. 2003. Disulfide bond isomerization in basic pancreatic trypsin inhibitor: Multisite chemical exchange quantified by CPMG relaxation dispersion and chemical shift modeling. J. Am. Chem. Soc. 125: 14324–14335.[CrossRef][Medline]

Grey, M.J., Tang, Y., Alexov, E., McKnight, C.J., Raleigh, D.P., and Palmer III, A.G. 2006. Characterizing a partially folded intermediate of the villin headpiece domain under non-denaturing conditions: Contribution of His41 to the pH-dependent stability of the N-terminal subdomain. J. Mol. Biol. 355: 1078–1094.[CrossRef][Medline]

Hollien, J. and Marqusee, S. 1999a. Structural distribution of stability in a thermophilic enzyme. Proc. Natl. Acad. Sci. 96: 13674–13678.[Abstract/Free Full Text]

Hollien, J. and Marqusee, S. 1999b. A thermodynamic comparison of mesophilic and thermophilic ribonucleases H. Biochemistry 38: 3831–3836.[CrossRef][Medline]

Ishikawa, K., Kimura, S., Kanaya, S., Morikawa, K., and Nakamura, H. 1993a. Structural study of mutants of Escherichia coli ribonuclease HI with enhanced thermostability. Protein Eng. 6: 85–91.[Abstract/Free Full Text]

Ishikawa, K., Nakamura, H., Morikawa, K., and Kanaya, S. 1993b. Stabilization of Escherichia coli ribonuclease HI by cavity-filling mutations within a hydrophobic core. Biochemistry 32: 6171–6178.[CrossRef][Medline]

Ishikawa, K., Nakamura, H., Morikawa, K., Kimura, S., and Kanaya, S. 1993c. Cooperative stabilization of Escherichia coli ribonuclease HI by insertion of Gly-80b and Gly-77 Ala substitution. Biochemistry 32: 7136–7142.[CrossRef][Medline]

Ishikawa, K., Okumura, M., Katayanagi, K., Kimura, S., Kanaya, S., Nakamura, H., and Morikawa, K. 1993d. Crystal structure of ribonuclease H from Thermus thermophilus HB8 refined at 2.8 Å resolution. J. Mol. Biol. 230: 529–542.[CrossRef][Medline]

Iwai, S., Wakasa, M., Ohtsuka, E., Kanaya, S., Kidera, A., and Nakamura, H. 1996. Interaction of the basic protrusion of Escherichia coli ribonuclease HI with its substrate. J. Mol. Biol. 263: 699–706.[CrossRef][Medline]

Jaenicke, R. and Böhm, G. 1998. The stability of proteins in extreme environments. Curr. Opin. Struct. Biol. 8: 738–748.[CrossRef][Medline]

Kanaya, S. and Itaya, M. 1992. Expression, purification, and characterization of a recombinant ribonuclease H from Thermus thermophilus HB8. J. Biol. Chem. 267: 10184–10192.[Abstract/Free Full Text]

Kanaya, S., Kohara, A., Miura, Y., Sekiguchi, A., Iwai, S., Inoue, H., Ohtsuka, E., and Ikehara, M. 1990. Identification of the amino acid residues involved in an active site of Escherichia coli ribonuclease H by site-directed mutagenesis. J. Biol. Chem. 265: 4615–4621.[Abstract/Free Full Text]

Kanaya, S., Katsuda-Nakai, C., and Ikehara, M. 1991. Importance of the positive charge cluster in Escherichia coli ribonuclease HI for the effective binding of the substrate. J. Biol. Chem. 266: 11621–11627.[Abstract/Free Full Text]

Kay, L.E. 2005. NMR studies of protein structure and dynamics. J. Magn. Reson. 173: 193–207.[CrossRef][Medline]

Keck, J.L. and Marqusee, S. 1995. Substitution of a highly basic helix/loop sequence into the RNase H domain of human immunodeficiency virus reverse transcriptase restores its Mn²⁺-dependent RNase H activity. Proc. Natl. Acad. Sci. 92: 2740–2744.[Abstract/Free Full Text]

Kern, D. and Zuiderweg, E.R. 2003. The role of dynamics in allosteric regulation. Curr. Opin. Struct. Biol. 13: 748–757.[CrossRef][Medline]

Kimura, S., Kanaya, S., and Nakamura, H. 1992a. Thermostabilization of Escherichia coli ribonuclease HI by replacing left-handed helical Lys95 with Gly or Asn. J. Biol. Chem. 267: 22014–22017.[Abstract/Free Full Text]

Kimura, S., Nakamura, H., Hashimoto, T., Oobatake, M., and Kanaya, S. 1992b. Stabilization of Escherichia coli ribonuclease HI by strategic replacement of amino acid residues with those from the thermophilic counterpart. J. Biol. Chem. 267: 21535–21542.[Abstract/Free Full Text]

Kohen, A. and Klinman, J.P. 2000. Protein flexibility correlates with degree of hydrogen tunneling in thermophilic and mesophilic alcohol dehydrogenases. J. Am. Chem. Soc. 122: 10738–10739.[CrossRef]

Koradi, R., Billeter, M., and Wüthrich, K. 1996. MOLMOL: A program for display and analysis of macromolecular structures. J. Mol. Graph. 14: 51–55: 29–32.

Kovrigin, E.L. and Loria, J.P. 2006. Enzyme dynamics along the reaction coordinate: Critical role of a conserved residue. Biochemistry 45: 2636–2647.[CrossRef][Medline]

Kroenke, C.D., Loria, J.P., Lee, L.K., Rance, M., and Palmer III, A.G. 1998. Longitudinal and transverse ¹H–¹⁵N dipolar/¹⁵N chemical shift anisotropy relaxation interference: Unambiguous determination of rotational diffusion tensors and chemical exchange effects in biological macromolecules. J. Am. Chem. Soc. 120: 7905–7915.[CrossRef]

Kumar, S., Tsai, C.J., and Nussinov, R. 2001. Thermodynamic differences among homologous thermophilic and mesophilic proteins. Biochemistry 40: 14152–14165.[CrossRef][Medline]

Mandel, A.M., Akke, M., and Palmer III, A.G. 1995. Backbone dynamics of Escherichia coli ribonuclease HI: Correlations with structure and function in an active enzyme. J. Mol. Biol. 246: 144–163.[CrossRef][Medline]

Mandel, A.M., Akke, M., and Palmer III, A.G. 1996. Dynamics of ribonuclease H: Temperature dependence of motions on multiple time scales. Biochemistry 35: 16009–16023.[CrossRef][Medline]

Massi, F., Johnson, E., Wang, C., Rance, M., and Palmer III, A.G. 2004. NMR R ₁ rotating-frame relaxation with weak radio frequency fields.