Using a strategy based on the concept of convergent evolution to identify residue substitutions responsible for thermal adaptation

Factors that are related to thermostability of proteins have been extensively studied in recent years, especially by comparing thermophiles and mesophiles. However, most of them are global characters. It is still not clear how to identify specific residues or fragments which may be more relevant to protein thermostability. Moreover, some of the differences among the thermophiles and mesophiles may be due to phylogenetic differences instead of thermal adaptation. To resolve these problems, I adopted a strategy to identify residue substitutions evolved convergently in thermophiles or mesophiles. These residues may therefore be responsible for thermal adaptation. Four classes of genomes were utilized in this study, including thermophilic archaea, mesophilic archaea, thermophilic bacteria, and mesophilic bacteria. For most clusters of orthologous groups (COGs) with sequences from all of these four classes of genomes, I can identify specific residues or fragments that may potentially be responsible for thermal adaptation. Functional or structural constraints (represented as sequence conservation) were suggested to have higher impact on thermal adaptation than secondary structure or solvent accessibility does. I further compared thermophilic archaea and mesophilic bacteria, and found that the most diverged fragments may not necessarily correspond to the thermostability-determining ones. The usual approach to compare thermophiles and mesophiles without considering phylogenetic relationships may roughly identify sequence features contributing to thermostability; however, to specifically identify residue substitutions responsible for thermal adaptation, one should take sequence evolution into consideration.     

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.     

Subunit interfaces of oligomeric hyperthermophilic enzymes display enhanced compactness

Enzymes from thermophilic and, particularly, from hyperthermophilic organisms are surprisingly stable. Understanding the molecular origin of protein thermostability and thermoactivity attracted the interest of many scientists both for the perspective comprehension of the principles of protein structure and for the possible biotechnological applications through protein engineering. Comparative studies at sequence and structure levels were aimed at detecting significant differences of structural parameters related to protein stability between thermophilic and hyperthermophilic proteins and their mesophilic homologs. In a recent work, we focused attention on structural adaptation occurring at the subunit interface of oligomeric hyper- and thermostable enzymes. A set of structural and chemico-physical parameters were compared to those observed at the corresponding interfaces of homologous mesophilic proteins. Among the most significant variations, a general increase of interface apolarity and packing density in hyperthermophilic enzymes were found. This work was therefore aimed at elucidating whether the increased packing observed is reached also through the reduction of interface cavity number and volume. The results indicate that number of cavities tends to be relatively constant while cavity volume tends to decrease in the hyperthermophilic interfaces. The cavity apolarity increases in thermophiles but, apparently, not in hyperthermophiles. Moreover, interface hot spot residues of the mesophilic interfaces tend to be conserved in the extremophilic counterparts.     

How do thermophilic proteins resist aggregation?

Aggregation is an ancient threat that must be overcome by proteins from all organisms to maintain their native functional states. This is essential for the maintenance of metabolic flux and viability of their cellular machineries. Here, we compare the aggregation-resistance strategies adapted by the thermophilic proteins and their mesophilic homologs using a dataset of 373 protein families. Like their mesophilic homologs, the thermophilic protein sequences also contain potential aggregation prone regions (APRs), capable of forming cross-β motif and amyloid-like fibrils. Tetrapeptide and hexapeptide amyloid-like fibril forming sequence patterns and experimentally proven amyloid-like fibril forming peptide sequences were also detected in the thermophilic proteins. Both the thermophilic and the mesophilic proteins use similar strategies to resist aggregation. However, the thermophilic proteins show superior utilization of these strategies. The thermophilic protein monomers show greater ability to "stow away" the APRs in the hydrophobic cores to protect them from solvent exposure. The thermophilic proteins are also better at gatekeeping the APRs by surrounding them with charged residues (Asp, Glu, Lys, and Arg) and Pro to a greater extent. While thermophilic and mesophilic proteins in our dataset are highly homologous and show strong overall sequence conservation, the APRs are not conserved between the homologs. These findings indicate that evolution is working to avoid amyloidogenic regions in proteins. Our results are also consistent with the observation that thermophilic cells often accumulate small molecule osmolytes capable of stabilizing their proteins and other macromolecules. This study has important implications for rational design and formulation of therapeutic proteins and antibodies.     

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.     

The thermostability of DNA-binding protein HU from mesophilic, thermophilic, and extreme thermophilic bacteria

Based on primary structure comparison between four highly homologous DNA-binding proteins (HUs) displaying differential thermostability, we have employed in vitro site-directed mutagenesis to decipher their thermostability mechanism at the molecular level. The contribution of the 11 amino acids that differ between the thermophilic HUBst from Bacillus stearothermophilus (Tm = 61.6 degrees C) and the mesophilic HUBsu from Bacillus subtilis (Tm = 39.7 degrees C) was evaluated by replacing these amino acids in HUBst with their mesophilic counterparts. Among 11 amino acids, three residues, Gly-15, Glu-34, and Val-42, which are highly conserved in the thermophilic HUs, have been found to be responsible for the thermostability of HUBst. These amino acids in combination (HUBst-G15E/E34D/V42I) reduce the thermostability of the protein (Tm = 45.1 degrees C) at the level of its mesophilic homologue HUBsu. By replacing these amino acids in HUBsu with their thermophilic counterparts, the HUBsu-E15G/D34E/142V mutant was generated with thermostability (Tm = 57.8 degrees C) at the level of thermophilic HUBst. Employing the same strategy, we generated several mutants in the extremely thermophilic HUTmar from Thermotoga maritima (Tm = 80.5 degrees C), and obtained data consistent with the previous results. The triplet mutant HUTmar-G15E/E34D/V421 (Tm = 35.9 degrees C) converted the extremely thermophilic protein HUTmar to mesophilic. The various forms of HU proteins were overproduced in Escherichia coli, highly purified, and the thermostability of the mutants confirmed by circular dichroism spectroscopy. The results presented here were elucidated on the basis of the X-ray structure of HUBst and HUTmar (our unpublished results), and their mechanism was proposed at the molecular level. The results clearly show that three individual local interactions located at the helix-turn-helix part of the protein are responsible for the stability of HU proteins by acting cooperatively in a common mechanism for thermostability.     

蛋白质空间结构属性与全基因组微生物耐热性的关系

 溶剂接触表面积、空腔个数和体积、紧密度、疏水性以及温度因子是影响蛋白质耐热的主要三维结构参数.挑选NCBI COG数据库中具有全基因组的单细胞微生物,选择其中三维结构已知的蛋白质作为研究对象,分析这些因素对细菌类和古细菌类微生物耐热性的影响.结果表明：（1）古细菌类蛋白质的空腔个数和体积与耐热性无关,极性面积和表面残基个数随耐热性增加而降低;（2）超高温细菌类蛋白质的分子量比较小,空腔个数和体积都小于常温蛋白质,而且空腔个数对稳定性的贡献大于空腔体积;（3）无论是古细菌还是细菌类蛋白质,疏水性和紧密度都不随耐热性变化,但暴露残基个数越多,蛋白质的耐热性越差;（4）两类蛋白质侧链的温度因子都高于主链,这与侧链的运动性（柔性）一般比主链高的实验结果一致;另外,超高温细菌类蛋白质的温度因子明显高于常温蛋白质.     

Role of loop dynamics in thermal stability of mesophilic and thermophilic adenylosuccinate synthetase: a molecular dynamics and normal mode analysis study

Enzymes from thermophiles are poorly active at temperatures at which their mesophilic homologs exhibit high activity and attain corresponding active states at high temperatures. In this study, comparative molecular dynamics (MD) simulations, supplemented by normal mode analysis, have been performed on an enzyme Adenylosuccinate synthetase (AdSS) from E. coli (mesophilic) and P. horikoshii (thermophilic) systems to understand the effects of loop dynamics on thermal stability of AdSS. In mesophilic AdSS, both ligand binding and catalysis are facilitated through the coordinated movement of five loops on the protein. The simulation results suggest that thermophilic P. horikoshii preserves structure and catalytic function at high temperatures by using the movement of only a subset of loops (two out of five) for ligand binding and catalysis unlike its mesophilic counterpart in E. coli. The pre-arrangement of the catalytic residues in P. horikoshii is well-preserved and salt bridges remain stable at high temperature (363K). The simulations suggest a general mechanism (including pre-arrangement of catalytic residues, increased polar residue content, stable salt bridges, increased rigidity, and fewer loop movements) used by thermophilic enzymes to preserve structure and be catalytically active at elevated temperatures.     

Using motif-based methods in multiple genome analyses: a case study comparing orthologous mesophilic and thermophilic proteins

Protein motifs represent highly conserved regions within protein families and are generally accepted to describe critical regions required for protein stability and/or function. In this comprehensive analysis, we present a robust, unique approach to identify and compare corresponding mesophilic and thermophilic sequence motifs between all orthologous proteins within 44 microbial genomes. Motif similarity is determined through global sequence alignment of mesophilic and thermophilic motif pairs, which are identified by a greedy algorithm. Our results reveal only modest correlation between motif and overall sequence similarity, highlighting the rationale of motif-based approaches in comprehensive multigenome comparisons. Conserved mutations reflect previously suggested physiochemical principles for conferring thermostability. Additionally, comparisons between corresponding mesophilic and thermophilic motif pairs provide key biochemical insights related to thermostability and can be used to test the evolutionary robustness of individual structural comparisons. We demonstrate the ability of our unique approach to provide key insights in two examples: the TATA-box binding protein and glutamate dehydrogenase families. In the latter example, conserved mutations hint at novel origins leading to structural stability differences within the hexamer structures. Additionally, we present amino acid composition data and average protein length comparisons for all 44 microbial genomes.     

Comparative molecular dynamics simulation studies for determining factors contributing to the thermostability of chemotaxis protein “CheY”

Comparative molecular dynamics simulations of chemotaxis protein “CheY” from thermophilic origin Thermotoga maritima and its mesophilic counterpart Salmonella enterica have been performed for 10?ns each at 300 and 350?K, and 20?ns each at 400 and 450?K. The trajectories were analyzed in terms of different factors like root-mean-square deviation, root-mean-square fluctuation, radius of gyration, solvent accessible surface area, H-bonds, salt bridge content, and protein–solvent interactions which indicate distinct differences between the two of them. The two proteins also follow dissimilar unfolding pathways. The overall flexibility calculated by the trace of the diagonalized covariance matrix displays similar flexibility of both the proteins near their optimum growth temperatures. However, at higher temperatures mesophilic protein shows increased overall flexibility than its thermophilic counterpart. Principal component analysis also indicates that the essential subspaces explored by the simulations of two proteins at different temperatures are nonoverlapping and they show significantly different directions of motion. However, there are significant overlaps within the trajectories and similar direction of motions are observed for both proteins at 300?K. Overall, the mesophilic protein leads to increased conformational sampling of the phase space than its thermophilic counterpart. This is the first ever study of thermostability of CheY protein homologs by using protein dynamism as a main impact. Our study might be used as a model for studying the molecular basis of thermostability of two homologous proteins from two organisms living at different temperatures with less visible differences.     

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.     

Density discriminates between thermophilic and mesophilic proteins

Despite an intense interest and a remarkable number of studies on the subject, the relationships between thermostability and (primary, secondary and tertiary) structure of proteins are still not fully understood. Here, comparing the protein density – defined by the ratio between the residue number and protein excluded volume – for a set of thermophilic/mesophilic pairs, we provide evidence that this property is connected to the optimal growth temperature. In particular, our results indicate that thermophilic proteins have – in general – a lower density with respect to the mesophilic counterparts, being such a correlation more pronounced for optimal growth temperature differences greater than 40°C. The effect of the protein thermostability changes on the molecular shape is also presented.     

Distance-dependent statistical potentials for discriminating thermophilic and mesophilic proteins

Identification of the characteristic structural patterns responsible for protein thermostability is theoretically important and practically useful but largely remains an open problem. These patterns may be revealed through comparative study on thermophilic and mesophilic proteins that have distinct thermostability. In this study, we constructed several distance-dependant potentials from thermophilic and mesophilic proteins. These potentials were then used to evaluate the structural difference between thermophilic and mesophilic proteins. We found that using the subtraction or division of the potentials derived from thermophilic and mesophilic proteins can dramatically increase the discriminatory ability. This approach revealed that the ability to distinct the subtle structural features responsible for protein thermostability may be effectively enhanced through rationally designed comparative study.     

"Hot cores" in proteins: Comparative analysis of the apolar contact area in structures from hyper/thermophilic and mesophilic organisms

Background  

A wide variety of stabilizing factors have been invoked so far to elucidate the structural basis of protein thermostability. These include, amongst the others, a higher number of ion-pairs interactions and hydrogen bonds, together with a better packing of hydrophobic residues. It has been frequently observed that packing of hydrophobic side chains is improved in hyperthermophilic proteins, when compared to their mesophilic counterparts. In this work, protein crystal structures from hyper/thermophilic organisms and their mesophilic homologs have been compared, in order to quantify the difference of apolar contact area and to assess the role played by the hydrophobic contacts in the stabilization of the protein core, at high temperatures.     

Amino acid coupling patterns in thermophilic proteins

Structural analysis is useful in elucidating structural features responsible for enhanced thermal stability of proteins. However, due to the rapid increase of sequenced genomic data, there are far more protein sequences than the corresponding three-dimensional (3D) structures. The usual sequence-based amino acid composition analysis provides useful but simplified clues about the amino acid types related to thermal stability of proteins. In this work, we developed a statistical approach to identify the significant amino acid coupling sequence patterns in thermophilic proteins. The amino acid coupling sequence pattern is defined as any 2 types of amino acids separated by 1 or more amino acids. Using this approach, we construct the rho profiles for the coupling patterns. The rho value gives a measure of the relative occurrence of a coupling pattern in thermophiles compared with mesophiles. We found that thermophiles and mesophiles exhibit significant bias in their amino acid coupling patterns. We showed that such bias is mainly due to temperature adaptation instead of species or GC content variations. Though no single outstanding coupling pattern can adequately account for protein thermostability, we can use a group of amino acid coupling patterns having strong statistical significance (p values < 10(-7)) to distinguish between thermophilic and mesophilic proteins. We found a good correlation between the optimal growth temperatures of the genomes and the occurrences of the coupling patterns (the correlation coefficient is 0.89). Furthermore, we can separate the thermophilic proteins from their mesophilic orthologs using the amino acid coupling patterns. These results may be useful in the study of the enhanced stability of proteins from thermophiles-especially when structural information is scarce. Proteins 2005. (c) 2005 Wiley-Liss, Inc.     

氨基酸组成对蛋白质耐热性的影响

为了研究一级结构对蛋白质耐热性的影响,利用软件DNAMAN对16个家族32种蛋白质序列进行了氨基酸含量分析,并统计分析了氨基酸组成对蛋白质耐热性的影响。通过比较同一家族的高低温蛋白质序列及16个家族中所有高温和低温蛋白质序列中氨基酸含量的变化可以推断(从低温到高温)：Ser、Cys．含量降低显著,Arg、Ile、Pro含量升高显著。由此可知高温蛋白质倾向于含有疏水性氨基酸而避免亲水性氨基酸。     

Thermostability of membrane protein helix-helix interaction elucidated by statistical analysis

A prerequisite for the survival of (micro)organisms at high temperatures is an adaptation of protein stability to extreme environmental conditions. In contrast to soluble proteins, where many factors have already been identified, the mechanisms by which the thermostability of membrane proteins is enhanced are almost unknown. The hydrophobic membrane environment constrains possible stabilizing factors for transmembrane domains, so that a difference might be expected between soluble and membrane proteins. Here we present sequence analysis of predicted transmembrane helices of the genomes from eight thermophilic and 12 mesophilic organisms. A comparison of the amino acid compositions indicates that more polar residues can be found in the transmembrane helices of thermophilic organisms. Particularly, the amino acids aspartic acid and glutamic acid replace the corresponding amides. Cysteine residues are found to be significantly decreased by about 70% in thermophilic membrane domains suggesting a non-specific function of most cysteine residues in transmembrane domains of mesophilic organisms. By a pair-motif analysis of the two sets of transmembrane helices, we found that the small residues glycine and serine contribute more to transmembrane helix-helix interactions in thermophilic organisms. This may result in a tighter packing of the helices allowing more hydrogen bond formation.     

Robust determinants of thermostability highlighted by a codon frequency index capable of discriminating thermophilic from mesophilic genomes

Can genome analysis tell us about the lifestyle of an organism? We ask this question considering a thorough cross comparison of thermophilic and mesophilic genomes, since presently the number of available genomes is enough to ensure statistical significance of the results. We analyze, by means of principal component analysis (PCA), the codon composition of a database comprising 116 genomes, selected so as to include one species for each genus and show that a cross genomic approach can allow the extraction of common determinants of thermostability at the genome level. The results of our analysis indicate that all the known features of thermostability can be found in the 64 component loadings of the second principal axis of PCA. By this, we develop an index of thermostability whose discriminative power between mesophiles and thermophiles scores with 98% accuracy at the genome level and with 95% accuracy at the protein sequence level. We also prove that these results are not due to phylogenetic differences between archaea and bacteria.     

嗜热与嗜常温微生物的蛋白质氨基酸组成比较

嗜热微生物的嗜热特性与其蛋白质的高度热稳定性紧密相关。为了探索嗜热蛋白质的热稳定机制,比较嗜热和嗜常温微生物的蛋白质在氨基酸组成上的差别,收集110对分别来自嗜热和嗜常温微生物的同源蛋白质序列,比较两组蛋白质各种氨基酸含量以及疏水性氨基酸组成、疏水性指数和荷电氨基酸组成的差别,结果两者在多种氨基酸含量上存在微小但统计学上显著的差别,嗜热蛋白质比嗜常温蛋白质具有较高的平均疏水性和荷电氨基酸组成。对两组蛋白质的“脂肪族氨基酸指数”进行分析,证明嗜热蛋白质之所以具有较高的脂肪族氨基酸指数是由于其亮氨酸含量较高,与影响该指数的其它几种氨基酸无关;从而认为该指数的意义值得怀疑。通过对大量同源嗜热蛋白质和嗜常温蛋白质氨基酸组成的比较,能够揭示一些有关蛋白质热稳定性的普遍规律。