A search for structural factors responsible for the stability of proteins from thermophilic organisms

A database was designed to include 392 pairs of homologous proteins from thermophilic and mesophilic organisms. Proteins from thermophilic organisms proved to contain more atom-atom contacts per residue as compared with their mesophilic homologs. Solvent-accessible exterior amino acid residues contribute to the increase in the number of contacts. The amino acid composition was analyzed for internal (solvent-inaccessible) and exterior amino acid residues of thermophilic and mesophilic proteins. The exterior residues of thermophils have higher contents of Lys, Arg, and Glu and lower contents of Ala, Asp, Asn, Gln, Ser, and Thr as compared with mesophilic proteins. Interior protein regions did not differ in amino acid composition.     

Comparative thermodynamic elucidation of the structural stability of thermophilic proteins

Differential scanning calorimetry, circular dichroism, and visible absorption spectrophotometry were employed to elucidate the structural stability of thermophilic phycocyanin derived from Cyanidium caldarium, a eucaryotic organism which contains a nucleus, grown in acidic conditions (pH 3.4) at 54°C. The obtained results were compared with those previously reported for thermophilic phycocyanin derived from Synechococcus lividus, a procaryote containing no organized nucleus, grown in alkaline conditions (pH 8.5) at 52°C. The temperature of thermal unfolding (t_d) was found to be comparable between C. caldarium (73°C) and S. lividus (74°C) phycocyanins. The apparent free energy of unfolding (ΔG_[urea]=0) at zero denaturant (urea) concentration was also comparable: 9.1 and 8.7 kcal/mole for unfolding the chromophore part of the protein, and 5.0 and 4.3 kcal/mole for unfolding the apoprotein part of the protein, respectively. These values of t_d and ΔG_[urea]=0 were significantly higher than those previously reported for mesophilic Phormidium luridum phycocyanin (grown at 25°C). These findings revealed that relatively higher values of t_d and ΔG_[urea]=0 were characteristics of thermophilic proteins. In contrast, the enthalpies of completed unfolding (ΔH_d) and the half-completed unfolding (ΔH_d)1/2 for C. caldarium phycocyanin were much lower than those for S. lividus protein (89 versus 180 kcal/mole and 62 versus 115 kcal/mole, respectively). Factors contributing to a lower ΔH_d in C. caldarium protein and the role of charged groups in enhancing the stability of thermophilic proteins were discusse.     

Elucidation of factors responsible for enhanced thermal stability of proteins: a structural genomics based study

Understanding the molecular basis for the enhanced stability of proteins from thermophiles has been hindered by a lack of structural data for homologous pairs of proteins from thermophiles and mesophiles. To overcome this difficulty, complete genome sequences from 9 thermophilic and 21 mesophilic bacterial genomes were aligned with protein sequences with known structures from the protein data bank. Sequences with high homology to proteins with known structures were chosen for further analysis. High quality models of these chosen sequences were obtained using homology modeling. The current study is based on a data set of models of 900 mesophilic and 300 thermophilic protein single chains and also includes 178 templates of known structure. Structural comparisons of models of homologous proteins allowed several factors responsible for enhanced thermostability to be identified. Several statistically significant, specific amino acid substitutions that occur going from mesophiles to thermophiles are identified. Most of these are at solvent-exposed sites. Salt bridges occur significantly more often in thermophiles. The additional salt bridges in thermophiles are almost exclusively in solvent-exposed regions, and 35% are in the same element of secondary structure. Helices in thermophiles are stabilized by intrahelical salt bridges and by an increase in negative charge at the N-terminus. There is an approximate decrease of 1% in the overall loop content and a corresponding increase in helical content in thermophiles. Previously overlooked cation-pi interactions, estimated to be twice as strong as ion-pairs, are significantly enriched in thermophiles. At buried sites, statistically significant hydrophobic amino acid substitutions are typically consistent with decreased side chain conformational entropy.     

Pressure enhances thermal stability of DNA polymerase from three thermophilic organisms

DNA polymerases derived from three thermophilic microorganisms, Pyrococcus strain ES4, Pyrococcus furiosus, and Thermus aquaticus, were stabilized in vitro by hydrostatic pressure at denaturing temperatures of 111°C, 107.5°C, and 100°C (respectively). Inactivation rates, as determined by enzyme activity measurements, were measured at 3, 45, and 89 MPa. Half-lives of P. strain ES4, P. furiosus, and T. aquaticus DNA polymerases increased from 5.0, 6.9, and 5.2 minutes (respectively) at 3 MPa to 12, 36, and 13 minutes (respectively) at 45 MPa. A pressure of 89 MPa further increased the half-lives of P. strain ES4 and T. aquaticus DNA polymerases to 26 and 39 minutes, while the half-life of P. furiosus DNA polymerase did not increase significantly from that at 45 MPa. The decay constant for P. strain ES4 and T. aquaticus polymerases decreased exponentially with increasing pressure, reflecting an observed change in volume for enzyme inactivation of 61 and 73 cm³/mol, respectively. Stabilization by pressure may result from pressure effects on thermal unfolding or pressure retardation of unimolecular inactivation of the unfolded state. Regardless of the mechanism, pressure stabilization of proteins could explain the previously observed extension of the maximum temperature for survival of P. strain ES4 and increase the survival of thermophiles in thermally variable deep-sea environments such as hydrothermal vents. Received: September 12, 1997 / Accepted: February 24, 1998     

Comparison of energy components of proteins from thermophilic and mesophilic organisms

In order to infer the energetic determinants of thermophilic proteins, molecular mechanics calculations were applied to five proteins from thermophilic eubacteria and their mesophilic homologs. The energy function includes a hydration term as well as the electrostatic contribution from the solvent in addition to the usual conformational energy terms. We calculated energy values for three different states of each protein: the native, near-native, and unfolded structures. The energy difference and its components between pairs of these states were compared. The hypothetical near-native structures have almost the same backbone conformation as the native structure but with largely distorted side-chain packing, thus enabling us to extract the energy components important for stabilizing the native backbone topology itself, irrespective of structural details. It was found that the sum of the electrostatic and hydration energies, although of large positive values, were consistently lower for the thermophilic proteins than for their mesophilic counterparts. This trend was observed in the energy difference not only between the native and unfolded structures, but also between the near-native and unfolded structures. In contrast, the energy components regarding side-chain packing did not show any clear tendency. These results suggest that the thermophilic proteins are stabilized so that the precise packing of the native structure does not significantly affect the stability. Implications of this conclusion are also discussed.     

Search for structural similarity in proteins

MOTIVATION: The expanding protein sequence and structure databases await methods allowing rapid similarity search. Geometric parameters-dihedral angle between two sequential peptide bond planes (V) and radius of curvature (R) as they appear in pentapeptide fragments in polypeptide chains-are proposed for use in evaluating structural similarity in proteins (VeaR). The parabolic (empirical) function expressing the radius of curvature's dependence on the V-angle in model polypeptides is altered in real proteins in a form characteristic for a particular protein. This can be used as a criterion for judging similarity. RESULTS: A structural comparison of proteins representing a wide spectrum of structures was assessed versus sequence similarity analysis based on the genetic semihomology algorithm. The term 'consensus structure', analogous to 'consensus sequence', was introduced for the serpine family. AVAILABILITY: Semihom-sequence comparison freely available on request from J. Leluk. VeaR-structural comparison freely available on request from I. Roterman.     

Thermodynamic stability and folding of proteins from hyperthermophilic organisms

Life grows almost everywhere on earth, including in extreme environments and under harsh conditions. Organisms adapted to high temperatures are called thermophiles (growth temperature 45-75 degrees C) and hyperthermophiles (growth temperature >or= 80 degrees C). Proteins from such organisms usually show extreme thermal stability, despite having folded structures very similar to their mesostable counterparts. Here, we summarize the current data on thermodynamic and kinetic folding/unfolding behaviors of proteins from hyperthermophilic microorganisms. In contrast to thermostable proteins, rather few (i.e. less than 20) hyperthermostable proteins have been thoroughly characterized in terms of their in vitro folding processes and their thermodynamic stability profiles. Examples that will be discussed include co-chaperonin proteins, iron-sulfur-cluster proteins, and DNA-binding proteins from hyperthermophilic bacteria (i.e. Aquifex and Theromotoga) and archea (e.g. Pyrococcus, Thermococcus, Methanothermus and Sulfolobus). Despite the small set of studied systems, it is clear that super-slow protein unfolding is a dominant strategy to allow these proteins to function at extreme temperatures.     

Important inter-residue contacts for enhancing the thermal stability of thermophilic proteins

Proteins from thermophilic organisms exhibit high thermal stability, but have structures that are very similar to their mesophilic homologues. In order to gain insight into the basis of thermostability, we have analyzed the medium- and long-range contacts in mesophilic and thermophilic proteins of 16 different families. We found that the thermophiles prefer to have contacts between residues with hydrogen-bond-forming capability. Apart from hydrophobic contacts, more contacts are observed between polar and non-polar residues in thermophiles than mesophiles. Residue-wise analysis showed that Tyr has good contacts with several other residues, and Cys has considerably higher long-range contacts in thermophiles compared with mesophiles. Furthermore, the residues occurring in the range of 31-34 residues apart in the sequence contribute significant long-range contacts to the stability of thermophilic proteins.     

Thermodynamics elucidation of the structural stability of a thermophilic protein

The structural stability of the protein, phycocyanin isolated from two strains of cyanophyta, Synechococcus lividus (thermophile) and Phormidium luridum (mesophile), are investigated by comparative thermal and denaturant unfolding, using differential scanning calorimetry, visible absorption spectrophotometry, and circular dichroism. The thermophilic protein exhibits a much higher temperature and enthalpy of unfolding from the native to the denatured state. The concentration of urea at half-completion of thermal unfolding is essentially the same between the thermophilic and mesophilic proteins; in contrast, the corresponding temperature and the enthalpy of thermal unfolding are much higher for the thermophilic protein. In addition, the concentration of urea at which the non-thermal (denaturant) unfolding of protein is half-completed, as detected by either circular dichroism or absorption spectroscopy, is significantly higher in the thermophilic protein, while the apparent free energy of unfolding only shows a moderate difference between the two proteins. The distinct differences in the enthalpy of thermal unfolding and the free energy of denaturant unfolding are interpreted in terms of a significant entropy change associated with the unfolding of these proteins. This entropy contribution is much higher in the thermophilic protein, and may be derived from its more rigid overall structure that possesses higher internal hydrophobicity and stronger internal packing.     

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.     

Role of cation-pi interactions to the stability of thermophilic proteins

Elucidating the factors responsible for exhibiting extreme thermal stability of thermophilic proteins is very important for an understanding of the mechanism of protein stability, as well as to design stable proteins. In this work, we have analyzed the influence of cation-pi interactions to enhance the stability from mesophilic to thermophilic proteins. The favorable residue pairs forming such a system of interactions have been brought out. We found that the Tyr has a greater number of such interactions with Lys in thermophilic proteins. Specifically, the same Lys would experience a greater number of cation-pi interactions with several Tyr residues in thermophiles. On the other hand, the influence of Phe in making cation-pi interactions is higher in mesophiles than in thermophiles. Further, a network of cation-pi interactions are maintained by Lys in thermophiles, whereas Arg plays a major role in mesophilic proteins. Moreover, atoms that have a substantial positive charge in both Lys and Arg make a more significant contribution for cation-pi interactions than do cationic group atoms.     

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.     

Explanation of the stability of thermophilic proteins based on unique micromorphology

Two mesophilic/thermophilic variants of the G-domain of the elongation factor Tu were studied via molecular dynamics simulations. By analyzing the simulation data via the Voronoi space tessellation, we have found that the two proteins have the same macromolecular packing, while the water-exposed surface area is larger for the thermophile. A larger coordination with water is probably due to a peculiar corrugation of the exposed surface of this species. From an enthalpic point of view, the thermophile shows a larger number of intramolecular hydrogen bonds, stronger electrostatic interactions, and a flatter free-energy landscape. Overall, the data suggest that the specific hydration state enhances macromolecular fluctuations but, at the same time, increases thermal stability.     

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.     

Packing-based difference of structural features between thermophilic and mesophilic proteins

Twenty pairs of thermophilic and mesophilic proteins were compared in terms of residue packing distribution to obtain structural features related to protein thermostability. Based on residue packing concept, structural features of residues such as residue packing distribution, inner/outer position, secondary structure and water solvation were investigated. The statistical tests revealed that higher frequency in well-packed state of residues, lower frequency in exposed state and higher frequency in well-packed state of inner positioned residues, and higher frequency in well-packed state of 3/10 helix residues could be general structural features thermophilic proteins have.     

Importance of main-chain hydrophobic free energy to the stability of thermophilic proteins

Living organisms are found in the most unexpected places, including deep-sea vents at 100 degrees C and several hundred bars pressure, in hot springs. Needless to say, the proteins found in thermophilic species are much more stable than their mesophilic counterparts. There are no obvious reasons to say that one would be more stable than others. Even examination of the amino acids and comparison of structural features of thermophiles with mesophilies cannot bring satisfactory explanation for the thermal stability of such proteins. In order to bring out the hidden information behind the thermal stabilization of such proteins in terms of energy factors and their combinations, analysis were made on good resolution structures of thermophilic and their mesophilic homologous from 23 different families. From the structural coordinates, free energy contributions due to hydrophobic, electrostatic, hydrogen bonding, disulfide bonding and van der Waals interactions are computed. In this analysis, a vast majority of thermophilic proteins adopt slightly lower free energy contribution in each energy terms than its mesophilic counterparts. The major observation noted from this study is the lower hydrophobic free energy contribution due to carbon atoms and main-chain nitrogen atoms in all the thermophilic proteins. The possible combination of different free energy terms shows majority of the thermophilic proteins have lower free energy strategy than their mesophilic homologous. The derived results show that the hydrophobic free energy due to carbon and nitrogen atoms and such combinations of free energy components play a vital role in the thermostablisation of such proteins.     

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.     

In vitro replication slippage by DNA polymerases from thermophilic organisms

Replication slippage of DNA polymerases is a potential source of spontaneous genetic rearrangements in prokaryotic and eukaryotic cells. Here we show that different thermostable DNA polymerases undergo replication slippage in vitro, during single-round replication of a single-stranded DNA template carrying a hairpin structure. Low-fidelity polymerases, such as Thermus aquaticus (Taq), high-fidelity polymerases, such as Pyrococcus furiosus (Pfu) and a highly thermostable polymerase from Pyrococcus abyssi (Pyra exo(-)) undergo slippage. Thermococcus litoralis DNA polymerase (Vent) is also able to slip; however, slippage can be inhibited when its strand-displacement activity is induced. Moreover, DNA polymerases that have a constitutive strand-displacement activity, such as Bacillus stearothermophilus DNA polymerase (Bst), do not slip. Polymerases that slip during single-round replication generate hairpin deletions during PCR amplification, with the exception of Vent polymerase because its strand-displacement activity is induced under these conditions. We show that these hairpin deletions occurring during PCR are due to replication slippage, and not to a previously proposed process involving polymerization across the hairpin base.