Environment specific substitution tables for thermophilic proteins

In silico characterization of differentially expressed short-read nucleotide sequences identified in dieback stress-induced transcriptomic analysis reveals their role as antimicrobial peptides

We investigated the in silico characterization of short-length nucleotide sequences that were differentially expressed in dieback stress-induced transcriptomic analysis. They displayed homology with C-terminal flanking peptides and defensins-like proteins, revealing their antimicrobial activity. Their predicted fingerprints displayed protein signatures related to antimicrobial peptides. These short-length RGAs have been shown to possess structural motifs such as APLT P-type ATPase, casein kinase II (CK2), protein kinase 3, protein kinase C (PKC), and N-glycosylation site that are the attributes of disease resistance genes. The prediction of arginine and lysine residues in active binding sites in ligand docking analysis prophesied them as antimicrobial peptides due to their strong relation with antimicrobial activity. The in silico structural-functional characterization has predicted their role in resistance against microbial pathogens. Moreover, the predicted antimicrobial peptide regions showed their homology with the signature domain of PR-5-like protein and AMP family Thaumatin.

Structural features, temperature adaptation and industrial applications of microbial lipases from psychrophilic, mesophilic and thermophilic origins

Microbial lipases are very prominent biocatalysts because of their ability to catalyze a wide variety of reactions in aqueous and non-aqueous media. Here microbial lipases from different origins (psychrophiles, mesophiles, and thermophiles) have been reviewed. This review emphasizes an update of structural diversity in temperature adaptation and industrial applications, of psychrophilic, mesophilic, and thermophilic lipases. The microbial origins of lipases are logically dynamic, proficient, and also have an extensive range of industrial uses with the manufacturing of altered molecules. It is therefore of interest to understand the molecular mechanisms of adaptation to temperature in occurring lipases. However, lipases from extremophiles (psychrophiles, and thermophiles) are widely used to design biotransformation reactions with higher yields, fewer byproducts, or useful side products and have been predicted to catalyze those reactions also, which otherwise are not possible with the mesophilic lipases. Lipases as a multipurpose biological catalyst have given a favorable vision in meeting the needs of several industries such as biodiesel, foods, and drinks, leather, textile, detergents, pharmaceuticals, and medicals.

Accurate Prediction of Protein Thermodynamic Stability Changes upon Residue Mutation using Free Energy Perturbation

This work describes the application of a physics-based computational approach to predict the relative thermodynamic stability of protein variants, and evaluates the quantitative accuracy of those predictions compared to experimental data obtained from a diverse set of protein systems assayed at variable pH conditions. Physical stability is a key determinant of the clinical and commercial success of biological therapeutics, vaccines, diagnostics, enzymes and other protein-based products. Although experimental techniques for measuring the impact of amino acid residue mutation on the stability of proteins exist, they tend to be time consuming and costly, hence the need for accurate prediction methods. In contrast to many of the commonly available computational methods for stability prediction, the Free Energy Perturbation approach applied in this paper explicitly accounts for solvent effects and samples conformational dynamics using a rigorous molecular dynamics simulation process. On the entire validation dataset, consisting of 328 single point mutations spread across 14 distinct protein structures, our results show good overall correlation with experiment with an R² of 0.65 and a low mean unsigned error of 0.95 kcal/mol. Application of the FEP approach in conjunction with experimental assessment techniques offers opportunities to lower the time and expense of product development and reduce the risk of costly late-stage failures.

Parallel molecular mechanisms for enzyme temperature adaptation

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. We identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

Exploring mutable conserved sites and fatal non-conserved sites by random mutation of esterase from Sulfolobus tokodaii and subtilisin from Thermococcus kodakarensis

Homologous proteins differ in their amino acid sequences at several positions. Generally, conserved sites are recognized as not suitable for amino acid substitution, and thus in evolutionary protein engineering, non-conserved sites are often selected as mutation sites. However, there have also been reports of possible mutations in conserved sites. In this study, we explored mutable conserved sites and immutable non-conserved sites by testing random mutations of two thermostable proteins, an esterase from Sulfolobus tokodaii (Sto-Est) and a subtilisin from Thermococcus kodakarensis (Tko-Sub). The subtilisin domain of Tko-Sub needs Ca²⁺ ions and the propeptide domain for stability, folding and maturation. The results from the two proteins showed that about one-third of the mutable sites were detected in conserved sites and some non-conserved sites lost enzymatic activity at high temperatures due to mutation. Of the conserved sites in Sto-Est, the sites on the loop, on the surface, and far from the active site are more resistant to mutation. In Tko-Sub, the sites flanking Ca²⁺-binding sites and propeptide were undesirable for mutation. The results presented here serve as an index for selecting mutation sites and contribute to the expansion of available sequence range by introducing mutations at conserved sites.

Evolution of Protein Structure and Stability in Global Warming

This review focuses on the molecular signatures of protein structures in relation to evolution and survival in global warming. It is based on the premise that the power of evolutionary selection may lead to thermotolerant organisms that will repopulate the planet and continue life in general, but perhaps with different kinds of flora and fauna. Our focus is on molecular mechanisms, whereby known examples of thermoresistance and their physicochemical characteristics were noted. A comparison of interactions of diverse residues in proteins from thermophilic and mesophilic organisms, as well as reverse genetic studies, revealed a set of imprecise molecular signatures that pointed to major roles of hydrophobicity, solvent accessibility, disulfide bonds, hydrogen bonds, ionic and π-electron interactions, and an overall condensed packing of the higher-order structure, especially in the hydrophobic regions. Regardless of mutations, specialized protein chaperones may play a cardinal role. In evolutionary terms, thermoresistance to global warming will likely occur in stepwise mutational changes, conforming to the molecular signatures, such that each "intermediate" fits a temporary niche through punctuated equilibrium, while maintaining protein functionality. Finally, the population response of different species to global warming may vary substantially, and, as such, some may evolve while others will undergo catastrophic mass extinction.

Investigation into archaeal extremophilic lifestyles through comparative proteogenomic analysis

Archaea are a group of primary life forms on Earth and could thrive in many unique environments. Their successful colonization of extreme niches requires corresponding adaptations at proteogenomic level in order to maintain stable cellular structures and active physiological functions. Although some studies have already investigated the extremophilic lifestyles of archaeal species based on genomic features and protein structures, there is a lack of comparative proteogenomic analysis in a large scale. In this study, we explored 686 high-quality archaeal genomes (proteomes) sourced from the Pathosystems Resource Integration Center (PATRIC) database. General patterns of genomic features such as genome size, coding capacity (coding genes and non-coding regions), and G + C contents were re-confirmed. Protein domain distribution patterns were then identified across archaeal species. Domains with unknown functions (DUFs) and mini proteins were investigated in terms of their distributions due to their importance in archaeal physiological functions. In addition, physicochemical properties of protein sequences, such as stability, hydrophobicity, isoelectric point, aromaticity and amino acid compositions in corresponding archaeal groups were compared. Unique features associated with extremophilic lifestyles were observed, which suggested that evolutionary adaptations to different extreme environments had intrinsic impacts on archaeal protein features. Taken together, this systematic study facilitates a better understanding of the mechanisms behind the extremophilic lifestyles of archaeal species, which will further contribute to the evolutionary explorations of archaeal adaptations both experimentally and theoretically in the future studies.Communicated by Ramaswamy H. Sarma.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

Novel Sequence Feature of SecA Translocase Protein Unique to the Thermophilic Bacteria: Bioinformatics Analyses to Investigate Their Potential Roles

SecA is an evolutionarily conserved protein that plays an indispensable role in the secretion of proteins across the bacterial cell membrane. Comparative analyses of SecA homologs have identified two large conserved signature inserts (CSIs) that are unique characteristics of thermophilic bacteria. A 50 aa conserved insert in SecA is exclusively present in the SecA homologs from the orders Thermotogales and Aquificales, while a 76 aa insert in SecA is specific for the order Thermales and Hydrogenibacillus schlegelii. Phylogenetic analyses on SecA sequences show that the shared presence of these CSIs in unrelated groups of thermophiles is not due to lateral gene transfers, but instead these large CSIs have likely originated independently in these lineages due to their advantageous function. Both of these CSIs are located in SecA protein in a surface exposed region within the ATPase domain. To gain insights into the functional significance of the 50 aa CSI in SecA, molecular dynamics (MD) simulations were performed at two different temperatures using ADP-bound SecA from Thermotoga maritima. These analyses have identified a conserved network of water molecules near the 50 aa insert in which the Glu185 residue from the CSI is found to play a key role towards stabilizing these interactions. The results provide evidence for the possible role of the 50 aa CSI in stabilizing the binding interaction of ADP/ATP, which is required for SecA function. Additionally, the surface-exposed CSIs in SecA, due to their potential to make novel protein-protein interactions, could also contribute to the thermostability of SecA from thermophilic bacteria.

GC bias lead to increased small amino acids and random coils of proteins in cold-water fishes

BACKGROUND

Temperature adaptation of biological molecules is fundamental in evolutionary studies but remains unsolved. Fishes living in cold water are adapted to low temperatures through adaptive modification of their biological molecules, which enables their functioning in extreme cold. To study nucleotide and amino acid preference in cold-water fishes, we investigated the substitution asymmetry of codons and amino acids in protein-coding DNA sequences between cold-water fishes and tropical fishes., The former includes two Antarctic fishes, Dissostichus mawsoni (Antarctic toothfish), Gymnodraco acuticeps (Antarctic dragonfish), and two temperate fishes, Gadus morhua (Atlantic cod) and Gasterosteus aculeatus (stickleback), and the latter includes three tropical fishes, including Danio rerio (zebrafish), Oreochromis niloticus (Nile tilapia) and Xiphophorus maculatus (Platyfish).

RESULTS

Cold-water fishes showed preference for Guanines and cytosines (GCs) in both synonymous and nonsynonymous codon substitution when compared with tropical fishes. Amino acids coded by GC-rich codons are favored in the temperate fishes, while those coded by AT-rich codons are disfavored. Similar trends were discovered in Antarctic fishes but were statistically weaker. The preference of GC rich codons in nonsynonymous substitution tends to increase ratio of small amino acid in proteins, which was demonstrated by biased small amino acid substitutions in the cold-water species when compared with the tropical species, especially in the temperate species. Prediction and comparison of secondary structure of the proteomes showed that frequency of random coils are significantly larger in the cold-water fish proteomes than those of the tropical fishes.

CONCLUSIONS

Our results suggested that natural selection in cold temperature might favor biased GC content in the coding DNA sequences, which lead to increased frequency of small amino acids and consequently increased random coils in the proteomes of cold-water fishes.

Characterization of the Dihydroorotase from Methanococcus jannaschii

The gene that codes for the putative dihydroorotase in the hyperthermophilic archaeon Methanococcus jannaschii was subcloned in pET-21a and expressed in Escherichia coli. A purification protocol was devised. The purity of the protein was evaluated by SDS-PAGE and the protein was confirmed by sequencing using LC-MS. The calculated molecular mass is 48104 Da. SEC-LS suggested that the protein is a monomer in solution. ICP-MS showed that there are two Zn ions per monomer. Kinetic analysis of the recombinant protein gave hyperbolic kinetics with Vmax = 12.2 µmol/min/mg and Km = 0.14 mM at 25 °C. Furthermore the activity of the protein increased with temperature consistent with the hyperthermophilic nature of the organism. A homology model was constructed using the mesophilic Bacillus anthracis protein as the template. Residues known to be critical for Zn and substrate binding were conserved. The activity of the enzyme at 85 and 90 °C was found to be relatively constant over 160 min and this correlates with the temperature of optimal growth of the organism of 85 °C. The amino acid sequences and structures of the two proteins were compared and this gave insight into some of the factors that may confer thermostability-more Lys and Ile, fewer Ala, Thr, Gln and Gly residues, and shorter N- and C-termini. Additional and better insight into the thermostabilization strategies adopted by this enzyme will be provided when its crystal structure is determined.

The Code of Silence: Widespread Associations Between Synonymous Codon Biases and Gene Function

Some mutations in gene coding regions exchange one synonymous codon for another, and thus do not alter the amino acid sequence of the encoded protein. Even though they are often called 'silent,' these mutations may exhibit a plethora of effects on the living cell. Therefore, they are often selected during evolution, causing synonymous codon usage biases in genomes. Comparative analyses of bacterial, archaeal, fungal, and human cancer genomes have found many links between a gene's biological role and the accrual of synonymous mutations during evolution. In particular, highly expressed genes in certain functional categories are enriched with optimal codons, which are decoded by the abundant tRNAs, thus enhancing the speed and accuracy of the translating ribosome. The set of genes exhibiting codon adaptation differs between genomes, and these differences show robust associations to organismal phenotypes. In addition to selection for translation efficiency, other distinct codon bias patterns have been found in: amino acid starvation genes, cyclically expressed genes, tissue-specific genes in animals and plants, oxidative stress response genes, cellular differentiation genes, and oncogenes. In addition, genomes of organisms harboring tRNA modifications exhibit particular codon preferences. The evolutionary trace of codon bias patterns across orthologous genes may be examined to learn about a gene's relevance to various phenotypes, or, more generally, its function in the cell.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

Overlapping genes: a new strategy of thermophilic stress tolerance in prokaryotes

Overlapping genes (OGs) draw the focus of recent day's research. However, the significance of OGs in prokaryotic genomes remained unexplored. As an adaptation to high temperature, thermophiles were shown to eliminate their intergenic regions. Therefore, it could be possible that prokaryotes would increase their OG content to adapt to high temperature. To test this hypothesis, we carried out a comparative study on OG frequency of 256 prokaryotic genomes comprising both thermophiles and non-thermophiles. It was found that thermophiles exhibit higher frequency of overlapping genes than non-thermophiles. Moreover, overlap frequency was found to correlate with optimal growth temperature (OGT) in prokaryotes. Long overlap frequency was found to hold a positive correlation with OGT resulting in an abundance of long overlaps in thermophiles compared to non-thermophiles. On the other hand, short overlap (1-4 nucleotides) frequency (SOF) did not yield any direct correlation with OGT. However, the correlation of SOF with CAIavg (extent of variation of codon usage bias measured as the mean of codon adaptation index of all genes in a given genome) and IG% (proportion of intergenic regions) indicate that they might upregulate the aforementioned factors (CAIavg and IG%) which are already known to be vital forces for thermophilic adaptation. From these evidences, we propose that the OG content bears a strong link to thermophily. Long overlaps are important for their genome compaction and short overlaps are important to uphold high CAIavg. Our findings will surely help in better understanding of the significance of overlapping gene content in prokaryotic genomes.

Reconstructed ancestral Myo-inositol-3-phosphate synthases indicate that ancestors of the Thermococcales and Thermotoga species were more thermophilic than their descendants

The bacterial genomes of Thermotoga species show evidence of significant interdomain horizontal gene transfer from the Archaea. Members of this genus acquired many genes from the Thermococcales, which grow at higher temperatures than Thermotoga species. In order to study the functional history of an interdomain horizontally acquired gene we used ancestral sequence reconstruction to examine the thermal characteristics of reconstructed ancestral proteins of the Thermotoga lineage and its archaeal donors. Several ancestral sequence reconstruction methods were used to determine the possible sequences of the ancestral Thermotoga and Archaea myo-inositol-3-phosphate synthase (MIPS). These sequences were predicted to be more thermostable than the extant proteins using an established sequence composition method. We verified these computational predictions by measuring the activities and thermostabilities of purified proteins from the Thermotoga and the Thermococcales species, and eight ancestral reconstructed proteins. We found that the ancestral proteins from both the archaeal donor and the Thermotoga most recent common ancestor recipient were more thermostable than their descendants. We show that there is a correlation between the thermostability of MIPS protein and the optimal growth temperature (OGT) of its host, which suggests that the OGT of the ancestors of these species of Archaea and the Thermotoga grew at higher OGTs than their descendants.

Identification and analysis of conserved pockets on protein surfaces

Background

The interaction between proteins and ligands occurs at pockets that are often lined by conserved amino acids. These pockets can represent the targets for low molecular weight drugs. In order to make the research for new medicines as productive as possible, it is necessary to exploit "in silico" techniques, high throughput and fragment-based screenings that require the identification of druggable pockets on the surface of proteins, which may or may not correspond to active sites.

Results

We developed a tool to evaluate the conservation of each pocket detected on the protein surface by CastP. This tool was named DrosteP because it recursively searches for optimal input sequences to be used to calculate conservation. DrosteP uses a descriptor of statistical significance, Poisson p-value, as a target to optimize the choice of input sequences. To benchmark DrosteP we used monomeric or homodimer human proteins with known 3D-structure whose active site had been annotated in UniProt. DrosteP is able to detect the active site with high accuracy because in 81% of the cases it coincides with the most conserved pocket. Comparing DrosteP with analogous programs is difficult because the outputs are different. Nonetheless we could assess the efficacy of the recursive algorithm in the identification of active site pockets by calculating conservation with the same input sequences used by other programs.

We analyzed the amino-acid composition of conserved pockets identified by DrosteP and we found that it differs significantly from the amino-acid composition of non conserved pockets.

Conclusions

Several methods for predicting ligand binding sites on protein surfaces, that combine 3D-structure and evolutionary sequence conservation, have been proposed. Any method relying on conservation mainly depends on the choice of the input sequences. DrosteP chooses how deeply distant homologs must be collected to evaluate conservation and thus optimizes the identification of active site pockets. Moreover it recognizes conserved pockets other than those coinciding with the sites annotated in UniProt that might represent useful druggable sites. The distinctive amino-acid composition of conserved pockets provides useful hints on the fundamental principles underlying protein-ligand interaction.

Availability

http://www.icb.cnr.it/project/drosteppy/

Environment specific substitution tables improve membrane protein alignment

MOTIVATION

Membrane proteins are both abundant and important in cells, but the small number of solved structures restricts our understanding of them. Here we consider whether membrane proteins undergo different substitutions from their soluble counterparts and whether these can be used to improve membrane protein alignments, and therefore improve prediction of their structure.

RESULTS

We construct substitution tables for different environments within membrane proteins. As data is scarce, we develop a general metric to assess the quality of these asymmetric tables. Membrane proteins show markedly different substitution preferences from soluble proteins. For example, substitution preferences in lipid tail-contacting parts of membrane proteins are found to be distinct from all environments in soluble proteins, including buried residues. A principal component analysis of the tables identifies the greatest variation in substitution preferences to be due to changes in hydrophobicity; the second largest variation relates to secondary structure. We demonstrate the use of our tables in pairwise sequence-to-structure alignments (also known as 'threading') of membrane proteins using the FUGUE alignment program. On average, in the 10-25% sequence identity range, alignments are improved by 28 correctly aligned residues compared with alignments made using FUGUE's default substitution tables. Our alignments also lead to improved structural models.

AVAILABILITY

Substitution tables are available at: http://www.stats.ox.ac.uk/proteins/resources.

Comparison of the structural basis for thermal stability between archaeal and bacterial proteins

In this study, the structural basis for thermal stability in archaeal and bacterial proteins was investigated. There were many common factors that confer resistance to high temperature in both archaeal and bacterial proteins. These factors include increases in the Lys content, the bends and blanks of secondary structure, the Glu content of salt bridge; decreases in the number of main-side chain hydrogen bond and exposed surface area, and changes in the bends and blanks of amino acids. Certainly, the utilization of charged amino acids to form salt bridges is a primary factor. In both heat-resistant archaeal and bacterial proteins, most Glu and Asp participate in the formation of salt bridges. Other factors may influence either archaeal or bacterial protein thermostability, which includes the more frequent occurrence of shorter 3(10)-helices and increased hydrophobicity in heat-resistant archaeal proteins. However, there were increases in average helix length, the Glu content in salt bridges, temperature factors and decreases in the number of main-side chain hydrogen bonds, uncharged-uncharged hydrogen bonds, hydrophobicity, and buried and exposed polar surface area in heat-resistant bacterial proteins. Evidently, there are few similarities and many disparities between the heat-resistant mechanisms of archaeal and bacterial proteins.

Using knowledge of protein structural constraints to predict the evolution of HIV-1

The high levels of sequence diversity and rapid rates of evolution of HIV-1 represent the main challenges for developing effective therapies. However, there are constraints imposed by the three-dimensional protein structure that affect the sequence space accessible to the evolution of HIV-1. Here, we present a strategy for predicting the set of possible amino acid replacements in HIV. Our approach is based on the identification of likely amino acid changes in the context of these structural constraints using environment-specific substitution matrices as well as considering the physical constraints imposed by local structure. Assessment of the power of various published algorithms in predicting the evolution of HIV-1 Gag P17 shows that it is possible to use these methods to make accurate predictions of the sequence diversity. Our own method, SubFit, uses knowledge of local structural constraints; it achieves similar prediction success with the best-performing methods. We also show that erroneous predictions are largely due to infrequently occurring amino acids that will probably have severe fitness costs for the protein. Future improvements; for example, incorporating covariation and immunological constraints will permit more reliable prediction of viral evolution.

A structural dissection of amino acid substitutions in helical transmembrane proteins

The evolution of protein folds is under strong constraints from their surrounding environment. Although folding in water-soluble proteins is driven primarily by hydrophobic forces, the nature of the forces that determine the folding and stability of transmembrane proteins are still not fully understood. Furthermore, the chemically heterogeneous lipid bilayer has a non-uniform effect on protein structure. In this article, we attempt to get an insight into the nature of this effect by examining the impact of various types of local structure environment on amino acid substitution, based on alignments of high-resolution structures of polytopic helical transmembrane proteins combined with sequences of close homologs. Compared to globular proteins, burying amino acid sidechains, especially hydrophilic ones, led to a lower increase in conservation in both the lipid-water interface region and the hydrocarbon core region. This observation is due to surface residues in HTM proteins especially in the HC region being relatively highly conserved, suggesting higher evolutionary constraints from their specific interactions with the surrounding lipid molecules. Polar and small residues, particularly Pro and Gly, show a noticeable increase in conservation as they are positioned more towards the centre of the membrane, which is consistent with their recognized key roles in structural stability. In addition, the examination of hydrogen bonds in the membrane environment identified some exposed hydrophilic residues being better conserved when not hydrogen-bonded to other residues, supporting the importance of lipid-protein sidechain interactions. The conclusions presented in this study highlight the distinct features of substitution matrices that take into account the membrane environment, and their potential role in improving sequence-structure alignments of transmembrane proteins.

Amino acid composition in endothermic vertebrates is biased in the same direction as in thermophilic prokaryotes

Background

Among bacteria and archaea, amino acid usage is correlated with habitat temperatures. In particular, protein surfaces in species thriving at higher temperatures appear to be enriched in amino acids that stabilize protein structure and depleted in amino acids that decrease thermostability. Does this observation reflect a causal relationship, or could the apparent trend be caused by phylogenetic relatedness among sampled organisms living at different temperatures? And do proteins from endothermic and exothermic vertebrates show similar differences?

Results

We find that the observed correlations between the frequencies of individual amino acids and prokaryotic habitat temperature are strongly influenced by evolutionary relatedness between the species analysed; however, a proteome-wide bias towards increased thermostability remains after controlling for phylogeny. Do eukaryotes show similar effects of thermal adaptation? A small shift of amino acid usage in the expected direction is observed in endothermic ('warm-blooded') mammals and chicken compared to ectothermic ('cold-blooded') vertebrates with lower body temperatures; this shift is not simply explained by nucleotide usage biases.

Conclusion

Protein homologs operating at different temperatures have different amino acid composition, both in prokaryotes and in vertebrates. Thus, during the transition from ectothermic to endothermic life styles, the ancestors of mammals and of birds may have experienced weak genome-wide positive selection to increase the thermostability of their proteins.

Reduction in structural disorder and functional complexity in the thermal adaptation of prokaryotes

Genomic correlates of evolutionary adaptation to very low or very high optimal growth temperature (OGT) values have been the subject of many studies. Whereas these provided a protein-structural rationale of the activity and stability of globular proteins/enzymes, the point has been neglected that adaptation to extreme temperatures could also have resulted from an increased use of intrinsically disordered proteins (IDPs), which are resistant to these conditions in vitro. Contrary to these expectations, we found a conspicuously low level of structural disorder in bacteria of very high (and very low) OGT values. This paucity of disorder does not reflect phylogenetic relatedness, i.e. it is a result of genuine adaptation to extreme conditions. Because intrinsic disorder correlates with important regulatory functions, we asked how these bacteria could exist without IDPs by studying transcription factors, known to harbor a lot of function-related intrinsic disorder. Hyperthermophiles have much less transcription factors, which have reduced disorder compared to their mesophilic counterparts. On the other hand, we found by systematic categorization of proteins with long disordered regions that there are certain functions, such as translation and ribosome biogenesis that depend on structural disorder even in hyperthermophiles. In all, our observations suggest that adaptation to extreme conditions is achieved by a significant functional simplification, apparent at both the level of the genome and individual genes/proteins.

Temperature adaptation at homologous sites in proteins from nine thermophile-mesophile species pairs

Whether particular amino acids are favored by selection at high temperatures over others has long been an open question in protein evolution. One way to approach this question is to compare homologous sites in proteins from one thermophile and a closely related mesophile; asymmetrical substitution patterns have been taken as evidence for selection favoring certain amino acids over others. However, most pairs of prokaryotic species that differ in optimum temperature also differ in genome-wide GC content, and amino acid content is known to be associated with GC content. Here, I compare homologous sites in nine thermophilic prokaryotes and their mesophilic relatives, all with complete published genome sequences. After adjusting for the effects of differing GC content with logistic regression, 139 of the 190 pairs of amino acids show significant substitutional asymmetry, evidence of widespread adaptive amino acid substitution. The patterns are fairly consistent across the nine pairs of species (after taking the effects of differing GC content into account), suggesting that much of the asymmetry results from adaptation to temperature. Some amino acids in some species pairs deviate from the overall pattern in ways indicating that adaptation to other environmental or physiological differences between the species may also play a role. The property that is best correlated with the patterns of substitutional asymmetry is transfer free energy, a measure of hydrophobicity, with more hydrophobic amino acids favored at higher temperatures. The correlation of asymmetry and hydrophobicity is fairly weak, suggesting that other properties may also be important.

Evolution and thermodynamics of the slow unfolding of hyperstable monomeric proteins

BACKGROUND

The unfolding speed of some hyperthermophilic proteins is dramatically lower than that of their mesostable homologs. Ribonuclease HII from the hyperthermophilic archaeon Thermococcus kodakaraensis (Tk-RNase HII) is stabilized by its remarkably slow unfolding rate, whereas RNase HI from the thermophilic bacterium Thermus thermophilus (Tt-RNase HI) unfolds rapidly, comparable with to that of RNase HI from Escherichia coli (Ec-RNase HI).

RESULTS

To clarify whether the difference in the unfolding rate is due to differences in the types of RNase H or differences in proteins from archaea and bacteria, we examined the equilibrium stability and unfolding reaction of RNases HII from the hyperthermophilic bacteria Thermotoga maritima (Tm-RNase HII) and Aquifex aeolicus (Aa-RNase HII) and RNase HI from the hyperthermophilic archaeon Sulfolobus tokodaii (Sto-RNase HI). These proteins from hyperthermophiles are more stable than Ec-RNase HI over all the temperature ranges examined. The observed unfolding speeds of all hyperstable proteins at the different denaturant concentrations studied are much lower than those of Ec-RNase HI, which is in accordance with the familiar slow unfolding of hyperstable proteins. However, the unfolding rate constants of these RNases H in water are dispersed, and the unfolding rate constant of thermophilic archaeal proteins is lower than that of thermophilic bacterial proteins.

CONCLUSIONS

These results suggest that the nature of slow unfolding of thermophilic proteins is determined by the evolutionary history of the organisms involved. The unfolding rate constants in water are related to the amount of buried hydrophobic residues in the tertiary structure.

Translational selection is ubiquitous in prokaryotes

Codon usage bias in prokaryotic genomes is largely a consequence of background substitution patterns in DNA, but highly expressed genes may show a preference towards codons that enable more efficient and/or accurate translation. We introduce a novel approach based on supervised machine learning that detects effects of translational selection on genes, while controlling for local variation in nucleotide substitution patterns represented as sequence composition of intergenic DNA. A cornerstone of our method is a Random Forest classifier that outperformed previous distance measure-based approaches, such as the codon adaptation index, in the task of discerning the (highly expressed) ribosomal protein genes by their codon frequencies. Unlike previous reports, we show evidence that translational selection in prokaryotes is practically universal: in 460 of 461 examined microbial genomes, we find that a subset of genes shows a higher codon usage similarity to the ribosomal proteins than would be expected from the local sequence composition. These genes constitute a substantial part of the genome—between 5% and 33%, depending on genome size—while also exhibiting higher experimentally measured mRNA abundances and tending toward codons that match tRNA anticodons by canonical base pairing. Certain gene functional categories are generally enriched with, or depleted of codon-optimized genes, the trends of enrichment/depletion being conserved between Archaea and Bacteria. Prominent exceptions from these trends might indicate genes with alternative physiological roles; we speculate on specific examples related to detoxication of oxygen radicals and ammonia and to possible misannotations of asparaginyl–tRNA synthetases. Since the presence of codon optimizations on genes is a valid proxy for expression levels in fully sequenced genomes, we provide an example of an “adaptome” by highlighting gene functions with expression levels elevated specifically in thermophilic Bacteria and Archaea.

Synonymous codons are not equally common in genomes. The main causes of unequal codon usage are varying nucleotide substitution patterns, as manifested in the wide range of genomic nucleotide compositions. However, since the first E. coli and yeast genes were sequenced, it became evident that there was also a bias towards codons that can be translated to protein faster and more accurately. This bias was stronger in highly expressed genes, and its driving force was termed translational selection. Researchers sought for effects of translational selection in microbial genomes as they became available, employing a flurry of mathematical approaches which sometimes led to contradictory conclusions. We introduce a sensitive and accurate machine learning-based methodology and find that highly expressed genes have a recognizable codon usage pattern in almost every bacterial and archaeal genome analyzed, even after accounting for large differences in background nucleotide composition. We also show that the gene functional category has a great bearing on whether that gene is subject to translational selection. Since presence of codon optimizations can be used as a purely sequence-derived proxy for expression levels, we can delineate “adaptomes” by relating predicted gene activity to organisms' phenotypes, which we demonstrate on genomes of temperature-resistant Bacteria and Archaea.

Gaining and losing the thermophilic adaptation in prokaryotes

We studied the evolution of thermophily in prokaryotes using the phylogenetic relationships between 279 bacteria and archaea and their thermophilic amino acid composition signature. Our findings suggest several examples in which the capacity of thermophilic adaptation has been gained or lost over relatively short evolutionary periods throughout the evolution of prokaryotes.