A novel strategy for stabilization of Escherichia coli ribonuclease HI involving a screen for an intragenic suppressor of carboxyl-terminal deletions

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape-all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs-questioning the differences between them and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the consensus protein derived from our full Ribonuclease H sequence alignment is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted set of sequences is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order sequence correlations using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape - all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs - questioning the differences and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the overall consensus protein is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted region is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order couplings using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Engineering a temperature sensitive tobacco etch virus protease

Since tobacco etch virus protease (TEVp) has a high specificity and efficiency in cleaving its target substrates, many groups have attempted to engineer conditional control of its activity. Temperature induction is widely used for modulating gene function because it has fast temporal response, good penetrability and applicability to many model organisms. Here, we engineered a temperature sensitive TEVp (tsTEVp) by using N-terminal truncations to TEVp that achieved efficient proteolysis on a timescale of 4 h after 30°C induction, while remaining relatively inactive at 37°C. As demonstration, tsTEVp was used to generate temperature-induced biological responses for protein translocation, protein degradation and Ca2+-mediated cellular blebbing. Lastly, tsTEVp and their engineered target substrates could find applications in engineered synthetic biological systems.

Simple models for nonpolar solvation: Parameterization and testing

Implicit solvent models are important for many biomolecular simulations. The polarity of aqueous solvent is essential and qualitatively captured by continuum electrostatics methods like Generalized Born (GB). However, GB does not account for the solvent-induced interactions between exposed hydrophobic sidechains or solute-solvent dispersion interactions. These "nonpolar" effects are often modeled through surface area (SA) energy terms, which lack realism, create mathematical singularities, and have a many-body character. We have explored an alternate, Lazaridis-Karplus (LK) gaussian energy density for nonpolar effects and a dispersion (DI) energy term proposed earlier, associated with GB electrostatics. We parameterized several combinations of GB, SA, LK, and DI energy terms, to reproduce 62 small molecule solvation free energies, 387 protein stability changes due to point mutations, and the structures of 8 protein loops. With optimized parameters, the models all gave similar results, with GBLK and GBDILK giving no performance loss compared to GBSA, and mean errors of 1.7 kcal/mol for the stability changes and 2 Å deviations for the loop conformations. The optimized GBLK model gave poor results in MD of the Trpcage mini-protein, but parameters optimized specifically for MD performed well for Trpcage and three other small proteins. Overall, the LK and DI nonpolar terms are valid alternatives to SA treatments for a range of applications. © 2017 Wiley Periodicals, Inc.

Conformational preferences underlying reduced activity of a thermophilic ribonuclease H

The conformational basis for reduced activity of the thermophilic ribonuclease HI enzyme from Thermus thermophilus, compared to its mesophilic homolog from Escherichia coli, is elucidated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. Explicit-solvent all-atom MD simulations of the two wild-type proteins and an E. coli mutant in which a glycine residue is inserted after position 80 to mimic the T. thermophilus protein reproduce the differences in conformational dynamics determined from (15)N spin-relaxation NMR spectroscopy of three loop regions that surround the active site and contain functionally important residues: the glycine-rich region, the handle region, and the β5/αE loop. Examination of the MD trajectories indicates that the thermophilic protein samples conformations productive for substrate binding and activity less frequently than the mesophilic enzyme, although these differences may manifest as either increased or decreased relative flexibility of the different regions. Additional MD simulations indicate that mutations increasing activity of the T. thermophilus enzyme at mesophilic temperatures do so by reconfiguring the local environments of the mutated sites to more closely resemble active conformations. Taken together, the results show that both locally increased and decreased flexibility contribute to an overall reduction in activity of T. thermophilus ribonuclease H compared to its mesophilic E. coli homolog.

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.

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.

Underexposed polar residues and protein stabilization

Increasing protein stability is interesting for practical reasons and because it tests our understanding of protein energetics. We explore here the feasibility of stabilizing proteins by replacing underexposed polar residues by apolar ones of similar size and shape. We have compared the stability of wild-type apoflavodoxin with that of a few carefully selected mutants carrying Y → F, Q → L, T → V or K → M replacements. Although a clear inverse correlation between native solvent exposures of replaced polar residues and stability of mutants is observed, most mutations fail to stabilize the protein. The promising exceptions are the two Q → L mutations tested, which characteristically combine the greatest reduction in polar burial with the greatest increase in apolar burial relative to wild type. Analysis of published stability data corresponding to a variety of mutant proteins confirms that, unlike Y → F or T → V replacements, Q → L mutations tend to be stabilizing, and it suggests that N → L mutations might be stabilizing as well. On the other hand, we show that the stability changes associated to the apoflavodoxin mutations can be rationalized in terms of differential polar and apolar burials upon folding plus a generic destabilizing penalty term. Simple equations combining these contributions predict stability changes in a large data set of 113 mutants (Y → F, Q → L or T → V) similarly well as more complex algorithms available on the Internet.

Inferring stabilizing mutations from protein phylogenies: application to influenza hemagglutinin

One selection pressure shaping sequence evolution is the requirement that a protein fold with sufficient stability to perform its biological functions. We present a conceptual framework that explains how this requirement causes the probability that a particular amino acid mutation is fixed during evolution to depend on its effect on protein stability. We mathematically formalize this framework to develop a Bayesian approach for inferring the stability effects of individual mutations from homologous protein sequences of known phylogeny. This approach is able to predict published experimentally measured mutational stability effects (ΔΔG values) with an accuracy that exceeds both a state-of-the-art physicochemical modeling program and the sequence-based consensus approach. As a further test, we use our phylogenetic inference approach to predict stabilizing mutations to influenza hemagglutinin. We introduce these mutations into a temperature-sensitive influenza virus with a defect in its hemagglutinin gene and experimentally demonstrate that some of the mutations allow the virus to grow at higher temperatures. Our work therefore describes a powerful new approach for predicting stabilizing mutations that can be successfully applied even to large, complex proteins such as hemagglutinin. This approach also makes a mathematical link between phylogenetics and experimentally measurable protein properties, potentially paving the way for more accurate analyses of molecular evolution.

Mutating a protein frequently causes a change in its stability. As scientists, we often care about these changes because we would like to engineer a protein's stability or understand how its stability is impacted by a naturally occurring mutation. Evolution also cares about mutational stability changes, because a basic evolutionary requirement is that proteins remain sufficiently stable to perform their biological functions. Our work is based on the idea that it should be possible to use the fact that evolution selects for stability to infer from related proteins the effects of specific mutations. We show that we can indeed use protein evolutionary histories to computationally predict previously measured mutational stability changes more accurately than methods based on either of the two main existing strategies. We then test whether we can predict mutations that increase the stability of hemagglutinin, an influenza protein whose rapid evolution is partly responsible for the ability of this virus to cause yearly epidemics. We experimentally create viruses carrying predicted stabilizing mutations and find that several do in fact improve the virus's ability to grow at higher temperatures. Our computational approach may therefore be of use in understanding the evolution of this medically important virus.

Identification of the gene encoding a type 1 RNase H with an N-terminal double-stranded RNA binding domain from a psychrotrophic bacterium

The gene encoding a bacterial type 1 RNase H, termed RBD-RNase HI, was cloned from the psychrotrophic bacterium Shewanella sp. SIB1, overproduced in Escherichia coli, and the recombinant protein was purified and biochemically characterized. SIB1 RBD-RNase HI consists of 262 amino acid residues and shows amino acid sequence identities of 26% to SIB1 RNase HI, 17% to E. coli RNase HI, and 32% to human RNase H1. SIB1 RBD-RNase HI has a double-stranded RNA binding domain (RBD) at the N-terminus, which is commonly present at the N-termini of eukaryotic type 1 RNases H. Gel mobility shift assay indicated that this domain binds to an RNA/DNA hybrid in an isolated form, suggesting that this domain is involved in substrate binding. SIB1 RBD-RNase HI exhibited the enzymatic activity both in vitro and in vivo. Its optimum pH and metal ion requirement were similar to those of SIB1 RNase HI, E. coli RNase HI, and human RNase H1. The specific activity of SIB1 RBD-RNase HI was comparable to that of E. coli RNase HI and was much higher than those of SIB1 RNase HI and human RNase H1. SIB1 RBD-RNase HI showed poor cleavage-site specificity for oligomeric substrates. SIB1 RBD-RNase HI was less stable than E. coli RNase HI but was as stable as human RNase H1. Database searches indicate that several bacteria and archaea contain an RBD-RNase HI. This is the first report on the biochemical characterization of RBD-RNase HI.

The screening of the second-site suppressor mutations of the common p53 mutants

Second-site suppressor (SSS) mutations in p53 found by random mutagenesis have shown to restore the inactivated function of some tumor-derived p53. To screen novel SSS mutations against common mutant p53s, intragenic second-site (SS) mutations were introduced into mutant p53 cDNA in a comprehensive manner by using a p53 missense mutation library. The resulting mutant p53s with background and SS mutations were assayed for their ability to restore the p53 transactivation function in both yeast and human cell systems. We identified 12 novel SSS mutations including H178Y against a common mutation G245S. Surprisingly, the G245S phenotype is rescued when coexpressed with p53 bearing the H178Y mutation. This result indicated that there is a possibility that intragenic suppressor mutations might restore the protein function in an intermolecular manner. The intermolecular mechanism may lead to novel strategies for restoring inactivated p53 function and tumor suppression in cancer treatment.

Identification of RNase HII from psychrotrophic bacterium, Shewanella sp. SIB1 as a high-activity type RNase H

The gene encoding RNase HII from the psychrotrophic bacterium, Shewanella sp. SIB1 was cloned, overexpressed in Escherichia coli, and the recombinant protein was purified and biochemically characterized. SIB1 RNase HII is a monomeric protein with 212 amino acid residues and shows an amino acid sequence identity of 64% to E. coli RNase HII. The enzymatic properties of SIB1 RNase HII, such as metal ion preference, pH optimum, and cleavage mode of substrate, were similar to those of E. coli RNase HII. SIB1 RNase HII was less stable than E. coli RNase HII, but the difference was marginal. The half-lives of SIB1 and E. coli RNases HII at 30 degrees C were approximately 30 and 45 min, respectively. The midpoint of the urea denaturation curve and optimum temperature of SIB1 RNase HII were lower than those of E. coli RNase HII by approximately 0.2 M and approximately 5 degrees C, respectively. However, SIB1 RNase HII was much more active than E. coli RNase HII at all temperatures studied. The specific activity of SIB1 RNase HII at 30 degrees C was 20 times that of E. coli RNase HII. Because SIB1 RNase HII was also much more active than SIB1 RNase HI, RNases HI and HII represent low- and high-activity type RNases H, respectively, in SIB1. In contrast, RNases HI and HII represent high- and low-activity type RNases H, respectively, in E. coli. We propose that bacterial cells usually contain low- and high-activity type RNases H, but these types are not correlated with RNase H families.

Stabilization of E. coli Ribonuclease HI by the ‘stability profile of mutant protein’ (SPMP)-inspired random and non-random mutagenesis

The change in the structural stability of Escherichia coli ribonuclease HI (RNase HI) due to single amino acid substitutions has been estimated computationally by the stability profile of mutant protein (SPMP) [Ota, M., Kanaya, S. Nishikawa, K., 1995. Desk-top analysis of the structural stability of various point mutations introduced into ribonuclease H. J. Mol. Biol. 248, 733-738]. As well, an effective strategy using random mutagenesis and genetic selection has been developed to obtain E. coli RNase HI mutants with enhanced thermostability [Haruki, M., Noguchi, E., Akasako, A., Oobatake, M., Itaya, M., Kanaya, S., 1994. A novel strategy for stabilization of Escherichia coli ribonuclease HI involving a screen for an intragenic suppressor of carboxyl-terminal deletions. J. Biol. Chem. 269, 26904-26911]. In this study, both methods were combined: random mutations were individually introduced to Lys99-Val101 on the N-terminus of the alpha-helix IV and the preceding beta-turn, where substitutions of other amino acid residues were expected to significantly increase the stability from SPMP, and then followed by genetic selection. Val101 to Ala, Gln, and Arg mutations were selected by genetic selection. The Val101-->Ala mutation increased the thermal stability of E. coli RNase HI by 2.0 degrees C in Tm at pH 5.5, whereas the Val101-->Gln and Val101-->Arg mutations decreased the thermostability. Separately, the Lys99-->Pro and Asn100-->Gly mutations were also introduced directly. The Lys99-->Pro mutation increased the thermostability of E. coli RNase HI by 1.8 degrees C in Tm at pH 5.5, whereas the Asn100-->Gly mutation decreased the thermostability by 17 degrees C. In addition, the Lys99-->Pro mutation altered the dependence of the enzymatic activity on divalent metal ions.

Structural perturbation and compensation by directed evolution at physiological temperature leads to thermostabilization of beta-lactamase

The choice of protein for use in technical and medical applications is limited by stability issues, making understanding and engineering of stability key. Here, enzyme destabilization by truncation was combined with directed evolution to create stable variants of TEM-1 beta-lactamase. This enzyme was chosen because of its implication in prodrug activation therapy, pathogen resistance to lactam antibiotics, and reporter enzyme bioassays. Removal of five N-terminal residues generated a mutant which did not confer antibiotic resistance at 37 degrees C. Accordingly, the half-life time in vitro was only 7 s at 40 degrees C. However, three cycles comprising random mutagenesis, DNA shuffling, and metabolic selection at 37 degrees C yielded mutants providing resistance levels significantly higher than that of the wild type. These mutants demonstrated increased thermoactivity and thermostability in time-resolved kinetics at various temperatures. Chemical denaturation revealed improved thermodynamic stabilities of a three-state unfolding pathway exceeding wild-type construct stability. Elongation of one optimized deletion mutant to full length increased its stability even further. Compared to that of the wild type, the temperature optimum was shifted from 35 to 50 degrees C, and the beginning of heat inactivation increased by 20 degrees C while full activity at low temperatures was maintained. We attribute these effects mainly to two independently acting boundary interface residue exchanges (M182T and A224V). Structural perturbation by terminal truncation, evolutionary compensation at physiological temperatures, and elongation is an efficient way to analyze and improve thermostability without the need for high-temperature selection, structural information, or homologous proteins.

Identification of the first archaeal Type 1 RNase H gene from Halobacterium sp. NRC-1: archaeal RNase HI can cleave an RNA-DNA junction

All the archaeal genomes sequenced to date contain a single Type 2 RNase H gene. We found that the genome of a halophilic archaeon, Halobacterium sp. NRC-1, contains an open reading frame with similarity to Type 1 RNase H. The protein encoded by the Vng0255c gene, possessed amino acid sequence identities of 33% with Escherichia coli RNase HI and 34% with a Bacillus subtilis RNase HI homologue. The B. subtilis RNase HI homologue, however, lacks amino acid sequences corresponding to a basic protrusion region of the E. coli RNase HI, and the Vng0255c has the similar deletion. As this deletion apparently conferred a complete loss of RNase H activity on the B. subtilis RNase HI homologue protein, the Vng0255c product was expected to exhibit no RNase H activity. However, the purified recombinant Vng0255c protein specifically cleaved an RNA strand of the RNA/DNA hybrid in vitro, and when the Vng0255c gene was expressed in an E. coli strain MIC2067 it could suppress the temperature-sensitive growth defect associated with the loss of RNase H enzymes of this strain. These results in vitro and in vivo strongly indicate that the Halobacterium Vng0255c is the first archaeal Type 1 RNase H. This enzyme, unlike other Type 1 RNases H, was able to cleave an Okazaki fragment-like substrate at the junction between the 3'-side of ribonucleotide and 5'-side of deoxyribonucleotide. It is likely that the archaeal Type 1 RNase H plays a role in the removal of the last ribonucleotide of the RNA primer from the Okazaki fragment during DNA replication.

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.

Multiple ribonuclease H-encoding genes in the Caenorhabditis elegans genome contrasts with the two typical ribonuclease H-encoding genes in the human genome

Database searches of the Caenorhabditis elegans and human genomic DNA sequences revealed genes encoding ribonuclease H1 (RNase H1) and RNase H2 in each genome. The human genome contains a single copy of each gene, whereas C. elegans has four genes encoding RNase H1-related proteins and one gene for RNase H2. By analyzing the mRNAs produced from the C. elegans genes, examining the amino acid sequence of the predicted protein, and expressing the proteins in Esherichia coli we have identified two active RNase H1-like proteins. One is similar to other eukaryotic RNases H1, whereas the second RNase H (rnh-1.1) is unique. The rnh-1.0 gene is transcribed as a dicistronic message with three dsRNA-binding domains; the mature mRNA is transspliced with SL2 splice leader and contains only one dsRNA-binding domain. Formation of RNase H1 is further regulated by differential cis-splicing events. A single rnh-2 gene, encoding a protein similar to several other eukaryotic RNase H2L's, also has been examined. The diversity and enzymatic properties of RNase H homologues are other examples of expansion of protein families in C. elegans. The presence of two RNases H1 in C. elegans suggests that two enzymes are required in this rather simple organism to perform the functions that are accomplished by a single enzyme in more complex organisms. Phylogenetic analysis indicates that the active C. elegans RNases H1 are distantly related to one another and that the C. elegans RNase H1 is more closely related to the human RNase H1. The database searches also suggest that RNase H domains of LTR-retrotransposons in C. elegans are quite unrelated to cellular RNases H1, but numerous RNase H domains of human endogenous retroviruses are more closely related to cellular RNases H.

Enhancement of the enzymatic activity of ribonuclease HI from Thermus thermophilus HB8 with a suppressor mutation method

A genetic method for isolating a mutant enzyme of ribonuclease HI (RNase HI) from Thermus thermophilus HB8 with enhanced activity at moderate temperatures was developed. T. thermophilus RNase HI has an ability to complement the RNase H-dependent temperature-sensitive (ts) growth phenotype of Escherichia coli MIC3001. However, this complementation ability was greatly reduced by replacing Asp(134), which is one of the active site residues, with His, probably due to a reduction in the catalytic activity. Random mutagenesis of the gene encoding the resultant D134H enzyme, followed by screening for second-site revertants, allowed us to isolate three single mutations (Ala(12) --> Ser, Lys(75) --> Met, and Ala(77) --> Pro) that restore the normal complementation ability to the D134H enzyme. These mutations were individually or simultaneously introduced into the wild-type enzyme, and the kinetic parameters of the resultant mutant enzymes for the hydrolysis of a DNA-RNA-DNA/DNA substrate were determined at 30 degrees C. Each mutation increased the k(cat)/K(m) value of the wild-type enzyme by 2.1-4.8-fold. The effects of the mutations on the enzymatic activity were roughly cumulative, and the combination of these three mutations increased the k(cat)/K(m) value of the wild-type enzyme by 40-fold (5.5-fold in k(cat)). Measurement of thermal stability of the mutant enzymes with circular dichroism spectroscopy in the presence of 1 M guanidine hydrochloride and 1 mM dithiothreitol showed that the T(m) value of the triple mutant enzyme, in which all three mutations were combined, was comparable to that of the wild-type enzyme (75.0 vs 77.4 degrees C). These results demonstrate that the activity of a thermophilic enzyme can be improved without a cost of protein stability.

Isolation of RNase H genes that are essential for growth of Bacillus subtilis 168

Two genes encoding functional RNase H (EC 3.1.26.4) were isolated from a gram-positive bacterium, Bacillus subtilis 168. Two DNA clones exhibiting RNase H activities both in vivo and in vitro were obtained from a B. subtilis DNA library. One (28.2 kDa) revealed high similarity to Escherichia coli RNase HII, encoded by the rnhB gene. The other (33.9 kDa) was designated rnhC and encodes B. subtilis RNase HIII. The B. subtilis genome has an rnhA homologue, the product of which has not yet shown RNase H activity. Analyses of all three B. subtilis genes revealed that rnhB and rnhC cannot be simultaneously inactivated. This observation indicated that in B. subtilis both the rnhB and rnhC products are involved in certain essential cellular processes that are different from those suggested by E. coli rnh mutation studies. Sequence conservation between the rnhB and rnhC genes implies that both originated from a single ancestral RNase H gene. The roles of bacterial RNase H may be indicated by the single rnhC homologue in the small genome of Mycoplasma species.

Identification of the genes encoding Mn2+-dependent RNase HII and Mg2+-dependent RNase HIII from Bacillus subtilis: classification of RNases H into three families

Database searches indicated that the genome of Bacillus subtilis contains three different genes encoding RNase H homologues. The ypdQ gene encodes an RNase HI homologue with 132 amino acid residues, whereas the rnh and ysgB genes encode RNase HII homologues with 255 and 313 amino acid residues, respectively. RNases HI and HII show no significant sequence similarity. These genes were individually expressed in Escherichia coli; the recombinant proteins were purified, and their enzymatic properties were compared with those of E. coli RNases HI and HII. We found that the ypdQ gene product showed no RNase H activity. The 2.2 kb pair genomic DNA containing this gene did not suppress the RNase H deficiency of an E. coli rnhA mutant, indicating that this gene product shows no RNase H activity in vivo as well. In contrast, the rnh (rnhB) gene product (RNase HII) showed a preference for Mn2+, as did E. coli RNase HII, whereas the ysgB (rnhC) gene product (RNase HIII) exhibited a Mg2+-dependent RNase H activity. Oligomeric substrates digested with these enzymes indicate similar recognition of these substrates by B. subtilis and E. coli RNases HII. Likewise, B. subtilis RNase HIII and E. coli RNase HI have generated similar products. These results suggest that B. subtilis RNases HII and HIII may be functionally similar to E. coli RNases HII and HI, respectively. We propose that Mn2+-dependent RNase HII is universally present in various organisms and Mg2+-dependent RNase HIII, which may have evolved from RNase HII, functions as a substitute for RNase HI.

Pressure-denatured state of Escherichia coli ribonuclease HI as monitored by Fourier transform infrared and NMR spectroscopy

Pressure denaturation of Escherichia coli ribonuclease HI (RNase HI) was studied by Fourier transform infrared (FTIR) and two-dimensional NMR spectroscopy at pD* 3.0 and 25 degrees C. A reversible transition in the pressure range of 0.1-1090 MPa was observed with second-derivative FTIR experiments. A cooperative and gradual denaturation, involving both the secondary and tertiary structures, was observed between 240 and 450 MPa. The two peaks at 1629 and 1652 cm(-1), due to beta-strands and alpha-helices, respectively, did not fully disappear after the denaturation, and are different from the spectra of the random coil peptides. The hydrogen-deuterium exchange rates of the individual backbone amide protons were determined by heteronuclear NMR combined with the pressure-jump technique at 500, 650, and 850 MPa. Although most of the amides protected in the native structure are also highly protected in the pressure-denatured state, the rate constants (0.048 +/- 0.007 min(-1)) for the amide protons at 500 MPa are similar regardless of their locations, which is an indication of the EX1 mechanism of hydrogen-deuterium exchange. The pressure-denatured state of RNase HI at 500 MPa represents a novel denatured state, which is different from a typical molten globule state at atmospheric pressure (0.1 MPa), from the viewpoint of the homogeneous rate constants. The observations at 650 MPa are essentially the same as those at 500 MPa. However, at 850 MPa, the amide exchange rates for the highly hydrophobic C-terminal half of alpha-helix I are significantly slower than those for the other part of the protein, which can be interpreted as a hydrophobic collapse centered at the C-terminal half of alpha-helix I.

Stabilization of ribonuclease HI from Thermus thermophilus HB8 by the spontaneous formation of an intramolecular disulfide bond

To identify factors that contribute to the thermal stability of ribonuclease HI (RNase HI) from Thermus thermophilus HB8, protein variants with a series of carboxyl-terminal truncations and Cys --> Ala mutations were constructed, and their thermal denaturations were analyzed by CD. The results indicate that Cys41 and Cys149 contribute to the protein stability, probably through the formation of a disulfide bond. Peptide mapping analysis for the mutant protein with only two cysteine residues, at positions 41 and 149, indicated that this disulfide bond is partially formed in a protein purified from Escherichia coli in the absence of a reducing reagent but is fully formed in a thermally denatured protein. These results suggest that the thermal stability of T. thermophilus RNase HI, determined in the absence of a reducing reagent, reflects that of an oxidized form of the protein. Comparison of the thermal stabilities and the enzymatic activities of the wild-type and truncated proteins, determined in the presence and absence of a reducing reagent, indicates that the formation of this disulfide bond increases the thermal stability of the protein by 6-7 degreesC in Tm and approximately 3 kcal/mol in DeltaG without seriously affecting the enzymatic activity. Since T. thermophilus RNase HI is present in a reducing environment in cells, this disulfide bond probably is not formed in vivo but is spontaneously formed in vitro in the absence of a reducing reagent.

Structural stability and internal motions of Escherichia coli ribonuclease HI: 15N relaxation and hydrogen-deuterium exchange analyses

The relationship between the structural stability and the internal motions of proteins was investigated through measurements of 15N relaxation and hydrogen-deuterium exchange rates of ribonuclease HI from Escherichia coli and its thermostable quintuple mutant (Gly23-->Ala, His62-->Pro, Val74-->Leu, Lys95-->Gly, and Asp134-->His), which has a higher melting temperature by 20.2 degreesC. For most of the residues, the generalized order parameters (S2) obtained from 15N relaxation analyses as well as the localized hydrogen-bond-breaking motions (local breathing) observed as fast H-D exchange rates were largely unaffected by the mutations, indicating no global mutational effect on the internal motions. Several local mutational effects were observed for residues close to the mutation sites as follows. The S2 value significantly increased for Lys96 and Val98, which indicated that motions on the pico- to nanosecond time-scale became restricted within a protruding region including the Lys95-->Gly mutation site. In contrast, slight decreases in S2, and drastic increases in the chemical exchange motion on the micro- to millisecond time-scale (Deltaex), were observed for residues located in the joining region between the protrusion and the major domain of the protein. These changes may be caused by the elimination of the bulky Lys95 side-chain at the center of the protrusion. Deltaex observed for residues in alpha-helix I of the wild-type protein was reduced for the mutant, probably because a cavity in the hydrophobic core is filled by the Val74-->Leu mutation. The local breathing at position 134 was restricted by the Asp134-->His mutation, probably because the reduction of the negative charge repulsion contributes to the stability of the native major conformation relative to the breathing conformations around position 134.

Directed evolution of enzyme catalysts

Directed enzyme evolution has emerged in the past few years as a powerful alternative to rational approaches for engineering biocatalysts. Prerequisites for successful directed evolution are functional expression in a suitable microbial host, a rapid screen for the desired feature(s) and a well-thought-out working strategy for navigating protein landscapes. The rapidly growing body of literature on enzyme evolution in vitro includes techniques for creating and searching combinatorial enzyme libraries, as well as several successful examples of different evolutionary strategies being used.

Conformational stabilities of Escherichia coli RNase HI variants with a series of amino acid substitutions at a cavity within the hydrophobic core

Escherichia coli ribonuclease HI has a cavity within the hydrophobic core. Two core residues, Ala52 and Val74, resided at both ends of this cavity. We have constructed a series of single mutant proteins at Ala52, and double mutant proteins, in which Ala52 was replaced by Gly, Val, Ile, Leu, or Phe, and Val74 was replaced by Ala or Leu. All of these mutant proteins, except for A52W, A52R, and A52G/V74A, were overproduced and purified. Measurement of the thermal denaturations of the proteins at pH 3.2 by CD suggests that the cavity is large enough to accommodate three methyl or methylene groups without creating serious strains. A correlation was observed between the protein stability and the hydrophobicity of the substituted residue. As a result, a number of the mutant proteins were more stable than the wild-type protein. The stabilities of the mutant proteins with charged or extremely bulky residues at the cavity were lower than those expected from the hydrophobicities of the substituted residues, suggesting that considerable strains are created at the mutation sites in these mutant proteins. However, examination of the far- and near-UV CD spectra and the enzymatic activities suggest that all of the mutant proteins have structures similar to that of the wild-type protein. These results suggest that the cavity in the hydrophobic core of E. coli RNase HI is conformationally fairly stable. This may be the reason why the cavity-filling mutations effectively increase the thermal stability of this protein.

Importance of the C-terminal helix to the stability and enzymatic activity of Escherichia coli ribonuclease H

The ribonuclease H (RNase H) family of enzymes selectively degrades the RNA strand of RNA-DNA hybrids. This activity is essential for retroviruses such as HIV and resides in a domain of the larger reverse transcriptase molecule. RNase H from Escherichia coli is the best-characterized member of the family and serves as a model for structure/function studies. Despite having almost identical alpha + beta folds, the isolated domain from HIV is inactive and much less stable than the E. coli homolog. The HIV domain also shows increased disorder in its C-terminal regions (E-helix and His-containing loop). We investigated the importance of this region by studying a variant of E. coli RNase H lacking these elements (RNHdeltaE). Despite the elimination of 33 of 155 residues (including a complete helix), this C-terminal deletion mutant folds cooperatively as a subdomain. Surprisingly, this protein lacking residues near the active site retains weak Mn2+-dependent activity. A peptide corresponding to the deleted E-helix is helical in isolation and stimulates the activity of the deletion mutant in vitro. These results have implications for the catalytic mechanism of RNase H and drug design targeted to HIV reverse transcriptase.

The putative substrate recognition loop of Escherichia coli ribonuclease H is not essential for activity

The RNase H family of enzymes catalyzes the hydrolysis of RNA from RNA DNA hybrids in a divalent metal-dependent fashion. To date, structure/function studies have focused on two members of this family: Escherichia coli RNase HI, a small monomeric protein; and human immunodeficiency virus, type I (HIV) RNase H, a domain of HIV reverse transcriptase. The isolated RNase H domain from HIV reverse transcriptase can be expressed independently and shares significant structural homology with its E. coli homologue; however, unlike the bacterial protein, it is inactive. The most notable difference between the inactive domain from HIV and the active E. coli protein is a basic helix/loop sequence, present in E. coli but absent from the HIV homologue. Substitution of this basic region into the HIV domain partially restores its activity and increases its thermodynamic stability. By deleting the basic helix/loop region, we have modeled the structural difference between these two polypeptides onto the E. coli homologue. Surprisingly, the resulting mutant protein is active in Mn2+-dependent fashion. Therefore, the basic helix/loop is not required for RNase H activity.

Characterization of the internal motions of Escherichia coli ribonuclease HI by a combination of 15N-NMR relaxation analysis and molecular dynamics simulation: examination of dynamic models

The backbone dynamics of Escherichia coli ribonuclease HI (RNase HI) in the picosecond to nanosecond time scale were characterized by a combination of measurements of 15N-NMR relaxation (T1, T2, and NOE), analyzed by a model-free approach, and molecular dynamics (MD) simulation in water. The MD simulations in water were carried out with long-range Coulomb interactions to avoid the artificial fluctuation caused by the cutoff approximation. The model-free analysis of the 15N-NMR relaxation indicated that RNase HI has a rotational correlation time of 10.9 ns at 27 degrees C. The generalized order parameter (S2) for the internal motions varied from 0.15 to 1.0, with an average value of 0.85, which is much larger than that of the RNase H domain of HIV-1 reverse transcriptase (0.78). Large internal motions (small order parameters) were observed in the N-terminal region (Leu2-Lys3), the loop between beta-strands A and B (Cys13-Gly15), the turn between alpha-helix I and beta-strand D (Glu61, His62), the loop between beta-strand D and alpha-helix II (Asp70-Tyr71), the loop between alpha-helices III and IV (Ala93-Lys96), the loop between beta-strand E and alpha-helix V (Gly123-His127), and the C-terminal region (Gln152-Val155). The effective correlation time observed in these regions varied from 0.45 ns (Glu61, Lys96) to 2.2 ns (Leu14). The order parameters calculated from the MD agreed well with those from the NMR experiment, with a few exceptions. The distributions of most of the backbone N-H vectors obtained by MD are approximately consistent with the diffusion-in-a-cone model. These distributions, however, were elliptic, with a long axis perpendicular to the plane defined by the N-H and N-C alpha vectors. Distributions supporting the axial fluctuation model or the jump-between-two-cones model were also observed in the MD simulation.