Sequence and structural parameters enhancing adaptation of proteins to low temperatures

A systematic analysis compared sequence and structural parameters distributions between 13 pairs of psychrophilic and mesophilic proteins for elucidating the cold adaptation parameters. The results of statistical test (t-test) revealed that helical content, tight turn content, disulfide bonds and hydrogen bonds do not show significant difference between psychrophilic and mesophilic proteins. However, it was demonstrated in this study that a larger proportion of open beta-turn in psychrophilic proteins is an effective parameter in specific activity at low temperature. In addition, substitution of amino acids of charged and aliphatic groups with amino acids of tiny and small groups in protein chains, tight turns and alpha-helices in the direction from mesophilic to psychrophilic proteins is one of the mechanisms of low temperature adaptation. Such sequence and structural parameter differences would help to develop a strategy for designing cold-adapted proteins.     

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.     

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.     

Nonlinear methods in the analysis of protein sequences: a case study in rubredoxins

Two computational methods widely used in time series analysis were applied to protein sequences, and their ability to derive structural information not directly accessible through classical sequence comparisons methods was assessed. The primary structures of 19 rubredoxins of both mesophilic and thermophilic bacteria, coded with hydrophobicity values of amino acid residues, were considered as time series and were analyzed by 1) recurrence quantification analysis and 2) spectral analysis of the sequence major eigenfunctions. The results of the two methods agreed to a large extent and generated a classification consistent with known 3D structural characteristics of the studied proteins. This classification separated in a clearcut manner a thermophilic protein from mesophilic proteins. The classification of primary structures given by the two dynamical methods was demonstrated to be basically different from classification stemming from classical sequence homology metrics. Moreover, on a more detailed scale, the method was able to discriminate between thermophilic and mesophilic proteins from a set of chimeric sequences generated from the mixing of a mesophilic (Rubr Clopa) and a thermophilic (Rubr Pyrfu) protein. Overall, our results point to a new way of looking at protein sequence comparisons.     

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.     

Contribution of electrostatic interactions, compactness and quaternary structure to protein thermostability: lessons from structural genomics of Thermotoga maritima

Studies of the structural basis of protein thermostability have produced a confusing picture. Small sets of proteins have been analyzed from a variety of thermophilic species, suggesting different structural features as responsible for protein thermostability. Taking advantage of the recent advances in structural genomics, we have compiled a relatively large protein structure dataset, which was constructed very carefully and selectively; that is, the dataset contains only experimentally determined structures of proteins from one specific organism, the hyperthermophilic bacterium Thermotoga maritima, and those of close homologs from mesophilic bacteria. In contrast to the conclusions of previous studies, our analyses show that oligomerization order, hydrogen bonds, and secondary structure play minor roles in adaptation to hyperthermophily in bacteria. On the other hand, the data exhibit very significant increases in the density of salt-bridges and in compactness for proteins from T.maritima. The latter effect can be measured by contact order or solvent accessibility, and network analysis shows a specific increase in highly connected residues in this thermophile. These features account for changes in 96% of the protein pairs studied. Our results provide a clear picture of protein thermostability in one species, and a framework for future studies of thermal adaptation.     

Intermolecular Interactions Drive Protein Adaptive and Coadaptive Evolution at Both Species and Population Levels

Proteins are the building blocks for almost all the functions in cells. Understanding the molecular evolution of proteins and the forces that shape protein evolution is essential in understanding the basis of function and evolution. Previous studies have shown that adaptation frequently occurs at the protein surface, such as in genes involved in host–pathogen interactions. However, it remains unclear whether adaptive sites are distributed randomly or at regions associated with particular structural or functional characteristics across the genome, since many proteins lack structural or functional annotations. Here, we seek to tackle this question by combining large-scale bioinformatic prediction, structural analysis, phylogenetic inference, and population genomic analysis of Drosophila protein-coding genes. We found that protein sequence adaptation is more relevant to function-related rather than structure-related properties. Interestingly, intermolecular interactions contribute significantly to protein adaptation. We further showed that intermolecular interactions, such as physical interactions, may play a role in the coadaptation of fast-adaptive proteins. We found that strongly differentiated amino acids across geographic regions in protein-coding genes are mostly adaptive, which may contribute to the long-term adaptive evolution. This strongly indicates that a number of adaptive sites tend to be repeatedly mutated and selected throughout evolution in the past, present, and maybe future. Our results highlight the important roles of intermolecular interactions and coadaptation in the adaptive evolution of proteins both at the species and population levels.     

Dynamic fluorescence studies of beta-glycosidase mutants from Sulfolobus solfataricus: effects of single mutations on protein thermostability

Multiple sequence alignment on 73 proteins belonging to glycosyl hydrolase family 1 reveals the occurrence of a segment (83-124) in the enzyme sequences from hyperthermophilic archaea bacteria, which is absent in all the mesophilic members of the family. The alignment of the known three-dimensional structures of hyperthermophilic glycosidases with the known ones from mesophilic organisms shows a similar spatial organizations of beta-glycosidases except for this sequence segment whose structure is located on the external surface of each of four identical subunits, where it overlaps two alpha-helices. Site-directed mutagenesis substituting N97 or S101 with a cysteine residue in the sequence of beta-glycosidase from hyperthermophilic archaeon Sulfolobus solfataricus caused some changes in the structural and dynamic properties as observed by circular dichroism in far- and near-UV light, as well as by frequency domain fluorometry, with a simultaneous loss of thermostability. The results led us to hypothesize an important role of the sequence segment present only in hyperthermophilic beta-glycosidases, in the thermal adaptation of archaea beta-glycosidases. The thermostabilization mechanism could occur as a consequence of numerous favorable ionic interactions of the 83-124 sequence with the other part of protein matrix that becomes more rigid and less accessible to the insult of thermal-activated solvent molecules.     

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.     

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.     

Transproteomic evidence of a loop-deletion mechanism for enhancing protein thermostability.

Understanding the molecular determinants of protein thermostability is of theoretical and practical importance. While numerous determinants have been suggested, no molecular feature has been judged of paramount importance, with the possible exception of ion-pair networks. The difficulty in identifying the main determinants may have been the limited structural information available on the thermostable proteins. Recently the complete genomes for mesophilic, thermophilic and hyperthermophilic organisms have been sequenced, vastly improving the potential for uncovering general trends in sequence and structure evolution related to thermostability and, thus, for isolating the more important determinants. From a comparative analysis of 20 complete genomes, we find a trend towards shortened thermophilic proteins relative to their mesophilic homologs. Moreover, sequence alignments to proteins of known structure indicate that thermophilic sequences are more likely than their mesophilic homologs to have deletions in exposed loop regions. The new genomes offer enough comparable sequences to compute meaningful statistics that point to loop deletion as a general evolutionary strategy for increasing thermostability.     

Different unfolding pathways for mesophilic and thermophilic homologues of serine hydroxymethyltransferase

To determine how much information can be transferred from folding and unfolding studies of one protein to another member of the same family or between the mesophilic and thermophilic homologues of a protein, we have characterized the equilibrium unfolding process of the dimeric enzyme serine hydroxymethyltransferase (SHMT) from two sources, Bacillus subtilis (bsSHMT) and Bacillus stearothermophilus (bstSHMT). Although the sequences of the two enzymes are highly identical ( approximately 77%) and homologous (89%), bstSHMT shows a significantly higher stability against both thermal and urea denaturation than bsSHMT. The GdmCl-induced unfolding of bsSHMT was found to be a two-step process with dissociation of the native dimer, resulting in stabilization of a monomeric species, followed by the unfolding of the monomeric species. A similar unfolding pathway has been reported for Escherichia coli aspartate aminotransferase, a member of the type I fold family of PLP binding enzymes such as SHMT, the sequence of which is only slightly identical ( approximately 14%) with that of SHMT. In contrast, for bstSHMT, a highly cooperative unfolding without stabilization of any monomeric intermediate was observed. These studies suggest that mesophilic proteins of the same structural family even sharing a low level of sequence identity may follow a common unfolding mechanism, whereas the mesophilic and thermophilic homologues of the same protein despite having a high degree of sequence identity may follow significantly different unfolding mechanisms.     

Oligomeric integrity--the structural key to thermal stability in bacterial alcohol dehydrogenases.

Principles of protein thermostability have been studied by comparing structures of thermostable proteins with mesophilic counterparts that have a high degree of sequence identity. Two tetrameric NADP(H)-dependent alcohol dehydrogenases, one from Clostridium beijerinckii (CBADH) and the other from Thermoanaerobacter brockii (TBADH), having exceptionally high (75%) sequence identity, differ by 30 degrees in their melting temperatures. The crystal structures of CBADH and TBADH in their holo-enzyme form have been determined at a resolution of 2.05 and 2.5 A, respectively. Comparison of these two very similar structures (RMS difference in Calpha = 0.8 A) revealed several features that can account for the higher thermal stability of TBADH. These include additional ion pairs, "charged-neutral" hydrogen bonds, and prolines as well as improved stability of alpha-helices and tighter molecular packing. However, a deeper structural insight, based on the location of stabilizing elements, suggests that enhanced thermal stability of TBADH is due mainly to the strategic placement of structural determinants at positions that strengthen the interface between its subunits. This is also supported by mutational analysis of structural elements at critical locations. Thus, it is the reinforcement of the quaternary structure that is most likely to be a primary factor in preserving enzymatic activity of this oligomeric bacterial ADH at elevated temperatures.     

A comparative infrared spectroscopic study of glycoside hydrolases from extremophilic archaea revealed different molecular mechanisms of adaptation to high temperatures

The identification of the determinants of protein thermal stabilization is often pursued by comparing enzymes from hyperthermophiles with their mesophilic counterparts while direct structural comparisons among proteins and enzymes from hyperthermophiles are rather uncommon. Here, oligomeric beta-glycosidases from the hyperthermophilic archaea Sulfolobus solfataricus (Ss beta-gly), Thermosphaera aggregans (Ta beta-gly), and Pyrococcus furiosus (Pf beta-gly), have been compared. Studies of FTIR spectroscopy and kinetics of thermal inactivation showed that the three enzymes had similar secondary structure composition, but Ss beta-gly and Ta beta-gly (temperatures of melting 98.1 and 98.4 degrees C, respectively) were less stable than Pf beta-gly, which maintained its secondary structure even at 99.5 degrees C. The thermal denaturation of Pf beta-gly, followed in the presence of SDS, suggested that this enzyme is stabilized by hydrophobic interactions. A detailed inspection of the 3D-structures of these enzymes supported the experimental results: Ss beta-gly and Ta beta-gly are stabilized by a combination of ion-pairs networks and intrasubunit S-S bridges while the increased stability of Pf beta-gly resides in a more compact protein core. The different strategies of protein stabilization give experimental support to recent theories on thermophilic adaptation and suggest that different stabilization strategies could have been adopted among archaea.     

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.     

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.     

The thermal adaptation of the nitrogenase Fe protein from thermophilic Methanobacter thermoautotrophicus

The nitrogenase Fe protein is a key component of the biochemical machinery responsible for the process of biological nitrogen fixation. The Fe protein is a member of a class of nucleotide-binding proteins that couple the binding and hydrolysis of nucleoside triphosphates to conformational changes. The nucleotide-dependent conformational changes modulate the formation of a macromolecular complex, and some members of the class include Galpha, EF-Tu, and myosin. The members of this class are highly interesting model systems for the analysis of aspects of thermal adaptability, since their mechanisms involve protein conformational change and protein-protein interactions. In this study, we have used our extensive knowledge of the structure of the Azotobacter vinelandii nitrogenase Fe protein in multiple structural conformations, and standard homology modeling approaches have been used to generate reliable models of the Fe protein from thermophilic Methanobacter thermoautotrophicus in the analogous structural conformations. The resulting structural comparison reveals that thermal adaptation of the M. thermoautotrophicus Fe protein is conferred by a number of factors, including increased structural rigidity that results from various structural changes within the protein interior. The analysis of hypothetical docking models and nitrogenase complex structures provides insights into the thermal adaptation of the protein-protein interactions that support macromolecular complex formation and catalysis at higher temperatures.     

HU histone-like DNA-binding protein from <Emphasis Type="Italic">Thermus thermophilus</Emphasis>: structural and evolutionary analyses

The histone-like DNA-binding proteins (HU) serve as model molecules for protein thermostability studies, as they function in different bacteria that grow in a wide range of temperatures and show sequence diversity under a common fold. In this work, we report the cloning of the hutth gene from Thermus thermophilus, the purification and crystallization of the recombinant HUTth protein, as well as its X-ray structure determination at 1.7 Å. Detailed structural and thermodynamic analyses were performed towards the understanding of the thermostability mechanism. The interaction of HUTth protein with plasmid DNA in solution has been determined for the first time with MST. Sequence conservation of an exclusively thermophilic order like Thermales, when compared to a predominantly mesophilic order (Deinococcales), should be subject, to some extent, to thermostability-related evolutionary pressure. This hypothesis was used to guide our bioinformatics and evolutionary studies. We discuss the impact of thermostability adaptation on the structure of HU proteins, based on the detailed evolutionary analysis of the Deinococcus–Thermus phylum, where HUTth belongs. Furthermore, we propose a novel method of engineering thermostable proteins, by combining consensus-based design with ancestral sequence reconstruction. Finally, through the structure of HUTth, we are able to examine the validity of these predictions. Our approach represents a significant advancement, as it explores for the first time the potential of ancestral sequence reconstruction in the divergence between a thermophilic and a mainly mesophilic taxon, combined with consensus-based engineering.