How do thermophilic proteins and proteomes withstand high temperature?

The C-terminus is essential for the stability of the mycobacterial channel protein MspA

Outer membrane proteins perform essential functions in uptake and secretion processes in bacteria. MspA is an octameric channel protein in the outer membrane of Mycobacterium smegmatis and is structurally distinct from any other known outer membrane protein. MspA is the founding member of a family with more than 3000 homologs and is one of the most widely used proteins in nanotechnological applications due to its advantageous pore structure and extraordinary stability. While a conserved C-terminal signal sequence is essential for folding and protein assembly in the outer membrane of Gram-negative bacteria, the molecular determinants of these processes are unknown for MspA. In this study, we show that mutation and deletion of methionine 183 in the highly conserved C-terminus of MspA and mutation of the conserved tryptophan 40 lead to a complete loss of protein in heat extracts of M. smegmatis. Swapping these residues partially restores the heat stability of MspA indicating that methionine 183 and tryptophan 40 form a conserved sulfur-π electron interaction, which stabilizes the MspA monomer. Flow cytometry showed that all MspA mutants are surface-accessible demonstrating that oligomerization and membrane integration in M. smegmatis are not affected. Thus, the conserved C-terminus of MspA is essential for its thermal stability, but it is not required for protein assembly in its native membrane, indicating that this process is mediated by a mechanism distinct from that in Gram-negative bacteria. These findings will benefit the rational design of MspA-like pores to tailor their properties in current and future applications.

Revisiting structural organization of proteins at high temperature from a network perspective

Interactions between distantly placed amino acids in the primary chain (long-range) play a very crucial role in the formation and stabilization of the tertiary structure of a protein, while interactions between closely placed amino acids in the primary chain (short-range) mostly stabilize the secondary structures. Every protein needs to maintain marginal stability in order to perform its physiological functions in its native environment. The requirements for this stability in mesophilic and thermophilic proteins are different. Thermophilic proteins need to form more interactions as well as more stable interactions to survive in the extreme environment, they live in. Here, we aim to find out how the interacting amino acids in three-dimensional space are positioned in the primary chains in thermophilic and mesophilic. How does this arrangement help thermophiles to maintain their structural integrity at high temperatures? Working on a dataset of 1560 orthologous pairs we perceive that thermophiles are not only enriched with long-range interactions, they feature bigger connected clusters and higher network densities compared to their mesophilic orthologs, at higher interaction strengths between the amino acids. Moreover, we have observed the enrichment of different types of interactions at different secondary structural regions.

Surface exposed and charged residues drive thermostability in fungi

Fungi, though mesophilic, include thermophilic and thermostable species, as well. The thermostability of proteins observed in these fungi is most likely to be attributed to several molecular factors, such as the presence of salt bridges and hydrogen bond interactions between side chains. These factors cannot be generalized for all fungi. Factors impacting thermostability can guide how fungal thermophilic proteins gain thermostability. We curated a dataset of proteins for 14 thermophilic fungi and their evolutionarily closer mesophiles. Additionally, the proteome of Chaetomium thermophilum and its evolutionarily related mesophile Chaetomium globosum was analyzed. Using eggNOG, we categorized the proteomes into clusters of orthologous groups (COGs). While the individual count of proteins is over-represented in mesophiles (for COGs S, G, L, and Q), there are certain features that are significantly enriched in thermophiles (such as charged residues, exposed residues, polar residues, etc.). Since fungi are known to be cellulolytic and chitinolytic by nature, we selected 37 existing carbohydrate-active enzymes (CAZyme) families in Eurotiales, Mucorales, and Sordariales. We looked at closely similar sequences and their modeled structures for further comparison. Comparing solvent accessibilities of thermophilic and mesophilic proteins, exposed and intermediate residues are observed higher in thermophiles whereas buried residues are observed higher in mesophiles. For specific five CAZYme families (GH7, GH11, GH18, GH45, and CBM1) we looked at position-specific substitutions between thermophiles and mesophiles. We also found that there are relatively more intramolecular interactions in thermophiles compared to mesophiles. Thus, we found factors such as surface exposed residues and charged residues that are highly likely to impart thermostability in fungi, and this study sets the stage for further studies in the area of fungal thermostability.

Homologous Pairs of Low and High Temperature Originating Proteins Spanning the Known Prokaryotic Universe

Stability of proteins at high temperature has been a topic of interest for many years, as this attribute is favourable for applications ranging from therapeutics to industrial chemical manufacturing. Our current understanding and methods for designing high-temperature stability into target proteins are inadequate. To drive innovation in this space, we have curated a large dataset, learn2thermDB, of protein-temperature examples, totalling 24 million instances, and paired proteins across temperatures based on homology, yielding 69 million protein pairs - orders of magnitude larger than the current largest. This important step of pairing allows for study of high-temperature stability in a sequence-dependent manner in the big data era. The data pipeline is parameterized and open, allowing it to be tuned by downstream users. We further show that the data contains signal for deep learning. This data offers a new doorway towards thermal stability design models.

The Action of Chemical Denaturants: From Globular to Intrinsically Disordered Proteins

Proteins perform their many functions by adopting either a minimal number of strictly similar conformations, the native state, or a vast ensemble of highly flexible conformations. In both cases, their structural features are highly influenced by the chemical environment. Even though a plethora of experimental studies have demonstrated the impact of chemical denaturants on protein structure, the molecular mechanism underlying their action is still debated. In the present review, after a brief recapitulation of the main experimental data on protein denaturants, we survey both classical and more recent interpretations of the molecular basis of their action. In particular, we highlight the differences and similarities of the impact that denaturants have on different structural classes of proteins, i.e., globular, intrinsically disordered (IDP), and amyloid-like assemblies. Particular attention has been given to the IDPs, as recent studies are unraveling their fundamental importance in many physiological processes. The role that computation techniques are expected to play in the near future is illustrated.

Creatinase: Using Increased Entropy to Improve the Activity and Thermostability

Improving protein thermostability in mutagenesis-based enzyme engineering was often achieved by enhancing interresidue interactions via mutation to increase the enthalpy penalty of unfolding. However, this approach may trade off the functional activity due to the loss of structural flexibility of the biomolecule. Here, by performing X-ray crystallography, enzymatic kinetic experiments, neutron scattering, and thermodynamical measurements, we compared the structures, catalytic behaviors, dynamics, and thermostability between a wild-type creatinase and its four-point mutant. We found that the mutant is an entropy-driven thermostable protein with higher structural flexibility, i.e., higher conformational entropy, in the folded state compared to the wild type. The increased conformational entropy of the mutant in the folded state can reduce the entropy gain during unfolding and thus renders it greater thermostability. Moreover, the increased structural flexibility, particularly around the catalytic site, can broaden the mutant's working temperature range and considerably improve its activity at ambient conditions, which is crucial for its application in diagnosing kidney diseases. Complementary all-atom molecular dynamics simulations indicated that the four mutations replaced several of the strong interresidue interactions (electrostatic interactions and hydrogen bonds) with weak hydrophobic interactions. These substitutions not only release the structural flexibility to promote the thermostability and enzymatic activity of the protein but they also preserve the protein structure from collapsing. Our findings may pave a route for the entropy-driven strategy to design proteins with high thermostability and activity.

Structural Studies of Pif1 Helicases from Thermophilic Bacteria

Pif1 proteins are DNA helicases belonging to Superfamily 1, with 5' to 3' directionality. They are conserved from bacteria to human and have been shown to be particularly important in eukaryotes for replication and nuclear and mitochondrial genome stability. However, Pif1 functions in bacteria are less known. While most Pif1 from mesophilic bacteria consist of the helicase core with limited N-terminal and C-terminal extensions, some Pif1 from thermophilic bacteria exhibit a C-terminal WYL domain. We solved the crystal structures of Pif1 helicase cores from thermophilic bacteria Deferribacter desulfuricans and Sulfurihydrogenibium sp. in apo and nucleotide bound form. We show that the N-terminal part is important for ligand binding. The full-length Pif1 helicase was predicted based on the Alphafold algorithm and the nucleic acid binding on the Pif1 helicase core and the WYL domain was modelled based on known crystallographic structures. The model predicts that amino acids in the domains 1A, WYL, and linker between the Helicase core and WYL are important for nucleic acid binding. Therefore, the N-terminal and C-terminal extensions may be necessary to strengthen the binding of nucleic acid on these Pif1 helicases. This may be an adaptation to thermophilic conditions.

Diffusive Dynamics of Bacterial Proteome as a Proxy of Cell Death

Temperature variations have a big impact on bacterial metabolism and death, yet an exhaustive molecular picture of these processes is still missing. For instance, whether thermal death is determined by the deterioration of the whole or a specific part of the proteome is hotly debated. Here, by monitoring the proteome dynamics of E. coli, we clearly show that only a minor fraction of the proteome unfolds at the cell death. First, we prove that the dynamical state of the E. coli proteome is an excellent proxy for temperature-dependent bacterial metabolism and death. The proteome diffusive dynamics peaks at about the bacterial optimal growth temperature, then a dramatic dynamical slowdown is observed that starts just below the cell's death temperature. Next, we show that this slowdown is caused by the unfolding of just a small fraction of proteins that establish an entangling interprotein network, dominated by hydrophobic interactions, across the cytoplasm. Finally, the deduced progress of the proteome unfolding and its diffusive dynamics are both key to correctly reproduce the E. coli growth rate.

A proteome-scale analysis of vertebrate protein amino acid occurrence: Thermoadaptation shows a correlation with protein solvation but less so with dynamics

Despite differences in behaviors and living conditions, vertebrate organisms share the great majority of proteins, often with subtle differences in amino acid sequence. Here, we present a simple way to analyze the difference in amino acid occurrence by comparing highly homologous proteins on a subproteome level between several vertebrate model organisms. Specifically, we use this method to identify a pattern of amino acid conservation as well as a shift in amino acid occurrence between homeotherms (warm-blooded species) and poikilotherms (cold-blooded species). Importantly, this general analysis and a specific example further establish a broad correlation, if not likely connection between the thermal adaptation of protein sequences and two of their physical features: on average a change in their protein dynamics and, even more strongly, in their solvation. For poikilotherms, such as frog and fish, the lower body temperature is expected to increase the protein-protein interaction due to a decrease in protein internal dynamics. In order to counteract the tendency for enhanced binding caused by low temperatures, poikilotherms enhance the solvation of their proteins by favoring polar amino acids. This feature appears to dominate over possible changes in dynamics for some proteins. The results suggest that a general trend for amino acid choice is part of the mechanism for thermoadaptation of vertebrate organisms at the molecular level.

The methanogen core and pangenome: conservation and variability across biology's growth temperature extremes

Temperature is a key variable in biological processes. However, a complete understanding of biological temperature adaptation is lacking, in part because of the unique constraints among different evolutionary lineages and physiological groups. Here we compared the genomes of cultivated psychrotolerant and thermotolerant methanogens, which are physiologically related and span growth temperatures from -2.5°C to 122°C. Despite being phylogenetically distributed amongst three phyla in the archaea, the genomic core of cultivated methanogens comprises about one-third of a given genome, while the genome fraction shared by any two organisms decreases with increasing phylogenetic distance between them. Increased methanogenic growth temperature is associated with reduced genome size, and thermotolerant organisms-which are distributed across the archaeal tree-have larger core genome fractions, suggesting that genome size is governed by temperature rather than phylogeny. Thermotolerant methanogens are enriched in metal and other transporters, and psychrotolerant methanogens are enriched in proteins related to structure and motility. Observed amino acid compositional differences between temperature groups include proteome charge, polarity and unfolding entropy. Our results suggest that in the methanogens, shared physiology maintains a large, conserved genomic core even across large phylogenetic distances and biology's temperature extremes.

Learning deep representations of enzyme thermal adaptation

Temperature is a fundamental environmental factor that shapes the evolution of organisms. Learning thermal determinants of protein sequences in evolution thus has profound significance for basic biology, drug discovery, and protein engineering. Here, we use a data set of over 3 million BRENDA enzymes labeled with optimal growth temperatures (OGTs) of their source organisms to train a deep neural network model (DeepET). The protein-temperature representations learned by DeepET provide a temperature-related statistical summary of protein sequences and capture structural properties that affect thermal stability. For prediction of enzyme optimal catalytic temperatures and protein melting temperatures via a transfer learning approach, our DeepET model outperforms classical regression models trained on rationally designed features and other deep-learning-based representations. DeepET thus holds promise for understanding enzyme thermal adaptation and guiding the engineering of thermostable enzymes.

The Thermostability of Nucleoside Phosphorylases from Prokaryotes. I. The Role of the Primary Structure of the N-terminal fragment of the Protein in the Thermostability of Uridine Phosphorylases

Mutant uridine phosphorylase genes from Shewanella oneidensis MR-1 (S. oneidensis) were constructed by site-directed mutagenesis and strains-producers of the corresponding recombinant (F5I and F5G) proteins were obtained on the basis of Escherichia coli cells. The mutant proteins were purified and their physicochemical and enzymatic properties were studied. It was shown that the N-terminal fragment of uridine phosphorylase plays an important role in the thermal stabilization of the enzyme as a whole. The role of the aminoacid (a.a.) residue phenylalanine (F5) in the formation of thermotolerance of uridine phosphorylases from gamma-proteobacteria was revealed.

Adaptive Response of Thermophiles to Redox Stress and Their Role in the Process of dye Degradation From Textile Industry Wastewater

Release of dye-containing textile wastewater into the environment causes severe pollution with serious consequences on aquatic life. Bioremediation of dyes using thermophilic microorganisms has recently attracted attention over conventional treatment techniques. Thermophiles have the natural ability to survive under extreme environmental conditions, including high dye concentration, because they possess stress response adaptation and regulation mechanisms. Therefore, dye detoxification by thermophiles could offer enormous opportunities for bioremediation at elevated temperatures. In addition, the processes of degradation generate reactive oxygen species (ROS) and subject cells to oxidative stress. However, thermophiles exhibit better adaptation to resist the effects of oxidative stress. Some of the major adaptation mechanisms of thermophiles include macromolecule repair system; enzymes such as superoxide dismutase, catalase, and glutathione peroxidase; and non-enzymatic antioxidants like extracellular polymeric substance (EPSs), polyhydroxyalkanoates (PHAs), etc. In addition, different bacteria also possess enzymes that are directly involved in dye degradation such as azoreductase, laccase, and peroxidase. Therefore, through these processes, dyes are first degraded into smaller intermediate products finally releasing products that are non-toxic or of low toxicity. In this review, we discuss the sources of oxidative stress in thermophiles, the adaptive response of thermophiles to redox stress and their roles in dye removal, and the regulation and crosstalk between responses to oxidative stress.

Some Clues about Enzymes from Psychrophilic Microorganisms

Enzymes purified from psychrophilic microorganisms prove to be efficient catalysts at low temperatures and possess a great potential for biotechnological applications. The low-temperature catalytic activity has to come from specific structural fluctuations involving the active site region, however, the relationship between protein conformational stability and enzymatic activity is subtle. We provide a survey of the thermodynamic stability of globular proteins and their rationalization grounded in a theoretical approach devised by one of us. Furthermore, we provide a link between marginal conformational stability and protein flexibility grounded in the harmonic approximation of the vibrational degrees of freedom, emphasizing the occurrence of long-wavelength and excited vibrations in all globular proteins. Finally, we offer a close view of three enzymes: chloride-dependent α-amylase, citrate synthase, and β-galactosidase.

The Effects of Temperature on Cellular Physiology

Temperature impacts biological systems across all length and timescales. Cells and the enzymes that comprise them respond to temperature fluctuations on short timescales, and temperature can affect protein folding, the molecular composition of cells, and volume expansion. Entire ecosystems exhibit temperature-dependent behaviors, and global warming threatens to disrupt thermal homeostasis in microbes that are important for human and planetary health. Intriguingly, the growth rate of most species follows the Arrhenius law of equilibrium thermodynamics, with an activation energy similar to that of individual enzymes but with maximal growth rates and over temperature ranges that are species specific. In this review, we discuss how the temperature dependence of critical cellular processes, such as the central dogma and membrane fluidity, contributes to the temperature dependence of growth. We conclude with a discussion of adaptation to temperature shifts and the effects of temperature on evolution and on the properties of microbial ecosystems.

High-resolution Single-molecule Magnetic Tweezers

Single-molecule magnetic tweezers deliver magnetic force and torque to single target molecules, permitting the study of dynamic changes in biomolecular structures and their interactions. Because the magnetic tweezer setups can generate magnetic fields that vary slowly over tens of millimeters-far larger than the nanometer scale of the single molecule events being observed-this technique can maintain essentially constant force levels during biochemical experiments while generating a biologically meaningful force on the order of 1-100 pN. When using bead-tether constructs to pull on single molecules, smaller magnetic beads and shorter submicrometer tethers improve dynamic response times and measurement precision. In addition, employing high-speed cameras, stronger light sources, and a graphics programming unit permits true high-resolution single-molecule magnetic tweezers that can track nanometer changes in target molecules on a millisecond or even submillisecond time scale. The unique force-clamping capacity of the magnetic tweezer technique provides a way to conduct measurements under near-equilibrium conditions and directly map the energy landscapes underlying various molecular phenomena. High-resolution single-molecule magnetic tweezers can thus be used to monitor crucial conformational changes in single-protein molecules, including those involved in mechanotransduction and protein folding. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

A Rationalization of the Effect That TMAO, Glycine, and Betaine Exert on the Collapse of Elastin-like Polypeptides

Elastin-like polypeptides (ELPs) are soluble in water at low temperature, but, on increasing the temperature, they undergo a reversible and cooperative, coil-to-globule collapse transition. It has been shown that the addition to water of either trimethylamine N-oxide (TMAO), glycine, or betaine causes a significant decrease of T(collapse) in the case of a specific ELP. Traditional rationalizations of these phenomena do not work in the present case. We show that an alternative approach, grounded in the magnitude of the solvent-excluded volume effect and its temperature dependence (strictly linked to the translational entropy of solvent and co-solute molecules), is able to rationalize the occurrence of ELP collapse in water on raising the temperature, as well as the T(collapse) lowering caused by the addition to water of either TMAO, glycine, or betaine.

Novel Enzymes From the Red Sea Brine Pools: Current State and Potential

The Red Sea is a marine environment with unique chemical characteristics and physical topographies. Among the various habitats offered by the Red Sea, the deep-sea brine pools are the most extreme in terms of salinity, temperature and metal contents. Nonetheless, the brine pools host rich polyextremophilic bacterial and archaeal communities. These microbial communities are promising sources for various classes of enzymes adapted to harsh environments - extremozymes. Extremozymes are emerging as novel biocatalysts for biotechnological applications due to their ability to perform catalytic reactions under harsh biophysical conditions, such as those used in many industrial processes. In this review, we provide an overview of the extremozymes from different Red Sea brine pools and discuss the overall biotechnological potential of the Red Sea proteome.

Geothermal heating: Is it a boon or a bane for bioremediation?

There has been a worldwide interest in renewable energy technologies, as a means of decreasing reliance on fossil fuels, minimizing climate change effects, and reducing greenhouse emissions. One such technology is geothermal heating where the constant subsurface temperature is used to cool or heat building interiors via heat pumps. In Canada, the use of geothermal heating has become a popular option for heating and cooling buildings, and it is anticipated that, in the near term, most large buildings will include geothermal heating as part of their climate control strategy. However, little is known about the environmental impacts of geothermal heating on the subsurface environment. The present review will examine the effect of geothermal heating on groundwater flow and remediation efforts, whereby the heat generated by geothermal systems may help with urban pollution. "Geothermal Remediation" could leverage the subsurface heating resulting from geothermal systems to accelerate biodegradation of certain petroleum-based pollutants at brown-field sites, while providing building(s) with sustainable heating and cooling. This idea coincides with the rising momentum towards sustainable and green remediation in Europe and the United States. To ensure that Geothermal Remediation is achievable, the effect of heat on bioremediation needs to be examined. This review provides an insight into what we know about heat effects on bioremediation activities and subsurface transport.

Amylases from thermophilic bacteria: structure and function relationship

Amylases hydrolyze starch to diverse products including dextrins and progressively smaller polymers of glucose units. Thermally stable amylases account for nearly 25% of the enzyme market. This review highlights the structural attributes of the α-amylases from thermophilic bacteria. Heterologous expression of amylases in suitable hosts is discussed in detail. Further, specific value maximization approaches, such as protein engineering and immobilization of the amylases are discussed in order to improve its suitability for varied applications on a commercial scale. The review also takes into account of the immobilization of the amylases on nanomaterials to increase the stability and reusability of the enzymes. The function-based metagenomics would provide opportunities for searching amylases with novel characteristics. The review is expected to explore novel amylases for future potential applications.

Novel insights of waterborne human rotavirus A in Riyadh (Saudi Arabia) involving G2 predominance and emergence of a thermotolerant sequence

The routine evaluation of water environments is necessary to manage enteric virus-mediated fecal contamination and the possible emergence of novel variants. Here, we detected human rotavirus A (HRVA) circulating in two wastewater treatment plants, two lakes, irrigation water and a wastewater landfill located in Riyadh. VP7-derived surface protein sequences were assessed by phylogenetic analyses and inspection of thermotolerance-mediated secondary structure and seasonal variation. HRVA was most prevalent at An-Nazim wastewater landfill (AN-WWLF; 63.89%). Phylogenetic analyzes revealed the predominance of HRVA G2 lineage for the first time in Saudi Arabia. Moreover, a single HRVA sequence (2B64I-ANLF3/2018) was recovered at 45 °C from AN-WWLF; secondary structure prediction indicated that this sequence was thermotolerant with a high hydrophobicity, an absence of Ramachandran outliers, and a higher content of proline patches on the protein surface. Varied relationships were significantly observed between sampling areas influenced by temperature ranges (p < 0.05). HRVA prevalence was influenced by seasonal variations, favoring moderate temperatures in late autumn and early winter in all locations. However, a significant temperature impact was detected in Wadi-Hanifah Lake (p = 0.01). Our study extends the knowledge of currently circulating HRVA genotypes, and indicates the probable emergence of thermotolerant strains and seasonally mediated HRVA prevalence.

Trait dimensions in bacteria and archaea compared to vascular plants

Bacteria and archaea have very different ecology compared to plants. One similarity, though, is that much discussion of their ecological strategies has invoked concepts such as oligotrophy or stress tolerance. For plants, so-called 'trait ecology'-strategy description reframed along measurable trait dimensions-has made global syntheses possible. Among widely measured trait dimensions for bacteria and archaea three main axes are evident. Maximum growth rate in association with rRNA operon copy number expresses a rate-yield trade-off that is analogous to the acquisitive-conservative spectrum in plants, though underpinned by different trade-offs. Genome size in association with signal transduction expresses versatility. Cell size has influence on diffusive uptake and on relative wall costs. These trait dimensions, and potentially others, offer promise for interpreting ecology. At the same time, there are very substantial differences from plant trait ecology. Traits and their underpinning trade-offs are different. Also, bacteria and archaea use a variety of different substrates. Bacterial strategies can be viewed both through the facet of substrate-use pathways, and also through the facet of quantitative traits such as maximum growth rate. Preliminary evidence shows the quantitative traits vary widely within substrate-use pathways. This indicates they convey information complementary to substrate use.

Dynamic cross correlation analysis of Thermus thermophilus alkaline phosphatase and determinants of thermostability

BACKGROUND

Understanding the determinants of protein thermostability is very important both from the theoretical and applied perspective. One emerging view in thermostable enzymes seems to indicate that a salt bridge/charged residue network plays a fundamental role in their thermostability.

METHODS

The structure of alkaline phosphatase (AP) from Thermus thermophilus HB8 was solved by X-ray crystallography at 2.1 Å resolution. The obtained structure was further analyzed by molecular dynamics studies at different temperatures (303 K, 333 K and 363 K) and compared to homologous proteins from the cold-adapted organisms Shewanella sp. and Vibrio strain G15-21. To analyze differences in measures of dynamic variation, several data reduction techniques like principal component analysis (PCA), residue interaction network (RIN) analysis and rotamer analysis were used. Using hierarchical clustering, the obtained results were combined to determine residues showing high degree dynamical variations due to temperature jumps. Furthermore, dynamic cross correlation (DCC) analysis was carried out to characterize networks of charged residues.

RESULTS

Top clustered residues showed a higher propensity for thermostabilizing mutations, indicating evolutionary pressure acting on thermophilic organisms. The description of rotamer distributions by Gini coefficients and Kullback-Leibler (KL) divergence both revealed significant correlations with temperature. DCC analysis revealed a significant trend to de-correlation of the movement of charged residues at higher temperatures.

SIGNIFICANCE

The de-correlation of charged residues detected in Thermus thermophilus AP, highlights the importance of dynamic electrostatic network interactions for the thermostability of this enzyme.

Discovery of Novel Cyclic Salt Bridge in Thermophilic Bacterial Protease and Study of its Sequence and Structure

The plausible explanation behind the stability of thermophilic protein is still yet to be defined more clearly. Here, an in silico study has been undertaken by investigating the sequence and structure of protease from thermophilic (tPro) bacteria and mesophilic (mPro) bacteria. Results showed that charged and uncharged polar residues have higher abundance in tPro. In extreme environment, the tPro is stabilized by high number of isolated and network salt bridges. A novel cyclic salt bridge is also found in a structure of tPro. High number of metal ion-binding site also helps in protein stabilization of thermophilic protease. Aromatic-aromatic interactions also play a crucial role in tPro stabilization. Formation of long network aromatic-aromatic interactions also first time reported here. Finally, the present study provides a major insight with a newly identified cyclic salt bridge in the stability of the enzyme, which may be helpful for protein engineering. It is also used in industrial applications for human welfare.

The Molecular Basis for Life in Extreme Environments

Sampling and genomic efforts over the past decade have revealed an enormous quantity and diversity of life in Earth's extreme environments. This new knowledge of life on Earth poses the challenge of understandingits molecular basis in such inhospitable conditions, given that such conditions lead to loss of structure and of function in biomolecules from mesophiles. In this review, we discuss the physicochemical properties of extreme environments. We present the state of recent progress in extreme environmental genomics. We then present an overview of our current understanding of the biomolecular adaptation to extreme conditions. As our current and future understanding of biomolecular structure-function relationships in extremophiles requires methodologies adapted to extremes of pressure, temperature, and chemical composition, advances in instrumentation for probing biophysical properties under extreme conditions are presented. Finally, we briefly discuss possible future directions in extreme biophysics.

Elevated temperature increases genome-wide selection on de novo mutations

Adaptation in new environments depends on the amount of genetic variation available for evolution, and the efficacy by which natural selection discriminates among this variation. However, whether some ecological factors reveal more genetic variation, or impose stronger selection pressures than others, is typically not known. Here, we apply the enzyme kinetic theory to show that rising global temperatures are predicted to intensify natural selection throughout the genome by increasing the effects of DNA sequence variation on protein stability. We test this prediction by (i) estimating temperature-dependent fitness effects of induced mutations in seed beetles adapted to ancestral or elevated temperature, and (ii) calculate 100 paired selection estimates on mutations in benign versus stressful environments from unicellular and multicellular organisms. Environmental stress per se did not increase mean selection on de novo mutation, suggesting that the cost of adaptation does not generally increase in new ecological settings to which the organism is maladapted. However, elevated temperature increased the mean strength of selection on genome-wide polymorphism, signified by increases in both mutation load and mutational variance in fitness. These results have important implications for genetic diversity gradients and the rate and repeatability of evolution under climate change.

Integrative structure of a 10-megadalton eukaryotic pyruvate dehydrogenase complex from native cell extracts

The pyruvate dehydrogenase complex (PDHc) is a giant enzymatic assembly involved in pyruvate oxidation. PDHc components have been characterized in isolation, but the complex's quaternary structure has remained elusive due to sheer size, heterogeneity, and plasticity. Here, we identify fully assembled Chaetomium thermophilum α-keto acid dehydrogenase complexes in native cell extracts and characterize their domain arrangements utilizing mass spectrometry, activity assays, crosslinking, electron microscopy (EM), and computational modeling. We report the cryo-EM structure of the PDHc core and observe unique features of the previously unknown native state. The asymmetric reconstruction of the 10-MDa PDHc resolves spatial proximity of its components, agrees with stoichiometric data (60 E2p:12 E3BP:∼20 E1p: ≤ 12 E3), and proposes a minimum reaction path among component enzymes. The PDHc shows the presence of a dynamic pyruvate oxidation compartment, organized by core and peripheral protein species. Our data provide a framework for further understanding PDHc and α-keto acid dehydrogenase complex structure and function.

Bayesian genome scale modelling identifies thermal determinants of yeast metabolism

The molecular basis of how temperature affects cell metabolism has been a long-standing question in biology, where the main obstacles are the lack of high-quality data and methods to associate temperature effects on the function of individual proteins as well as to combine them at a systems level. Here we develop and apply a Bayesian modeling approach to resolve the temperature effects in genome scale metabolic models (GEM). The approach minimizes uncertainties in enzymatic thermal parameters and greatly improves the predictive strength of the GEMs. The resulting temperature constrained yeast GEM uncovers enzymes that limit growth at superoptimal temperatures, and squalene epoxidase (ERG1) is predicted to be the most rate limiting. By replacing this single key enzyme with an ortholog from a thermotolerant yeast strain, we obtain a thermotolerant strain that outgrows the wild type, demonstrating the critical role of sterol metabolism in yeast thermosensitivity. Therefore, apart from identifying thermal determinants of cell metabolism and enabling the design of thermotolerant strains, our Bayesian GEM approach facilitates modelling of complex biological systems in the absence of high-quality data and therefore shows promise for becoming a standard tool for genome scale modeling.

While temperature impacts the function of all cellular components, it’s hard to rule out how the temperature dependence of cell phenotypes emerged from the dependence of individual components. Here, the authors develop a Bayesian genome scale modelling approach to identify thermal determinants of yeast metabolism.

Hierarchical Conformational Dynamics Confers Thermal Adaptability to preQ1 RNA Riboswitches

Single-stranded noncoding regulatory RNAs, as exemplified by bacterial riboswitches, are highly dynamic. The conformational dynamics allow the riboswitch to reach maximum switching efficiency under appropriate conditions. Here we characterize the conformational dynamics of preQ₁ riboswitches from mesophilic and thermophilic bacterial species at various temperatures. With the integrative use of small-angle X-ray scattering, NMR, and molecular dynamics simulations, we model the ensemble-structures of the preQ₁ riboswitch aptamers without or with a ligand bound. We show that the preQ₁ riboswitch is sufficiently dynamic and fluctuating among multiple folding intermediates only near the physiological temperature of the microorganism. The hierarchical folding dynamics of the RNA involves the docking of 3'-tail to form a second RNA helix and the helical stacking to form an H-type pseudoknot structure. Further, we show that RNA secondary and tertiary dynamics can be modulated by temperature and by the length of an internal loop. The coupled equilibria between RNA folding intermediates are essential for preQ₁ binding, and a four-state exchange model can account for the change of ligand-triggered switching efficiency with temperature. Together, we have established a relationship between the hierarchical dynamics and riboswitch function, and illustrated how the RNA adapts to high temperature.

Why small proteins tend to have high denaturation temperatures

Data indicate that small globular proteins (consisting of less than about 70 residues) tend to have high denaturation temperatures. This finding is analysed by comparing experimental denaturation enthalpy and entropy changes of a selected set of small proteins with values calculated on the basis of average and common properties of globular proteins. The conclusion is that the denaturation entropy change is smaller than expected, leading to an increase in denaturation temperature. The proposed molecular rationalization considers the existence of long-wavelength, low-frequency vibrational modes in the native state of small proteins due to their large surface-to-interior ratio. The effect of decreasing the conformational entropy gain associated with denaturation on thermal stability is directly verified by means of an already devised theoretical model [G. Graziano, Phys. Chem. Chem. Phys. 2010, 12, 14245-14252; 2014, 16, 21755-21767].

Survival of the cheapest: how proteome cost minimization drives evolution

Darwin's theory of evolution emphasized that positive selection of functional proficiency provides the fitness that ultimately determines the structure of life, a view that has dominated biochemical thinking of enzymes as perfectly optimized for their specific functions. The 20th-century modern synthesis, structural biology, and the central dogma explained the machinery of evolution, and nearly neutral theory explained how selection competes with random fixation dynamics that produce molecular clocks essential e.g. for dating evolutionary histories. However, quantitative proteomics revealed that selection pressures not relating to optimal function play much larger roles than previously thought, acting perhaps most importantly via protein expression levels. This paper first summarizes recent progress in the 21st century toward recovering this universal selection pressure. Then, the paper argues that proteome cost minimization is the dominant, underlying 'non-function' selection pressure controlling most of the evolution of already functionally adapted living systems. A theory of proteome cost minimization is described and argued to have consequences for understanding evolutionary trade-offs, aging, cancer, and neurodegenerative protein-misfolding diseases.

How Do Cells Adapt? Stories Told in Landscapes

Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.

Common mechanism of thermostability in small α- and β-proteins studied by molecular dynamics

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (T_m ) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to T_m . Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH - TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.

Protein Stability in Titan's Subsurface Water Ocean

Models of Titan predict that there is a subsurface ocean of water and ammonia under a layer of ice. Such an ocean would be important in the search for extraterrestrial life since it provides a potentially habitable environment. To evaluate how Earth-based proteins would behave in Titan's subsurface ocean environment, we used molecular dynamics simulations to calculate the properties of proteins with the most common secondary structure types (alpha helix and beta sheet) in both Earth and Titan-like conditions. The Titan environment was simulated by using a temperature of 300 K, a pressure of 1000 bar, and a eutectic mixture of water and ammonia. We analyzed protein compactness, flexibility, and backbone dihedral distributions to identify differences between the two environments. Secondary structures in the Titan environment were found to be less long-lasting, less flexible, and had small differences in backbone dihedral preferences (e.g., in one instance a pi helix formed). These environment-driven differences could lead to changes in how these proteins interact with other biomolecules and therefore changes in how evolution would potentially shape proteins to function in subsurface ocean environments.

Protein thermal stability

Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.

Modeling and Exploiting Microbial Temperature Response

Temperature is an important parameter in bioprocesses, influencing the structure and functionality of almost every biomolecule, as well as affecting metabolic reaction rates. In industrial biotechnology, the temperature is usually tightly controlled at an optimum value. Smart variation of the temperature to optimize the performance of a bioprocess brings about multiple complex and interconnected metabolic changes and is so far only rarely applied. Mathematical descriptions and models facilitate a reduction in complexity, as well as an understanding, of these interconnections. Starting in the 19th century with the “primal” temperature model of Svante Arrhenius, a variety of models have evolved over time to describe growth and enzymatic reaction rates as functions of temperature. Data-driven empirical approaches, as well as complex mechanistic models based on thermodynamic knowledge of biomolecular behavior at different temperatures, have been developed. Even though underlying biological mechanisms and mathematical models have been well-described, temperature as a control variable is only scarcely applied in bioprocess engineering, and as a conclusion, an exploitation strategy merging both in context has not yet been established. In this review, the most important models for physiological, biochemical, and physical properties governed by temperature are presented and discussed, along with application perspectives. As such, this review provides a toolset for future exploitation perspectives of temperature in bioprocess engineering.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes

The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.

Unusually high α-proton acidity of prolyl residues in cyclic peptides

The acidity of the α-proton in peptides has an essential role in numerous biochemical reactions and underpins their stereochemical integrity, which is critical to their biological function. We report a detailed kinetic and computational study of the acidity of the α-proton in two cyclic peptide systems: diketopiperazine (DKP) and triketopiperazine (TKP). The kinetic acidity (protofugality) of the α-protons were determined though hydrogen deuterium exchange studies in aqueous solutions. The acidities of the α-proton in prolyl residues were increased by 3–89 fold relative to other amino acid residues (prolyl > glycyl ≫ alanyl > tyrosyl). Experimental and computational evidence for the stereoelectronic origins of this enhanced prolyl reactivity is presented. TKPs were 10⁶-fold more reactive than their DKP analogues towards deprotonation, which we attribute to the advanced development of aromaticity in the earlier transition state for proton transfer in these cases. A Brønsted linear free energy analysis of the reaction data was conducted to provide estimates of α-proton pK_as.

Kinetic and computational studies reveal that prolyl residues in cyclic peptides are substantially more acidic than other residues due to a stereoelectronic effect.

Finding the generalized molecular principles of protein thermal stability

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

Proteostasis collapse is a driver of cell aging and death

Cells and organisms grow old and die. We develop a biophysical model of the mechanism. Young cells are kept healthy by the positive processes of protein synthesis, degradation, and chaperoning (the activity of keeping proteins properly folded). But, with age, negative processes increase: Oxidative damage accumulates randomly in the cell’s proteins, healthy synthesis and degradation slow down, and—like overfilled garbage cans—chaperone capacity is exceeded. The chaperones are distracted trying to fold irreversibly damaged proteins, leading to accumulating misfolded and aggregated proteins in the cell. The tipping point to death happens when the negative overwhelms the positive. The model makes several quantitative predictions of the life span of the worm Caenorhabditis elegans.

What molecular processes drive cell aging and death? Here, we model how proteostasis—i.e., the folding, chaperoning, and maintenance of protein function—collapses with age from slowed translation and cumulative oxidative damage. Irreparably damaged proteins accumulate with age, increasingly distracting the chaperones from folding the healthy proteins the cell needs. The tipping point to death occurs when replenishing good proteins no longer keeps up with depletion from misfolding, aggregation, and damage. The model agrees with experiments in the worm Caenorhabditis elegans that show the following: Life span shortens nonlinearly with increased temperature or added oxidant concentration, and life span increases in mutants having more chaperones or proteasomes. It predicts observed increases in cellular oxidative damage with age and provides a mechanism for the Gompertz-like rise in mortality observed in humans and other organisms. Overall, the model shows how the instability of proteins sets the rate at which damage accumulates with age and upends a cell’s normal proteostasis balance.

Metal biorecovery and bioremediation: Whether or not thermophilic are better than mesophilic microorganisms

Metal mobilization and immobilization catalyzed by microbial action are key processes in environmental biotechnology. Metal mobilization from ores, mining wastes, or solid residues can be used for recovering metals and/or remediating polluted environments; furthermore, immobilization reduces the migration of metals; cleans up effluents plus ground- and surface water; and, moreover, can help to concentrate and recover metals. Usually these processes provide certain advantages over traditional technologies such as more efficient economical and environmentally sustainable results. Since elevated temperatures typically increase chemical kinetics, it could be expected that bioprocesses should also be enhanced by replacing mesophiles with thermophiles or hyperthermophiles. Nevertheless, other issues like process stability, flexibility, and thermophile-versus-mesophile resistance to acidity and/or metal toxicity should be carefully considered. This review critically analyzes and compares thermophilic and mesophilic microbial performances in recent and selected representative examples of metal bioremediation and biorecovery.

Structure and function of the Ts2631 endolysin of Thermus scotoductus phage vB_Tsc2631 with unique N-terminal extension used for peptidoglycan binding

To escape from hosts after completing their life cycle, bacteriophages often use endolysins, which degrade bacterial peptidoglycan. While mesophilic phages have been extensively studied, their thermophilic counterparts are not well characterized. Here, we present a detailed analysis of the structure and function of Ts2631 endolysin from thermophilic phage vB_Tsc2631, which is a zinc-dependent amidase. The active site of Ts2631 consists of His30, Tyr58, His131 and Cys139, which are involved in Zn²⁺ coordination and catalysis. We found that the active site residues are necessary for lysis yet not crucial for peptidoglycan binding. To elucidate residues involved in the enzyme interaction with peptidoglycan, we tested single-residue substitution variants and identified Tyr60 and Lys70 as essential residues. Moreover, substitution of Cys80, abrogating disulfide bridge formation, inactivates Ts2631, as do substitutions of His31, Thr32 and Asn85 residues. The endolysin contains a positively charged N-terminal extension of 20 residues that can protrude from the remainder of the enzyme and is crucial for peptidoglycan binding. We show that the deletion of 20 residues from the N-terminus abolished the bacteriolytic activity of the enzyme. Because Ts2631 exhibits intrinsic antibacterial activity and unusual thermal stability, it is perfectly suited as a scaffold for the development of antimicrobial agents.

Sequence charge decoration dictates coil-globule transition in intrinsically disordered proteins

We present an analytical theory to compute conformations of heteropolymers-applicable to describe disordered proteins-as a function of temperature and charge sequence. The theory describes coil-globule transition for a given protein sequence when temperature is varied and has been benchmarked against the all-atom Monte Carlo simulation (using CAMPARI) of intrinsically disordered proteins (IDPs). In addition, the model quantitatively shows how subtle alterations of charge placement in the primary sequence-while maintaining the same charge composition-can lead to significant changes in conformation, even as drastic as a coil (swelled above a purely random coil) to globule (collapsed below a random coil) and vice versa. The theory provides insights on how to control (enhance or suppress) these changes by tuning the temperature (or solution condition) and charge decoration. As an application, we predict the distribution of conformations (at room temperature) of all naturally occurring IDPs in the DisProt database and notice significant size variation even among IDPs with a similar composition of positive and negative charges. Based on this, we provide a new diagram-of-states delineating the sequence-conformation relation for proteins in the DisProt database. Next, we study the effect of post-translational modification, e.g., phosphorylation, on IDP conformations. Modifications as little as two-site phosphorylation can significantly alter the size of an IDP with everything else being constant (temperature, salt concentration, etc.). However, not all possible modification sites have the same effect on protein conformations; there are certain "hot spots" that can cause maximal change in conformation. The location of these "hot spots" in the parent sequence can readily be identified by using a sequence charge decoration metric originally introduced by Sawle and Ghosh. The ability of our model to predict conformations (both expanded and collapsed states) of IDPs at a high-throughput level can provide valuable insights into the different mechanisms by which phosphorylation/charge mutation controls IDP function.