Multiple time scale backbone dynamics of homologous thermophilic and mesophilic ribonuclease HI enzymes

Backbone conformational fluctuations on multiple time scales in a cysteine-free Thermus thermophilus ribonuclease HI mutant (ttRNH(*)) are quantified using (15)N nuclear magnetic spin relaxation. Laboratory-frame relaxation data acquired at 310 K and at static magnetic field strengths of 11.7, 14.1 and 18.8 T are analysed using reduced spectral density mapping and model-free approaches. Chemical exchange line broadening is characterized using Hahn-echo transverse and multiple quantum relaxation data acquired over a temperature range of 290-320 K and at a static magnetic field strength of 14.1 T. Results for ttRNH(*) are compared to previously published data for a mesophilic homologue, Escherichia coli ribonuclease HI (ecRNH). Intramolecular conformational fluctuations on the picosecond-to-nanosecond time scale generally are similar for ttRNH(*) and ecRNH. beta-Strands 3 and 5 and the glycine-rich region are more rigid while the substrate-binding handle region and C-terminal tail are more flexible in ttRNH(*) than in ecRNH. Rigidity in the two beta-strands and the glycine-rich region, located along the periphery of the central beta-sheet, may be associated with the increased thermodynamic stability of the thermophilic enzyme. Chemical exchange line broadening, reflecting microsecond-to-millisecond time scale conformational changes, is more pronounced in ttRNH(*) than in ecRNH, particularly for residues in the handle and surrounding the catalytic site. The temperature dependence of chemical exchange show an increase of approximately 15 kJ/mol in the apparent activation energies for ttRNH(*) residues in the handle compared to ecRNH. Increased activation barriers, coupled with motion between alpha-helices B and C not present in ecRNH, may be associated with the reduced catalytic activity of the thermophilic enzyme at 310 K.     

Autonomously Folding Protein Fragments Reveal Differences in the Energy Landscapes of Homologous RNases H

An important approach to understanding how a protein sequence encodes its energy landscape is to compare proteins with different sequences that fold to the same general native structure. In this work, we compare E. coli and T. thermophilus homologs of the protein RNase H. Using protein fragments, we create equilibrium mimics of two different potential partially-folded intermediates (I_core and I_core+1) hypothesized to be present on the energy landscapes of these two proteins. We observe that both T. thermophilus RNase H (ttRNH) fragments are folded and have distinct stabilities, indicating that both regions are capable of autonomous folding and that both intermediates are present as local minima on the ttRNH energy landscape. In contrast, the two E. coli RNase H (ecRNH) fragments have very similar stabilities, suggesting that the presence of additional residues in the I_core+1 fragment does not affect the folding or structure as compared to I_core. NMR experiments provide additional evidence that only the I_core intermediate is populated by ecRNH. This is one of the biggest differences that has been observed between the energy landscapes of these two proteins. Additionally, we used a FRET experiment in the background of full-length ttRNH to specifically monitor the formation of the I_core+1 intermediate. We determine that the ttRNH I_core+1 intermediate is likely the intermediate populated prior to the rate-limiting barrier to global folding, in contrast to E. coli RNase H for which I_core is the folding intermediate. This result provides new insight into the nature of the rate-limiting barrier for the folding of RNase H.     

Thermodynamic System Drift in Protein Evolution

Proteins from thermophiles are generally more thermostable than their mesophilic homologs, but little is known about the evolutionary process driving these differences. Here we attempt to understand how the diverse thermostabilities of bacterial ribonuclease H1 (RNH) proteins evolved. RNH proteins from Thermus thermophilus (ttRNH) and Escherichia coli (ecRNH) share similar structures but differ in melting temperature (T_m) by 20°C. ttRNH''s greater stability is caused in part by the presence of residual structure in the unfolded state, which results in a low heat capacity of unfolding (ΔC_p) relative to ecRNH. We first characterized RNH proteins from a variety of extant bacteria and found that T_m correlates with the species'' growth temperatures, consistent with environmental selection for stability. We then used ancestral sequence reconstruction to statistically infer evolutionary intermediates along lineages leading to ecRNH and ttRNH from their common ancestor, which existed approximately 3 billion years ago. Finally, we synthesized and experimentally characterized these intermediates. The shared ancestor has a melting temperature between those of ttRNH and ecRNH; the T_ms of intermediate ancestors along the ttRNH lineage increased gradually over time, while the ecRNH lineage exhibited an abrupt drop in T_m followed by relatively little change. To determine whether the underlying mechanisms for thermostability correlate with the changes in T_m, we measured the thermodynamic basis for stabilization—ΔC_p and other thermodynamic parameters—for each of the ancestors. We observed that, while the T_m changes smoothly, the mechanistic basis for stability fluctuates over evolutionary time. Thus, even while overall stability appears to be strongly driven by selection, the proteins explored a wide variety of mechanisms of stabilization, a phenomenon we call “thermodynamic system drift.” This suggests that even on lineages with strong selection to increase stability, proteins have wide latitude to explore sequence space, generating biophysical diversity and potentially opening new evolutionary pathways.     

Hydrogen-exchange stabilities of RNase T1 and variants with buried and solvent-exposed Ala --> Gly mutations in the helix

Hydrogen-exchange rates were measured for RNase T1 and three variants with Ala --> Gly substitutions at a solvent-exposed (residue 21) and a buried (residue 23) position in the helix: A21G, G23A, and A21G + G23A. These results were used to measure the stabilities of the proteins. The hydrogen-exchange stabilities (DeltaG(HX)) for the most stable residues in each variant agree with the equilibrium conformational stability measured by urea denaturation (DeltaG(U)), if the effects of D(2)O and proline isomerization are included [Huyghues-Despointes, B. M. P., Scholtz, J. M., and Pace, C. N. (1999) Nat. Struct. Biol. 6, 210-212]. These residues also show similar changes in DeltaG(HX) upon Ala --> Gly mutations (DeltaDeltaG(HX)) as compared to equilibrium measurements (DeltaDeltaG(U)), indicating that the most stable residues are exchanging from the globally unfolded ensemble. Alanine is stabilizing compared to glycine by 1 kcal/mol at a solvent-exposed site 21 as seen by other methods for the RNase T1 protein and peptide helix [Myers, J. K., Pace, C. N., and Scholtz, J. M. (1997) Proc. Natl. Acad. Sci. U.S.A. 94, 3833-2837], while it is destabilizing at the buried site 23 by the same amount. For the A21G variant, only local NMR chemical shift perturbations are observed compared to RNase T1. For the G23A variant, large chemical shift changes are seen throughout the sequence, although X-ray crystal structures of the variant and RNase T1 are nearly superimposable. Ala --> Gly mutations in the helix of RNase T1 at both helical positions alter the native-state hydrogen-exchange stabilities of residues throughout the sequence.     

Factors influencing the stability of alpha-helices and beta-strands in thermophilic ribonuclease H.

Understanding the influence of structural parameters is crucial to enhance the thermal stability of proteins. In this work, the stability (deltaG) of residues in different secondary structures of Ribonuclease H (RNase H) has been analyzed with 48 amino acid properties. The properties reflecting hydrophobicity show a good correlation with stability. Further, the linear distribution of surrounding hydrophobicity in alpha-helices, obtained from the three dimensional structure of thermophilic RNase H, agrees well with experimental deltaG values. Moreover, the stability parameters correlate better in alpha-helices than those did in beta-strand segments. Multiple regression analysis, incorporating combinations of three properties from among all possible combinations of the 48 properties, increased the correlation coefficient to 0.77.     

Flexibility of NS5 Methyltransferase-Polymerase Linker Region Is Essential for Dengue Virus Replication

We examined the function of the conserved Val/Ile residue within the dengue virus NS5 interdomain linker (residues 263 to 272) by site-directed mutagenesis. Gly substitution or Gly/Pro insertion after the conserved residue increased the linker flexibility and created slightly attenuated viruses. In contrast, Pro substitution abolished virus replication by imposing rigidity in the linker and restricting NS5''s conformational plasticity. Our biochemical and reverse genetics experiments demonstrate that NS5 utilizes conformational regulation to achieve optimum viral replication.     

Effect of the disease-causing mutations identified in human ribonuclease (RNase) H2 on the activities and stabilities of yeast RNase H2 and archaeal RNase HII

Eukaryotic ribonuclease (RNase) H2 consists of one catalytic and two accessory subunits. Several single mutations in any one of these subunits of human RNase H2 cause Aicardi-Goutières syndrome. To examine whether these mutations affect the complex stability and activity of RNase H2, three mutant proteins of His-tagged Saccharomyces cerevisiae RNase H2 (Sc-RNase H2*) were constructed. Sc-G42S*, Sc-L52R*, and Sc-K46W* contain single mutations in Sc-Rnh2Ap*, Sc-Rnh2Bp*, and Sc-Rnh2Cp*, respectively. The genes encoding the three subunits were coexpressed in Escherichia coli, and Sc-RNase H2* and its derivatives were purified in a heterotrimeric form. All of these mutant proteins exhibited enzymatic activity. However, only the enzymatic activity of Sc-G42S* was greatly reduced compared to that of the wild-type protein. Gly42 is conserved as Gly10 in Thermococcus kodakareansis RNase HII. To analyze the role of this residue, four mutant proteins, Tk-G10S, Tk-G10A, Tk-G10L, and Tk-G10P, were constructed. All mutant proteins were less stable than the wild-type protein by 2.9-7.6 degrees C in T(m). A comparison of their enzymatic activities, substrate binding affinities, and CD spectra suggests that the introduction of a bulky side chain into this position induces a local conformational change, which is unfavorable for both activity and substrate binding. These results indicate that Gly10 is required to make the protein fully active and stable.     

Thermal Adaptation of Conformational Dynamics in Ribonuclease H

The relationship between inherent internal conformational processes and enzymatic activity or thermodynamic stability of proteins has proven difficult to characterize. The study of homologous proteins with differing thermostabilities offers an especially useful approach for understanding the functional aspects of conformational dynamics. In particular, ribonuclease HI (RNase H), an 18 kD globular protein that hydrolyzes the RNA strand of RNA:DNA hybrid substrates, has been extensively studied by NMR spectroscopy to characterize the differences in dynamics between homologs from the mesophilic organism E. coli and the thermophilic organism T. thermophilus. Herein, molecular dynamics simulations are reported for five homologous RNase H proteins of varying thermostabilities and enzymatic activities from organisms of markedly different preferred growth temperatures. For the E. coli and T. thermophilus proteins, strong agreement is obtained between simulated and experimental values for NMR order parameters and for dynamically averaged chemical shifts, suggesting that these simulations can be a productive platform for predicting the effects of individual amino acid residues on dynamic behavior. Analyses of the simulations reveal that a single residue differentiates between two different and otherwise conserved dynamic processes in a region of the protein known to form part of the substrate-binding interface. Additional key residues within these two categories are identified through the temperature-dependence of these conformational processes.     

Ribonuclease A suggests how proteins self-chaperone against amyloid fiber formation

Genomic analyses have identified segments with high fiber-forming propensity in many proteins not known to form amyloid. Proteins are often protected from entering the amyloid state by molecular chaperones that permit them to fold in isolation from identical molecules; but, how do proteins self-chaperone their folding in the absence of chaperones? Here, we explore this question with the stable protein ribonuclease A (RNase A). We previously identified fiber-forming segments of amyloid-related proteins and demonstrated that insertion of these segments into the C-terminal hinge loop of nonfiber-forming RNase A can convert RNase A into the amyloid state through three-dimensional domain-swapping, where the inserted fiber-forming segments interact to create a steric zipper spine. In this study, we convert RNase A into amyloid-like fibers by increasing the loop length and hence conformational freedom of an endogenous fiber-forming segment, SSTSAASS, in the N-terminal hinge loop. This is accomplished by sandwiching SSTSAASS between inserted Gly residues. With these inserts, SSTSAASS is now able to form the steric zipper spine, allowing RNase A to form amyloid-like fibers. We show that these fibers contain RNase A molecules retaining their enzymatic activity and therefore native-like structure. Thus, RNase A appears to prevent fiber formation by limiting the conformational freedom of this fiber-forming segment from entering a steric zipper. Our observations suggest that proteins have evolved to self-chaperone by using similar protective mechanisms.     

Structural equilibrium fluctuations in mesophilic and thermophilic alpha-amylase

By comparing a mesophilic alpha-amylase with its thermophilic homolog, we investigated the relationship between thermal stability and internal equilibrium fluctuations. Fourier transform infrared spectroscopy monitoring hydrogen/deuterium (H/D) exchange kinetics and incoherent neutron scattering measuring picosecond dynamics were used to study dynamic features of the folded state at room temperature. Fairly similar rates of slowly exchanging amide protons indicate about the same free energy of stabilization DeltaG(stab) for both enzymes at room temperature. With respect to motions on shorter time scales, the thermophilic enzyme is characterized by an unexpected higher structural flexibility as compared to the mesophilic counterpart. In particular, the picosecond dynamics revealed a higher degree of conformational freedom for the thermophilic alpha-amylase. The mechanism proposed for increasing thermal stability in the present case is characterized by entropic stabilization and by flattening of the curvature of DeltaG(stab) as a function of temperature.     

Comparison of the folding processes of T. thermophilus and E. coli ribonucleases H

In order to examine how the stabilization of thermophilic proteins affects their folding, we have characterized the folding process of Thermus thermophilus ribonuclease H using circular dichroism, fluorescence, and pulse-labeling hydrogen exchange. Like its homolog from Escherichia coli, this thermophilic protein populates a partially folded kinetic intermediate within the first few milliseconds of folding. The structure of this intermediate is similar to that of E.coli RNase H and corresponds remarkably well to a partially folded form that is populated at low levels in the native state of the protein. Proline isomerization appears to partly limit the folding of the thermophilic but not the mesophilic protein. Lastly, unlike other thermophilic proteins, which unfold much more slowly than their mesophilic counterparts, T.thermophilus RNase H folds and unfolds with overall rates similar to those of E.coli RNase H.     

High pressure NMR reveals that apomyoglobin is an equilibrium mixture from the native to the unfolded

Pressure-induced reversible conformational changes of sperm whale apomyoglobin have been studied between 30 bar and 3000 bar on individual residue basis by utilizing 1H/15N hetero nuclear single-quantum coherence two-dimensional NMR spectroscopy at pH 6.0 and 35 degrees C. Apomyoglobin showed a series of pressure-dependent NMR spectra as a function of pressure, assignable to the native (N), intermediates (I), molten globule (MG) and unfolded (U) conformers. At 30 bar, the native fold (N) shows disorder only in the F helix. Between 500 bar and 1200 bar, a series of locally disordered conformers I are produced, in which local disorder occurs in the C helix, the CD loop, the G helix and part of the H helix. At 2000 bar, most cross-peaks exhibit severe line-broadening, suggesting the formation of a molten globule, but at 3000 bar all the cross-peaks reappear, showing that the molten globule turns into a well-hydrated, mobile unfolded conformation U. Since all the spectral changes were reversible with pressure, apomyoglobin is considered to exist as an equilibrium mixture of the N, I, MG and U conformers at all pressures. MG is situated at 2.4+/-(0.1) kcal/mol above N at 1 bar and the unfolding transition from the combined N-I state to MG is accompanied by a loss of partial molar volume by 75+/-(3) ml/mol. On the basis of these observations, we postulate a theorem that the partial molar volume of a protein decreases in parallel with the loss of its conformational order.     

Characterization of internal motions of Escherichia coli ribonuclease H by Monte Carlo simulation

The backbone dynamics of Escherichia coli ribonuclease H (RNase H) is studied by a recently developed off-lattice Monte Carlo/Metropolis simulation technique. A low-resolution model (virtual-bond model) is used together with knowledge-based potentials. The calculated mean-square fluctuations in alpha carbons are in good agreement with crystallographic temperature factors. The conformations generated around the native state are analyzed by time-dependent orientational and conformational correlation functions to study the internal motions of RNase at different time windows. A correlation between the free-energy changes for native-state hydrogen exchange (HX) and the extent of the autocorrelation in the rotations of the virtual bonds at long times has been observed. Cross-correlations between the rotations of the bonds, which are near-neighbor in the sequence, are effective in all time windows and help the secondary structures to preserve their kinetic stability. On the other hand, the existence of cross-correlations at long times help the tertiary contacts be maintained. The order parameter of NH bond vector for each residue has been calculated and compared with 15N-NMR relaxation measurements.     

Native-state energetics of a thermostabilized variant of ribonuclease HI.

Escherichia coli RNase HI is a well-characterized model system for protein folding and stability. Controlling protein stability is critical for both natural proteins and for the development of engineered proteins that function under extreme conditions. We have used native-state hydrogen exchange on a variant containing the stabilizing mutation Asp10 to alanine in order to determine its residue-specific stabilities. On average, the DeltaG(unf) value for each residue was increased by 2-3 kcal/mol, resulting in a lower relative population of partially unfolded forms. Though increased in stability by a uniform factor, D10A shows a distribution of stabilities in its secondary structural units that is similar to that of E. coli RNase H, but not the closely related protein from Thermus thermophilus. Hence, the simple mutation used to stabilize the enzyme does not recreate the balance of conformational flexibility evolved in the thermophilic protein.     

Solution structure of the cytotoxic RNase 4 from oocytes of bullfrog Rana catesbeiana

Cytotoxic ribonucleases with antitumor activity are mainly found in the oocytes and early embryos of frogs. Native RC-RNase 4 (RNase 4), consisting of 106 residues linked with four disulfide bridges, is a cytotoxic ribonuclease isolated from oocytes of bullfrog Rana catesbeiana. RNase 4 belongs to the bovine pancreatic ribonuclease (RNase A) superfamily. Recombinant RC-RNase 4 (rRNase 4), which contains an additional Met residue and glutamine instead of pyroglutamate at the N terminus, was found to possess less catalytic and cytotoxic activities than RNase 4. Equilibrium thermal and guanidine-HCl denaturation CD measurements revealed that RNase 4 is more thermally and chemically stable than rRNase 4. However, CD and NMR data showed that there is no gross conformational change between native and recombinant RNase 4. The NMR solution structure of rRNase 4 was determined to comprise three alpha-helices and two sets of antiparallel beta-sheets. Superimposition of each structure with the mean structure yielded an average root mean square deviation (RMSD) of 0.72(+/-0.14)A for the backbone atoms, and 1.42(+/-0.19)A for the heavy atoms in residues 3-105. A comparison of the 3D structure of rRNase 4 with the structurally and functionally related cytotoxic ribonuclease, onconase (ONC), showed that the two H-bonds in the N-terminal pyroglutamate of ONC were not present at the corresponding glutamine residue of rRNase 4. We suggest that the loss of these two H-bonds is one of the key factors responsible for the reductions of the conformational stability, catalytic and cytotoxic activities in rRNase 4. Furthermore, the differences of side-chain conformations of subsite residues among RNase A, ONC and rRNase 4 are related to their distinct catalytic activities and base preferences.     

Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK(a)s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (epsilon(prot)) of four and by explicit conformational changes. The pK(a) of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and >90% have errors of <1 pH units whereas only 3% have errors >2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (epsilon(prot) 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK(a).     

Contributions of folding cores to the thermostabilities of two ribonucleases H

To investigate the contribution of the folding cores to the thermodynamic stability of RNases H, we used rational design to create two chimeras composed of parts of a thermophilic and a mesophilic RNase H. Each chimera combines the folding core from one parent protein and the remaining parts of the other. Both chimeras form active, well-folded RNases H. Stability curves, based on CD-monitored chemical denaturations, show that the chimera with the thermophilic core is more stable, has a higher midpoint of thermal denaturation, and a lower change in heat capacity (DeltaCp) upon unfolding than the chimera with the mesophilic core. A possible explanation for the low DeltaCp of both the parent thermophilic RNase H and the chimera with the thermophilic core is the residual structure in the denatured state. On the basis of the studied parameters, the chimera with the thermophilic core resembles a true thermophilic protein. Our results suggest that the folding core plays an essential role in conferring thermodynamic parameters to RNases H.     

Role of the switch II region in the conformational transition of activation of Ha-ras-p21

The role of the switch II region in the conformational transition of activation of Ha-ras-p21 has been investigated by mutating residues predicted to act as hinges for the conformational transition of this loop (Ala59, Gly60, and Gly75) (Díaz JF, Wroblowski B, Schlitter J, Engelborghs Y, 1997, Proteins 28:434-451), as well as mutating the catalytic residue Gln61. The proposed mutations of the hinge residues decrease the rate of the conformational transition of activation as measured by the binding of BeF3- to the GDP-p21 complex. Also, the thermodynamic parameters of the binding reaction are altered by a factor between three and five, depending on the temperature. (Due to changes in activation and reaction enthalpies, partially compensated by entropy changes.) The control mutation Q61H in which only the catalytic residue is changed has only a limited effect on the kinetic rate constants of the conformational transition and on the thermodynamic parameters of the reaction. The fact that mutations of the hinge residues of the switch II region affect both the binding of the phosphate analog and the conformational transition of activation indicates that the switch II is implicated both in the early and the late states of the transition.