key: cord-0825091-fxvq5xqy
authors: Barik, Sailen
title: Evolution of Protein Structure and Stability in Global Warming
date: 2020-12-18
journal: Int J Mol Sci
DOI: 10.3390/ijms21249662
sha: 5580eceec6dc0ba4de62192247a4cfd3d5016c0e
doc_id: 825091
cord_uid: fxvq5xqy

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.

The rapidly accelerating pace of climate change is most pronounced as a rapid rise in global temperature and the serious harm it is inflicting on the environment and life, threatening our very existence on Earth [1] [2] [3] [4] [5] [6] . The "global warming", as the temperature change is commonly called, has drawn worldwide attention for decades and is rapidly approaching the "tipping point" [7] [8] [9] ; unfortunately, however, multiple attempts to generate a universal consensus to combat it have also failed, mainly for short-term economic and political priorities [10] [11] [12] . Most experts now believe that global warming will soon become irreversible, or may already have, which means that it is, or soon will be, too late for preventive measures to reverse the process. In either case, it is safe to assume that global warming is here to stay, and, therefore, we must find ways to survive it. Considering this premise, it is important to understand how high temperatures affect living beings, and how they cope with it.

The majority of organisms, including all vertebrates, such as humans, grows in moderate ambient temperatures and are called "mesophilic" (less often, "neutrophilic"). The heat-tolerant organisms are classified as either "thermoresistant' or 'thermophilic". The thermoresistant organisms require ambient temperatures for optimal growth, but can survive higher temperatures to some extent. The thermophilic organisms, in contrast, have an obligatory requirement of high temperatures for growth; however, they may retain viability at lower temperatures. Some investigators also use the terms "extreme thermophile" and/or "hyperthermophile" in describing those that grow at significantly higher temperatures, some approaching that of boiling water. In general, albeit somewhat arbitrary, the guideline range of optimal growth temperature is as follows: mesophile, 20-45 • C; simple thermophile, 50-64 • C; extreme thermophile, 65-79 • C; and hyperthermophile, ≥80 • C [13, 14] . All high-temperature organisms belong to the larger family of "extremophiles", defined as organisms residing in environmental conditions that are extreme in some respect, such as extreme cold (Arctic and Antarctic), alkaline or acidic pH (hot springs and mudpots), high salinity (oceans and various lakes) or pressure (deep sea), and desiccation (deserts). Reciprocally, the residents of extreme cold, the psychrophiles, as they are called, are defined as living below 15 • C; while they present an interesting contrast to the thermophiles in terms of both habitat and molecular mechanism, they are beyond the scope of this review. An organism can also be a crossover extremophile or "polyextremophile", surviving and growing under multiple extreme conditions, such as cold and pressure, or heat and acidity [15] . In fact, oftentimes, the habitats of the extreme thermophiles combine high temperature with some of the other extreme conditions, both physical and chemical, aerobic or anaerobic [16] . The different classes of extremophiles may share some molecular mechanisms, such as those governing protein folding, but this review focusses specifically on high heat, which is most pertinent to global warming and its effect on terrestrial species, including humans. Lastly, unless otherwise mentioned, all types of thermoresistant and thermophilic organisms are addressed synonymously, or simply referred to as thermophiles, for the sake of brevity.

As mentioned above, unicellular protobacteria and archaea constitute a large fraction of thermophiles, but the hyperthermophilic organisms are mostly archaea [17] . Very few eukaryotes are thermophiles; however, some fungi grow at 50-55 • C [18] . Nevertheless, a significant amount of information is available on the heat shock response of eukaryotes. In fact, it is fair to assume that every organism is exposed to higher-than-ambient temperature at some point in its life, and, therefore, the response to heat is likely to be universal, and subsets of the response genes are phylogenetically conserved. Early evidence documented the survival benefit of heat shock to heat tolerance, i.e., when the heat shock response genes are induced by high but nonlethal temperature, the cells can survive a subsequent exposure to lethal temperature, which can be viewed as "stress immunization". The technical advantage of the heat shock response is that it can be applied to mesophilic organisms that are easily grown in the laboratory and then shifted to desired higher temperature for defined periods, providing a detailed kinetics of gene induction and potentially, a glimpse of the regulatory circuitry that connects them. Natural thermophiles, on the other hand, are often polyextremophiles that must be grown not only at high temperature but also in specially reconstituted media and gas mixture. Ergo, the more extreme the natural habitat of the organism is, the harder is it to culture. The Archaea are particularly difficult to use due to their relatively exotic habitats and rudimentary genetics, and possibly unique genomic features (see later).

Heat denatures a protein by destroying the noncovalent bonds and altering the protein's higher-order structure [19] . Evidently, an understanding of the effect of heat on a protein belongs to the realm of protein folding, which is a highly complex process and still not fully retractable-even less predictable-despite remarkable theoretical and technological achievements [19] [20] [21] [22] . Although folding and refolding can be achieved with purified proteins in vitro, the role of the intracellular environment and interacting proteins, such as chaperons [23, 24] , which are relevant for a real-life survival, cannot be ignored. Many proteins also interact with RNA, both large and small, and a large variety of small molecules can act as regulators and modifiers [25, 26] . Given these caveats, we should at least focus on studies comparing orthologous proteins in closely related species. For example, a major survival mechanism adopted by some bacteria (such as Bacillus sp.) is the formation of spores at high temperature, which occurs via a unique but well-researched gene expression program [27] . However, sporulation is a monopoly of select organisms and is of little use to others. In the higher eukaryotes, such as humans, live animal experiments are ruled out, and, therefore, results are obtained from in vitro studies of cell lines and proteins. In this review, relevant organisms and proteins are adhered to as much as possible, although the molecular principles of protein folding, which are fundamental to all life, were gathered from diverse organisms and summarized.

Early studies showed that enzymes, purified from thermophilic organisms, exhibit high temperature optima, comparable to the optimal growth temperature of the parent organism [28, 29] . Moreover, a majority of the hyperthermophilic enzymes, recombinantly expressed in mesophilic hosts, retain all of their native properties, such as catalytic parameters [30] [31] [32] and three-dimensional structure [33] [34] [35] [36] . These results reinforced the fundamental tenet that thermophilia, likely a function of the higher-order structure, is encoded in the primary structure (amino acid sequence) of a protein, and opened the door for analysis of thermoresistance at the protein level. In parallel, much of our knowledge of thermosurvival also derived from studies of heat shock response (often synonymously called "stress response") [37] [38] [39] , mentioned earlier, in which organisms and cultured cells, shifted to higher temperature, induce the expression of a large set of genes for survival. Evidently, these genes can play a critical role in protection from global warming, as well, and, therefore, they are explored in this review. Other areas include mutational analysis of thermostability of specific proteins acquired through natural selection, site-directed mutagenesis and reiteration, and comparison between orthologs of differing thermotolerance. When multiple orthologs are available, multiple sequence alignment points to specific amino acids that may have a role; these are then interrogated by mutational analysis to identify the residues required for thermophilia, since some differences may be due to evolutionary selection for other functions. As suggested previously [40, 41] , all mutations found in a thermoresistant protein (compared to its mesophilic ortholog) are unlikely to have been selected for thermoresistance, i.e., some may be unintended consequence of other selections [42, 43] . Proteins with assayable function, such as enzymes, are useful in this regard, as they can allow distinction of functional residues from thermotolerant ones, and also ensure that functionality is not compromised to achieve thermal stability in reverse genetic mutational studies [41, 42] .

Even though global warming looms in the future, geobiological records from the past should facilitate our understanding of the mechanism and evolution of thermostability. A technical problem, however, is the scarcity of ancient biological material that is well-preserved, undegraded, and uncontaminated. Counter-intuitively, this is also true of cold-preserved material, such as the permafrost proxies in the Arctic that remain the least characterized [44] . An alternative, bioinformatic approach, named "ancestral sequence reconstruction" (ASR), has received significant attention [45] [46] [47] [48] [49] . It is based on the well-founded notion that the earth was extremely hot early in creation and gradually cooled down to allow the development of life. By inference, the earliest life forms must have been hyperthermophiles [50, 51] , which then evolved into simple thermophiles or mesophiles with the gradual cooling of the climate. This "hot-start" hypothesis [52] is supported by environmental evidence [53] , as well as by the location of thermophilic organisms near the roots of the universal tree of life [54] [55] [56] , such as the bacteria Thermotoga, Aquifex, and Hydrogenbacter, and many Archaea, such as Thermoproteus and Thermococcus. Thus, it is conceivable that extrapolation of an extant protein sequence, a few gigayears into the past, using a combination of molecular phylogenetic strategy, computer simulations, and maximum-likelihood methods, will lead to its ancient thermophilic or hyperthermophilic parent that is otherwise intractable, which is the essence of ASR [45] [46] [47] . The beauty of ASR is that it is testable, since the ancestral amino acids can be introduced into the extant proteins by recombinant means, and the resultant protein can be studied. In several such studies, an ASR-derived enzyme indeed exhibited higher thermal stability than the wild-type parent. In a representative example, the ASR-derived sequence of lignin peroxidase, a heme-containing fungal enzyme that has been extensively studied in Phanerochaete chrysosporium, was found to differ from the extant enzyme by eleven residues [57] . Eleven mutant enzymes containing each of these residues were then constructed and expressed in Escherichia coli; when the purified enzymes were assayed in vitro, several of them showed increased thermal stabilities and increased specific activities [57] . After 30 min incubation at 37 • C, the wild-type enzyme retained only 17% of its activity, whereas the single ASR mutants Ala110Gln, Gly198Ser, and Phe254 retained respectively 28.0, 33.5, and 18.4% of their activity. A triple mutant, in which three neighboring residues were all changed to ASR residues to generate the His239Phe/Thr240Leu/Ile241Leu enzyme, showed the highest thermoresistance, retaining 86.2% of its activity. The triple mutant also showed higher activity at its optimum temperature, i.e., specific activity of 68 units/mg at 40 • C, compared to the wild-type enzyme's activity of 23 units/mg at 29 • C. When the authors examined the location of the residues on the crystal structure [57] , residues 110, 254, and the 239/240/241 trio ( Figure 1 ) were all found in helical regions, whereas only 198 was in a flexible loop. As we see later (Section 3), these locations agree with a role of helical stability in thermoresistance; they probably also improved the packing of the hydrophobic core of the peroxidase. The facts that these residues were far from the catalytic pocket, and that a few were located in non-helical regions, suggest that thermostability of a protein is dictated by the overall structure and not by any one residue. each of these residues were then constructed and expressed in Escherichia coli; when the purified enzymes were assayed in vitro, several of them showed increased thermal stabilities and increased specific activities [57] . After 30 min incubation at 37 °C, the wild-type enzyme retained only 17% of its activity, whereas the single ASR mutants Ala110Gln, Gly198Ser, and Phe254 retained respectively 28.0, 33.5, and 18.4% of their activity. A triple mutant, in which three neighboring residues were all changed to ASR residues to generate the His239Phe/Thr240Leu/Ile241Leu enzyme, showed the highest thermoresistance, retaining 86.2% of its activity. The triple mutant also showed higher activity at its optimum temperature, i.e., specific activity of 68 units/mg at 40 °C, compared to the wild-type enzyme's activity of 23 units/mg at 29 °C. When the authors examined the location of the residues on the crystal structure [57] , residues 110, 254, and the 239/240/241 trio ( Figure 1 ) were all found in helical regions, whereas only 198 was in a flexible loop. As we see later (Section 3), these locations agree with a role of helical stability in thermoresistance; they probably also improved the packing of the hydrophobic core of the peroxidase. The facts that these residues were far from the catalytic pocket, and that a few were located in non-helical regions, suggest that thermostability of a protein is dictated by the overall structure and not by any one residue. [57] , and (B) isocitrate dehydrogenase (ICDH) Caldococcus noboribetus, a thermophilic archaeon [58] .

In both examples, only selected ASR-predicted mutations that led to a substantial increase in thermotolerance are shown, for brevity. The residues are colored, as well as named, and the peptide backbone is green. The model is based on the available crystal structures of the respective close homologs, viz lignin peroxidase (PDB 1B85) of Phanerochaete chrysosporium and ICDH of Escherichia coli (PDB 1CW1). Note that several but not all residues are located on helices. The directions of the amino (N) and carboxy (C) termini of the peptide backbone are indicated.

In another study [58] , the isocitrate dehydrogenase (ICDH) enzyme of Caldococcus noboribetus, a thermophilic archaeon that is often found in hot springs and hyperthermal vents, was subjected to ASR. Sequence comparison suggested a role of multiple mutations in thermostability, four of which were Tyr309Ile/Ile310Leu/Ile321Leu/Gly326Ser ( Figure 1 ). These residues were located in multiple structures, including in a helix and in a flexible loop (Figure 1 ), documenting that diverse regions of a protein can contribute to thermostability. The half-inactivation temperature of the mutant was found to be ~93 °C, whereas that of the parent enzyme was 87 °C. Thus, use of ASR pointed to progenotes that were even more thermophilic than the Archaea, consistent with the high heat of the primordial environment. Similar studies with many other enzymes have further confirmed the conclusion that early organisms were indeed endowed with superior thermoresistance, adding confidence to ASR [50, [59] [60] [61] [62] [63] . To the molecular biologist, the ASR studies provide an opportunity to understand the roles of specific amino acids in thermophilia; however, this requires analyses of [57] , and (B) isocitrate dehydrogenase (ICDH) Caldococcus noboribetus, a thermophilic archaeon [58] . In both examples, only selected ASR-predicted mutations that led to a substantial increase in thermotolerance are shown, for brevity. The residues are colored, as well as named, and the peptide backbone is green. The model is based on the available crystal structures of the respective close homologs, viz lignin peroxidase (PDB 1B85) of Phanerochaete chrysosporium and ICDH of Escherichia coli (PDB 1CW1). Note that several but not all residues are located on helices. The directions of the amino (N) and carboxy (C) termini of the peptide backbone are indicated.

In another study [58] , the isocitrate dehydrogenase (ICDH) enzyme of Caldococcus noboribetus, a thermophilic archaeon that is often found in hot springs and hyperthermal vents, was subjected to ASR. Sequence comparison suggested a role of multiple mutations in thermostability, four of which were Tyr309Ile/Ile310Leu/Ile321Leu/Gly326Ser ( Figure 1 ). These residues were located in multiple structures, including in a helix and in a flexible loop (Figure 1 ), documenting that diverse regions of a protein can contribute to thermostability. The half-inactivation temperature of the mutant was found to be~93 • C, whereas that of the parent enzyme was 87 • C. Thus, use of ASR pointed to progenotes that were even more thermophilic than the Archaea, consistent with the high heat of the primordial environment. Similar studies with many other enzymes have further confirmed the conclusion that early organisms were indeed endowed with superior thermoresistance, adding confidence to ASR [50, [59] [60] [61] [62] [63] . To the molecular biologist, the ASR studies provide an opportunity to understand the roles of specific amino acids in thermophilia; however, this requires analyses of additional proteins and determination of their structures. To the technologist, ASR offers a convenient and frugal route to design hyperthermophilic enzymes for industrial use (also see Section 4.2).

As stated, the current review interrogates structural features and molecular interactions of proteins in the perspective of biological survival in global warming. Readers interested in the fundamental physicochemical principles of such interactions may consult several excellent reviews and articles published over the years [64] [65] [66] [67] [68] [69] [70] .

A large number of studies have compared amino acid compositions of the proteome or a specific set of proteins in mesophiles and thermophiles [64, [71] [72] [73] . The observed trend is that the thermophilic proteins favor charged, as well as large, hydrophobic residues, and also aromatic residues, but they disfavor uncharged polar ones. In a study spanning the temperature range from −10 to 110 • C, the amino acid set of Ile, Val, Tyr, Trp, Arg, Glu, and Leu (IVYWREL) was found to have a high correlation coefficient with optimal growth temperature [71] . The same study found a correlation of temperature with the total purine (A + G) content of the coding regions, somewhat matching with the IVYWRL composition.

The altered amino acid composition of thermophiles appears to be related to the altered overall nucleotide composition of the genomic DNA, which co-evolved with the translational machinery to prevent melting of the double helix at the higher temperature [14, 71] . Higher proline content in proteins of bacterial thermophiles, for example, coincides with higher percentage of GC in the coding sequences, since CCN is the codon for Pro. Similarly, the hyperthermophilic archaeal genome is AG-rich, which matches with its codon bias, such as abundance of Thr (ACN codons) and Ser (AGU/C codons) [70] . The amino acids Asn, Gln, Met, and Cys are thermolabile for chemical reasons: While the amides (Asn and Gln) undergo deamidation in high heat, the sulfur-containing amino acids (Met and Cys) suffer from oxidation [74] . In thermophilic proteins, therefore, these amino acids are relatively sparse or internally buried, as we discuss later in the paper. Other studies showed that the archaeal thermophiles use an amino acid pattern slightly different from that of the eubacterial orthologs [75] ; specifically, the archaeal thermostability was achieved by substituting non-charged polar amino acids (such as Gln) with Glu and Lys, and non-polar amino acids with Ile on the protein surface, which signals caution in comparing thermophilia across kingdoms, as mentioned earlier. Overall, codon bias and related translational components affect translational kinetics and in turn, the folding of the nascent proteins; this is a complex, multifactorial process, and thus, its correlation with thermophilia may not be straightforward [71, 76] .

Many of the thermophilia-determining residues are in critical locations in the secondary structure of the protein, several being located inside α-helices [71] . The thermophilic proteins tend to have a reduced α-helix content and deleted loops, while being richer in β-sheets [74, 77, 78] ; however, this rule does not appear to apply to membrane proteins [71] . In general, unless specifically mentioned, the thermophilic signatures reported in published studies and discussed in most reviews [70] refer to soluble globular proteins. It is possible that incorporation of the helices of the membrane proteins into the lipid bilayer makes the mechanism of its thermostabilization different from that of the soluble ones [79, 80] . Finally, for several amino acids, the hyperthermophilic signature differs from that of the simple thermophiles [70] , which can only be attributed to the higher temperature (20-30 • C higher) to which the hyperthermophile cells are subjected to, requiring a higher molecular safeguard to prevent denaturation of the DNA helix and protein structures.

The severity of the thermal effects on a protein hinges on the nature of the constituent amino acids. Several investigators have, therefore, searched for consensus amino acids in thermophilic proteins, and have indeed noted various signatures, in which certain classes of amino acids are either favored or disfavored for a variety of mechanisms [64, 70, 72, 73, [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] , the predominant ones among which are summarized here. It should be mentioned that these observations are generalized, and that exceptions often exist, some of which have been noted.

(a) Charged amino acids (acidic and basic): The charged residues are hydrophilic, and in contrast to the hydrophobic ones, which are located to the center of the protein, are on the surface of the protein, where they promote interaction with the solvent and salt bridges with other moieties. The solvation and ionic interactions facilitate solubility and prevent heat denaturation and precipitation. In fact, an increase in surface ion pairs and a more strongly hydrophobic interior were concluded to be the two strongest contributors to thermostability [70, 78, 87, [91] [92] [93] [94] .

(b) Bulky hydrophobic amino acids (Ile, Leu, Met, and Val): While the charged side chains are exposed to the exterior, it is the opposite with the hydrophobic ones, which are typically located to the interior of proteins through hydrophobic interaction, so that they are shielded from the solvent. This allows for tighter packing of the structure of a globular protein, making it more resistant to unfolding [95] [96] [97] . However, they are more prevalent in the hyperthermophiles than in the simple thermophiles, and their beneficial effect is more pronounced on larger proteins [98] .

(c) Aromatic amino acids (Phe, Tyr, Trp, and His): The electron cloud above and below the plane of the cyclic structure sets the aromatic side chains apart from those of all other residues. They stabilize the core of the protein through interactions with other π electron clouds (another aromatic side chain within atomic distance) and with cations, which thwarts unfolding. However, to what extent this interaction contributes to thermophilia has been a matter of debate [70] . It appears that the simple abundance of aromatic amino acids is not a signature of thermophilia [99] and that the location is also critical; with Trp, for example, much of the observed difference relates to an increased number of Trp residues in α-helices and near helix termini (4.4% in mesophiles vs. 7.9% for thermophile protein orthologs) [84] .

(d) Tryptophan (also see above): Although it is an "aromatic" amino acid, Trp is a unique amino acid in many respects [100] , and it deserves special mention. It is the largest and rarest amino acid in the cell, has the largest π electron cloud in its side chain, and is produced by the most complex and energetically most expensive biosynthetic pathway. Besides Met (AUG), it is the only amino acid that is coded for by a single codon (UGG) in the standard codon repertoire. Thus, it stands to reason that Nature would use Trp sparingly, in only the strategically essential locations.

To test if the helix-terminal location of Trp is unique to thermophilic proteins, I analyzed the helices in a large number of helix-rich domains, namely the tetratricopeptide repeat (TPR) and the pentatricopeptide repeat (PPR) [101] [102] [103] [104] . TPR and PPR are, respectively, 34-and 35-amino acid bi-helical motifs, found in diverse proteins that overwhelmingly belong to mesophiles [105] [106] [107] . Both repeats offer an excellent opportunity for this analysis because of the large number of helices they collectively offer. Moreover, they can be identified, and their secondary structure predicted with confidence, due to the conservation of signature amino acid residues [80] . A total of~23,000 TPR and PPR sequences, representing >48,000 α-helices, were analyzed, and the largest concentration of Trp residues were found near the termini of the helices (Figure 2 ), confirming and extending early observations with different sets [102, 107] . Thus, protection by helix-terminal locations are not unique to constitutive thermophiles, but can also help the mesophiles when temperature rises. It would be interesting to investigate the interacting partners of Trp, which may play critical roles in heat-stability. [102, 103] were visually analyzed for the location of Trp residues, and the numbers were plotted, using Excel. The two classic helices of the tetratricopeptide repeat (TPR) (Top) and pentatricopeptide repeat (PPR) (Bottom) motif are shown as Helix A and Helix B, and their amino acid position numbers are presented on the horizontal axis. The two most abundant Trp residues (Trp4/3 and Trp17/16, in Helix A and Helix B, respectively) are indicated as red bars; note that both are near the N-terminal end of the helices.

(e) Acid amides (Asn and Gln): At physiological pH, deamidation of Asn and Gln is the most common nonenzymatic change in a protein that is promoted by elevated temperatures [98] . Deamidation increases the negative charge of the residue, which could cause structural change, unfolding and irreversible inactivation [108] [109] [110] . In fact, site-directed mutation of a strategic Asn to Asp has been shown to eliminate heat-mediated deamidation and inactivation of a fungal glucoamylase [111] . Asn and Gln are, therefore, highly disfavored in thermophiles, but more so in hyperthermophiles, since deamidation is more pronounced at higher temperatures.

(f) Aliphatic hydroxy-amino acids (Ser and Thr): A negative correlation between Ser/Thr content and optimum growth temperature has been noted [84] , i.e., these amino acids are disfavored in thermophiles. Again, this is better correlated in hyperthermophiles than in simple thermophiles, suggesting that the negative effect is marginal when the temperature rise is small. It is also less pronounced in hyperthermophilic archaea [70] . (e) Acid amides (Asn and Gln): At physiological pH, deamidation of Asn and Gln is the most common nonenzymatic change in a protein that is promoted by elevated temperatures [98] . Deamidation increases the negative charge of the residue, which could cause structural change, unfolding and irreversible inactivation [108] [109] [110] . In fact, site-directed mutation of a strategic Asn to Asp has been shown to eliminate heat-mediated deamidation and inactivation of a fungal glucoamylase [111] . Asn and Gln are, therefore, highly disfavored in thermophiles, but more so in hyperthermophiles, since deamidation is more pronounced at higher temperatures.

(f) Aliphatic hydroxy-amino acids (Ser and Thr): A negative correlation between Ser/Thr content and optimum growth temperature has been noted [84] , i.e., these amino acids are disfavored in thermophiles. Again, this is better correlated in hyperthermophiles than in simple thermophiles, suggesting that the negative effect is marginal when the temperature rise is small. It is also less pronounced in hyperthermophilic archaea [70] .

(g) Proline: Pro is the only residue with an alicyclic secondary amine (-NH) in the side chain, to which the α carbon is a direct substituent. The NH loses the H in forming the peptide band, so that no hydrogen is available to donate a hydrogen bond, and, therefore, Pro disrupts an α-helix. On the other hand, the unique cyclic structure of the Pro side chain imparts conformational rigidity to the protein backbone, beneficial for thermotolerance. In other words, the introduction of Pro residues increases the thermal stability of a protein due to decrease in the entropy of the unfolded state [64, 70, 81, 112] . It is, therefore, highly favored in hyperthermophilic organisms. The higher content of Pro, a helix-breaker, is in conformity with the results of an early survey of available 3D structures of thermophilic proteins in the PDB database, which revealed a significant increase in β-strands, compared to their mesophilic homologs [78] .

(h) Disulfide and hydrogen bonds: The Cys-Cys disulfide bonds are covalent and hold a protein structure together, thus boosting thermal stability of proteins in which they occur [113] . For instance, the intramolecular disulfide linkages of the capsid protein of betanodavirus, a deadly pathogen of young marine fish, plays an important role, not only in the assembly of the virion, but also in its thermal stability [114] . However, not all Cys residues in a protein engage in disulfide bonds, since it requires spatial proximity of two Cys residues, and, therefore, occurrence of Cys alone does not predicate thermostability. Indeed, attempts to stabilize proteins by the introduction of extra Cys-Cys disulfide bonds through site-directed mutagenesis of recombinant proteins have yielded inconsistent results, improving the stability and temperature optimum (Tm) for some locations, but not in others [115, 116] . In a representative study [116] , when seven residue-pairs in the alkaline α-amylase from Alkalimonas amylolytica were replaced with cysteines, only three of them (C35-C426, C116-C120, and C436-C480) showed significant enhancement of thermostability (Figure 3) , while the remaining four actually decreased it. (g) Proline: Pro is the only residue with an alicyclic secondary amine (-NH) in the side chain, to which the α carbon is a direct substituent. The NH loses the H in forming the peptide band, so that no hydrogen is available to donate a hydrogen bond, and, therefore, Pro disrupts an α-helix. On the other hand, the unique cyclic structure of the Pro side chain imparts conformational rigidity to the protein backbone, beneficial for thermotolerance. In other words, the introduction of Pro residues increases the thermal stability of a protein due to decrease in the entropy of the unfolded state [64, 70, 81, 112] . It is, therefore, highly favored in hyperthermophilic organisms. The higher content of Pro, a helix-breaker, is in conformity with the results of an early survey of available 3D structures of thermophilic proteins in the PDB database, which revealed a significant increase in β-strands, compared to their mesophilic homologs [78] .

(h) Disulfide and hydrogen bonds: The Cys-Cys disulfide bonds are covalent and hold a protein structure together, thus boosting thermal stability of proteins in which they occur [113] . For instance, the intramolecular disulfide linkages of the capsid protein of betanodavirus, a deadly pathogen of young marine fish, plays an important role, not only in the assembly of the virion, but also in its thermal stability [114] . However, not all Cys residues in a protein engage in disulfide bonds, since it requires spatial proximity of two Cys residues, and, therefore, occurrence of Cys alone does not predicate thermostability. Indeed, attempts to stabilize proteins by the introduction of extra Cys-Cys disulfide bonds through site-directed mutagenesis of recombinant proteins have yielded inconsistent results, improving the stability and temperature optimum (Tm) for some locations, but not in others [115, 116] . In a representative study [116] , when seven residue-pairs in the alkaline α-amylase from Alkalimonas amylolytica were replaced with cysteines, only three of them (C35-C426, C116-C120, and C436-C480) showed significant enhancement of thermostability (Figure 3) , while the remaining four actually decreased it. Site-directed mutagenesis to create Cys-Cys disulfide bonds that generated thermoresistance. The 3D structure of the triple-pair Cys mutant of A. amylolytica α-amylase was generated in silico, as described in Reference [116] , and then subjected to homology-based I-TASSER modeling [117] . The modeling selected α-amylase (3BC9) of Halothermothrix orenii as a high "C-score" homology match among eight closest PDB templates (the other seven being 3BCF, 4JCL, 2DIJ, 2GJP, 6GXV, 4UZU, and 4JCL), independently confirming the previous studies [116] . The PDB structure, thus obtained, is displayed through PyMol (https://pymol.org, Schrödinger, LLC, USA). Each Cys pair, engaged in disulfide bonding and marked with a box, is shown as stick side chains and has a distinctive color: Cys35-Cys426 (blue), Cys116-Cys120 (orange), and Cys436-Cys480 (magenta). For the sake of clarity, portions of the structure and the disulfide bonds are not shown; moreover, no attempt was made to depict the proper rotamer of the side chains, and, therefore, the S atoms of a pair may not be directly facing each other in this representation. The N-and C-terminal directions are also indicated. The structural elements are color-coded: helix, red; strand, yellow; loop, green. A more detailed structure can be found in the original study [116] . Site-directed mutagenesis to create Cys-Cys disulfide bonds that generated thermoresistance. The 3D structure of the triple-pair Cys mutant of A. amylolytica α-amylase was generated in silico, as described in Reference [116] , and then subjected to homology-based I-TASSER modeling [117] . The modeling selected α-amylase (3BC9) of Halothermothrix orenii as a high "C-score" homology match among eight closest PDB templates (the other seven being 3BCF, 4JCL, 2DIJ, 2GJP, 6GXV, 4UZU, and 4JCL), independently confirming the previous studies [116] . The PDB structure, thus obtained, is displayed through PyMol (https://pymol.org, Schrödinger, LLC, USA). Each Cys pair, engaged in disulfide bonding and marked with a box, is shown as stick side chains and has a distinctive color: Cys35-Cys426 (blue), Cys116-Cys120 (orange), and Cys436-Cys480 (magenta). For the sake of clarity, portions of the structure and the disulfide bonds are not shown; moreover, no attempt was made to depict the proper rotamer of the side chains, and, therefore, the S atoms of a pair may not be directly facing each other in this representation. The N-and C-terminal directions are also indicated. The structural elements are color-coded: helix, red; strand, yellow; loop, green. A more detailed structure can be found in the original study [116] .

The Tm of C116-C120 amylase activity was~5 • C higher than that of the wild type (50 to 55 • C), where some other pairs (not shown) had a lower Tm. A triple mutant, C35-C426/C116-C120/C436-C480, showed a six-fold increase in half-life at 60 • C and a 5.2 • C increase in Tm over the wild type. Note that the last pair, viz C436-C480, introduced alone, had little effect on the Tm.

These results clearly demonstrate that rational design of Cys disulfide bonds to improve thermotolerance is possible, but the effect is contextual and likely influenced by other intramolecular interactions; thus, the properties of the mutants need to be experimentally tested. Lastly, the disulfide effect should be most pronounced in oxidizing environments that prevent reduction of the disulfide bond.

In contrast to disulfide bonds, the hydrogen bonds (H-bonds) are weak, but they form between a greater variety of amino acid side chains, as well as with water, and thus provide significant strength to the overall structure. As a result, they are slightly more abundant in thermostable proteins; however, like disulfide bonds, their effect cannot be generalized [70, 78, 118] . Moreover, introduction of a new H-bond at one site may have a distant effect elsewhere in the structure, which is hard to predict.

In conclusion, there is no amino acid consensus that is universally and equally applicable to all thermotolerant proteins. Indeed, it is tempting to speculate that even the different thermotolerant mutant proteins, independently arising from the same mesophilic cell, will not contain the same mutations in the same genes [119] . Nevertheless, the broad and generalized features of thermotolerance are summarized here ( Figure 4) and should serve as a useful reference.

The Tm of C116-C120 amylase activity was ~5 °C higher than that of the wild type (50 to 55 °C), where some other pairs (not shown) had a lower Tm. A triple mutant, C35-C426/C116-C120/C436-C480, showed a six-fold increase in half-life at 60 °C and a 5.2 °C increase in Tm over the wild type. Note that the last pair, viz C436-C480, introduced alone, had little effect on the Tm.

These results clearly demonstrate that rational design of Cys disulfide bonds to improve thermotolerance is possible, but the effect is contextual and likely influenced by other intramolecular interactions; thus, the properties of the mutants need to be experimentally tested. Lastly, the disulfide effect should be most pronounced in oxidizing environments that prevent reduction of the disulfide bond.

In contrast to disulfide bonds, the hydrogen bonds (H-bonds) are weak, but they form between a greater variety of amino acid side chains, as well as with water, and thus provide significant strength to the overall structure. As a result, they are slightly more abundant in thermostable proteins; however, like disulfide bonds, their effect cannot be generalized [70, 78, 118] . Moreover, introduction of a new H-bond at one site may have a distant effect elsewhere in the structure, which is hard to predict.

In conclusion, there is no amino acid consensus that is universally and equally applicable to all thermotolerant proteins. Indeed, it is tempting to speculate that even the different thermotolerant mutant proteins, independently arising from the same mesophilic cell, will not contain the same mutations in the same genes [119] . Nevertheless, the broad and generalized features of thermotolerance are summarized here ( Figure 4) and should serve as a useful reference. 

A futuristic and evolutionary question is whether we can harness the molecular knowledge of thermophilic mutations in designing proteins or creating organisms that will survive in global warming. In order to answer this query, we need to first determine how many mutations are needed to achieve useful thermotolerance. From the molecular differences described earlier, it is clear that stable and optimal thermophilia requires multiple mutational changes; however, there are a few examples where altering one or two amino acid(s) in a protein sometimes makes a difference, as described later. Another problem in designing a functional thermostable protein, particularly a large one, stems from the fact that the majority of missense mutations affects protein function [120] . This is because a single missense amino acid substitution can cause a free energy change of 0.5 to 5 kcal/mol, which most proteins can barely tolerate without losing stability [44, [120] [121] [122] . The native functional state of a protein is thus marginally stabilized by the balance of the opposing forces of folding and unfolding [121, 123] . Lastly, as alluded to earlier, structural rigidity improves thermal stability, but proteins also need to flexible in order to function and evolve. Thus, there is a trade-off between stability and flexibility, and directed evolution towards thermostability must balance between the two [124, 125] . Designing of a thermostable organism from a mesophilic one is a quantum jump in complexity and requires many mutations of the proteome that also need to be mutually compatible; this may be fodder for science fiction but not a currently achievable goal.

On the other hand, spontaneous generation of a phenotype in simple and rapidly replicating unicellular organisms happens under all environmental challenges, and the same is true of temperature rise, as seen time and again in Nature and in laboratory conditions [126, 127] . In an interesting early experiment [128] , 3-isopropylmalate dehydrogenase (LeuB) of the extreme thermophile, Thermus thermophilus, was replaced with a temperature-sensitive chimeric LeuB. A thermostable mutant was then selected by culturing the transformant to higher temperature, and was found to contain Leu93 in the place of Ile93, which coincided with the Ile93 in the natural T. thermophilus LeuB. The mutant enzyme was confirmed to be more stable than the original chimeric enzyme. In another example [129] , a mesophilic protein was converted to a highly thermostable form by changing the Glu3 and Glu66 to Arg and Leu, respectively. These were surface residues of the protein, which exemplify a simple and powerful approach for increasing the thermostability of a protein in small steps, on the way to higher levels of tolerance through multiple additional changes [70, 130] . In a more realistic, tour-de-force experiment [131] , mesophilic E. coli, grown in long-term culture, was subjected to increasing temperature, and the sequentially collected survivors were subjected to whole-genome sequencing. A correlative analysis of genotype and growth phenotype of the mutants identified the glpF and fab genes as being mutated in the high-temperature isolates. The glpF gene product is involved in glycerol transport and is essential for glycerol utilization; fab codes for a fatty acid desaturase/isomerase that adds double bonds to fatty acid moieties. While the glpF mutation was a stop codon at position three and, therefore, essentially did not make any protein, the glpF mutant had a single missense mutation that significantly lowered its activity. Thus, two metabolically important enzymes played the most critical role in organismic thermal resistance. These results provide strong indication that mesophilic organisms will be successful in producing temperature-resistant mutants under challenge. However, these studies also provided direct evidence that only a small range of temperature rise is tolerated, above which optimal metabolism and life cycle cannot be sustained. In other words, only a few degrees separate life from death, and this is why a rapid rise in global temperature may lead to extinction of many species, a process that has already begun for many nonhuman species, accompanied by a cascade effect [132] . Genetic changes in evolution take a generation to be inherited, which for larger animals with long life spans and slow reproductive maturity can be many years, making them particularly vulnerable.

Although spontaneous natural mutations followed by survival of the thermotolerant ones will perpetuate some form of life on the planet, our molecular understanding of thermophilia and thermoresistance, boosted by the attention to global warming, has become profoundly useful in a multitude of industrial applications, such as high-temperature environmental cleaning, biodegradation, and process engineering. In recent years, thermophilic organisms and enzymes have been optimized for hydrolysis of various biopolymers [16] [17] [18] 116, 133, 134] . The initial hurdle of the paucity of molecular genetic tools is being overcome [135] , and the expanding list of thermophiles used for recombinant gene manipulation now includes, but is not limited to, several genera of bacteria (Caldicellulosiruptor and Thermotoga) and archaea (Sulfolobus, Thermococcus, and Pyrococcus).

Long-chain polysaccharides such as cellulose, starch, and pectin, do not enter cells, and, therefore, their utilization must precede hydrolysis to smaller units by extracellular hydrolases. Use of the cognate enzymes at high temperature allows access to bonds buried in the polysaccharide branches and tertiary structures. The enzymatic processing of starch into various sugars constitutes an important area in food, pharmaceutical, paper pulp, and textile industries. In the example shown earlier (Figure 3 ), the thermostability of alkaline α-amylases was successfully improved through rational in silico design to add Cys-Cys disulfide bridges in the catalytic domain [116] . Reciprocally, the naturally thermophilic chloramphenicol acetyltransferase enzyme from Clostridium thermocellum, a thermophilic anaerobic bacterium, was engineered to enable it to convert isobutanol into isobutyl acetate through site-directed mutation of the conserved Phe97 in the substrate-binding pocket to Trp, which allowed the mutant bacteria to produce isobutyl acetate directly from cellulose [133] . Structure-based mutagenesis like these can be applied to other enzymes, as well to improve either thermotolerance or catalytic parameters. In a different example, thermophiles are being engineered to degrade plastics, a perilous pollutant of the environment. Polyethylene terephthalate (PET), a major component of solid plastic waste, including microplastics, which cannot be easily degraded, partly due to its solid state [136] . However, effective enzymatic hydrolysis of PET can occur at a glass transition temperature >70 • C. In one study [134] , C. thermocellum was genetically engineered to express extracellular thermophilic cutinase (LCC), isolated from a plant compost metagenome, which could degrade PET at up to 70 • C. In a pilot test, the resultant bacteria degraded~60% commercial PET in 14 h of incubation at 60 • C (the optimum growth temperature of C. thermocellum), which is much more efficient than mesophilic bacteria or algae, and produced soluble monomer feedstocks [134] .

Industrial fermentation at elevated temperature also reduces the chances of contamination or bacteriophage infection, which is common with mesophilic hosts such as E. coli or Salmonella [16] . Nevertheless, recombinant thermophilic enzymes can be expressed in a mesophile that is easier to grow in a fermenter and which has several technical advantages [44] : The recombinant enzymes can be easily purified by using a heat-treatment step, whereby most of the host bacterial proteins will denature and precipitate out; oftentimes, a thermophilic enzyme is also tolerant to chemical denaturants, such as guanidinium hydrochloride, which is useful if the recombinant needs to be extracted from inclusion bodies.

Natural selection, by definition, engender a fitter variant, which is, therefore, different from the original, and this would also apply to the survivors of global warming. The effect of warming on a variety of organisms across phylogeny has been documented; however, the underlying genetic and molecular determinants have remain largely unraveled, clearly due to the multiplicity and complexity of the genes involved. Most organisms that have survived so far have used mitigation strategies, generally changing their habitats and habits such as foraging places and periods [1, 2, 137] . It is common knowledge, for example, that most animals rest in the shades at mid-noon on the hot summer days, and forage or hunt in cooler periods. Unfortunately, despite its obvious life-saving value, mitigation without acquiring mutational resistance to heat is only a temporary solution that shifts the problem to another day for eventual extinction. Nonetheless, some of these phenotypic changes in thermal adaptation may be accompanied by genetic alterations that have not been fully delineated or interrogated for their long-term value in anthropogenic climate change [138] . A few examples, in which there is some gleaning of the molecular mechanism, are presented below.

Insects exhibit a complex phenotype of heat tolerance [139] . Perhaps counter-intuitively, insects often shrink in size and other morphological measurements as temperatures rises. Studies in the thermophilic ant [140] , Aphaenogaster senilis, showed lack of larval development at the suboptimal temperature of 20 • C. At higher temperatures, developmental speed increased, but adults were smaller. At the same time, thermal resistance tests confirmed that ants raised at 28 and 32 • C had higher half-lethal temperatures than those raised at 24 • C. These results and others [141] show that the ants can use phenotypic plasticity to adjust to heat before genetic mutational selection occurs, and it is independent of morphological dimensions. This trait is obviously useful for the survival of the insect in rapid climate warming, but it currently lacks a mechanistic explanation. It would be interesting to test how the molecular features of Aphaenogaster senilis proteins compare with their homologs in other insects of lower or higher optimal growth temperatures.

Much more is known about chemical changes in the insect exoskeleton induced by heat. In a very recent study [142] , three model but diverse ant species, namely Atta sexdens, Ectatomma brunneum, and Odontomachus bauri, were subjected to different temperatures for 5 h, and their fatty acids and cuticle hydrocarbon (CHC) contents were analyzed. None of the three was affected by temperatures below 30 • C, but each was differentially sensitive to higher temperatures. The capacity to survive was unrelated to fatty acids, but the CHCs showed significant quantitative and qualitative variation with temperature, having longer chain length of linear alkanes above 30 • C. These results are in agreement with the role of the cuticular exoskeleton of insects as a desiccation barrier, which is critical in preventing water loss at higher temperatures. The biosynthetic pathway of CHCs from fatty acids involves a complex network of "very long chain" fatty acid synthases, elongases, and desaturases, reductases, and finally, oxidative decarbonylases [143] . How the expression or activity of these enzymes is regulated by short-term heat exposure can be a highly rewarding area of research. It has been suggested that in the fatty acid synthase complex, the length of the fatty acid chain (usually 16-carbon-long palmitic acid) is sensed by a "molecular caliper" mechanism, in which the size of the hydrophobic pocket, where the synthase sits, determines the length, and that in the case of very long chain synthase, this pocket is larger [144] . With this model in mind, it is tempting to speculate that increased temperature may further enlarge this pocket, allowing a longer alkene chain to be formed. It is also worth testing whether the heat shock proteins play an as yet unknown role in this regulation, perhaps stabilizing the enlarged synthase.

At odds with our expectations, there are a few cases where a specific organism has actually benefited from global warming, at least in an isolated, short-term consideration. Higher underwater temperatures as a rule increase a plant's demand for carbon, but at the same time, the accompanying floodwaters increase water turbidity and lower the light levels, compromising photosynthesis, and hence, creating carbon deficiency. The presence of multiple other stress factors accompanying a temperature shift creates a synergistic, lethal effect; however, for some seagrasses, warmer waters appear beneficial, resulting in faster reproduction, the basis of which remains to be studied [1] . Due to the complexity of multifactorial phenotypic changes and strain variations, it has been cautioned that climatic results without a mechanism may not have the predictive value expected of them [145] [146] [147] . It has been conjectured that the most extreme organisms may be found in exoplanets and atmospheres in outer space [148] [149] [150] . If indeed found, these alien organisms will most certainly reveal new molecular and metabolic features, but their very exotic nature will also make them nearly impossible to study in terrestrial laboratories.

It is important to realize that biological evolution in global warming is indeed a "global" phenomenon, in which selection of the fittest operates not only at the level of molecular, cellular, or individual species level, but also at the population level, which includes all types reproduction and genetic recombination, and perhaps also interaction between heterologous populations (e.g., "host-pathogen interactions", Section 4.3.2). Thus, long-term survival is dictated by a concerted and complex interplay between genomic changes, animal behavior, and mitigation by migration. Nevertheless, as migration is the most immediate response of animals to temperature change, it has received significant attention. Habitat range shifts and range contractions have played particularly major roles in the history of all species; however, the molecular consequences of these processes have been investigated only recently. One population study [151] employed computer simulations to study the patterns of molecular diversity both within and between populations, and for multiple types of range contractions and shifts. The results revealed that changes in range contractions affect genetic diversity but the speed of change have contrasting effects; specifically, fast range contractions lead to higher diversity and lower genetic differentiation than slow contractions. This contrasts with the effect of range shifts: Fast shifts lead to lower levels of diversity, when compared to slow shifts. Overall, the results suggested that shifts and contractions of range interact in a complex but tractable manner and that the levels of diversity preserved after a climate change both within and between refuge areas will not only depend on the dispersal abilities of a species but also on the speed of the change. It also implies that, in any climatic change, species with different generation times will be impacted differently.

For many pathogens that are transmitted by vectors, such as tick-borne parasites, a visible change is their gradual migration to higher altitudes as the normal habitat of their vector gets warmer, which is again a mitigating response [12] ; as mentioned before, a more universal adaptation to temperature can occur only by genetic changes in the host or the pathogen. The obligatory parasites, such as viruses, must evolve in co-ordination with the host cellular machinery, which they co-opt. The interaction of their mutually evolving protein domains is currently an interesting area of structural and biochemical analysis, as shown recently [152] . Infectious pathogens that are also free-living, such as the disease-causing bacteria, yeast and fungi, and some unicellular protozoa (e.g., Trypanosomes), can evolve to encode heat-tolerant mutations independent of their hosts. However, for optimal infection, they too need to fit to the pathways of the thermotolerant host, in which the mutations to heat-resistance may have produced subtle changes in their interaction with the cognate parasite function. Important molecular mechanisms in this area have been gleaned from viral pathogens, in part because host-virus interactions lend themselves to tractable genetic and biochemical assays. Viruses are by far the most abundant pathogens in the world's oceans that are estimated to contain 10 30 virions that infect diverse marine organisms [153] . The large pool not only allows for host-virus co-evolution, but also for the generation of novel virus strains by genetic assortment, such as the recent discovery of apparently novel coronaviruses in beluga whale and dolphins [154] [155] [156] . A case of apparent benefit of the host has been noted in Emiliania huxleyi, a ubiquitous marine phytoplankton, responsible for dense annual blooms that have a profound impact on planetary carbon and sulfur cycles [157] . The blooms are regularly terminated by dsDNA viruses known as E. huxleyi Viruses or EhVs. In an interesting recent study, E. huxleyi was found to develop EhV-resistance when the temperature increased by only 3 • C; this also resulted in increased production of dimethyl sulfide, a major plankton metabolite that forms "cloud condensation nuclei", a global rain-maker. The authors proposed that alterations to the E. huxleyi surface receptors are responsible for the temperature-induced resistance.

Finally, there are many other peripheral areas of global warming that will also affect biological life, such as possible alterations of light intensity and circadian rhythm, and change of the environment from aerobic to anaerobic [158] [159] [160] [161] , which will add to the burden of selection.

It should be noted that in spite of the most careful and methodical analyses, the complex nature of the protein folding ensemble, entropy landscape, and kinetic parameters makes it nearly impossible to derive an invariable universal set of rules for thermotolerance in all organisms or all proteins. As we have seen in the review, practically all molecular parameters derived from thermotolerance studies have some exceptions [70] . Nonetheless, these "nearly universal" molecular signature(s) will form a useful framework against which essentially all evolutionary selections in climate change can be interrogated. A simplified summary of interactive molecular changes is schematically depicted here ( Figure 5 ).

It should be noted that in spite of the most careful and methodical analyses, the complex nature of the protein folding ensemble, entropy landscape, and kinetic parameters makes it nearly impossible to derive an invariable universal set of rules for thermotolerance in all organisms or all proteins. As we have seen in the review, practically all molecular parameters derived from thermotolerance studies have some exceptions [70] . Nonetheless, these "nearly universal" molecular signature(s) will form a useful framework against which essentially all evolutionary selections in climate change can be interrogated. A simplified summary of interactive molecular changes is schematically depicted here ( Figure 5 ). Figure 5 . Selection for thermostable mutations, as well as mitigation by organismal behavior and habitat (range) change, in global warming. A pathogen (the smaller structure with flagella) and its host (the larger structure) are shown as a pair, to indicate that they will mutate together, with increasing temperature, as time passes. The changes are indicated by the altered structure of the protein (shown in three colors for the domains, to match Figure 3 ) and possibly by altered morphology. Simultaneously, the organisms may migrate to higher habitats, for a cooler atmosphere.

As indicated in the diagram, mutation and mitigation will both work concurrently, the first as a genetic change and the second as an immediate behavioral reaction. In persistent warming, various forms of life will prevail, but radical mutational changes will create highly modified, newer organisms, many of which will be morphologically different, with altered community structures [162] . Invasive and predatory species may also take over as privileged survivors. Finally, in terms of novel mechanisms, it is tempting to postulate that molecular changes may also result from the non-genetic mechanisms (i.e., without genomic mutations) that include but are not limited to the following: (a) transcriptional or translation induction of a protein by a thermo-sensing mechanism, not unlike heatshock response; (b) processing of a pro-enzyme into its active form; (c) alternative splicing, promoted by a heat-induced splicing factor; (d) altered RNA fold; (e) post-translational modifications, such as phosphorylation or methylation; (f) change of conformation of a protein, accompanying temperature changes [163] , for which dramatic examples are known, e.g., lymphotactin [164] ; and (g) epigenetic changes in the genome. These will constitute novel areas of future investigation of molecular biological strategies adopted by living beings in their effort to survive global warming. Figure 5 . Selection for thermostable mutations, as well as mitigation by organismal behavior and habitat (range) change, in global warming. A pathogen (the smaller structure with flagella) and its host (the larger structure) are shown as a pair, to indicate that they will mutate together, with increasing temperature, as time passes. The changes are indicated by the altered structure of the protein (shown in three colors for the domains, to match Figure 3 ) and possibly by altered morphology. Simultaneously, the organisms may migrate to higher habitats, for a cooler atmosphere.

As indicated in the diagram, mutation and mitigation will both work concurrently, the first as a genetic change and the second as an immediate behavioral reaction. In persistent warming, various forms of life will prevail, but radical mutational changes will create highly modified, newer organisms, many of which will be morphologically different, with altered community structures [162] . Invasive and predatory species may also take over as privileged survivors. Finally, in terms of novel mechanisms, it is tempting to postulate that molecular changes may also result from the non-genetic mechanisms (i.e., without genomic mutations) that include but are not limited to the following: (a) transcriptional or translation induction of a protein by a thermo-sensing mechanism, not unlike heat-shock response; (b) processing of a pro-enzyme into its active form; (c) alternative splicing, promoted by a heat-induced splicing factor; (d) altered RNA fold; (e) post-translational modifications, such as phosphorylation or methylation; (f) change of conformation of a protein, accompanying temperature changes [163] , for which dramatic examples are known, e.g., lymphotactin [164] ; and (g) epigenetic changes in the genome. These will constitute novel areas of future investigation of molecular biological strategies adopted by living beings in their effort to survive global warming.

Funding: This review received no external funding. The article-processing charge was paid for by the author's personal funds.

The author declares no conflict of interest.

Fingerprints of global warming on wild animals and plants

Fingerprints of global warming on wild animals and plants

Climate change, habitat loss, protected areas and the climate adaptation potential of species in Mediterranean ecosystems worldwide

Biological consequences of global warming: Is the signal already apparent?

Hyperthermal-driven mass extinctions: Killing models during the Permian-Triassic mass extinction

Ocean warming alters cellular metabolism and induces mortality in fish early life stages: A proteomic approach

Investigating cellular stress response to heat stress in the seagrass Posidonia oceanica in a global change scenario

Climate change and an invasive, tropical milkweed: An ecological trap for monarch butterflies

Plant functional trait change across a warming tundra biome

Past and future global transformation of terrestrial ecosystems under climate change

The economics of the Kyoto protocol

The brave new world of carbon trading

Global warming and its health impact

History of discovery of the first hyperthermophiles

Molecular mechanisms of adaptation emerging from the physics and evolution of nucleic acids and proteins

The more adaptive to change, the more likely you are to survive: Protein adaptation in extremophiles

Extremely thermophilic microorganisms as metabolic engineering platforms for production of fuels and industrial chemicals

Potential and utilization of thermophiles and thermostable enzymes in biorefining. Microb. Cell Fact

Probing the folding-unfolding transition of a thermophilic protein

Is there a unifying mechanism for protein folding?

Phi-value analysis and the nature of protein-folding transition states

Identification and characterization of protein folding intermediates

A novel class of dual-family immunophilins

Chaperones and proteostasis: Role in Parkinson's Disease. Diseases 2020, 8, 24

A comparison of structural and evolutionary attributes of Escherichia coli and Thermus thermophilus small ribosomal subunits: Signatures of thermal adaptation

Inhibitor-induced conformational stabilization and structural alteration of a mip-like peptidyl prolyl cis-trans isomerase and its C-terminal domain

Bacillus subtilis builds structurally and functionally different spores in response to the temperature of growth

Recombinant phosphoglycerate kinase from the hyperthermophilic bacterium Thermotoga maritima: Catalytic, spectral and thermodynamic properties

The family 1b-glucosidases from Pyrococcus furiosus and Agrobacterium faecalis share a common catalytic mechanism

xylA cloning and sequencing and biochemical characterization of xylose from Thermotoga neapolitana

Glyceraldehyde-3-phosphate dehydrogenase from the hyperthermophilic archaebacterium Pyrococcus woesei: Characterization of the enzyme, cloning and sequencing of the gene, and expression in Escherichia coli

Lactate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima: The crystal structure at 2.1 A resolution reveals strategies for intrinsic protein stabilization

Crystal structure of the beta-glycosidase from the hyperthermophile Thermosphaera aggregans: Insights into its activity and thermostability

Crystal structure of methionine aminopeptidase from hyperthermophile, Pyrococcus furiosus

The crystal structure of citrate synthase from the hyperthermophilic archaeon Pyrococcus furiosus at 1.9A resolution

Cells in stress: Transcriptional activation of heat shock genes

The heat shock response: Life on the verge of death

Hormetic heat stress and HSF-1 induce autophagy to improve survival and proteostasis in C. elegans

Punctuated equilibrium comes of age

Protein stability promotes evolvability

A relationship between protein stability and protein function

How protein stability and new functions trade off

Ancient DNA and human history

Analyzing protein structure and function using ancestral gene reconstruction

Maximum-likelihood phylogenetic inference with selection on protein folding stability

Characterization of reconstructed ancestral proteins suggests a change in temperature of the ancient biosphere

ProtASR: An evolutionary framework for ancestral protein reconstruction with selection on folding stability

How to resurrect ancestral proteins as proxies for ancient biogeochemistry. Free Radic

Experimental evidence for the thermophilicity of ancestral life

The physiology and habitat of the last universal common ancestor

Evolutionary drivers of thermoadaptation in enzyme catalysis

Neoarchean deep marine paleotemperature: Evidence from turbidite successions

Rooting the archaebacterial tree: The pivotal role of Thermococcus celer in archaebacterial evolution

Toward automatic reconstruction of a highly resolved tree of life

On the role, ecology, phylogeny, and structure of dual-family immunophilins

Ancestral amino acid substitution improves the thermal stability of recombinant lignin-peroxidase from white-rot fungi, Phanerochaete chrysosporium strain UAMH 3641

Thermostability of ancestral mutants of Caldococcus noboribetus isocitrate dehydrogenase

Ancestral residues stabilizing 3-isopropylmalate dehydrogenase of an extreme thermophile: Experimental evidence supporting the thermophilic common ancestor hypothesis

Palaeotemperature trend for Precambrian life inferred from resurrected proteins

Single-molecule paleoenzymology probes the chemistry of resurrected enzymes

Epistasis effects of multiple ancestral-consensus amino acid substitutions on the thermal stability of glycerol kinase from Cellulomonas sp. NT3060

Ancestral sequence reconstruction produces thermally stable enzymes with mesophilic enzyme-like catalytic properties

Hyperthermophilic enzymes: Sources, uses, and molecular mechanisms for thermostability. Microbiol

Searching for folded proteins in vitro and in silico

Solvent accessibility, protein surfaces, and protein folding

Helix to helix packing in proteins

Structural mechanisms for domain movements in proteins

What has de novo protein design taught us about protein folding and biophysics?

Finding the generalized molecular principles of protein thermal stability

Protein and DNA sequence determinants of thermophilic adaptation

Compositional and structural features related to thermal stability in the archaea SRP19 and SRP54 signal recognition particle proteins

Protein thermostability: Structure-based difference of amino acid between thermophilic and mesophilic proteins

Factors enhancing protein thermostability

Environment specific substitution tables for thermophilic proteins

Deciphering the rationale behind specific codon usage pattern in extremophiles

Insights into thermal resistance of proteins from the intrinsic stability of their α-helices

Structural differences between mesophilic, moderately thermophilic and extremely thermophilic protein subunits: Results of a comprehensive survey

Solving the membrane protein folding problem

Energetics of membrane protein folding and stability

Enhanced protein thermostability from site-directed mutations that decrease the entropy of unfolding

Effective factors in thermostability of thermophilic proteins

Differences in amino acids composition and coupling patterns between mesophilic and thermophilic proteins

Stability and solubility of proteins from extremophiles

Thermal adaptation analyzed by comparison of protein sequences from mesophilic and extremely thermophilic Methanococcus species

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

Discrimination of thermophilic and mesophilic proteins

Key role of proximal water in regulating thermostable proteins

Comparison of protein-water interactions in psychrophilic, mesophilic, and thermophilic Fe-SOD

Mechanisms of thermal adaptation revealed from the genomes of the Antarctic archaea Methanogenium frigidum and Methanococcoides burtonii

Ion pairs and the thermotolerance of proteins from hyperthermophiles: A "Traffic Rule" for hot roads

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

Protein stability governed by its structural plasticity is inferred by physicochemical factors and salt bridges

Consensus sequence design as a general strategy to create hyperstable, biologically active proteins

Evolution and thermodynamics of the slow unfolding of hyperstable monomeric proteins

Hydrophobic environment is a key factor for the stability of thermophilic proteins

Evolvability of thermophilic proteins from archaea and bacteria

Asparagine deamidation dependence on buffer type, pH, and temperature

Amino acid contacts in proteins adapted to different temperatures: Hydrophobic interactions and surface charges play a key role

The uniqueness of tryptophan in biology: Properties, metabolism, interactions and localization in proteins

Design of stable alpha-helical arrays from an idealized TPR motif

Protein tetratricopeptide repeat and the companion non-tetratricopeptide repeat helices: Bioinformatic analysis of interhelical interactions

The nature and arrangement of pentatricopeptide domains and the linker sequences between them

Tandem repeats in proteins: From sequence to structure

Pentatricopeptide repeat proteins in plants

The tetratricopeptide-repeat motif is a versatile platform that enables diverse modes of molecular recognition

All repeats are not equal: A module-based approach to guide repeat protein design

Mechanisms of irreversible thermal inactivation of Bacillus alpha-amylases

Glutamine deamidation destabilizes human gammaD-crystallin and lowers the kinetic barrier to unfolding

Does deamidation cause protein unfolding? A top-down tandem mass spectrometry study

Substitution of asparagine residues in Aspergillus awamori glucoamylase by site-directed mutagenesis to eliminate N-glycosylation and inactivation by deamidation

Molecular determinants of xylose isomerase thermal stability and activity: Analysis of thermozymes by site-directed mutagenesis

The crystal structure of 5 -deoxy-5 -methylthioadenosine phosphorylase II from Sulfolobus solfataricus, a thermophilic enzyme stabilized by intramolecular disulfide bonds

Roles of cysteines Cys115 and Cys201 in the assembly and thermostability of grouper betanodavirus particles

The effect of additional disulfide bonds on the stability and folding of ribonuclease A

In silico rational design and systems engineering of disulfide bridges in the catalytic domain of an alkaline α-amylase from Alkalimonas amylolytica to improve thermostability

The I-TASSER suite: Protein structure and function prediction

Thermophilic proteins as versatile scaffolds for protein engineering

Structural basis of thermostability in hyperthermophilic proteins, or "there's more than one way to skin a cat

Genetic constraints on protein evolution

Lessons in stability from thermophilic proteins

Missense meanderings in sequence space: A biophysical view of protein evolution

How quickly do proteins fold and unfold, and what structural parameters correlate with these values? Biomolecules

How proteins adapt: Lessons from directed evolution

Rigidity versus flexibility: The dilemma of understanding protein thermal stability

Genetic architecture of thermal adaptation in Escherichia coli

Genetic variation in clones of Pseudomonas pseudoalcaligenes after ten months of selection in different thermal environments in the laboratory

Screening of stable proteins in an extreme thermophile, Thermus thermophilus

Two exposed amino acid residues confer thermostability on a cold shock protein

Serial increase in the thermal stability of 3-isopropylmalate dehydrogenase from Bacillus subtilis by experimental evolution

Experimental evolution of a facultative thermophile from a mesophilic ancestor

Trophic downgrading of planet Earth

Single mutation at a highly conserved region of chloramphenicol acetyltransferase enables isobutyl acetate production directly from cellulose by Clostridium thermocellum at elevated temperatures

Thermophilic whole-cell degradation of polyethylene terephthalate using engineered Clostridium thermocellum

Developing a genetic manipulation system for the Antarctic archaeon, Halorubrum lacusprofundi: Investigating acetamidase gene function

Microplastic pollution in water and sediment in a textile industrial area

Emergent sustainability in open property regimes

Genomic changes associated with adaptation to arid environments in cactophilic Drosophila species

Complex body size differences in thermal tolerance among army ant workers (Eciton burchellii parvispinum)

Is phenotypic plasticity a key mechanism for responding to thermal stress in ants

Phenotypic plasticity of cuticular hydrocarbon profiles in insects

Effect of temperature on survival and cuticular composition of three different ant species

An insect-specific P450 oxidative decarbonylase for cuticular hydrocarbon biosynthesis

A molecular caliper mechanism for determining very long-chain fatty acid length

Biophysics, physiological ecology, and climate change: Does mechanism matter?

Thermal tolerance, acclimatory capacity and vulnerability to global climate change

Long-term dynamics among Wolbachia strains during thermal adaptation of their Drosophila melanogaster hosts

Microbial extremophiles at the limits of life

Exoplanet biosignatures: A review of remotely detectable signs of life

The adaptability of life on Earth and the diversity of planetary habitats

Consequences of range contractions and range shifts on molecular diversity

Temporal dynamics of protein complex formation and dissociation during human cytomegalovirus infection

Marine viruses-major players in the global ecosystem

Identification of a novel coronavirus from a beluga whale by using a panviral microarray

Discovery of a novel bottlenose dolphin coronavirus reveals a distinct species of marine mammal coronavirus in Gammacoronavirus

Genus-specific pattern of intrinsically disordered central regions in the nucleocapsid protein of coronaviruses

Temperature-induced viral resistance in Emiliania huxleyi (Prymnesiophyceae)

Climate warming and disease risks for terrestrial and marine biota

Influence of irradiance and temperature on the virus MpoV-45T infecting the arctic picophytoplankter Micromonas polaris. Viruses

Molecular interactions between pathogens and the circadian clock

Regulation of NF-kappa B and cell death by bacterial gingipains

Temperature increase altered Daphnia community structure in artificially heated lakes: A potential scenario for a warmer future

Proteins that switch folds

Interconversion between two unrelated protein folds in the lymphotactin native state