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.     

Exploiting the Link between Protein Rigidity and Thermostability for Data‐Driven Protein Engineering

Understanding and exploiting the relationship between microscopic structure and macroscopic stability is important for developing strategies to improve protein stability at high temperatures. The thermostability of proteins has been repeatedly linked to an enhanced structural rigidity of the folded native state. In the current study, the rigidity of protein structures from mesophilic and thermophilic organisms along a thermal unfolding trajectory is directly probed. In order to perform this, protein structures were modeled as constraint networks, and the rigidity in these networks was quantified using the Floppy Inclusion and Rigid Substructure Topography (FIRST) method. During the thermal unfolding, a phase transition was observed that defines the rigidity percolation threshold and corresponds to the folded‐unfolded transition in protein folding. Using concepts from percolation theory and network science, a higher phase transition temperature was observed for ca. two‐thirds of the proteins from thermophilic organisms compared to their mesophilic counterparts, when applied to a data set of 20 pairs of homologues. From both the analysis of the microstructure of the constraint networks and monitoring the macroscopic behavior during the thermal unfolding, direct evidence was found for the “corresponding states” concept, which states that mesophilic and thermophilic enzymes are in corresponding states of similar flexibility at their respective optimal temperature. Finally, the current approach facilitated the identification of structural features from which a destabilization of the structure originates upon thermal unfolding. These predictions show a good agreement with the experimental data. Therefore, the information might be exploited in data‐driven protein engineering by pointing to residues that should be varied to obtain a protein with higher thermostability.     

Electrostatic contributions to the stability of a thermophilic cold shock protein

The thermophilic Bacillus caldolyticus cold shock protein (Bc-Csp) differs from the mesophilic Bacillus subtilis cold shock protein B (Bs-CspB) in 11 of the 66 residues. Stability measurements of Schmid and co-workers have implicated contributions of electrostatic interactions to the thermostability. To further elucidate the physical basis of the difference in stability, previously developed theoretical methods that treat electrostatic effects in both the folded and the unfolded states were used in this paper to study the effects of mutations, ionic strength, and temperature. For 27 mutations that narrow the difference in sequence between Bc-Csp and Bs-CspB, calculated changes in unfolding free energy (Delta G) and experimental results have a correlation coefficient of 0.98. Bc-Csp appears to use destabilization of the unfolded state by unfavorable charge-charge interactions as a mechanism for increasing stability. Accounting for the effects of ionic strength and temperature on the electrostatic free energies in both the folded and the unfolded states, explanations for two important experimental observations are presented. The disparate ionic strength dependences of Delta G for Bc-Csp and Bs-CspB were attributed to the difference in the total charges (-2e and -6e, respectively). A main contribution to the much higher unfolding entropy of Bs-CspB was found to come from the less favorable electrostatic interactions in the folded state. These results should provide insight for understanding the thermostability of other thermophilic proteins.     

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.     

Aromatic clusters: a determinant of thermal stability of thermophilic proteins

A number of factors have been elucidated as responsible for the thermal stability of thermophilic proteins. However, the contribution of aromatic interactions to thermal stability has not been systematically studied. In the present investigation we used a graph spectral method to identify aromatic clusters in a dataset of 24 protein families for which the crystal structures of both the thermophilic and their mesophilic homologues are known. Our analysis shows a presence of additional aromatic clusters or enlarged aromatic networks in 17 different thermophilic protein families, which are absent in the corresponding mesophilic homologue. The additional aromatic clusters identified in the thermophiles are smaller in size and are largely found on the protein surface. The aromatic clusters are found to be relatively rigid regions of the surface and often the additional aromatic cluster is located close to the active site of the thermophilic enzyme. The residues in the additional aromatic clusters are preferably mutated to Leu, Ser or Ile in the mesophilic homologue. An analysis of the packing geometry of the pairwise aromatic interaction in the additional aromatic clusters shows a preference for a T-shaped orthogonal packing geometry. The present study also provides new insights for protein engineers to design thermostable and thermophilic proteins.     

Dynamic arrangement of ion pairs and individual contributions to the thermal stability of the cofactor-binding domain of glutamate dehydrogenase from Thermotoga maritima

The dynamics of a hyperthermophilic protein fragment in a water environment, as studied by performing molecular dynamics (MD) simulations at various temperatures, is compared to the dynamical behavior of a homologous mesophilic protein simulated under identical conditions. The effects on the stability of the spatial arrangement and mobility of the charged residues in solution were quantified by calculating free energy changes upon salt bridge formation in these proteins. Electrostatic free energy terms derived from a thermodynamic cycle were obtained by solving the linearized Poisson-Boltzmann equation for a series of protein conformations generated by MD simulations and placed subsequently in a continuum solvent medium. Our results show that the ion pairs are electrostatically stabilizing in most of the cases, but their individual contributions vary significantly. The greater contribution of the charged residues to the stability of the hyperthermophilic protein as compared with the mesophilic counterpart was evidenced only by the calculations that included conformations sampled at 343 and 373 K. The "dynamic" structure of the hyperthermophilic protein fragment simulated at elevated temperatures reveals an optimum placement of the ionizable residues within the protein structure as well as the role of their cooperative interactions in promoting thermal stability. The thermodynamic properties such as electrostatic free energy differences, configurational entropies, and specific heat capacities calculated in the dynamic context of the protein structure provided new insight into the mechanism of protein thermostabilization.     

Toward the physical basis of thermophilic proteins: linking of enriched polar interactions and reduced heat capacity of unfolding

The enrichment of salt bridges and hydrogen bonding in thermophilic proteins has long been recognized. Another tendency, featuring lower heat capacity of unfolding (DeltaC(p)) than found in mesophilic proteins, is emerging from the recent literature. Here we present a simple electrostatic model to illustrate that formation of a salt-bridge or hydrogen-bonding network around an ionized group in the folded state leads to increased folding stability and decreased DeltaC(p). We thus suggest that the reduced DeltaC(p) of thermophilic proteins could partly be attributed to enriched polar interactions. A reduced DeltaC(p) might serve as an indicator for the contribution of polar interactions to folding stability.     

Hydrophobic environment is a key factor for the stability of thermophilic proteins

The stability of thermophilic proteins has been viewed from different perspectives and there is yet no unified principle to understand this stability. It would be valuable to reveal the most important interactions for designing thermostable proteins for such applications as industrial protein engineering. In this work, we have systematically analyzed the importance of various interactions by computing different parameters such as surrounding hydrophobicity, inter‐residue interactions, ion‐pairs and hydrogen bonds. The importance of each interaction has been determined by its predicted relative contribution in thermophiles versus the same contribution in mesophilic homologues based on a dataset of 373 protein families. We predict that hydrophobic environment is the major factor for the stability of thermophilic proteins and found that 80% of thermophilic proteins analyzed showed higher hydrophobicity than their mesophilic counterparts. Ion pairs, hydrogen bonds, and interaction energy are also important and favored in 68%, 50%, and 62% of thermophilic proteins, respectively. Interestingly, thermophilic proteins with decreased hydrophobic environments display a greater number of hydrogen bonds and/or ion pairs. The systematic elimination of mesophilic proteins based on surrounding hydrophobicity, interaction energy, and ion pairs/hydrogen bonds, led to correctly identifying 95% of the thermophilic proteins in our analyses. Our analysis was also applied to another, more refined set of 102 thermophilic–mesophilic pairs, which again identified hydrophobicity as a dominant property in 71% of the thermophilic proteins. Further, the notion of surrounding hydrophobicity, which characterizes the hydrophobic behavior of residues in a protein environment, has been applied to the three‐dimensional structures of elongation factor‐Tu proteins and we found that the thermophilic proteins are enriched with a hydrophobic environment. The results obtained in this work highlight the importance of hydrophobicity as the dominating characteristic in the stability of thermophilic proteins, and we anticipate this will be useful in our attempts to engineering thermostable proteins. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Interface Matters: The Stiffness Route to Stability of a Thermophilic Tetrameric Malate Dehydrogenase

In this work we investigate by computational means the behavior of two orthologous bacterial proteins, a mesophilic and a thermophilic tetrameric malate dehydrogenase (MalDH), at different temperatures. Namely, we quantify how protein mechanical rigidity at different length- and time-scales correlates to protein thermophilicity as commonly believed. In particular by using a clustering analysis strategy to explore the conformational space of the folded proteins, we show that at ambient conditions and at the molecular length-scale the thermophilic variant is indeed more rigid that the mesophilic one. This rigidification is the result of more efficient inter-domain interactions, the strength of which is further quantified via ad hoc free energy calculations. When considered isolated, the thermophilic domain is indeed more flexible than the respective mesophilic one. Upon oligomerization, the induced stiffening of the thermophilic protein propagates from the interface to the active site where the loop, controlling the access to the catalytic pocket, anchors down via an extended network of ion-pairs. On the contrary in the mesophilic tetramer the loop is highly mobile. Simulations at high temperature, could not re-activate the mobility of the loop in the thermophile. This finding opens questions on the similarities of the binding processes for these two homologues at their optimal working temperature and suggests for the thermophilic variant a possible cooperative role of cofactor/substrate.     

Two exposed amino acid residues confer thermostability on a cold shock protein

Thermophilic organisms produce proteins of exceptional stability. To understand protein thermostability at the molecular level we studied a pair of cold shock proteins, one of mesophilic and one of thermophilic origin, by systematic mutagenesis. Although the two proteins differ in sequence at 12 positions, two surface-exposed residues are responsible for the increase in stability of the thermophilic protein (by 15.8 kJ mol-1 at 70 degrees C). 11.5 kJ mol-1 originate from a predominantly electrostatic contribution of Arg 3 and 5.2 kJ mol-1 from hydrophobic interactions of Leu 66 at the carboxy terminus. The mesophilic protein could be converted to a highly thermostable form by changing the Glu residues at positions 3 and 66 to Arg and Leu, respectively. The variation of surface residues may thus provide a simple and powerful approach for increasing the thermostability of a protein.     

A comparative molecular dynamics study of thermophilic and mesophilic ribonuclease HI enzymes

We studied a pair of homologous thermophilic and mesophilic ribonuclease HI enzymes by molecular dynamics simulations. Each protein was subjected to three 5 ns simulations in explicit water at both 310 K and 340 K. The thermophilic enzyme showed larger overall positional fluctuations at both temperatures, while only the mesophilic enzyme at the higher temperature showed significant instability. When the temperature is changed, the relative flexibility of different local segments on the two proteins changed differently. Principal component analysis showed that the simulations of the two proteins explored largely overlapping regions in the conformational space. However, at 340 K, the collective structure variations of the thermophilic protein are different from those of the mesophilic protein. Our results, although not in accordance with the view that hyperthermostability of proteins may originate from their conformational rigidity, are consistent with several recent experimental and simulation studies which showed that thermophilic proteins may be conformationally more flexible than their mesophilic counterparts. The decorrelation between conformational rigidity and hyperthermostability may be attributed to the temperature dependence and long range nature of electrostatic interactions that play more important roles in the structural stability of thermophilic proteins.     

Differences in amino acids composition and coupling patterns between mesophilic and thermophilic proteins

Summary. Thermophilic proteins show substantially higher intrinsic thermal stability than their mesophilic counterparts. Amino acid composition is believed to alter the intrinsic stability of proteins. Several investigations and mutagenesis experiment have been carried out to understand the amino acid composition for the thermostability of proteins. This review presents some generalized features of amino acid composition found in thermophilic proteins, including an increase in residue hydrophobicity, a decrease in uncharged polar residues, an increase in charged residues, an increase in aromatic residues, certain amino acid coupling patterns and amino acid preferences for thermophilic proteins. The differences of amino acids composition between thermophilic and mesophilic proteins are related to some properties of amino acids. These features provide guidelines for engineering mesophilic protein to thermophilic protein. Authors’ addresses: Yuan-Jiang Pan, Institute of Chemical Biology and Pharmaceutical Chemistry, Zhejiang University, Zhejiang University Road 38, Hangzhou 310027, China; Wei-Fen Li, Microbiology Division, College of Animal Science, Zhejiang University, Hangzhou 310029, China     

Solution structure and thermal stability of ribosomal protein L30e from hyperthermophilic archaeon Thermococcus celer

To understand the structural basis of thermostability, we have determined the solution structure of a thermophilic ribosomal protein L30e from Thermococcus celer by NMR spectroscopy. The conformational stability of T. celer L30e was measured by guanidine and thermal-induced denaturation, and compared with that obtained for yeast L30e, a mesophilic homolog. The melting temperature of T. celer L30e was 94 degrees C, whereas the yeast protein denatured irreversibly at temperatures >45 degrees C. The two homologous proteins also differ greatly in their stability at 25 degrees C: the free energy of unfolding was 45 kJ/mole for T. celer L30e and 14 kJ/mole for the yeast homolog. The solution structure of T. celer L30e was compared with that of the yeast homolog. Although the two homologous proteins do not differ significantly in their number of hydrogen bonds and the amount of solvent accessible surface area buried with folding, the thermophilic T. celer L30e was found to have more long-range ion pairs, more proline residues in loops, and better helix capping residues in helix-1 and helix-4. A K9A variant of T. celer L30e was created by site-directed mutagenesis to examine the role of electrostatic interactions on protein stability. Although the melting temperatures of the K9A variant is approximately 8 degrees C lower than that of the wild-type L30e, their difference in T(m) is narrowed to approximately 4.2 degrees C at 0.5 M NaCl. This salt-dependency of melting temperatures strongly suggests that electrostatic interactions contribute to the thermostability of T. celer L30e.     

A Comparative Molecular Dynamics Study of Thermophilic and Mesophilic Ribonuclease HI Enzymes

Abstract

We studied a pair of homologous thermophilic and mesophilic ribonuclease HI enzymes by molecular dynamics simulations. Each protein was subjected to three 5 ns simulations in explicit water at both 310 K and 340 K. The thermophilic enzyme showed larger overall positional fluctuations at both temperatures, while only the mesophilic enzyme at the higher temperature showed significant instability. When the temperature is changed, the relative flexibility of different local segments on the two proteins changed differently. Principal component analysis showed that the simulations of the two proteins explored largely overlapping regions in the conformational space. However, at 340 K, the collective structure variations of the thermophilic protein are different from those of the mesophilic protein. Our results, although not in accordance with the view that hyperthermostability of proteins may originate from their conformational rigidity, are consistent with several recent experimental and simulation studies which showed that thermophilic proteins may be conformationally more flexible than their mesophilic counterparts. The decorrelation between conformational rigidity and hyperthermostability may be attributed to the temperature dependence and long range nature of electrostatic interactions that play more important roles in the structural stability of thermophilic proteins.     

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.     

Effective factors in thermostability of thermophilic proteins

Thermostability of proteins in general and especially thermophilic proteins has been subject of a wide variety of studies based on theoretical and experimental investigation. Thermostability seems to be a property obtained through many minor structural modifications rather than certain amino acids substitution. In comparison with its mesophile homologue in a thermostable protein, usually a number of amino acids are exchanged. A wide variety of theoretical studies are based on comparative investigation of thermophilic proteins characteristics with their mesophilic counterparts in order to reveal their sequences, structural differences and consequently, to relate these observed differences to the thermostability properties. In this work we have compared a dataset of thermophilic proteins with their mesophilic homologues and furthermore, a mesophilic proteins dataset was also compared with its mesophilic homologue. This strategy enabled us first, to eliminate noise or background differences from signals and moreover, the important factors which were related to the thermostability were recognized too. Our results reveal that thermophilic and mesophilic proteins have both similar polar and nonpolar contribution to the surface area and compactness. On the other hand, salt bridges and main chain hydrogen bonds show an increase in the majority of thermophilic proteins in comparison to their mesophilic homologues. In addition, in thermophilic proteins hydrophobic residues are significantly more frequent, while polar residues are less. These findings indicate that thermostable proteins through evolution adopt several different strategies to withstand high temperature environments.     

Electrostatic contribution of surface charge residues to the stability of a thermophilic protein: benchmarking experimental and predicted pKa values

Optimization of the surface charges is a promising strategy for increasing thermostability of proteins. Electrostatic contribution of ionizable groups to the protein stability can be estimated from the differences between the pKa values in the folded and unfolded states of a protein. Using this pKa-shift approach, we experimentally measured the electrostatic contribution of all aspartate and glutamate residues to the stability of a thermophilic ribosomal protein L30e from Thermococcus celer. The pKa values in the unfolded state were found to be similar to model compound pKas. The pKa values in both the folded and unfolded states obtained at 298 and 333 K were similar, suggesting that electrostatic contribution of ionizable groups to the protein stability were insensitive to temperature changes. The experimental pKa values for the L30e protein in the folded state were used as a benchmark to test the robustness of pKa prediction by various computational methods such as H++, MCCE, MEAD, pKD, PropKa, and UHBD. Although the predicted pKa values were affected by crystal contacts that may alter the side-chain conformation of surface charged residues, most computational methods performed well, with correlation coefficients between experimental and calculated pKa values ranging from 0.49 to 0.91 (p<0.01). The changes in protein stability derived from the experimental pKa-shift approach correlate well (r = 0.81) with those obtained from stability measurements of charge-to-alanine substituted variants of the L30e protein. Our results demonstrate that the knowledge of the pKa values in the folded state provides sufficient rationale for the redesign of protein surface charges leading to improved protein stability.     

Protein thermostability: structure-based difference of residual properties between thermophilic and mesophilic proteins

Structure-based differences of residual properties between 20 pairs of thermophilic and mesophilic proteins were statistically analyzed to elucidate the factors governing protein thermostability. This study analyzed the distributions of outer residues, inner residues, flexible residues, rigid residues, hydrogen bonds, salt bridges, cation–pi interactions, and disulfide bonds in each protein in terms of residual structural states, which were determined as five kinds of states under residual packing value. Their structural patterns found in thermophilic protein groups were compared with those of mesophilic protein groups for showing distinctive difference of residual properties. The results of statistical tests (t-test) revealed that flexible residues in fully-exposed state and boundary state, salt bridges in exposed state, and hydrogen bonds in well-buried state could be critical factors related with protein thermostability. Such structure-based differences of residual properties would help to develop a strategy for enhancing protein thermostability.