Conferring thermostability to mesophilic proteins through optimized electrostatic surfaces

Recently, there have been several experimental reports of proteins displaying appreciable stability gains through mutation of one or two amino acid residues. Here, we employ a simple theoretical model to quickly screen mutant structures for increased thermostability through optimization of the protein's electrostatic surface. Our results are able to reproduce the experimental observation that elimination of like-charge repulsions and creation of opposite-charge attractions on the protein surface is an efficient method to confer thermostability to a mesophilic protein. Using Poisson-Boltzmann electrostatics, we calculate relative protein stabilities for the exhaustive surface mutagenesis of the cold shock, RNase T1, and CheY proteins. Comparison with 25 experimentally characterized cold shock protein mutants reveals an average correlation of 0.86. The model is also quantitatively accurate when reproducing the experimental D49A and D49H mutant stabilities of RNase T1. This work represents the first comprehensive in silico screening of mutant candidates likely to confer thermostability to mesophilic proteins through optimization of surface electrostatics. Systematic single mutant, followed by double mutant, screening yields a limited number of mutant structures displaying significant stability gains suitable for subsequent experimental characterization.     

A thermodynamic comparison of mesophilic and thermophilic ribonucleases H

The mechanisms by which thermophilic proteins attain their increased thermostability remain unclear, as usually the sequence and structure of these proteins are very similar to those of their mesophilic homologues. To gain insight into the basis of thermostability, we have determined protein stability curves describing the temperature dependence of the free energy of unfolding for two ribonucleases H, one from the mesophile Escherichia coli and one from the thermophile Thermus thermophilus. The circular dichroism signal was monitored as a function of temperature and guanidinium chloride concentration, and the resulting free energies of unfolding were fit to the Gibbs-Helmholtz equation to obtain a set of thermodynamic parameters for these proteins. Although the maximal stabilities for these proteins occur at similar temperatures, the heat capacity of unfolding for T. thermophilus RNase H is lower, resulting in a smaller temperature dependence of the free energy of unfolding and therefore a higher thermal melting temperature. In addition, the stabilities of these proteins are similar at the optimal growth temperatures for their respective organisms, suggesting that a balance of thermodynamic stability and flexibility is important for function.     

Increasing protein stability: Importance of ΔCp and the denatured state

Increasing the conformational stability of proteins is an important goal for both basic research and industrial applications. In vitro selection has been used successfully to increase protein stability, but more often site‐directed mutagenesis is used to optimize the various forces that contribute to protein stability. In previous studies, we showed that improving electrostatic interactions on the protein surface and improving the β‐turn sequences were good general strategies for increasing protein stability, and used them to increase the stability of RNase Sa. By incorporating seven of these mutations in RNase Sa, we increased the stability by 5.3 kcal/mol. Adding one more mutation, D79F, gave a total increase in stability of 7.7 kcal/mol, and a melting temperature 28°C higher than the wild‐type enzyme. Surprisingly, the D79F mutation lowers the change in heat capacity for folding, ΔC_p, by 0.6 kcal/mol/K. This suggests that this mutation stabilizes structure in the denatured state ensemble. We made other mutants that give some insight into the structure present in the denatured state. Finally, the thermodynamics of folding of these stabilized variants of RNase Sa are compared with those observed for proteins from thermophiles.     

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

Thermus thermophilus ribonuclease H is exceptionally stable against thermal and guanidine hydrochloride denaturations as compared to Escherichia coli ribonuclease HI (Kanaya, S., and Itaya, M. (1992) J. Biol. Chem. 267, 10184-10192). The identity in the amino acid sequences of these enzymes is 52%. As an initial step to elucidate the stabilization mechanism of the thermophilic RNase H, we examined whether certain regions in its amino acid sequence confer the thermostability. A variety of mutant proteins of E. coli RNase HI were constructed and analyzed for protein stability. In these mutant proteins, amino acid sequences in loops or terminal regions were systematically replaced with the corresponding sequences from T. thermophilus RNase H. Of the nine regions examined, replacement of the amino acid sequence in each of four regions (R4-R7) resulted in an increase in protein stability. Simultaneous replacements of these amino acid sequences revealed that the effect of each replacement on protein stability is independent of each other and cumulative. Replacement of all four regions (R4-R7) gave the most stable mutant protein. The temperature of the midpoint of the transition in the thermal unfolding curve and the free energy change of unfolding in the absence of denaturant of this mutant protein were increased by 16.7 degrees C and 3.66 kcal/mol, respectively, as compared to those of E. coli RNase HI. These results suggest that individual local interactions contribute to the stability of thermophilic proteins in an independent manner, rather than in a cooperative manner.     

Acetylation regulates the stability of a bacterial protein: growth stage-dependent modification of RNase R

RNase R, an Escherichia coli exoribonuclease important for degradation of structured RNAs, increases 3-?to 10-fold under certain stress conditions, due to an increased half-life for this usually unstable protein. Components of the trans-translation machinery, tmRNA, and its associated protein, SmpB, are essential for RNase R instability. However, it is not understood why exponential phase RNase R is unstable or how it becomes stabilized in stationary phase. Here, we show that these phenomena are regulated by acetylation catalyzed by YfiQ protein. One residue, Lys544, is acetylated in exponential phase RNase R, but not in the stationary phase protein, resulting in tighter binding of tmRNA-SmpB to the C-terminal region of exponential phase RNase R and subsequent proteolytic degradation. Removal of the positive charge at Lys544 or a negative charge in the C-terminal region likely disrupts their interaction, facilitating tmRNA-SmpB binding. These findings indicate that acetylation can regulate the stability of a bacterial protein.     

Remarkable stabilization of a psychrotrophic RNase HI by a combination of thermostabilizing mutations identified by the suppressor mutation method

Ribonuclease HI from the psychrotrophic bacterium Shewanella oneidensis MR-1 (So-RNase HI) is much less stable than Escherichia coli RNase HI (Ec-RNase HI) by 22.4 degrees C in T m and 12.5 kJ mol (-1) in Delta G(H 2O), despite their high degrees of structural and functional similarity. To examine whether the stability of So-RNase HI increases to a level similar to that of Ec-RNase HI via introduction of several mutations, the mutations that stabilize So-RNase HI were identified by the suppressor mutation method and combined. So-RNase HI and its variant with a C-terminal four-residue truncation (154-RNase HI) complemented the RNase H-dependent temperature-sensitive (ts) growth phenotype of E. coli strain MIC3001, while 153-RNase HI with a five-residue truncation could not. Analyses of the activity and stability of these truncated proteins suggest that 153-RNase HI is nonfunctional in vivo because of a great decrease in stability. Random mutagenesis of 153-RNase HI using error-prone PCR, followed by screening for the revertants, allowed us to identify six single suppressor mutations that make 153-RNase HI functional in vivo. Four of them markedly increased the stability of the wild-type protein by 3.6-6.7 degrees C in T m and 1.7-5.2 kJ mol (-1) in Delta G(H 2O). The effects of these mutations were nearly additive, and combination of these mutations increased protein stability by 18.7 degrees C in T m and 12.2 kJ mol (-1) in Delta G(H 2O). These results suggest that several residues are not optimal for the stability of So-RNase HI, and their replacement with other residues strikingly increases it to a level similar to that of the mesophilic counterpart.     

Complementation analysis of the cold-sensitive phenotype of the Escherichia coli csdA deletion strain

The cold shock response of Escherichia coli is elicited by downshift of temperature from 37 degrees C to 15 degrees C and is characterized by induction of several cold shock proteins, including CsdA, during the acclimation phase. CsdA, a DEAD-box protein, has been proposed to participate in a variety of processes, such as ribosome biogenesis, mRNA decay, translation initiation, and gene regulation. It is not clear which of the functions of CsdA play a role in its essential cold shock function or whether all do, and so far no protein has been shown to complement its function in vivo. Our screening of an E. coli genomic library for an in vivo counterpart of CsdA that can compensate for its absence at low temperature revealed only one protein, RhlE, another DEAD-box RNA helicase. We also observed that although not detected in our genetic screening, two cold shock-inducible proteins, namely, CspA, an RNA chaperone, and RNase R, an exonuclease, can also complement the cold shock function of CsdA. Interestingly, the absence of CsdA and RNase R leads to increased sensitivity of the cells to even moderate temperature downshifts. The correlation between the helicase activity of CsdA and the stability of mRNAs of cold-inducible genes was shown using cspA mRNA, which was significantly stabilized in the DeltacsdA cells, an effect counteracted by overexpression of wild-type CsdA or RNase R but not by that of the helicase-deficient mutant of CsdA. These results suggest that the primary role of CsdA in cold acclimation of cells is in mRNA decay and that its helicase activity is pivotal for promoting degradation of mRNAs stabilized at low temperature.     

Intraspecies regulation of ribonucleolytic activity

The evolutionary rate of proteins involved in obligate protein-protein interactions is slower and the degree of coevolution higher than that for nonobligate protein-protein interactions. The coevolution of the proteins involved in certain nonobligate interactions is, however, essential to cell survival. To gain insight into the coevolution of one such nonobligate protein pair, the cytosolic ribonuclease inhibitor (RI) proteins and secretory pancreatic-type ribonucleases from cow (Bos taurus) and human (Homo sapiens) were produced in Escherichia coli and purified, and their physicochemical properties were analyzed. The two intraspecies complexes were found to be extremely tight (bovine Kd = 0.69 fM; human Kd = 0.34 fM). Human RI binds to its cognate ribonuclease (RNase 1) with 100-fold greater affinity than to the bovine homologue (RNase A). In contrast, bovine RI binds to RNase 1 and RNase A with nearly equal affinity. This broader specificity is consistent with there being more pancreatic-type ribonucleases in cows (20) than humans (13). Human RI (32 cysteine residues) also has 4-fold less resistance to oxidation by hydrogen peroxide than does bovine RI (29 cysteine residues). This decreased oxidative stability of human RI, which is caused largely by Cys74, implies a larger role for human RI as an antioxidant. The conformational and oxidative stabilities of both RIs increase upon complex formation with ribonucleases. Thus, RI has evolved to maintain its inhibition of invading ribonucleases, even when confronted with extreme environmental stress. That role appears to take precedence over its role in mediating oxidative damage.     

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.     

The C-terminal half of RNase E, which organizes the Escherichia coli degradosome, participates in mRNA degradation but not rRNA processing in vivo

RNase E is an essential Escherichia coli endonuclease, which controls both 5S rRNA maturation and bulk mRNA decay. While the C-terminal half of this 1061-residue protein associates with polynucleotide phosphorylase (PNPase) and several other enzymes into a 'degradosome', only the N-terminal half, which carries the catalytic activity, is required for growth. We characterize here a mutation (rne131 ) that yields a metabolically stable polypeptide lacking the last 477 residues of RNAse E. This mutation resembles the N-terminal conditional mutation rne1 in stabilizing mRNAs, both in bulk and individually, but differs from it in leaving rRNA processing and cell growth unaffected. Another mutation (rne105 ) removing the last 469 residues behaves similarly. Thus, the C-terminal half of RNase E is instrumental in degrading mRNAs, but dispensable for processing rRNA. A plausible interpretation is that the former activity requires that RNase E associates with other degradosome proteins; however, PNPase is not essential, as RNase E remains fully active towards mRNAs in rne+pnp mutants. All mRNAs are not stabilized equally by the rne131 mutation: the greater their susceptibility to RNase E, the larger the stabilization. Artificial mRNAs generated by E. coli expression systems based on T7 RNA polymerase can be genuinely unstable, and we show that the mutation can improve the yield of such systems without compromising cell growth.     

RNase P of the Cyanophora paradoxa cyanelle: a plastid ribozyme

Ribonuclease P (RNase P) is a ribonucleoprotein enzyme that generates the mature 5' ends of tRNAs. Ubiquitous across all three kingdoms of life, the composition and functional contributions of the RNA and protein components of RNase P differ between the kingdoms. RNA-alone catalytic activity has been reported throughout bacteria, but only for some archaea, and only as trace activity for eukarya. Available information for RNase P from photosynthetic organelles points to large differences to bacterial as well as to eukaryotic RNase P: for spinach chloroplasts, protein-alone activity has been discussed; for RNase P from the cyanelle of the glaucophyte Cyanophora paradoxa, a type of organelle sharing properties of both cyanobacteria and chloroplasts, the proportion of protein was found to be around 80% rather than the usual 10% in bacteria. Furthermore, the latter RNase P was previously found catalytically inactive in the absence of protein under a variety of conditions; however, the RNA could be activated by a cyanobacterial protein, but not by the bacterial RNase P protein from Escherichia coli. Here we demonstrate that, under very high enzyme concentrations, the RNase P RNA from the cyanelle of C. paradoxa displays RNA-alone activity well above the detection level. Moreover, the RNA can be complemented to a functional holoenzyme by the E. coli RNase P protein, further supporting its overall bacterial-like architecture. Mutational analysis and domain swaps revealed that this A,U-rich cyanelle RNase P RNA is globally optimized but conformationally unstable, since changes as little as a single point mutation or a base pair identity switch at positions that are not part of the universally conserved catalytic core led to a complete loss of RNA-alone activity. Likely related to this low robustness, extensive structural changes towards an E. coli-type P5-7/P15-17 subdomain as a canonical interaction site for tRNA 3'-CCA termini could not be coaxed into increased ribozyme activity.     

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

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

Structural and thermodynamic analyses of Escherichia coli RNase HI variant with quintuple thermostabilizing mutations

A combination of five thermostabilizing mutations, Gly23-->Ala, His62-->Pro, Val74-->Leu, Lys95-->Gly, and Asp134-->His, has been shown to additively enhance the thermostability of Escherichia coli RNase HI [Akasako A, Haruki M, Oobatake M & Kanaya S (1995) Biochemistry34, 8115-8122]. In this study, we determined the crystal structure of the protein with these mutations (5H-RNase HI) to analyze the effects of the mutations on the structure in detail. The structures of the mutation sites were almost identical to those of the mutant proteins to which the mutations were individually introduced, except for G23A, for which the structure of the single mutant protein is not available. Moreover, only slight changes in the backbone conformation of the protein were observed, and the interactions of the side chains were almost conserved. These results indicate that these mutations almost independently affect the protein structure, and are consistent with the fact that the thermostabiling effects of the mutations are cumulative. We also determined the protein stability curve describing the temperature dependence of the free energy of unfolding of 5H-RNase HI to elucidate the thermostabilization mechanism. The maximal stability for 5H-RNase HI was as high as that for the cysteine-free variant of Thermus thermophilus RNase HI. In contrast, the heat capacity of unfolding for 5H-RNase H was similar to that for E. coli RNase HI, which is considerably higher than that for T. thermophilus RNase HI. These results suggest that 5H-RNase HI is stabilized, in part, by the thermostabilization mechanism adopted by T. thermophilus RNase HI.     

Divalent metal cofactor binding in the kinetic folding trajectory of Escherichia coli ribonuclease HI

Proteins often require cofactors to perform their biological functions and must fold in the presence of their cognate ligands. Using circular dichroism spectroscopy. we investigated the effects of divalent metal binding upon the folding pathway of Escherichia coli RNase HI. This enzyme binds divalent metal in its active site, which is proximal to the folding core of RNase HI as defined by hydrogen/deuterium exchange studies. Metal binding increases the apparent stability of native RNase HI chiefly by reducing the unfolding rate. As with the apo-form of the protein, refolding from high denaturant concentrations in the presence of Mg2+ follows three-state kinetics: formation of a rapid burst phase followed by measurable single exponential kinetics. Therefore, the overall folding pathway of RNase HI is minimally perturbed by the presence of metal ions. Our results indicate that the metal cofactor enters the active site pocket only after the enzyme reaches its native fold, and therefore, divalent metal binding stabilizes the protein by decreasing its unfolding rate. Furthermore, the binding of the cofactor is dependent upon a carboxylate critical for activity (Asp10). A mutation in this residue (D10A) alters the folding kinetics in the absence of metal ions such that they are similar to those observed for the unaltered enzyme in the presence of metal.     

The rate of decay of Rhodobacter capsulatus-specific puf mRNA segments is differentially affected by RNase E activity in Escherichia coli.

In Rhodobacter capsulatus the puf operon encodes proteins of the photosynthetic apparatus. The polycistronic puf mRNA is comprised of segments that show differential stability. Here, we show that the rate of decay of the 2.7-kb pufBALMX mRNA species in Escherichia coli depends on the activity of ribonuclease E (RNase E), whereas the degradation of the 0.5-kb pufBA mRNA segment is not affected by a mutation in the rne gene. The RNase E-promoted decay of the pufLMX mRNA depends on the presence of a 1.4-kb pufLM mRNA segment, in which rate-limiting endonucleolytic cleavage was postulated to occur in R. capsulatus. The insertion of 185 bp of this 1.4-kb segment into pufB results in an RNase E-dependent decay of the modified pufBA mRNA segment in E. coli. Our findings suggest that in R. capsulatus an RNase E-like activity is responsible for the rate-limiting endonucleolytic cleavage occurring within the pufLM mRNA segment, whereas the 0.5-kb pufBA mRNA segment is degraded by a different RNase E-independent decay mechanism.     

Protein activation of a ribozyme: the role of bacterial RNase P protein

Bacterial ribonuclease P (RNase P) belongs to a class of enzymes that utilize both RNAs and proteins to perform essential cellular functions. The bacterial RNase P protein is required to activate bacterial RNase P RNA in vivo, but previous studies have yielded contradictory conclusions regarding its specific functions. Here, we use biochemical and biophysical techniques to examine all of the proposed functions of the protein in both Escherichia coli and Bacillus subtilis RNase P. We demonstrate that the E. coli protein, but not the B. subtilis protein, stabilizes the global structure of RNase P RNA, although both proteins influence holoenzyme dimer formation and precursor tRNA recognition to different extents. By comparing each protein in complex with its cognate and noncognate RNA, we show that differences between the two types of holoenzymes reside primarily in the RNA and not the protein components of each. Our results reconcile previous contradictory conclusions regarding the role of the protein and support a model where the protein activates local RNA structures that manifest multiple holoenzyme properties.     

RNase ES of Streptomyces coelicolor A3(2) can complement the rne and rng mutations in Escherichia coli

The Streptomyces coelicolor gene SCC88.10c encodes a protein (RNase ES) which is homologous to endoribonucleases in the RNase E/G family. We expressed S. coelicolor RNase ES as a 6 x His-tagged protein in an Escherichia coli mutant carrying a rng (which encodes RNase G) or a rne (which encodes RNase E) mutation to study whether S. coelicolor RNase ES is able to complement these mutations in host E. coli cells. The results clearly indicated that the S. coelicolor RNase ES can partially abrogate either the rng::cat or rne-1 mutation, as measured by the ability to suppress the several aberrant phenotypes resulting from the rng or rne mutation. Thus, S. coelicolor RNase ES appears to have the dual ability to supplant the functions of both RNase G and RNase E in E. coli.     

Heat labile ribonuclease HI from a psychrotrophic bacterium: gene cloning, characterization and site-directed mutagenesis.

The rnhA gene encoding RNase HI from a psychrotrophic bacterium, Shewanella sp. SIB1, was cloned, sequenced and overexpressed in an rnh mutant strain of Escherichia coli. SIB1 RNase HI is composed of 157 amino acid residues and shows 63% amino acid sequence identity to E.coli RNase HI. Upon induction, the recombinant protein accumulated in the cells in an insoluble form. This protein was solubilized and purified in the presence of 7 M urea and refolded by removing urea. Determination of the enzymatic activity using M13 DNA-RNA hybrid as a substrate revealed that the enzymatic properties of SIB1 RNase HI, such as divalent cation requirement, pH optimum and cleavage mode of a substrate, are similar to those of E.coli RNase HI. However, SIB1 RNase HI was much less stable than E.coli RNase HI and the temperature (T(1/2)) at which the enzyme loses half of its activity upon incubation for 10 min was approximately 25 degrees C for SIB1 RNase HI and approximately 60 degrees C for E.coli RNase HI. The optimum temperature for the SIB1 RNase HI activity was also shifted downward by 20 degrees C compared with that of E.coli RNase HI. Nevertheless, SIB1 RNase HI was less active than E.coli RNase HI even at low temperatures. The specific activity determined at 10 degrees C was 0.29 units/mg for SIB1 RNase HI and 1.3 units/mg for E.coli RNase HI. Site-directed mutagenesis studies suggest that the amino acid substitution in the middle of the alphaI-helix (Pro52 for SIB1 RNase HI and Ala52 for E.coli RNase HI) partly accounts for the difference in the stability and activity between SIB1 and E.coli RNases HI.