Stabilization of Escherichia coli ribonuclease HI by strategic replacement of amino acid residues with those from the thermophilic counterpart

Reductionism in the study of enzyme adaptation

One of the principal goals of comparative biology is the elucidation of mechanisms by which organisms adapt to different environments. The study of enzyme structure, function, and stability has contributed significantly to this effort, by revealing adaptation at a molecular level. Comparative biochemistry, including enzymology, necessarily pursues a reductionist approach in describing the function and structure of biomolecules, allowing more straightforward study of molecular systems by removing much of the complexity of their biological milieu. Although this reductionism has allowed a remarkable series of discoveries linking chemical processes to metabolism and to whole-organism function in the context of the environment, it also has the potential to mislead when careful consideration is not made of the simplifying assumptions inherent to such research. In this review, a brief history of the growth of enzymology, its reliance on a reductionist philosophy, and its contributions to our understanding of biological systems is given. Examples then are provided of research techniques, based on a reductionist approach, that have advanced our knowledge about enzyme adaptation to environmental stresses, including stability assays, enzyme kinetics, and the impact of solute composition on enzyme function. In each case, the benefits of the reductionist nature of the approach is emphasized, notable advances are described, but potential drawbacks due to inherent oversimplification of the study system are also identified.

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.

Redirection of the Reaction Specificity of a Thermophilic Acetolactate Synthase toward Acetaldehyde Formation

Acetolactate synthase and pyruvate decarboxylase are thiamine pyrophosphate-dependent enzymes that convert pyruvate into acetolactate and acetaldehyde, respectively. Although the former are encoded in the genomes of many thermophiles and hyperthermophiles, the latter has been found only in mesophilic organisms. In this study, the reaction specificity of acetolactate synthase from Thermus thermophilus was redirected to catalyze acetaldehyde formation to develop a thermophilic pyruvate decarboxylase. Error-prone PCR and mutant library screening led to the identification of a quadruple mutant with 3.1-fold higher acetaldehyde-forming activity than the wild-type. Site-directed mutagenesis experiments revealed that the increased activity of the mutant was due to H474R amino acid substitution, which likely generated two new hydrogen bonds near the thiamine pyrophosphate-binding site. These hydrogen bonds might result in the better accessibility of H⁺ to the substrate-cofactor-enzyme intermediate and a shift in the reaction specificity of the enzyme.

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.

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.

Redesign the α/β fold to enhance the stability of mannanase Man23 from Bacillus subtilis

In this work, we engineered the α/β fold of mannanase Man23 based on its molecular structure analysis to obtain more stable variants. By introducing 31 single-site mutations in the α/β fold and shuffling them, the incorporation of four mutations (K178R, K207R, N340R, and S354R) displayed a good balance between high activity and stability at higher temperature and broader pH. This quartet variant was characterized by an almost threefold increased activity and a sevenfold increased stability compared to native mannanase Man23. Our results suggest that such work is safe to increase our target protein stability with no loss of activity.

Structure, stability, and folding of ribonuclease H1 from the moderately thermophilic Chlorobium tepidum: comparison with thermophilic and mesophilic homologues

Proteins from thermophilic organisms are able to function under conditions that render a typical mesophilic protein inactive. Pairwise comparisons of homologous mesophilic and thermophilic proteins can help to identify the energetic features of a protein's energy landscape that lead to such thermostability. Previous studies of bacterial ribonucleases H (RNases H) from the thermophile Thermus thermophilus and the mesophile Escherichia coli revealed that the thermostability arises in part from an unusually low change in heat capacity upon unfolding (DeltaC(p)) for the thermophilic protein [Hollien, J., and Marqusee, S. (1999) Biochemistry 38, 3831-3836]. Here, we have further examined how nearly identical proteins can adapt to different thermal constraints by adding a moderately thermophilic homologue to the previously characterized mesophilic and thermophilic pair. We identified a putative RNase H from Chlorobium. tepidum and demonstrated that it is an active RNase H and adopts the RNase H fold. The moderately thermophilic protein has a melting temperature (T(m)) similar to that of the mesophilic homologue yet also has a surprisingly low DeltaC(p), like the thermophilic homologue. This new RNase H folds through a pathway similar to that of the previously studied RNases H. These results suggest that lowering the DeltaC(p) may be a general strategy for achieving thermophilicity for some protein families and implicate the folding core as the major contributor to this effect. It should now be possible to design RNases H that display the desired thermophilic or mesophilic properties, as defined by their DeltaC(p) values, and therefore fine-tune the energy landscape in a predictable fashion.

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

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

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.

Relationship between local structural entropy and protein thermostability

We developed a technique to compute structural entropy directly from protein sequences. We explored the possibility of using structural entropy to identify residues involved in thermal stabilization of various protein families. Examples include methanococcal adenylate kinase, Ribonuclease HI and holocytochrome c(551). Our results show that the positions of the largest structural entropy differences between wild type and mutant usually coincide with the residues relevant to thermostability. We also observed a good linear relationship between the average structural entropy and the melting temperatures for adenylate kinase and its chimeric constructs. To validate this linear relationship, we compiled a large dataset comprised of 1153 sequences and found that most protein families still display similar linear relationships. Our results suggest that the multitude of interactions involved in thermal stabilization may be generalized into the tendency of proteins to maintain local structural conservation. The linear relationship between structural entropy and protein thermostability should be useful in the study of protein thermal stabilization.

Conserved quantitative stability/flexibility relationships (QSFR) in an orthologous RNase H pair

Many reports qualitatively describe conserved stability and flexibility profiles across protein families, but biophysical modeling schemes have not been available to robustly quantify both. Here we investigate an orthologous RNase H pair by using a minimal distance constraint model (DCM). The DCM is an all atom microscopic model [Jacobs and Dallakyan, Biophys J 2005;88(2):903-915] that accurately reproduces heat capacity measurements [Livesay et al., FEBS Lett 2004;576(3):468-476], and is unique in its ability to harmoniously calculate thermodynamic stability and flexibility in practical computing times. Consequently, quantified stability/flexibility relationships (QSFR) can be determined using the DCM. For the first time, a comparative QSFR analysis is performed, serving as a paradigm study to illustrate the utility of a QSFR analysis for elucidating evolutionarily conserved stability and flexibility profiles. Despite global conservation of QSFR profiles, distinct enthalpy-entropy compensation mechanisms are identified between the RNase H pair. In both cases, local flexibility metrics parallel H/D exchange experiments by correctly identifying the folding core and several flexible regions. Remarkably, at appropriately shifted temperatures (e.g., melting temperature), these differences lead to a global conservation in Landau free energy landscapes, which directly relate thermodynamic stability to global flexibility. Using ensemble-based sampling within free energy basins, rigidly, and flexibly correlated regions are quantified through cooperativity correlation plots. Five conserved flexible regions are identified within the structures of the orthologous pair. Evolutionary conservation of these flexibly correlated regions is strongly suggestive of their catalytic importance. Conclusions made herein are demonstrated to be robust with respect to the DCM parameterization.

Thermal denaturation of a recombinant mouse amelogenin: circular dichroism and differential scanning calorimetric studies

Conformational analyses of a recombinant mouse tooth enamel amelogenin (rM179) were performed using circular dichroism (CD), fluorescence, differential scanning calorimetry, and sedimentation equilibrium studies. The results show that the far-UV CD spectra of rM179 at acidic pH and 10 degrees C are different from the spectra of random coil in 6 M GdnHCl. A near-UV CD spectrum of rM179 at 10 degrees C is similar to that of rM179 in 6 M GdnHCl, which indicates that aromatic residues of native structure are exposed to solvent and rotate freely. Far-UV CD values of rM179 at 80 degrees C are different from that of random-coil structure in 6 M GdnHCl, which suggests that rM179 at 80 degrees C has specific secondary structures. A gradual thermal transition was observed by far-UV CD, which is interpreted as a weak cooperative transition from specific secondary structures to other specific secondary structures. The fluorescence emission maximum for the spectrum due to Trp residues in rM179 at 10 degrees C shows the same fluorescence emission maximum as rM179 in 6 M GdnHCl and amino acid Trp, which indicates that the three Trp in rM179 are exposed to solvent. Deconvolution of differential scanning calorimetry curve gives the population of three states (A, I, and C states). These results indicate that three states (A, I, and C) have specific secondary structures, in which hydrophobic and Trp residues are exposed to the solvent. The thermodynamic characteristics of rM179 are unique and different from a typical globular protein, proline-rich peptides, and a molten globule state.

Stabilities of chimeras of hyperthermophilic and mesophilic glycerol kinases constructed by DNA shuffling

Glycerol kinases from Thermococcus kodakaraensis KOD1 (Tk-GK) and Escherichia coli (Ec-GK) greatly differ in thermostability. The temperature (T(1/2)) at which the enzymes lose half of their activity upon incubation for 20 min is 50-55 degrees C for Ec-GK and approximately 95 degrees C for Tk-GK. To examine whether the amino acid substitutions that make Tk-GK more stable than Ec-GK are localized in a limited region, the chimeras of two parental genes encoding Tk-GK and Ec-GK were constructed by DNA shuffling. E. coli cells were transformed with a plasmid library harboring these chimeras and screened for those tht produce chimeric enzymes which are more stable than Ec-GK. Four chimeric enzymes were isolated and purified, and their biochemical properties characterized. Replacement of 83 or 93 residues in the C-terminus of Ec-GK with the corresponding ones of Tk-GK increased the T(1/2) value of Ec-GK by 25-30 degrees C. In contrast, replacement of 85 residues in the N-terminus of Ec-GK with the corresponding ones of Tk-GK reduced the T(1/2) value by 5-10 degrees C. In addition, replacement of 10 residues in the C-terminus of Tk-GK with the corresponding ones of Ec-GK reduced the T(1/2) value ot Tk-GK by approximately 15 degrees C. Measurement of the far-UV CD spectra indicates that the three-dimensional structures of the chimeric enzymes, as well as those of the parent enzymes, are similar to one another. These results suggest that the amino acid substitutions responsible for the high stability of Tk-GK are largely localized in the C-terminal region.

Thermostability of ancestral mutants of Caldococcus noboribetus isocitrate dehydrogenase

We constructed mutant genes of Caldococcus noboribetus isocitrate dehydrogenase containing ancestral amino acid residues that were inferred using the maximal likelihood method and a composite phylogenetic tree of isocitrate dehydrogenase and 3-isopropylmalate dehydrogenase. The mutant genes were expressed in Escherichia coli and the protein products purified. Thermostabilities, reported as the half-inactivation temperatures, for the purified enzymes were determined and compared with that of the wild-type enzyme. Four of the five mutant enzymes have greater thermal stabilities than wild-type isocitrate dehydrogenase. The results are compatible with the hyperthermophilic universal ancestor (commonote) hypothesis. Incorporation of ancestral residues into a modern-day protein sequence can be used to improve protein thermostability.

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.

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.

Folding of bovine pancreatic trypsin inhibitor (BPTI) variants in which almost half the residues are alanine

Recent studies indicate that a fraction of the information contained in an amino acid sequence may be sufficient for specifying a native protein structure. An earlier alanine-scanning experiment conducted on bovine pancreatic trypsin inhibitor (BPTI; 58 residues) suggested that if cumulative mutations have additive effects on protein stability, a native protein structure could be built from BPTI sequences that contained many alanine residues distributed throughout the protein. To test this hypothesis, we designed and produced six BPTI mutants containing from 21 to 29 alanine residues. We found that the melting temperature of mutants containing up to 27 alanine residues (48 % of the total number of residues) could be predicted quite well by the sum of the change in melting temperature for the single mutations. Additionally, these same mutants folded into a native-like structure, as judged by their cooperative thermal denaturation curves and heteronuclear multiple quantum correlation (HMQC) NMR spectra. A BPTI mutant containing 22 alanine residues was further shown by 2D and 3D-NMR to fold into a structure very similar to that of native BPTI, and to be a functional trypsin inhibitor. These results provide insight into the extent to which native protein structure and function can be achieved with a highly simplified amino acid sequence.

Structural distribution of stability in a thermophilic enzyme

Stability parameters for individual residues in Thermus thermophilus cysteine-free RNase H were determined by native state hydrogen exchange, thus providing a unique comparison of regional thermodynamics between thermophilic and mesophilic homologues. The general distribution of stability in the thermophilic protein is similar to that of its mesophilic homologue, with a proportional increase in stability for almost all residues. As a consequence, the residue-specific stabilities of the two proteins are remarkably similar under conditions where their global stabilities are the same. These results indicate that T. thermophilus RNase H is stabilized in a delocalized fashion, preserving a finely tuned balance of stabilizing interactions throughout the structure. Therefore, although protein stability can be altered by single amino acid substitution, evolution for optimal function may require more subtle and delocalized mechanisms.

Important amino acid properties for enhanced thermostability from mesophilic to thermophilic proteins

Understanding the role of various interactions in enhancing the thermostability of proteins is important not only for clarifying the mechanism of protein stability but also for designing stable proteins. In this work, we have analyzed the thermostability of 16 different families by comparing mesophilic and thermophilic proteins with 48 various physicochemical, energetic and conformational properties. We found that the increase in shape, s (location of branch point in side chain) increases the thermostability, whereas, an opposite trend is observed for Gibbs free energy change of hydration for native proteins, GhN, in 14 families. A good correlation is observed between these two properties and the simultaneous increases of -GhN and s is necessary to enhance the thermostability from mesophile to thermophile. The increase in shape, which tends to increase with increasing number of carbon atoms both for polar and non-polar residues, may generate more packing and compactness, and the position of beta and higher order branches may be important for better packing. On the other hand, the increase in -GhN in thermophilic proteins increases the solubility of the proteins. This tendency counterbalances the increases in insolubility and unfolding heat capacity change due to the increase in the number of carbon atoms. Thus, the present results suggest that the stability of thermophilic proteins may be achieved by a balance between better packing and solubility.

Thermozymes

Enzymes synthesized by thermophiles (organisms with optimal growth temperatures > 60 degrees C) and hyperthermophiles (optimal growth temperatures > 80 degrees C) are typically thermostable (resistant to irreversible inactivation at high temperatures) and thermophilic (optimally active at high temperatures, i.e., > 60 degrees C). These enzymes, called thermozymes, share catalytic mechanisms with their mesophilic counterparts. When cloned and expressed in mesophilic hosts, thermozymes usually retain their thermal properties, suggesting that these properties are genetically encoded. Sequence alignments, amino acid content comparisons, and crystal structure comparisons indicate that thermozymes are, indeed, very similar to mesophilic enzymes. No obvious sequence or structural features account for enzyme thermostability and thermophilicity. Thermostability and thermophilicity molecular mechanisms are varied, differing from enzyme to enzyme. Thermostability and thermophilicity are usually caused by the accumulation of numerous subtle sequence differences. This review concentrates on the mechanisms involved in enzyme thermostability and thermophilicity. Their relationships with protein rigidity and flexibility and with protein folding and unfolding are discussed. Intrinsic stabilizing forces (e.g., salt bridges, hydrogen bonds, hydrophobic interactions) and extrinsic stabilizing factors are examined. Finally, thermozymes' potential as catalysts for industrial processes and specialty uses are discussed, and lines of development (through new applications, and protein engineering) are also proposed.

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.

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.

The refined crystal structure of Bacillus cereus oligo-1,6-glucosidase at 2.0 A resolution: structural characterization of proline-substitution sites for protein thermostabilization

The crystal structure of oligo-1,6-glucosidase (dextrin 6-alpha-glucanohydrolase, EC 3.2.1.10) from Bacillus cereus ATCC7064 has been refined to 2.0 A resolution with an R-factor of 19.6% for 43,328 reflections. The final model contains 4646 protein atoms and 221 ordered water molecules with respective root-mean-square deviations of 0.015 A for bond lengths and of 3.166 degrees for bond angles from the ideal values. The structure consists of three domains: the N-terminal domain (residues 1 to 104 and 175 to 480), the subdomain (residues 105 to 174) and the C-terminal domain (residues 481 to 558). The N-terminal domain takes a (beta/alpha)8-barrel structure with additions of an alpha-helix (N alpha6') between the sixth strand Nbeta6 and the sixth helix N alpha6, an alpha-helix (N alpha7') between the seventh strand Nbeta7 and the seventh helix N alpha7 and three alpha-helices (N alpha8', N alpha8" and N alpha8'" between the eighth strand Nbeta8 and the eighth helix N alpha8. The subdomain is composed of an alpha-helix, a three-stranded antiparallel beta-sheet, and long intervening loops. The C-terminal domain has a beta-barrel structure of eight antiparallel beta-strands folded in double Greek key motifs, which is distorted in the sixth strand Cbeta6. Three catalytic residues, Asp199, Glu255 and Asp329, are located at the bottom of a deep cleft formed by the subdomain and a cluster of the two additional alpha-helices N alpha8' and N alpha8" in the (beta/alpha)8-barrel. The refined structure reveals the locations of 21 proline-substitution sites that are expected to be critical to protein thermostabilization from a sequence comparison among three Bacillus oligo-1,6-glucosidases with different thermostability. These sites lie in loops, beta-turns and alpha-helices, in order of frequency, except for Cys515 in the fourth beta-strand Cbeta4 of the C-terminal domain. The residues in beta-turns (Lys121, Glu208, Pro257, Glu290, Pro443, Lys457 and Glu487) are all found at their second positions, and those in alpha-helices (Asn109, Glu175, Thr261 and Ile403) are present at their N1 positions of the first helical turns. Those residues in both secondary structures adopt phi and phi values favorable for proline substitution. Residues preceding the 21 sites are mostly conserved upon proline occurrence at these 21 sites in more thermostable Bacillus oligo-1,6-glucosidases. These structural features with respect to the 21 sites indicate that the sites in beta-turns and alpha-helices have more essential prerequisites for proline substitution to thermostabilize the protein than those in loops. This well supports the previous finding that the replacement at the appropriate positions in beta-turns or alpha-helices is the most effective for protein thermostabilization by proline substitution.

Thermal stability of Escherichia coli ribonuclease HI and its active site mutants in the presence and absence of the Mg2+ ion. Proposal of a novel catalytic role for Glu48

Escherichia coli ribonuclease HI, which requires divalent cations (Mg2+ or Mn2+) for activity, was thermostabilized by 2.6-3.0 kcal/mol in the presence of the Mg2+, Mn2+, or Ca2+ ion, probably because the negative charge repulsion around the active site was canceled upon the binding of these metal ions. The dissociation constants were determined to be 0.71 mM for Mg2+, 0.035 mM for Mn2+, and 0.16 mM for Ca2+. Likewise, various active site mutants at Asp10, Glu48, Asp70, or Asp134 were thermostabilized by 0.4-3.0 kcal/mol in the presence of the Mg2+ ion, suggesting that this ion binds to these mutant proteins as well. The dissociation constants of Mg2+ were determined to be 9.8 mM for D10N, 1.1 mM for E48Q, 18.8 mM for D70N, and 1.8 mM for D134N. Thus, the mutation of Asp10 or Asp70 to Asn considerably impairs the Mg2+ binding, whereas the mutation of Glu48 to Gln or Asp134 to Asn does not. Comparison of the thermal stability of the mutant proteins with that of the wild-type protein in the absence of the Mg2+ ion suggests that the negative charge repulsion between Asp10 and Asp70 is responsible for the binding of the metal cofactor. Glu48 may be required to anchor a water molecule, which functions as a general acid.

Reconstitution of Escherichia coli RNase HI from the N-fragment with high helicity and the C-fragment with a disordered structure

The Escherichia coli RNase HI variant with the Lys86-->Ala mutation is purified in two forms, as nicked and intact proteins. The nicked K86A protein, in which the N-fragment (Met1-Lys87) and the C-fragment (Arg88-Val155) remain associated, is enzymatically active. These N- and C-fragments were isolated and examined for reassociation. These peptides did not associate to form the nicked K86A protein at pH 3.0 in the absence of salt, but were associated, with a yield of 30-80%, when the pH was raised to 5.5 or when salt was added. Measurements of the CD spectra show that the alpha-helices are partially formed in the N-fragment at pH 3.0 in the absence of salt and are almost fully formed either at pH 5.5 or at pH 3.0 in the presence of 0.15 M NaCl. In contrast, the C-fragment remains almost fully disordered under these conditions. The N-fragment with this high (native-like) helicity shows the characteristics of a molten globule with respect to the content of the secondary and tertiary structures, the ability to bind a fluorescent probe (1-anilinonaphthalene-8-sulfonic acid), and the behavior on the thermal transition. These results suggest that the N-fragment contains an initial folding site, probably the alpha I-helix, and the completion of the folding in this site provides a surface that facilitates the folding of the C-fragment. This folding process may represent that of the intact RNase HI molecule.

Multiple proline substitutions cumulatively thermostabilize Bacillus cereus ATCC7064 oligo-1,6-glucosidase. Irrefragable proof supporting the proline rule

Nine residues of Bacillus cereus ATCC7064 oligo-1,6-glucosidase were replaced stepwise with proline residues. Of the nine residues, Lys121, Glu208 and Glu290 were at second sites of beta turns; Asn109, Glu175 and Thr261 were at N-caps of alpha helices; Glu216, Glu270 and Glu378 were in coils within loops. The replacements were carried out in the order, Lys121-->Pro, Glu175-->Pro, Glu290-->Pro, Glu208-->Pro, Glu270-->Pro, Glu378-->Pro, Thr261-->Pro, Glu216-->Pro and Asn109-->Pro. The resultant nine active mutant enzymes contained 1-9 more proline residues than B. cereus oligo-1,6-glucosidase. The thermostability of these mutants was additively enhanced with the increase in the number of proline residues introduced. The increase in the thermostability was most remarkable when proline residues were introduced at second sites of beta turns or at N-caps of alpha helices. The above results afforded irrefragable proof for the proline rule as an effective principle for increasing protein thermostability [Suzuki, Y., Oishi, K., Nakano, H. & Nagayama, T. (1987) Appl. Microbiol. Biotechnol. 26, 546-551].

Investigating the role of conserved residue Asp134 in Escherichia coli ribonuclease HI by site-directed random mutagenesis

The role of the conserved Asp134 residue in Escherichia coli ribonuclease HI, which is located at the center of the alpha V helix and lies close to the active site, was analyzed by means of site-directed random mutagenesis. Mutant rnhA genes encoding proteins with ribonuclease H activities were screened by their ability to suppress the ribonuclease-H-dependent, temperature-sensitive growth phenotype of E. coli strain MIC3001. Based on the DNA sequences, nine mutant proteins were predicted to have ribonuclease H activity in vivo. All of these mutant proteins were purified to homogeneity and examined for enzymic activity and protein stability. Among them, only the mutant proteins [D134H]RNase H and [D134N]RNase H were shown to have considerable ribonuclease H activities. Determination of the kinetic parameters revealed that replacement of Asp134 by amino acid residues other than asparagine and histidine dramatically decreased the enzymic activity without seriously affecting the substrate binding. Determination of the CD spectra indicated that none of the mutations seriously affected secondary and tertiary structure. The protein stability was determined from the thermal denaturation curves. All mutant proteins were more stable than the wild-type protein. Such stabilization effects would be a result of a reduction in the negative charge repulsion between Asp134 and the active-site residues, and/or an enhancement of the stability of the alpha V helix. These results strongly suggest that Asp134 does not contribute to the maintenance of the molecular architecture but the carboxyl oxygen at its delta 1 position impacts catalysis.

Protein engineering for unusual environments

The ability to use proteins in unusual or non-natural environments greatly expands their potential applications in biotechnology. Because natural selection has neither maximized the stability of proteins nor optimized them to function under unusual conditions, there is considerable room for their improvement by protein engineering. Significant advances reported within the past year include a dramatic demonstration of a protein's ability to tolerate changes in its amino acid sequence, large increases in protein stability, and the use of random mutagenesis to obtain novel enzymatic properties. Approaches using random or site-directed mutagenesis have been successful in generating proteins able to function in an extended range of environments.