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

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

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.     

Experiment-guided thermodynamic simulations on reversible two-state proteins: implications for protein thermostability

Here, we perform protein thermodynamic simulations within a set of boundary conditions, effectively blanketing the experimental data. The thermodynamic parameters, melting temperature (TG), enthalpy change at the melting temperature (DeltaHG) and heat capacity change (DeltaCp) were systematically varied over the experimentally observed ranges for small single domain reversible two-state proteins. Parameter sets that satisfy the Gibbs-Helmholtz equation and yield a temperature of maximal stability (TS) around room temperature were selected. The results were divided into three categories by arbitrarily chosen TG ranges. The TG ranges in these categories correspond to typical values of the melting temperatures observed for the majority of the proteins from mesophilic, thermophilic and hyperthermophilic organisms. As expected, DeltaCp values tend to be high in mesophiles and low in hyperthermophiles. An increase in TG is accompanied by an up-shift and broadening of the protein stability curves, however, with a large scatter. Furthermore, the simulations reveal that the average DeltaHG increases with TG up to approximately 360 K and becomes constant thereafter. DeltaCp decreases with TG with different rates before and after approximately 360 K. This provides further justification for the separate grouping of proteins into thermophiles and hyperthermophiles to assess their thermodynamic differences. This analysis of the Gibbs-Helmholtz equation has allowed us to study the interdependence of the thermodynamic parameters TG, DeltaHG and DeltaCp and their derivatives in a more rigorous way than possible by the limited experimental protein thermodynamics data available in the literature. The results provide new insights into protein thermostability and suggest potential strategies for its manipulation.     

A thermodynamic study of mesophilic, thermophilic, and hyperthermophilic L-arabinose isomerases: the effects of divalent metal ions on protein stability at elevated temperatures

To gain insight into the structural stability of homologous homo-tetrameric l-arabinose isomerases (AI), we have examined the isothermal guanidine hydrochloride (GdnHCl)-induced unfolding of AIs from mesophilic Bacillus halodurans (BHAI), thermophilic Geobacillus stearothermophilus (GSAI), and hyperthermophilic Thermotoga maritima (TMAI) using circular dichroism spectroscopy. The GdnHCl-induced unfolding of the AIs can be well described by a two-state reaction between native tetramers and unfolded monomers, which directly confirms the validity of the linear extrapolation method to obtain the intrinsic stabilities of these proteins. The resulting unfolding free energy (DeltaGU) values of the AIs as a function of temperature were fit to the Gibbs-Helmholtz equation to determine their thermodynamic parameters based on a two-state mechanism. Compared with the stability curves of BHAI in the presence and absence of Mn2+, those of holo GSAI and TMAI were more broadened than those of the apo enzymes at all temperatures, indicating increased melting temperatures (Tm) due to decreased heat capacity (DeltaGp). Moreover, the extent of difference in DeltaCp between the apo and holo thermophilic AIs is larger than that of BHAI. From these studies, we suggest that the metal dependence of the thermophilic AIs, resulting in the reduced DeltaCp, may play a significant role in structural stability compared to their mesophilic analogues, and that the extent of metal dependence of AI stability seems to be highly correlated to oligomerization.     

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

Mutations derived from the thermophilic polyhydroxyalkanoate synthase PhaC enhance the thermostability and activity of PhaC from Cupriavidus necator H16

The thermophile Cupriavidus sp. strain S-6 accumulated polyhydroxybutyrate (PHB) from glucose at 50°C. A 9.0-kbp EcoRI fragment cloned from the genomic DNA of Cupriavidus sp. S-6 enabled Escherichia coli XL1-Blue to synthesize PHB at 45°C. Nucleotide sequence analysis showed a pha locus in the clone. The thermophilic polyhydroxyalkanoate (PHA) synthase (PhaC(Csp)) shared 81% identity with mesophilic PhaC of Cupriavidus necator H16. The diversity between these two strains was found dominantly on their N and C termini, while the middle regions were highly homologous (92% identity). We constructed four chimeras of mesophilic and thermophilic phaC genes to explore the mutations related to its thermostability. Among the chimeras, only PhaC(H16β), which was PhaC(H16) bearing 30 point mutations derived from the middle region of PhaC(Csp), accumulated a high content of PHB (65% [dry weight]) at 45°C. The chimera phaC(H16)(β) and two parental PHA synthase genes were overexpressed in E. coli BLR(DE3) cells and purified. At 30°C, the specific activity of the chimera PhaC(H16β) (172 ± 17.8 U/mg) was 3.45-fold higher than that of the parental enzyme PhaC(H16) (50 ± 5.2 U/mg). At 45°C, the half-life of the chimera PhaC(H16β) (11.2 h) was 127-fold longer than that of PhaC(H16) (5.3 min). Furthermore, the chimera PhaC(H16β) accumulated 1.55-fold (59% [dry weight]) more PHA content than the parental enzyme PhaC(H16) (38% [dry weight]) at 37°C. This study reveals a limited number of point mutations which enhance not only thermostability but also PhaC(H16) activity. The highly thermostable and active PHA synthase will provide advantages for its promising applications to in vitro PHA synthesis and recombinant E. coli PHA fermentation.     

Studies on structural stability of thermophilic Sulfolobus acidocaldarius ribosomes

Structural stability of thermophilic archaeon Sulfolobus acidocaldarius ribosomes, with respect their susceptibility to pancreatic RNase A and stability to temperature (deltaTm), on treatment with various stabilizing (polyamines) and destabilizing (sulfhydryl and intercalating) agents were studied and compared with mesophilic E. coli ribosomes, to understand the structural differences between thermophilic and mesophilic ribosomes. Thermophilic archaeal ribosomes and their subunits were 10-times less susceptible to pancreatic RNase A, compared to mesophilic ribosomes, showing the presence of strong and compact structural organization in them. Thermophilic ribosomes treated with destabilizing agents, such as sulfhydryl reagents [5,5'-Dithio-bis-(2-nitrobenzoic acid), N-ethylmaleimide and p-hydroxymercurybenzoate) and intercalating agents (ethidium bromide, EtBr) showed higher stability to RNase A, compared to similarly treated mesophilic ribosomes, indicating the unavailability of thiol-reactive groups and the presence of strong solvent inaccessible inner core. Higher stability of thermophilic ribosomes compared to mesophilic ribosomes to unfolding agents like urea further supported the presence of strong inner core particle. Thermophilic ribosomes treated with intercalating agents, such as EtBr were less susceptible to RNase A, though they bound to more reagent, showing the rigidity or resilience of their macromolecular structure to alterations caused by destabilizing agents. Overall, these results indicated that factors such as presence of strong solvent inaccessible inner core and rigidity of ribosome macromolecular structure contributed stability of thermophilic ribosomes to RNase A and other destabilizing agents, when compared to mesophilic ribosomes.     

Different packing of external residues can explain differences in the thermostability of proteins from thermophilic and mesophilic organisms

MOTIVATION: Understanding the basis of protein stability in thermophilic organisms raises a general question: what structural properties of proteins are responsible for the higher thermostability of proteins from thermophilic organisms compared to proteins from mesophilic organisms? RESULTS: A unique database of 373 structurally well-aligned protein pairs from thermophilic and mesophilic organisms is constructed. Comparison of proteins from thermophilic and mesophilic organisms has shown that the external, water-accessible residues of the first group are more closely packed than those of the second. Packing of interior parts of proteins (residues inaccessible to water molecules) is the same in both cases. The analysis of amino acid composition of external residues of proteins from thermophilic organisms revealed an increased fraction of such amino acids as Lys, Arg and Glu, and a decreased fraction of Ala, Asp, Asn, Gln, Thr, Ser and His. Our theoretical investigation of folding/unfolding behavior confirms the experimental observations that the interactions that differ in thermophilic and mesophilic proteins form only after the passing of the transition state during folding. Thus, different packing of external residues can explain differences in thermostability of proteins from thermophilic and mesophilic organisms. AVAILABILITY: The database of 373 structurally well-aligned protein pairs is available at http://phys.protres.ru/resources/termo_meso_base.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

The mechanism of folding of pancreatic ribonucleases is independent of the presence of covalently linked carbohydrate

The role of asparagine-linked oligosaccharides for the mechanism of protein folding was investigated. We compared the stability and folding kinetics for two sets of pancreatic ribonucleases (RNases) with identical amino acid sequences and differences in glycosylation. First the folding of RNases A (carbohydrate free) and B (a single N-linked oligosaccharide) from bovine pancreas was investigated. The kinetics of refolding were identical under a wide range of conditions. The rate of unfolding by guanidinium chloride was decreased in RNase B. In further experiments the folding of porcine RNase (three carbohydrate chains at Asn-21, -34, and -76) was compared with the corresponding data for the deglycosylated protein. Even for this RNase with almost 40% carbohydrate content the mechanism of refolding is independent of glycosylation. Although the folding mechanism is conserved, the rates of individual steps in folding are decreased about 2-fold upon deglycosylation. We interpret this to originate from a slight destabilization of folding intermediates by carbohydrate depletion. In control experiments with nonglycosylated bovine RNase A it was ascertained that treatment with HF (as used for deglycosylation) did not affect the folding kinetics. The in vitro folding mechanism of glycosylated RNases apparently does not depend on the presence of N-linked oligosaccharide chains. The information for the folding of glycoproteins is contained exclusively in the protein moiety, i.e. in the amino acid sequence. Carbohydrate chains are attached at chain positions which remain solvent exposed. This ensures that the presence of oligosaccharides does not interfere with correct folding of the polypeptide chain.     

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.     

Influence of key residues on the heterologous extracellular production of fungal ribonuclease U2 in the yeast Pichia pastoris

Ribonuclease U2, secreted by the smut fungus Ustilago sphaerogena, is a cyclizing ribonuclease that displays a rather unusual specificity within the group of microbial extracellular RNases, best represented by RNase T1. Superposition of the three-dimensional structures of RNases T1 and U2 suggests that the RNase U2 His 101 would be the residue equivalent to the RNase T1 catalytically essential His 92. RNase U2 contains three disulfide bridges but only two of them are conserved among the family of fungal extracellular RNases. The non-conserved disulfide bond is established between Cys residues 1 and 54. Mispairing of the disulfide network due to the presence of two consecutive Cys residues (54 and 55) has been invoked to explain the presence of wrongly folded RNase U2 species when produced in Pichia pastoris. In order to study both hypotheses, the RNase U2 H101Q and C1/54S variants have been produced, purified, and characterized. The results obtained support the major conclusion that His 101 is required for proper protein folding when secreted by the yeast P. pastoris. On the other hand, substitution of the first Cys residue for Ser results in a mutant version which is more efficiently processed in terms of a more complete removal of the yeast α-factor signal peptide. In addition, it has been shown that elimination of the Cys 1–Cys 54 disulfide bridge does not interfere with RNase U2 proper folding, generating a natively folded but much less stable protein.     

Stepwise conversion of a mesophilic to a thermophilic ribozyme

A fundamental question in RNA folding is the mechanism of thermodynamic stability. We investigated the equilibrium folding of a series of sequence variants in which one to three motifs of a 255-nucleotide mesophilic ribozyme were substituted with the corresponding motifs from its thermophilic homologue. Substitution of three crucial motifs individually or in groups results in a continual increase in the stability and folding cooperativity in a stepwise fashion. We find an unexpected relationship between stability and folding cooperativity. Without changing the folding cooperativity, RNAs having a similar native structure can only achieve moderate change in stability and likewise, without changing stability, RNAs having a similar native structure can only achieve moderate change in folding cooperativity. This intricate relationship must be included in the predictions of tertiary RNA stability.     

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.     

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

Ribonuclease (RNase) HI from the psychrotrophic bacterium Shewanella oneidensis MR-1 was overproduced in Escherichia coli, purified, and structurally and biochemically characterized. The amino acid sequence of MR-1 RNase HI is 67% identical to that of E. coli RNase HI. The crystal structure of MR-1 RNase HI determined at 2.0 A resolution was highly similar to that of E. coli RNase HI, except that the number of intramolecular ion pairs and the fraction of polar surface area of MR-1 RNase HI were reduced compared to those of E. coli RNase HI. The enzymatic properties of MR-1 RNase HI were similar to those of E. coli RNase HI. However, MR-1 RNase HI was much less stable than E. coli RNase HI. The stability of MR-1 RNase HI against heat inactivation was lower than that of E. coli RNase HI by 19 degrees C. The conformational stability of MR-1 RNase HI was thermodynamically analyzed by monitoring the CD values at 220 nm. MR-1 RNase HI was less stable than E. coli RNase HI by 22.4 degrees C in Tm and 12.5 kJ/mol in DeltaG(H2O). The thermodynamic stability curve of MR-1 RNase HI was characterized by a downward shift and increased curvature, which results in an increased DeltaCp value, compared to that of E. coli RNase HI. Site-directed mutagenesis studies suggest that the difference in the number of intramolecular ion pairs partly accounts for the difference in stability between MR-1 and E. coli RNases HI.     

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.     

Confirmation of the hierarchical folding of RNase H: a protein engineering study.

The kinetic intermediate of RNase H is structured in a core region of the protein. To probe the role of this intermediate in the folding of RNase H, the folding kinetics of mutant proteins with altered native state stabilities were investigated. Mutations within the folding core destabilize the kinetic intermediate and slow refolding in a manner consistent with an obligatory intermediate model. Mutations outside of the folding core, however, do not affect the stability of the kinetic intermediate but do perturb the native state and transition state. These results indicate that interactions formed in the intermediate persist in the transition and native states and that RNase H folds through a hierarchical mechanism.     

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.     

Activity-stability relationships in extremophilic enzymes

Psychrophilic, mesophilic, and thermophilic alpha-amylases have been studied as regards their conformational stability, heat inactivation, irreversible unfolding, activation parameters of the reaction, properties of the enzyme in complex with a transition state analog, and structural permeability. These data allowed us to propose an energy landscape for a family of extremophilic enzymes based on the folding funnel model, integrating the main differences in conformational energy, cooperativity of protein unfolding, and temperature dependence of the activity. In particular, the shape of the funnel bottom, which depicts the stability of the native state ensemble, also accounts for the thermodynamic parameters of activation that characterize these extremophilic enzymes, therefore providing a rational basis for stability-activity relationships in protein adaptation to extreme temperatures.     

Structural basis for altering the stability of homologous RNAs from a mesophilic and a thermophilic bacterium

Tertiary RNA structures from thermophilic bacteria generally are more stable than their mesophilic homologs. To understand the structural basis of the increase in stability, we investigated equilibrium folding of the specificity domain (S-domain) of RNase P RNA from a mesophilic (Escherichia coli) and a thermophilic (Thermus thermophilus) bacterium. Equilibrium folding of both S-domains is described by a minimal, three-state folding scheme, U-to-I-to-N. In the I-to-N transition of the thermophilic S-domain, more structure forms and protections are stronger against T1 nuclease and hydroxyl radical reactions. Phylogenetic comparison in the context of the native structure reveals that among 39 nucleotide differences between these S-domains, 12 likely contribute to higher stability. These residues participate in extensive networks of hydrogen bonding, stacking, and metal ion coordination throughout the molecule. The thermophilic S-domain achieves higher stability by mutating strategic base pairs to G-C, decreasing surface accessibility of the native state, and increasing the amount of structure formation in the native folding transition. An E. coli S-domain mutant containing these 12 nt has the same stability and folding cooperativity as the T. thermophilus S-domain. E. coli S-domain mutants containing a subset of 4 or 6 nt have the same stability as the T. thermophilus S-domain but the same folding cooperativity as the E. coli S-domain. These results show that increasing stability can be accomplished by mutations within a local structure, but increasing folding cooperativity needs concerted changes among multiple structural units.