Structural, thermodynamic, and mutational analyses of a psychrotrophic RNase HI

Quantifying the Relationship between Conformational Dynamics and Enzymatic Activity in Ribonuclease HI Homologues

Ribonuclease HI (RNHI), a ubiquitous, non-sequence-specific endonuclease, cleaves the RNA strand in RNA/DNA hybrids. RNHI functions in replication and genome maintenance, and retroviral reverse transcriptases contain an essential ribonuclease H domain. Nuclear magnetic resonance (NMR) spectroscopy combined with molecular dynamics (MD) simulations suggests a model in which the extended handle region domain of Escherichia coli RNHI populates (substrate-binding-competent) "open" and (substrate-binding-incompetent) "closed" states, while the thermophilic Thermus thermophilus RNHI mainly populates the closed state at 300 K [Stafford, K. A., Robustelli, P., and Palmer, A. G., III (2013) PLoS Comput. Biol. 9, 1-10]. In addition, an in silico-designed mutant E. coli Val98Ala RNHI was predicted to populate primarily the closed state. The work presented here validates this model and confirms the predicted properties of the designed mutant. MD simulations suggest that the conformational preferences of the handle region correlate with the conformations of Trp85, Thr92, and Val101. NMR residual dipolar coupling constants, three-bond scalar coupling constants, and chemical shifts experimentally define the conformational states of these residues and hence of the handle domain. These NMR parameters correlate with the Michaelis constants for RNHI homologues, confirming the important role of the handle region in the modulation of substrate recognition and illustrating the power of NMR spectroscopy in dissecting the conformational preferences underlying enzyme function.

RNases H: Structure and mechanism

RNases H are a family of endonucleases that hydrolyze RNA residues in various nucleic acids. These enzymes are present in all branches of life, and their counterpart domains are also found in reverse transcriptases (RTs) from retroviruses and retroelements. RNases H are divided into two main classes (RNases H1 and H2 or type 1 and type 2 enzymes) with common structural features of the catalytic domain but different range of substrates for enzymatic cleavage. Additionally, a third class is found in some Archaea and bacteria. Besides distinct cellular functions specific for each type of RNases H, this family of proteins is generally involved in the maintenance of genome stability with overlapping and cooperative role in removal of R-loops thus preventing their accumulation. Extensive biochemical and structural studies of RNases H provided not only a comprehensive and complete picture of their mechanism but also revealed key basic principles of nucleic acid recognition and processing. RNase H1 is present in prokaryotes and eukaryotes and cleaves RNA in RNA/DNA hybrids. Its main function is hybrid removal, notably in the context of R-loops. RNase H2, which is also present in all branches of life, can play a similar role but it also has a specialized function in the cleavage of single ribonucleotides embedded in the DNA. RNase H3 is present in Archaea and bacteria and is closely related to RNase H2 in sequence and structure but has RNase H1-like biochemical properties. This review summarizes the mechanisms of substrate recognition and enzymatic cleavage by different classes of RNases H with particular insights into structural features of nucleic acid binding, specificity towards RNA and/or DNA strands and catalysis.

Thermostability of Enzymes from Molecular Dynamics Simulations

Thermodynamic stability is a central requirement for protein function, and one goal of protein engineering is improvement of stability, particularly for applications in biotechnology. Herein, molecular dynamics simulations are used to predict in vitro thermostability of members of the bacterial ribonuclease HI (RNase H) family of endonucleases. The temperature dependence of the generalized order parameter, S, for four RNase H homologues, from psychrotrophic, mesophilic, and thermophilic organisms, is highly correlated with experimentally determined melting temperatures and with calculated free energies of folding at the midpoint temperature of the simulations. This study provides an approach for in silico mutational screens to improve thermostability of biologically and industrially relevant enzymes.

Thermodynamic system drift in protein evolution

Tracking the evolution of thermostability in resurrected ancestors of a heat-tolerant extremophile protein and its less heat tolerant Escherichia coli homologue shows how thermostability has probably explored different mechanisms of protein stabilization over evolutionary time.

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.

The biophysical properties of proteins must adjust to accommodate environmental temperatures because of the narrow range over which any given protein sequence can remain folded and functional. We compared the evolution of homologous bacterial enzymes (ribonucleases H1) from two lineages: one from Escherichia coli, which live at moderate temperatures, the other from Thermus thermophilus, which live at extremely high temperatures. Our aim was to investigate how these structurally homologous proteins can have such different thermostabilities, unfolding at temperatures that are 20°C apart. We used bioinformatics to reconstruct the sequences of ancestral proteins along each lineage, synthesized the proteins in the lab, and experimentally traced the evolution of ribonuclease H1 stability. While thermostability appears to have been strongly shaped by selection, the biophysical mechanisms used to tune protein stability appear to have varied throughout evolutionary history; this suggests that proteins have wide latitude to explore different mechanisms of stabilization, generating biophysical diversity and opening up new evolutionary pathways.

Evidence from molecular dynamics simulations of conformational preorganization in the ribonuclease H active site

Ribonuclease H1 (RNase H) enzymes are well-conserved endonucleases that are present in all domains of life and are particularly important in the life cycle of retroviruses as domains within reverse transcriptase. Despite extensive study, especially of the E. coli homolog, the interaction of the highly negatively charged active site with catalytically required magnesium ions remains poorly understood. In this work, we describe molecular dynamics simulations of the E. coli homolog in complex with magnesium ions, as well as simulations of other homologs in their apo states. Collectively, these results suggest that the active site is highly rigid in the apo state of all homologs studied and is conformationally preorganized to favor the binding of a magnesium ion. Notably, representatives of bacterial, eukaryotic, and retroviral RNases H all exhibit similar active-site rigidity, suggesting that this dynamic feature is only subtly modulated by amino acid sequence and is primarily imposed by the distinctive RNase H protein fold.

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.

The relationship between enzymatic activity and protein stability has long been a difficult problem in the study of protein biochemistry. Enzymes may undergo structural changes in order to bind substrates, catalyze chemical reactions, and release products, but flexibility often is inversely correlated with thermodynamic stability. Proteins from organisms that are adapted to high temperature can be both more rigid and less active at ambient temperature than their homologs from organisms that grow at lower temperatures. For this reason, studying homologous pairs of proteins from organisms adapted to different thermal environments is a productive way to identify functionally important motions. In this work we perform comparative analyses of molecular dynamics simulations for five ribonuclease H proteins of varying thermal stabilities, isolated from organisms that grow in varying thermal environments. We identify two different mechanisms of motion in a region of the protein that interacts with substrate molecules, suggesting at least two forms of thermal adaptation in this protein family.

Stabilization by fusion to the C-terminus of hyperthermophile Sulfolobus tokodaii RNase HI: a possibility of protein stabilization tag

RNase HI from the hyperthermophile Sulfolobus tokodaii (Sto-RNase HI) is stabilized by its C-terminal residues. In this work, the stabilization effect of the Sto-RNase HI C-terminal residues was investigated in detail by thermodynamic measurements of the stability of variants lacking the disulfide bond (C58/145A), or the six C-terminal residues (ΔC6) and by structural analysis of ΔC6. The results showed that the C-terminal does not affect overall structure and stabilization is caused by local interactions of the C-terminal, suggesting that the C-terminal residues could be used as a “stabilization tag.” The Sto-RNase HI C-terminal residues (-IGCIILT) were introduced as a tag on three proteins. Each chimeric protein was more stable than its wild-type protein. These results suggested the possibility of a simple stabilization technique using a stabilization tag such as Sto-RNase HI C-terminal residues.

Evolution and thermodynamics of the slow unfolding of hyperstable monomeric proteins

BACKGROUND

The unfolding speed of some hyperthermophilic proteins is dramatically lower than that of their mesostable homologs. Ribonuclease HII from the hyperthermophilic archaeon Thermococcus kodakaraensis (Tk-RNase HII) is stabilized by its remarkably slow unfolding rate, whereas RNase HI from the thermophilic bacterium Thermus thermophilus (Tt-RNase HI) unfolds rapidly, comparable with to that of RNase HI from Escherichia coli (Ec-RNase HI).

RESULTS

To clarify whether the difference in the unfolding rate is due to differences in the types of RNase H or differences in proteins from archaea and bacteria, we examined the equilibrium stability and unfolding reaction of RNases HII from the hyperthermophilic bacteria Thermotoga maritima (Tm-RNase HII) and Aquifex aeolicus (Aa-RNase HII) and RNase HI from the hyperthermophilic archaeon Sulfolobus tokodaii (Sto-RNase HI). These proteins from hyperthermophiles are more stable than Ec-RNase HI over all the temperature ranges examined. The observed unfolding speeds of all hyperstable proteins at the different denaturant concentrations studied are much lower than those of Ec-RNase HI, which is in accordance with the familiar slow unfolding of hyperstable proteins. However, the unfolding rate constants of these RNases H in water are dispersed, and the unfolding rate constant of thermophilic archaeal proteins is lower than that of thermophilic bacterial proteins.

CONCLUSIONS

These results suggest that the nature of slow unfolding of thermophilic proteins is determined by the evolutionary history of the organisms involved. The unfolding rate constants in water are related to the amount of buried hydrophobic residues in the tertiary structure.

Sequence-based analysis of protein energy landscapes reveals nonuniform thermal adaptation within the proteome

Thermal adaptation of individual proteins is often achieved through modulating protein stability, with proteins that are adapted to extreme cold environments having increased conformational flexibility when brought to mesophilic conditions. Conversely, proteins adapted to higher temperatures appear less dynamic and are found to be much more stable against thermal denaturation than their mesophilic counterparts. According to the current paradigm, the adaptation of an organism for survival at higher or lower temperatures is facilitated by the adaptation of the component proteins. We note, however, that these observations have been carried out on relatively few proteins. The extent to which the conformational stabilities of all members of the proteome have been modulated for thermal adaptation remains unclear, with no direct experimental strategies to address this issue. Adapted extremophilies are likely to use a multitude of molecular and biophysical strategies for survival and, therefore, evolution of specific biophysical properties of proteins for optimal function may not be necessary for all proteins in the proteome. Using a sequence-based predictor of protein stability, eScape, an in silico examination of several extremophilic proteomes shows a correlation between the collective stability of the proteins and the thermal range of survival for the organism as expected. Unexpectedly, however, the analysis shows that protein thermostability is modified to different extents across the proteome and depends on the functional role for which the protein is involved. Identification of these differences provides unique opportunities to study interdependence within the proteome as well as the role that the proteome plays in the process of evolutionary thermal adaptation.

Ribonuclease H: molecular diversities, substrate binding domains, and catalytic mechanism of the prokaryotic enzymes

The prokaryotic genomes, for which complete nucleotide sequences are available, always contain at least one RNase H gene, indicating that RNase H is ubiquitous in all prokaryotic cells. Coupled with its unique substrate specificity, the enzyme has been expected to play crucial roles in the biochemical processes associated with DNA replication, gene expression and DNA repair. The physiological role of prokaryotic RNases H, especially of type 1 RNases H, has been extensively studied using Escherichia coli strains that are defective in RNase HI activity or overproduce RNase HI. However, it is not fully understood yet. By contrast, significant progress has been made in this decade in identifying novel RNases H with respect to their biochemical properties and structures, and elucidating catalytic mechanism and substrate recognition mechanism of RNase H. We review the results of these studies.

Destabilization of psychrotrophic RNase HI in a localized fashion as revealed by mutational and X-ray crystallographic analyses

The Arg97 --> Gly and Asp136 --> His mutations stabilized So-RNase HI from the psychrotrophic bacterium Shewanella oneidensis MR-1 by 5.4 and 9.7 degrees C, respectively, in T(m), and 3.5 and 6.1 kJ x mol(-1), respectively, in DeltaG(H2O). These mutations also stabilized the So-RNase HI derivative (4x-RNase HI) with quadruple thermostabilizing mutations in an additive manner. As a result, the resultant sextuple mutant protein (6x-RNase HI) was more stable than the wild-type protein by 28.8 degrees C in T(m) and 27.0 kJ x mol(-1) in DeltaG(H2O). To analyse the effects of the mutations on the protein structure, the crystal structure of the 6x-RNase HI protein was determined at 2.5 A resolution. The main chain fold and interactions of the side-chains of the 6x-RNase HI protein were basically identical to those of the wild-type protein, except for the mutation sites. These results indicate that all six mutations independently affect the protein structure, and are consistent with the fact that the thermostabilizing effects of the mutations are roughly additive. The introduction of favourable interactions and the elimination of unfavourable interactions by the mutations contribute to the stabilization of the 6x-RNase HI protein. We propose that So-RNase HI is destabilized when compared with its mesophilic and thermophilic counterparts in a localized fashion by increasing the number of amino acid residues unfavourable for protein stability.