Structure-dependent relationships between growth temperature of prokaryotes and the amino acid frequency in their proteins

Taxonomic identification and temperature stress tolerance mechanisms of Aequorivita marisscotiae sp. nov

The deep sea harbours microorganisms with unique life characteristics and activities due to adaptation to particular environmental conditions, but the limited sample collection and pure culture techniques available constrain the study of deep-sea microorganisms. In this study, strain Ant34-E75 was isolated from Antarctic deep-sea sediment samples and showed the highest 16 S rRNA gene sequence similarity (97.18%) with the strain Aequorivita viscosa 8-1bT. Strain Ant34-E75 is psychrotrophic and can effectively increase the cold tolerance of Chlamydomonas reinhardtii (a model organism). Subsequent transcriptome analysis revealed multiple mechanisms involved in the Ant34-E75 response to temperature stress, and weighted gene co-expression network analysis (WGCNA) showed that the peptidoglycan synthesis pathway was the key component. Overall, this study provides insights into the characteristics of a deep-sea microorganism and elucidates mechanisms of temperature adaptation at the molecular level.

Fever as an evolutionary agent to select immune complexes interfaces

We herein analyzed all available protein–protein interfaces of the immune complexes from the Protein Data Bank whose antigens belong to pathogens or cancers that are modulated by fever in mammalian hosts. We also included, for comparison, protein interfaces from immune complexes that are not significantly modulated by the fever response. We highlight the distribution of amino acids at these viral, bacterial, protozoan and cancer epitopes, and at their corresponding paratopes that belong strictly to monoclonal antibodies. We identify the “hotspots”, i.e. residues that are highly connected at such interfaces, and assess the structural, kinetic and thermodynamic parameters responsible for complex formation. We argue for an evolutionary pressure for the types of residues at these protein interfaces that may explain the role of fever as a selective force for optimizing antibody binding to antigens.

Morphological Phenotypes, Cell Division, and Gene Expression of Escherichia coli under High Concentration of Sodium Sulfate

Sodium and sulfate ions are among the suggested abundant ions on Europa, a moon of Jupiter. In order to investigate the potential habitability of Europa, we study the effects of sodium sulfate (Na₂SO₄) on a non-halophilic bacterium by subjecting Escherichia coli (E. coli) to a wide range of Na₂SO₄ concentrations (0–1.0 m). We discover that, as the concentration of sodium sulfate increases, the biomass doubling time increases and the cell growth is completely inhibited at 1.0 m Na₂SO₄. Furthermore, we find that E. coli exhibits three distinct morphological phenotypes—(i) shortened, (ii) normal, and (iii) elongated/filamented cells at 0.6 m and 0.8 m Na₂SO₄. We have examined the expression of different genes involved in sodium and sulfate transport (nhaA, nhaB, cysZ, sbp), osmotically driven transport of water (aqpZ), sulfate metabolism (cysN), fatty acid production (fabA), and a global transcriptional regulator (osmZ). Our results suggest that the expression of these genes is not affected significantly at high concentrations of sodium sulfate in the exponential growth phase. Using our experimental data and the existing data in the literature, we show that the osmotic pressure difference may play a major role in determining the growth inhibition of E. coli and B. subtilis at high concentrations of salt.

Building a tRNA thermometer to estimate microbial adaptation to temperature

Because ambient temperature affects biochemical reactions, organisms living in extreme temperature conditions adapt protein composition and structure to maintain biochemical functions. While it is not feasible to experimentally determine optimal growth temperature (OGT) for every known microbial species, organisms adapted to different temperatures have measurable differences in DNA, RNA and protein composition that allow OGT prediction from genome sequence alone. In this study, we built a 'tRNA thermometer' model using tRNA sequence to predict OGT. We used sequences from 100 archaea and 683 bacteria species as input to train two Convolutional Neural Network models. The first pairs individual tRNA sequences from different species to predict which comes from a more thermophilic organism, with accuracy ranging from 0.538 to 0.992. The second uses the complete set of tRNAs in a species to predict optimal growth temperature, achieving a maximum ${r^2}$ of 0.86; comparable with other prediction accuracies in the literature despite a significant reduction in the quantity of input data. This model improves on previous OGT prediction models by providing a model with minimum input data requirements, removing laborious feature extraction and data preprocessing steps and widening the scope of valid downstream analyses.

Predicting the optimal growth temperatures of prokaryotes using only genome derived features

MOTIVATION

Optimal growth temperature is a fundamental characteristic of all living organisms. Knowledge of this temperature is central to the study of a prokaryote, the thermal stability and temperature dependent activity of its genes, and the bioprospecting of its genome for thermally adapted proteins. While high throughput sequencing methods have dramatically increased the availability of genomic information, the growth temperatures of the source organisms are often unknown. This limits the study and technological application of these species and their genomes. Here, we present a novel method for the prediction of growth temperatures of prokaryotes using only genomic sequences.

RESULTS

By applying the reverse ecology principle that an organism's genome includes identifiable adaptations to its native environment, we can predict a species' optimal growth temperature with an accuracy of 5.17°C root-mean-square error and a coefficient of determination of 0.835. The accuracy can be further improved for specific taxonomic clades or by excluding psychrophiles. This method provides a valuable tool for the rapid calculation of organism growth temperature when only the genome sequence is known.

AVAILABILITY AND IMPLEMENTATION

Source code, genomes analyzed and features calculated are available at: https://github.com/DavidBSauer/OGT_prediction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

High-throughput quantification of protein structural change reveals potential mechanisms of temperature adaptation in Mytilus mussels

Background

Temperature exerts a strong influence on protein evolution: species living in thermally distinct environments often exhibit adaptive differences in protein structure and function. However, previous research on protein temperature adaptation has focused on small numbers of proteins and on proteins adapted to extreme temperatures. Consequently, less is known about the types and quantity of evolutionary change that occurs to proteins when organisms adapt to small shifts in environmental temperature. In this study, these uncertainties were addressed by developing software that enabled comparison of structural changes associated with temperature adaptation (hydrogen bonding, salt bridge formation, and amino acid use) among large numbers of proteins from warm- and cold-adapted species of marine mussels, Mytilus galloprovincialis and Mytilus trossulus, respectively.

Results

Small differences in habitat temperature that characterize the evolutionary history of Mytilus mussels were sufficient to cause protein structural changes consistent with temperature adaptation. Hydrogen bonds and salt bridges that increase stability and protect against heat-induced denaturation were more abundant in proteins from warm-adapted M. galloprovincialis compared with proteins from cold-adapted M. trossulus. These structural changes were related to deviations in the use of polar and charged amino acids that facilitate formation of hydrogen bonds and salt bridges within proteins, respectively. Enzymes, in particular those within antioxidant and cell death pathways, were over-represented among proteins with the most hydrogen bonds and salt bridges in warm-adapted M. galloprovincialis. Unlike extremophile proteins, temperature adaptation in Mytilus proteins did not involve substantial changes in the number of hydrophobic or large volume amino acids, nor in the content of glycine or proline.

Conclusions

Small shifts in organism temperature tolerance, such as that needed to cope with climate warming, may result from structural and functional changes to a small percentage of the proteome. Proteins in which function is dependent on large conformational change, notably enzymes, may be particularly sensitive to temperature perturbation and represent foci for natural selection. Protein temperature adaptation can occur through different types and frequencies of structural change, and adaptive mechanisms used to cope with small shifts in habitat temperature appear different from mechanisms used to retain protein function at temperature extremes.

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.

Extremozymes from metagenome: Potential applications in food processing

The long-established use of enzymes for food processing and product formulation has resulted in an increased enzyme market compounding to 7.0% annual growth rate. Advancements in molecular biology and recognition that enzymes with specific properties have application for industrial production of infant, baby and functional foods boosted research toward sourcing the genes of microorganisms for enzymes with distinctive properties. In this regard, functional metagenomics for extremozymes has gained attention on the premise that such enzymes can catalyze specific reactions. Hence, metagenomics that can isolate functional genes of unculturable extremophilic microorganisms has expanded attention as a promising tool. Developments in this field of research in relation to food sector are reviewed.

Comparative analysis of microbial diversity in two hot springs of Bakreshwar, West Bengal, India

Various aspects of hot springs at Bakreshwar (Lat. 23°52'48″N; Long. 87°22'40″E) in West Bengal, India have been investigated since the middle of 20th century, but comprehending the complete diversity and the complexity of the microbial population therein has been in the continuing process. Some of these microorganisms are found to have immense industrial importance. Microbes generally exist in milieus of varying complexities and diversities. Attempting the usually employed cultivation-based techniques in experimentation with those microbes had confronted various limitations. To overcome these limitations a strategy based on high-throughput sequencing of 16S rRNA gene amplicon analysis was employed for studying the differential diversity and the detailed nature of microbial population of the two hot springs of Bakreshwar (54 °C & 65 °C). Paired-end libraries of amplified V-3 hyper-variable 16S rDNA fragments from sets of samples that varied in their contents, ranging from a single bacterium to highly complex communities were sequenced. The comparison revealed the differential aspects in the two hot spring waters; the samples at 54 °C showed the bacterial phylum Firmicutes (65.85%) and Synergistetes (27.24%) predominating and those from hot spring water at 65 °C showed the abundance of the phyla Firmicutes (96.10%) and Proteobacteria (3.36%). The presence of Archaea in the hot springs could not be ascertained.

Comparative Metagenomics of Eight Geographically Remote Terrestrial Hot Springs

Hot springs are natural habitats for thermophilic Archaea and Bacteria. In this paper, we present the metagenomic analysis of eight globally distributed terrestrial hot springs from China, Iceland, Italy, Russia, and the USA with a temperature range between 61 and 92 (∘)C and pH between 1.8 and 7. A comparison of the biodiversity and community composition generally showed a decrease in biodiversity with increasing temperature and decreasing pH. Another important factor shaping microbial diversity of the studied sites was the abundance of organic substrates. Several species of the Crenarchaeal order Thermoprotei were detected, whereas no single bacterial species was found in all samples, suggesting a better adaptation of certain archaeal species to different thermophilic environments. Two hot springs show high abundance of Acidithiobacillus, supporting the idea of a true thermophilic Acidithiobacillus species that can thrive in hyperthermophilic environments. Depending on the sample, up to 58 % of sequencing reads could not be assigned to a known phylum, reinforcing the fact that a large number of microorganisms in nature, including those thriving in hot environments remain to be isolated and characterized.

The thermostability of two kinds of recombinant ∆6-fatty acid desaturase with different N-terminal sequence lengths in low temperature

Two recombinant Rhizopus stolonifer ∆6-fatty acid desaturase enzymes with different-length N-termini were cloned and expressed in Saccharomyces cerevisiae strain INVScl: LRsD6D begins with the sequence of the N-terminal of the R. stolonifer ∆6-fatty acid desaturase native, encoding a deduced polypeptide of 459 amino acids (M-S-T-L-D-R-Q-S-I-F-T-I-K-E-L-E-S-I-S-Q-R-I-H-D-G-D-E-E-A-M-K-F), whereas SRsD6D begins with the amino acid sequence of the predicted ORF, encoding a deduced polypeptide of 430 amino acids (M-K-F) and LRsD6D is longer than SRsD6D by 29 amino acids (M-S-T-L-D-R-Q-S-I-F-T-I-K-E-L-E-S-I-S-Q-R-I-H-D-G-D-E-E-A). Bioinformatic analysis characterized the two recombinant ∆6-fatty acid desaturase enzymes with different-length N-termini, including three conserved histidine-rich motifs, hydropathy profile, and a cytochrome b5-like domain in the N-terminus. When the coding sequence was expressed in S. cerevisiae strain INVScl, the coding produced ∆6-fatty acid desaturase activity exhibited by RsD6D, leading to a novel peak corresponding to γ-linolenic acid methyl ester standards, which was detected with the same retention time. The residual activity of LRsD6D was 74 % at 15 °C for 4 h and that of SRsD6D was 43 %. Purified recombinant LRsD6D was more stable than SRsD6D, indicating that the N-terminal extension, containing mostly hydrophobic residues, affected the overall stability of recombinant LRsD6D.

Psychrophilic enzymes: from folding to function and biotechnology

Psychrophiles thriving permanently at near-zero temperatures synthesize cold-active enzymes to sustain their cell cycle. Genome sequences, proteomic, and transcriptomic studies suggest various adaptive features to maintain adequate translation and proper protein folding under cold conditions. Most psychrophilic enzymes optimize a high activity at low temperature at the expense of substrate affinity, therefore reducing the free energy barrier of the transition state. Furthermore, a weak temperature dependence of activity ensures moderate reduction of the catalytic activity in the cold. In these naturally evolved enzymes, the optimization to low temperature activity is reached via destabilization of the structures bearing the active site or by destabilization of the whole molecule. This involves a reduction in the number and strength of all types of weak interactions or the disappearance of stability factors, resulting in improved dynamics of active site residues in the cold. These enzymes are already used in many biotechnological applications requiring high activity at mild temperatures or fast heat-inactivation rate. Several open questions in the field are also highlighted.

Optimization to low temperature activity in psychrophilic enzymes

Psychrophiles, i.e., organisms thriving permanently at near-zero temperatures, synthesize cold-active enzymes to sustain their cell cycle. These enzymes are already used in many biotechnological applications requiring high activity at mild temperatures or fast heat-inactivation rate. Most psychrophilic enzymes optimize a high activity at low temperature at the expense of substrate affinity, therefore reducing the free energy barrier of the transition state. Furthermore, a weak temperature dependence of activity ensures moderate reduction of the catalytic activity in the cold. In these naturally evolved enzymes, the optimization to low temperature activity is reached via destabilization of the structures bearing the active site or by destabilization of the whole molecule. This involves a reduction in the number and strength of all types of weak interactions or the disappearance of stability factors, resulting in improved dynamics of active site residues in the cold. Considering the subtle structural adjustments required for low temperature activity, directed evolution appears to be the most suitable methodology to engineer cold activity in biological catalysts.

Characterization of a zinc-containing alcohol dehydrogenase with stereoselectivity from the hyperthermophilic archaeon Thermococcus guaymasensis

An alcohol dehydrogenase (ADH) from hyperthermophilic archaeon Thermococcus guaymasensis was purified to homogeneity and was found to be a homotetramer with a subunit size of 40 ± 1 kDa. The gene encoding the enzyme was cloned and sequenced; this gene had 1,095 bp, corresponding to 365 amino acids, and showed high sequence homology to zinc-containing ADHs and l-threonine dehydrogenases with binding motifs of catalytic zinc and NADP(+). Metal analyses revealed that this NADP(+)-dependent enzyme contained 0.9 ± 0.03 g-atoms of zinc per subunit. It was a primary-secondary ADH and exhibited a substrate preference for secondary alcohols and corresponding ketones. Particularly, the enzyme with unusual stereoselectivity catalyzed an anti-Prelog reduction of racemic (R/S)-acetoin to (2R,3R)-2,3-butanediol and meso-2,3-butanediol. The optimal pH values for the oxidation and formation of alcohols were 10.5 and 7.5, respectively. Besides being hyperthermostable, the enzyme activity increased as the temperature was elevated up to 95°C. The enzyme was active in the presence of methanol up to 40% (vol/vol) in the assay mixture. The reduction of ketones underwent high efficiency by coupling with excess isopropanol to regenerate NADPH. The kinetic parameters of the enzyme showed that the apparent K(m) values and catalytic efficiency for NADPH were 40 times lower and 5 times higher than those for NADP(+), respectively. The physiological roles of the enzyme were proposed to be in the formation of alcohols such as ethanol or acetoin concomitant to the NADPH oxidation.

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.

Changes in non-core regions stabilise plastocyanin from the thermophilic cyanobacterium Phormidium laminosum

We report a theoretical investigation on the different stabilities of two plastocyanins. The first one belongs to the thermophilic cyanobacterium Phormidium laminosum and the second one belongs to its mesophilic relative Synechocystis sp. These proteins share the same topology and secondary-structure elements; however, the melting temperatures of their oxidised species differ by approximately 15 K. Long-time-scale molecular dynamics simulations, performed at different temperatures, show that the thermophilic protein optimises a set of intramolecular interactions (interstrand hydrogen bonding, salt bridging and hydrophobic clustering) within the region that comprises the strands beta 5 and beta 6, loop L5 and the helix. This region exhibits most of the differences in the primary sequence between the two proteins and, in addition, it is involved in the interaction with known physiological partners. Further work is in progress to unveil the specific structural features responsible for the different thermal stability of the two proteins.

Structural basis for thermostability revealed through the identification and characterization of a highly thermostable phosphotriesterase-like lactonase from Geobacillus stearothermophilus

A new enzyme homologous to phosphotriesterase was identified from the bacterium Geobacillus stearothermophilus (GsP). This enzyme belongs to the amidohydrolase family and possesses the ability to hydrolyze both lactone and organophosphate (OP) compounds, making it a phosphotriesterase-like lactonase (PLL). GsP possesses higher OP-degrading activity than recently characterized PLLs, and it is extremely thermostable. GsP is active up to 100 degrees C with an energy of activation of 8.0 kcal/mol towards ethyl paraoxon, and it can withstand an incubation temperature of 60 degrees C for two days. In an attempt to understand the thermostability of PLLs, the X-ray structure of GsP was determined and compared to those of existing PLLs. Based upon a comparative analysis, a new thermal advantage score and plot was developed and reveals that a number of different factors contribute to the thermostability of PLLs.

A genomic analysis of the archaeal system Ignicoccus hospitalis-Nanoarchaeum equitans

Sequencing of the complete genome of Ignicoccus hospitalis gives insight into its association with another species of Archaea, Nanoarchaeum equitans.

Background

The relationship between the hyperthermophiles Ignicoccus hospitalis and Nanoarchaeum equitans is the only known example of a specific association between two species of Archaea. Little is known about the mechanisms that enable this relationship.

Results

We sequenced the complete genome of I. hospitalis and found it to be the smallest among independent, free-living organisms. A comparative genomic reconstruction suggests that the I. hospitalis lineage has lost most of the genes associated with a heterotrophic metabolism that is characteristic of most of the Crenarchaeota. A streamlined genome is also suggested by a low frequency of paralogs and fragmentation of many operons. However, this process appears to be partially balanced by lateral gene transfer from archaeal and bacterial sources.

Conclusions

A combination of genomic and cellular features suggests highly efficient adaptation to the low energy yield of sulfur-hydrogen respiration and efficient inorganic carbon and nitrogen assimilation. Evidence of lateral gene exchange between N. equitans and I. hospitalis indicates that the relationship has impacted both genomes. This association is the simplest symbiotic system known to date and a unique model for studying mechanisms of interspecific relationships at the genomic and metabolic levels.

Teaching noncovalent interactions using protein molecular evolution

Noncovalent interactions and physicochemical properties of amino acids are important topics in biochemistry courses. Here, we present a computational laboratory where the capacity of each of the 20 amino acids to maintain different noncovalent interactions are used to investigate the stabilizing forces in a set of proteins coming from organisms adapted to different environments. Using protein sequence and structure information it is possible to evaluate the noncovalent contributions to the stabilization of a given protein fold. As a case study, we use the protein lumazine synthase from three different organisms adapted to live in extreme temperatures: one psychrophilic (optimal growth temperature, 0-20 °C), one mesophilic (optimal growth temperature, 20-50 °C), and one thermophilic (optimal growth temperature, 80-110 °C). We found that this computational laboratory for biochemistry and molecular biology courses enhances student amino acid noncovalent interaction understanding and how these interactions are involved in protein stability.

Gene cloning of cold-adapted isocitrate lyase from a psychrophilic bacterium, Colwellia psychrerythraea, and analysis of amino acid residues involved in cold adaptation of this enzyme

The gene (icl) encoding cold-adapted isocitrate lyase (ICL) of a psychrophilic bacterium, Colwellia psychrerythraea, was cloned and sequenced. Open reading frame of the gene was 1,587 bp in length and corresponded to a polypeptide composed of 528 amino acids. The deduced amino acid sequence showed high homology with that of cold-adapted ICL from other psychrophilic bacterium, C. maris (88% identity), but the sequential homology with that of the Escherichia coli ICL was low (28% identity). Primer extension analysis revealed that transcriptional start site for the C. psychrerythraea icl gene was guanine, located at 87 bases upstream of translational initiation codon. The expression of this gene in the cells of an E. coli mutant defective in ICL was induced by not only low temperature but also acetate. However, cis-acting elements for cold-inducible expression known in the several other bacterial genes were absent in the promoter region of the C. psychrerythraea icl gene. The substitution of Ala214 for Ser in the C. psychrerythraea ICL introduced by point mutation resulted in the increased thermostability and lowering of the specific activity at low temperature, indicating that Ala214 is important for psychrophilic properties of this enzyme.