Elucidation of factors responsible for enhanced thermal stability of proteins: a structural genomics based study

TemBERTure: advancing protein thermostability prediction with deep learning and attention mechanisms

Motivation

Understanding protein thermostability is essential for numerous biotechnological applications, but traditional experimental methods are time-consuming, expensive, and error-prone. Recently, deep learning (DL) techniques from natural language processing (NLP) was extended to the field of biology, since the primary sequence of proteins can be viewed as a string of amino acids that follow a physicochemical grammar.

Results

In this study, we developed TemBERTure, a DL framework that predicts thermostability class and melting temperature from protein sequences. Our findings emphasize the importance of data diversity for training robust models, especially by including sequences from a wider range of organisms. Additionally, we suggest using attention scores from Deep Learning models to gain deeper insights into protein thermostability. Analyzing these scores in conjunction with the 3D protein structure can enhance understanding of the complex interactions among amino acid properties, their positioning, and the surrounding microenvironment. By addressing the limitations of current prediction methods and introducing new exploration avenues, this research paves the way for more accurate and informative protein thermostability predictions, ultimately accelerating advancements in protein engineering.

Availability and implementation

TemBERTure model and the data are available at: https://github.com/ibmm-unibe-ch/TemBERTure.

Supra-biological performance of immobilized enzymes enabled by chaperone-like specific non-covalent interactions

Designing complex synthetic materials for enzyme immobilization could unlock the utility of biocatalysis in extreme environments. Inspired by biology, we investigate the use of random copolymer brushes as dynamic immobilization supports that enable supra-biological catalytic performance of immobilized enzymes. This is demonstrated by immobilizing Bacillus subtilis Lipase A on brushes doped with aromatic moieties, which can interact with the lipase through multiple non-covalent interactions. Incorporation of aromatic groups leads to a 50 °C increase in the optimal temperature of lipase, as well as a 50-fold enhancement in enzyme activity. Single-molecule FRET studies reveal that these supports act as biomimetic chaperones by promoting enzyme refolding and stabilizing the enzyme's folded and catalytically active state. This effect is diminished when aromatic residues are mutated out, suggesting the importance of π-stacking and π-cation interactions for stabilization. Our results underscore how unexplored enzyme-support interactions may enable uncharted opportunities for using enzymes in industrial biotransformations.

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

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.

Instability Challenges and Stabilization Strategies of Pharmaceutical Proteins

Maintaining the structure of protein and peptide drugs has become one of the most important goals of scientists in recent decades. Cold and thermal denaturation conditions, lyophilization and freeze drying, different pH conditions, concentrations, ionic strength, environmental agitation, the interaction between the surface of liquid and air as well as liquid and solid, and even the architectural structure of storage containers are among the factors that affect the stability of these therapeutic biomacromolecules. The use of genetic engineering, side-directed mutagenesis, fusion strategies, solvent engineering, the addition of various preservatives, surfactants, and additives are some of the solutions to overcome these problems. This article will discuss the types of stress that lead to instabilities of different proteins used in pharmaceutics including regulatory proteins, antibodies, and antibody-drug conjugates, and then all the methods for fighting these stresses will be reviewed. New and existing analytical methods that are used to detect the instabilities, mainly changes in their primary and higher order structures, are briefly summarized.

Structure of the hypothetical protein TTHA1873 from Thermus thermophilus

The crystal structure of an uncharacterized hypothetical protein, TTHA1873 from Thermus thermophilus, has been determined by X-ray crystallography to a resolution of 1.78 Å using the single-wavelength anomalous dispersion method. The protein crystallized as a dimer in two space groups: P43212 and P6122. Structural analysis of the hypothetical protein revealed that the overall fold of TTHA1873 has a β-sandwich jelly-roll topology with nine β-strands. TTHA1873 is a dimeric metal-binding protein that binds to two Ca2+ ions per chain, with one on the surface and the other stabilizing the dimeric interface of the two chains. A structural homology search indicates that the protein has moderate structural similarity to one domain of cell-surface proteins or agglutinin receptor proteins. Red blood cells showed visible agglutination at high concentrations of the hypothetical protein.

Fever as an evolutionary agent to select immune complexes interfaces

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

Insights on Rigidity and Flexibility at the Global and Local Levels of Protein Structures and Their Roles in Homologous Psychrophilic, Mesophilic, and Thermophilic Proteins: A Computational Study

The rigidity and flexibility of homologous psychrophilic (P), mesophilic (M), and thermophilic (T) proteins have been investigated at the global and local levels in terms of "packing factors" and "atomic fluctuations" obtained from B-factors. For comparison of atomic fluctuations, correction of errors by considering errors in B-factors from all sources in a consolidated manner and conversion of the fluctuations to the same temperature have been suggested and validated. The results indicate no differences in the global values like the average packing factor among the three classes of protein homologues, but at local levels there are differences. A comparison of homologous protein triplets show that the average atomic fluctuations at a given temperature mainly obey the order P > M > T. Packing factors and the atomic fluctuations are anti-correlated, suggesting that altering the rigidity of the active site might be a potential strategy to make tailor-made psychrophilic or thermophilic proteins from their mesophilic homologues. The computer codes developed and used in this work are available at the link https://github.com/Munna-Sarkar/proteins-rigidity-flexibility.git.

A novel sequence-based predictor for identifying and characterizing thermophilic proteins using estimated propensity scores of dipeptides

Owing to their ability to maintain a thermodynamically stable fold at extremely high temperatures, thermophilic proteins (TTPs) play a critical role in basic research and a variety of applications in the food industry. As a result, the development of computation models for rapidly and accurately identifying novel TTPs from a large number of uncharacterized protein sequences is desirable. In spite of existing computational models that have already been developed for characterizing thermophilic proteins, their performance and interpretability remain unsatisfactory. We present a novel sequence-based thermophilic protein predictor, termed SCMTPP, for improving model predictability and interpretability. First, an up-to-date and high-quality dataset consisting of 1853 TPPs and 3233 non-TPPs was compiled from published literature. Second, the SCMTPP predictor was created by combining the scoring card method (SCM) with estimated propensity scores of g-gap dipeptides. Benchmarking experiments revealed that SCMTPP had a cross-validation accuracy of 0.883, which was comparable to that of a support vector machine-based predictor (0.906-0.910) and 2-17% higher than that of commonly used machine learning models. Furthermore, SCMTPP outperformed the state-of-the-art approach (ThermoPred) on the independent test dataset, with accuracy and MCC of 0.865 and 0.731, respectively. Finally, the SCMTPP-derived propensity scores were used to elucidate the critical physicochemical properties for protein thermostability enhancement. In terms of interpretability and generalizability, comparative results showed that SCMTPP was effective for identifying and characterizing TPPs. We had implemented the proposed predictor as a user-friendly online web server at http://pmlabstack.pythonanywhere.com/SCMTPP in order to allow easy access to the model. SCMTPP is expected to be a powerful tool for facilitating community-wide efforts to identify TPPs on a large scale and guiding experimental characterization of TPPs.

Aspartate-phobia of thermophiles as a reaction to deleterious chemical transformations

Prokaryotes growing at high temperatures have a high proportion of charged residues in their proteins to stabilize their 3D structure. By mining 175 disparate bacterial and archaeal proteomes we found that, against the general trend for charged residues, the frequency of aspartic acid residues decreases strongly as natural growth temperature increases. In search of the explanation, we hypothesized that the reason for such unusual correlation is the deleterious consequences of spontaneous chemical transformations of aspartate at high temperatures. Our subsequent statistical analysis supported this hypothesis. This finding reveals that organisms have likely adapted to high temperatures by minimizing the harmful consequences of spontaneous chemical transformations.

Comprehensive mutagenesis to identify amino acid residues contributing to the difference in thermostability between two originally thermostable ancestral proteins

Further improvement of the thermostability of inherently thermostable proteins is an attractive challenge because more thermostable proteins are industrially more useful and serve as better scaffolds for protein engineering. To establish guidelines that can be applied for the rational design of hyperthermostable proteins, we compared the amino acid sequences of two ancestral nucleoside diphosphate kinases, Arc1 and Bac1, reconstructed in our previous study. Although Bac1 is a thermostable protein whose unfolding temperature is around 100°C, Arc1 is much more thermostable with an unfolding temperature of 114°C. However, only 12 out of 139 amino acids are different between the two sequences. In this study, one or a combination of amino acid(s) in Bac1 was/were substituted by a residue(s) found in Arc1 at the same position(s). The best mutant, which contained three amino acid substitutions (S108D, G116A and L120P substitutions), showed an unfolding temperature more than 10°C higher than that of Bac1. Furthermore, a combination of the other nine amino acid substitutions also led to improved thermostability of Bac1, although the effects of individual substitutions were small. Therefore, not only the sum of the contributions of individual amino acids, but also the synergistic effects of multiple amino acids are deeply involved in the stability of a hyperthermostable protein. Such insights will be helpful for future rational design of hyperthermostable proteins.

Overview of strategies for developing high thermostability industrial enzymes: Discovery, mechanism, modification and challenges

Biocatalysts such as enzymes are environmentally friendly and have substrate specificity, which are preferred in the production of various industrial products. However, the strict reaction conditions in industry including high temperature, organic solvents, strong acids and bases and other harsh environments often destabilize enzymes, and thus substantially compromise their catalytic functions, and greatly restrict their applications in food, pharmaceutical, textile, bio-refining and feed industries. Therefore, developing industrial enzymes with high thermostability becomes very important in industry as thermozymes have more advantages under high temperature. Discovering new thermostable enzymes using genome sequencing, metagenomics and sample isolation from extreme environments, or performing molecular modification of the existing enzymes with poor thermostability using emerging protein engineering technology have become an effective means of obtaining thermozymes. Based on the thermozymes as biocatalytic chips in industry, this review systematically analyzes the ways to discover thermostable enzymes from extreme environment, clarifies various interaction forces that will affect thermal stability of enzymes, and proposes different strategies to improve enzymes' thermostability. Furthermore, latest development in the thermal stability modification of industrial enzymes through rational design strategies is comprehensively introduced from structure-activity relationship point of view. Challenges and future research perspectives are put forward as well.

The Mutational Robustness of the Genetic Code and Codon Usage in Environmental Context: A Non-Extremophilic Preference?

The genetic code was evolved, to some extent, to minimize the effects of mutations. The effects of mutations depend on the amino acid repertoire, the structure of the genetic code and frequencies of amino acids in proteomes. The amino acid compositions of proteins and corresponding codon usages are still under selection, which allows us to ask what kind of environment the standard genetic code is adapted to. Using simple computational models and comprehensive datasets comprising genomic and environmental data from all three domains of Life, we estimate the expected severity of non-synonymous genomic mutations in proteins, measured by the change in amino acid physicochemical properties. We show that the fidelity in these physicochemical properties is expected to deteriorate with extremophilic codon usages, especially in thermophiles. These findings suggest that the genetic code performs better under non-extremophilic conditions, which not only explains the low substitution rates encountered in halophiles and thermophiles but the revealed relationship between the genetic code and habitat allows us to ponder on earlier phases in the history of Life.

Stability of Engineered Ferritin Nanovaccines Investigated by Combined Molecular Simulation and Experiments

Human ferritin is regarded as an attractive and promising vaccine platform because of its uniform structure, good plasticity, and desirable thermal and chemical stabilities. Besides, it is biocompatible and presumed safe when used as a vaccine carrier. However, there is a lack of knowledge of how different antigen insertion sites on the ferritin nanocage impact the resulting protein stability and performance. To address this question, we selected Epstein-Barr nuclear antigen 1 as a model epitope and fused it at the DNA level with different insertion sites, namely, the N- and C-termini of ferritin, to engineer proteins E1F1 and F1E1, respectively. Protein properties including hydrophobicity and thermal, pH, and chemical stability were investigated both by molecular dynamics (MD) simulation and by experiments. Both methods demonstrate that the insertion site plays an important role in protein properties. The C-terminus insertion (F1E1) leads to a less hydrophobic surface and more tolerance to the external influence of high temperature, pH, and high concentration of chemical denaturants compared to N-terminus insertion (E1F1). Simulated protein hydrophobicity and thermal stability by MD were in high accordance with experimental results. Thus, MD simulation can be used as a valuable tool to engineer nanovaccine candidates, cutting down costs by reducing the experimental effort and accelerating vaccine design.

Thermal and spectroscopic analysis of nitrated compounds and their break-down products using gas chromatography/vacuum UV spectroscopy (GC/VUV)

Gas chromatography/vacuum UV spectroscopy (GC/VUV) was utilized to study various explosives and pharmaceuticals in the nitrate ester and nitramine structural classes. In addition to generating specific VUV spectra for each compound, VUV was used to indicate the onset of thermal decomposition based upon the appearance of break-down products such as nitric oxide, carbon monoxide, formaldehyde, water, and molecular oxygen. The effect of temperature on decomposition could be fit to a logistical function where the fraction of intact compound remaining decreased as the transfer line/flow cell temperature was increased from 200 °C to 300 °C. Utilizing this relationship, the decomposition temperatures for the nitrate ester and nitramine compounds were determined to range between 244 °C and 277 °C. It was also discovered that the decomposition temperature was dependent on the GC carrier gas flow rate and, therefore, the residence time of the compounds in the transfer line/flow cell. For example, the measured decomposition temperature of nitroglycerine ranged from 222 °C to 253 °C across four flow rates. Tracking the appearance/disappearance of decomposition products across this temperature range indicated that NO, CO, and H₂CO are final decomposition products while O₂ and H₂O are intermediate products. The decomposition temperatures for all explosives were highly correlated to similar decomposition measurements taken by differential scanning calorimetry (DSC) (r = 0.91) and thermal gravimetric analysis (TGA) (r = 0.90-0.98). In addition, the decomposition temperatures for all explosives were negatively correlated to the heat of explosion at constant volume (r = -0.68) and strongly positively correlated to the oxygen balance (r = 0.92).

Investigation on the correlation between changes in water and texture properties during the processing of surimi from golden pompano (Trachinotus ovatus)

Water and texture properties are important indicators for evaluating surimi and surimi seafood quality. The present study assessed the impact of surimi processing on the water and texture properties and analyzed their correlations. The results showed that the moisture content exhibited significant positive correlations with T₂₂ and A₂₂ (P < 0.05). However, water-holding capacity (WHC) revealed significant negative correlation with T₂₁ and A₂₂ (P < 0.05). Moreover, correlations between water and texture properties were observed. Moisture content and A₂₂ displayed significant negative correlation with hardness (P < 0.05). WHC exhibited significant positive correlations with cohesiveness, springiness, and chewiness (P < 0.05), whereas T₂₁ exhibited significant negative correlation with those properties (P < 0.05). Springiness is an important indicator for texture properties of surimi gel. T₂₁ and WHC exhibited significant correlations with springiness (P < 0.05). We used two equations to assess the relationship between springiness and T₂₁ (y = 9.64 × x² - 33.79 × x + 29.62, R² = 0.995, P < 0.05) and between springiness and WHC (y = 0.57 + 9.05/{1+e^{[-29.29 × ( x - 83.90)]} }^0.01 , R² = 0.999, P < 0.05), respectively. Verification experiments proved that these two regression equations could be used to predict the change law among feature indicators during the processing of surimi seafood. The present study finds an easy-to-control method to monitor the quality of surimi production and processing for all the stages, and provides ideas for quality design of surimi products. PRACTICAL APPLICATION: The quality of surimi seafood is affected by the process used for surimi production. Early detection of the material quality can efficiently prevent resource wastage. The feature indicators and their correlations help in easy-to-control of all surimi processing stages. Low-field nuclear magnetic resonance (LF-NMR) is a fast, nondestructive method for monitoring food real-time, but the instrument is expensive; whereas water-holding capacity and texture profile analysis (TPA) are traditional methods but time-consuming and material-consuming. Surimi enterprises or researchers choose the right methods based on their needs.

Computation-aided engineering of starch-debranching pullulanase from Bacillus thermoleovorans for enhanced thermostability

Pullulanases are widely used in food, medicine, and other industries because they specifically hydrolyze α-1,6-glycosidic linkages in starch and oligosaccharides. In addition, high-temperature thermostable pullulanase has multiple advantages, including decreasing saccharification solution viscosity accompanied with enhanced mass transfer and reducing microbial contamination in starch hydrolysis. However, thermophilic pullulanase availability remains limited. Additionally, most do not meet starch-manufacturing requirements due to weak thermostability. Here, we developed a computation-aided strategy to engineer the thermophilic pullulanase from Bacillus thermoleovorans. First, three computational design predictors (FoldX, I-Mutant 3.0, and dDFIRE) were combined to predict stability changes introduced by mutations. After excluding conserved and catalytic sites, 17 mutants were identified. After further experimental verification, we confirmed six positive mutants. Among them, the G692M mutant had the highest thermostability improvement, with 3.8 °C increased T_m and 2.1-fold longer half-life than the wild type at 70 °C. We then characterized the mechanism underlying increased thermostability, such as rigidity enhancement, closer conformation, and strengthened motion correlation using root mean square fluctuation (RMSF), principal component analysis (PCA), dynamic cross-correlation map (DCCM), and free energy landscape (FEL) analysis. KEY POINTS: • A computation-aided strategy was developed to engineer pullulanase thermostability. • Seventeen mutants were identified by combining three computational design predictors. • The G692M mutant was obtained with increased T_mand half-life at 70 °C.

Uncovering the evolutionary origin of blue anthocyanins in cereal grains

Functional divergence after gene duplication plays a central role in plant evolution. Among cereals, only Hordeum vulgare (barley), Triticum aestivum (wheat) and Secale cereale (rye) accumulate delphinidin-derived (blue) anthocyanins in the aleurone layer of grains, whereas Oryza sativa (rice), Zea mays (maize) and Sorghum bicolor (sorghum) do not. The underlying genetic basis for this natural occurrence remains elusive. Here, we mapped the barley Blx1 locus involved in blue aleurone to an approximately 1.13 Mb genetic interval on chromosome 4HL, thus identifying a trigenic cluster named MbHF35 (containing HvMYB4H, HvMYC4H and HvF35H). Sequence and expression data supported the role of these genes in conferring blue-coloured (blue aleurone) grains. Synteny analyses across monocot species showed that MbHF35 has only evolved within distinct Triticeae lineages, as a result of dispersed gene duplication. Phylogeny analyses revealed a shared evolution pattern for MbHF35 in Triticeae, suggesting that these genes have co-evolved together. We also identified a Pooideae-specific flavonoid 3',5'-hydroxylase (F3'5'H) lineage, termed here Mo_F35H2, which has a higher amino acid similarity with eudicot F3'5'Hs, demonstrating a scenario of convergent evolution. Indeed, selection tests identified 13 amino acid residues in Mo_F35H2 that underwent positive selection, possibly driven by protein thermostablility selection. Furthermore, through the interrogation of barley germplasm there is evidence that HvMYB4H and HvMYC4H have undergone human selection. Collectively, our study favours blue aleurone as a recently evolved trait resulting from environmental adaptation. Our findings provide an evolutionary explanation for the absence of blue anthocyanins in other cereals and highlight the importance of gene functional divergence for plant diversity and environmental adaptation.

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

Background

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

Results

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

Conclusions

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

Parallel Strategy Increases the Thermostability and Activity of Glutamate Decarboxylase

Glutamate decarboxylase (GAD; EC 4.1.1.15) is a unique pyridoxal 5-phosphate (PLP)-dependent enzyme that specifically catalyzes the decarboxylation of L-glutamic acid to produce γ-aminobutyric acid (GABA), which exhibits several well-known physiological functions. However, glutamate decarboxylase from different sources has the common problem of poor thermostability that affects its application in industry. In this study, a parallel strategy comprising sequential analysis and free energy calculation was applied to identify critical amino acid sites affecting thermostability of GAD and select proper mutation contributing to improve structure rigidity of the enzyme. Two mutant enzymes, D203E and S325A, with higher thermostability were obtained, and their semi-inactivation temperature (T₅₀¹⁵) values were 2.3 °C and 1.4 °C higher than the corresponding value of the wild-type enzyme (WT), respectively. Moreover, the mutant, S325A, exhibited enhanced activity compared to the wild type, with a 1.67-fold increase. The parallel strategy presented in this work proved to be an efficient tool for the reinforcement of protein thermostability.

Extremophile - An Adaptive Strategy for Extreme Conditions and Applications

The concurrence of microorganisms in niches that are hostile like extremes of temperature, pH, salt concentration and high pressure depends upon novel molecular mechanisms to enhance the stability of their proteins, nucleic acids, lipids and cell membranes. The structural, physiological and genomic features of extremophiles that make them capable of withstanding extremely selective environmental conditions are particularly fascinating. Highly stable enzymes exhibiting several industrial and biotechnological properties are being isolated and purified from these extremophiles. Successful gene cloning of the purified extremozymes in the mesophilic hosts has already been done. Various extremozymes such as amylase, lipase, xylanase, cellulase and protease from thermophiles, halothermophiles and psychrophiles are of industrial interests due to their enhanced stability at forbidding conditions. In this review, we made an attempt to point out the unique features of extremophiles, particularly thermophiles and psychrophiles, at the structural, genomic and proteomic levels, which allow for functionality at harsh conditions focusing on the temperature tolerance by them.

Amino acid exchangeabilities vary across the tree of life

Different amino acid pairs have drastically different relative exchangeabilities (REs), and accounting for this variation is an important and common practice in inferring phylogenies, testing selection, and predicting mutational effects, among other analyses. In all such endeavors, REs have been generally considered invariant among species; this assumption, however, has not been scrutinized. Using maximum likelihood to analyze 180 genome sequences, we estimated REs from 90 clades representing all three domains of life, and found numerous instances of substantial between-clade differences in REs. REs show more differences between orthologous proteins of different clades than unrelated proteins of the same clade, suggesting that REs are genome-wide, clade-specific features, probably a result of proteome-wide evolutionary changes in the physicochemical environments of amino acid residues. The discovery of among-clade RE variations cautions against assuming constant REs in various analyses and demonstrates a higher-than-expected complexity in mechanisms of proteome evolution.

Structure of a hyperthermostable dimeric archaeal Rubisco fromHyperthermus butylicus

The crystal structure of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) from the hyperthermophilic archaeonHyperthermus butylicusis presented at 1.8 Å resolution. Previous structures of archaeal Rubisco have been found to assemble into decamers, and this oligomerization was thought to be required for a highly thermally stable enzyme. In the current study,H. butylicusRubisco is shown to exist as a dimer in solution, yet has a thermal denaturation midpoint of 114°C, suggesting that high thermal stability can be achieved without an increased oligomeric state. This increased thermal stability appears to be due to an increased number of electrostatic interactions within the monomeric subunit. As such,H. butylicusRubisco presents a well characterized system in which to investigate the role of assembly and thermal stability in enzyme function.

Investigation of the thermophilic mechanism in the genus Porphyrobacter by comparative genomic analysis

Background

Type strains of the genus Porphyrobacter belonging to the family Erythrobacteraceae and the class Alphaproteobacteria have been isolated from various environments, such as swimming pools, lake water and hot springs. P. cryptus DSM 12079^T and P. tepidarius DSM 10594^T out of all Erythrobacteraceae type strains, are two type strains that have been isolated from geothermal environments. Next-generation sequencing (NGS) technology offers a convenient approach for detecting situational types based on protein sequence differences between thermophiles and mesophiles; amino acid substitutions can lead to protein structural changes, improving the thermal stabilities of proteins. Comparative genomic studies have revealed that different thermal types exist in different taxa, and few studies have been focused on the class Alphaproteobacteria, especially the family Erythrobacteraceae. In this study, eight genomes of Porphyrobacter strains were compared to elucidate how Porphyrobacter thermophiles developed mechanisms to adapt to thermal environments.

Results

P. cryptus DSM 12079^T grew optimally at 50 °C, which was higher than the optimal growth temperature of other Porphyrobacter type strains. Phylogenomic analysis of the genus Porphyrobacter revealed that P. cryptus DSM 12079^T formed a distinct and independent clade. Comparative genomic studies uncovered that 1405 single-copy genes were shared by Porphyrobacter type strains. Alignments of single-copy proteins showed that various types of amino acid substitutions existed between P. cryptus DSM 12079^T and the other Porphyrobacter strains. The primary substitution types were changes from glycine/serine to alanine.

Conclusions

P. cryptus DSM 12079^T was the sole thermophile within the genus Porphyrobacter. Phylogenomic analysis and amino acid frequencies indicated that amino acid substitutions might play an important role in the thermophily of P. cryptus DSM 12079^T. Bioinformatic analysis revealed that major amino acid substitutional types, such as changes from glycine/serine to alanine, increase the frequency of α-helices in proteins, promoting protein thermostability in P. cryptus DSM 12079^T. Hence, comparative genomic analysis broadens our understanding of thermophilic mechanisms in the genus Porphyrobacter and may provide a useful insight in the design of thermophilic enzymes for agricultural, industrial and medical applications.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4789-4) contains supplementary material, which is available to authorized users.

Molecular dynamics perspective on the thermal stability of mandelate racemase

Mandelate racemase from Pseudomonas putida is a promising candidate for the dynamic kinetic resolution of α-hydroxy carboxylic acids. In the present study, the thermal stability of mandelate racemase was investigated through molecular dynamics simulations in the temperature range of 303-363 K, which can guide the design of mandelate racemase with higher stability. The basic features such as radius of gyration, surface accessibility, and secondary structure content suggested the instability of mandelate racemase at high temperatures. With increase in temperature, α-helix content reduced significantly, especially the α-helices exposed to the environment. At the simulation time scale considered, intra-protein hydrogen bonds, hydrogen bonds between protein and water decreased at 363 K, while the number of salt-bridges increased. The long-distance networks remarkably changed at 363 K. A considerable number of long-lived (percentage existence time higher than 90%) hydrogen bonds and Cα contacts were lost. Root mean square fluctuation analysis revealed regions with high fluctuation, which should be helpful in the reengineering of mandelate racemase for enhanced thermal stability.

Functional and cooperative stabilization of a two-metal (Ca, Zn) center in α-amylase derived from Flavobacteriaceae species

Mesophilic α-amylase from Flavobacteriaceae (FSA) is evolutionary closely related to thermophilic archaeal Pyrococcus furiosus α-amylase (PWA), but lacks the high thermostability, despite the conservation of most residues involved in the two-metal (Ca, Zn) binding center of PWA. In this study, a disulfide bond was introduced near the two-metal binding center of FSA (designated mutant EH-CC) and this modification resulted in a slight improvement in thermostability. As expected, E204G mutations in FSA and EH-CC led to the recovery of Ca²⁺-binding site. Interestingly, both Ca²⁺- and Zn²⁺-dependent thermostability were significantly enhanced; 153.1% or 50.8% activities was retained after a 30-min incubation period at 50 °C, in the presence of Ca²⁺ or Zn²⁺. The C214S mutation, which affects Zn²⁺-binding, also remarkably enhanced Zn²⁺- and Ca²⁺- dependent thermostability, indicating that Ca²⁺- and Zn²⁺-binding sites function cooperatively to maintain protein stability. Furthermore, an isothermal titration calorimetry (ITC) analysis revealed a novel Zn²⁺-binding site in mutant EH-CC-E204G. This metal ion cooperation provides a possible method for the generation of α-amylases with desired thermal properties by in silico rational design and systems engineering, to generate a Zn²⁺-binding site adjacent to the conserved Ca²⁺-binding site.

Structural signatures of thermal adaptation of bacterial ribosomal RNA, transfer RNA, and messenger RNA

Temperature adaptation of bacterial RNAs is a subject of both fundamental and practical interest because it will allow a better understanding of molecular mechanism of RNA folding with potential industrial application of functional thermophilic or psychrophilic RNAs. Here, we performed a comprehensive study of rRNA, tRNA, and mRNA of more than 200 bacterial species with optimal growth temperatures (OGT) ranging from 4°C to 95°C. We investigated temperature adaptation at primary, secondary and tertiary structure levels. We showed that unlike mRNA, tRNA and rRNA were optimized for their structures at compositional levels with significant tertiary structural features even for their corresponding randomly permutated sequences. tRNA and rRNA are more exposed to solvent but remain structured for hyperthermophiles with nearly OGT-independent fluctuation of solvent accessible surface area within a single RNA chain. mRNA in hyperthermophiles is essentially the same as random sequences without tertiary structures although many mRNA in mesophiles and psychrophiles have well-defined tertiary structures based on their low overall solvent exposure with clear separation of deeply buried from partly exposed bases as in tRNA and rRNA. These results provide new insight into temperature adaptation of different RNAs.

Histidine substitution in the most flexible fragments of firefly luciferase modifies its thermal stability

Molecular dynamics (MD) at two temperatures of 300 and 340 K identified two histidine residues, His461 and His489, in the most flexible regions of firefly luciferase, a light emitting enzyme. We therefore designed four protein mutants H461D, H489K, H489D and H489M to investigate their enzyme kinetic and thermodynamic stability changes. Substitution of His461 by aspartate (H461D) decreased ATP binding affinity, reduced the melting temperature of protein by around 25 °C and shifted its optimum temperature of activity to 10 °C. In line with the common feature of psychrophilic enzymes, the MD data showed that the overall flexibility of H461D was relatively high at low temperature, probably due to a decrease in the number of salt bridges around the mutation site. On the other hand, substitution of His489 by aspartate (H489D) introduced a new salt bridge between the C-terminal and N-terminal domains and increased protein rigidity but only slightly improved its thermal stability. Similar changes were observed for H489K and, to a lesser degree, H489M mutations. Based on our results we conclude that the MD simulation-based rational substitution of histidines by salt-bridge forming residues can modulate conformational dynamics in luciferase and shift its optimal temperature activity.

Engineering and introduction of de novo disulphide bridges in organophosphorus hydrolase enzyme for thermostability improvement

The organophosphorus hydrolase (OPH) has been used to degrade organophosphorus chemicals, as one of the most frequently used decontamination methods. Under chemical and thermal denaturing conditions, the enzyme has been shown to unfold. To utilize this enzyme in various applications, the thermal stability is of importance. The engineering of de novo disulphide bridges has been explored as a means to increase the thermal stability of enzymes in the rational method of protein engineering. In this study, Disulphide by Design software, homology modelling and molecular dynamics simulations were used to select appropriate amino acid pairs for the introduction of disulphide bridge to improve protein thermostability. The thermostability of the wild-type and three selected mutant enzymes were evaluated by half-life, delta G inactivation (ΔGi) and structural studies (fluorescence and far-UV CD analysis). Data analysis showed that half-life of A204C/T234C and T128C/E153C mutants were increased up to 4 and 24 min, respectively; however, for the G74C/A78C mutant, the half-life was decreased up to 9 min. For the T128C/E124C mutant, both thermal stability and Catalytic efficiency (kcat) were also increased. The half-life and ΔGi results were correlated to the obtained information from structural studies by circular dichroism (CD) spectrometry and extrinsic fluorescence experiments; as rigidity increased in A204C/T2234C and T128C/E153C mutants, half-life and ΔGi also increased. For G74C/A78C mutant, these parameters decreased due to its higher flexibility. The results were submitted a strong evidence for the possibility to improve the thermostability of OPH enzyme by introducing a disulphide bridge after bioinformatics design, even though this design would not be always successful.

Multifactorial level of extremostability of proteins: can they be exploited for protein engineering?

Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.

Establishing knowledge on the sequence arrangement pattern of nucleated protein folding

The heat-tolerance mechanisms of (hyper)thermophilic proteins provide a unique opportunity to investigate the unsolved protein folding problem. In an attempt to determine whether the interval between residues in sequence might play a role in determining thermostability, we constructed a sequence interval-dependent value function to calculate the residue pair frequency. Additionally, we identified a new sequence arrangement pattern, where like-charged residues tend to be adjacently assembled, while unlike-charged residues are distributed over longer intervals, using statistical analysis of a large sequence database. This finding indicated that increasing the intervals between unlike-charged residues can increase protein thermostability, with the arrangement patterns of these charged residues serving as thermodynamically favorable nucleation points for protein folding. Additionally, we identified that the residue pairs K-E, R-E, L-V and V-V involving long sequence intervals play important roles involving increased protein thermostability. This work demonstrated a novel approach for considering sequence intervals as keys to understanding protein folding. Our findings of novel relationships between residue arrangement and protein thermostability can be used in industry and academia to aid the design of thermostable proteins.

Improved insights into protein thermal stability: from the molecular to the structurome scale

Despite the intense efforts of the last decades to understand the thermal stability of proteins, the mechanisms responsible for its modulation still remain debated. In this investigation, we tackle this issue by showing how a multiscale perspective can yield new insights. With the help of temperature-dependent statistical potentials, we analysed some amino acid interactions at the molecular level, which are suggested to be relevant for the enhancement of thermal resistance. We then investigated the thermal stability at the protein level by quantifying its modification upon amino acid substitutions. Finally, a large scale analysis of protein stability-at the structurome level-contributed to the clarification of the relation between stability and natural evolution, thereby showing that the mutational profile of proteins differs according to their thermal properties. Some considerations on how the multiscale approach could help in unravelling the protein stability mechanisms are briefly discussed.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Improvement of the thermostability and catalytic efficiency of a highly active β-glucanase from Talaromyces leycettanus JCM12802 by optimizing residual charge-charge interactions

BACKGROUND

β-Glucanase is one of the most extensively used biocatalysts in biofuel, food and animal feed industries. However, the poor thermostability and low catalytic efficiency of most reported β-glucanases limit their applications. Currently, two strategies are used to overcome these bottlenecks, i.e., mining for novel enzymes from extremophiles and engineering existing enzymes.

RESULTS

A novel endo-β-1,3-1,4-glucanase of GH16 (Tlglu16A) from the thermophilic fungus Talaromyces leycettanus JCM12802 was produced in Pichia pastoris and characterized. For potential industrial applications, recombinant TlGlu16A exhibits favorable enzymatic properties over most reported glucanases, i.e., remarkable stability over a wide pH range from 1.0 to 10.0 and superior activity on glucan substrates (up to 15,197 U/mg). The only weakness of TlGlu16A is the thermolability at 65 °C and higher. To improve the thermostability, the enzyme thermal stability system was then used to engineer TlGlu16A through optimization of residual charge-charge interactions. Eleven mutants were constructed and compared to the wild-type TlGlu16A. Four mutants, H58D, E134R, D235G and D296K, showed longer half-life time at 80 °C (31, 7, 25, 22 vs. 0.5 min), and two mutants, D235G and D296K, had greater specific activities (158.2 and 122.2 %, respectively) and catalytic efficiencies (k cat/K m, 170 and 114 %, respectively).

CONCLUSIONS

The engineered TlGlu16A has great application potentials from the perspectives of enzyme yield and properties. Its thermostability and activity were apparently improved in the engineered enzymes through charge optimization. This study spans the genetic, functional and structural fields, and provides a combination of gene mining and protein engineering approaches for the systematic improvement of enzyme performance.

Predicting protein thermal stability changes upon point mutations using statistical potentials: Introducing HoTMuSiC

The accurate prediction of the impact of an amino acid substitution on the thermal stability of a protein is a central issue in protein science, and is of key relevance for the rational optimization of various bioprocesses that use enzymes in unusual conditions. Here we present one of the first computational tools to predict the change in melting temperature ΔT_m upon point mutations, given the protein structure and, when available, the melting temperature T_m of the wild-type protein. The key ingredients of our model structure are standard and temperature-dependent statistical potentials, which are combined with the help of an artificial neural network. The model structure was chosen on the basis of a detailed thermodynamic analysis of the system. The parameters of the model were identified on a set of more than 1,600 mutations with experimentally measured ΔT_m. The performance of our method was tested using a strict 5-fold cross-validation procedure, and was found to be significantly superior to that of competing methods. We obtained a root mean square deviation between predicted and experimental ΔT_m values of 4.2 °C that reduces to 2.9 °C when ten percent outliers are removed. A webserver-based tool is freely available for non-commercial use at soft.dezyme.com.

A Miniature Protein Stabilized by a Cation-π Interaction Network

The design of folded miniature proteins is predicated on establishing noncovalent interactions that direct the self-assembly of discrete thermostable tertiary structures. In this work, we describe how a network of cation-π interactions present in proteins containing "WSXWS motifs" can be emulated to stabilize the core of a miniature protein. This 19-residue protein sequence recapitulates a set of interdigitated arginine and tryptophan residues that stabilize a distinctive β-strand:loop:PPII-helix topology. Validation of the compact fold determined by NMR was carried out by mutagenesis of the cation-π network and by comparison to the corresponding disulfide-bridged structure. These results support the involvement of a coordinated set of cation-π interactions that stabilize the tertiary structure.

Origins and Early Evolution of the tRNA Molecule

Modern transfer RNAs (tRNAs) are composed of ~76 nucleotides and play an important role as "adaptor" molecules that mediate the translation of information from messenger RNAs (mRNAs). Many studies suggest that the contemporary full-length tRNA was formed by the ligation of half-sized hairpin-like RNAs. A minihelix (a coaxial stack of the acceptor stem on the T-stem of tRNA) can function both in aminoacylation by aminoacyl tRNA synthetases and in peptide bond formation on the ribosome, indicating that it may be a vestige of the ancestral tRNA. The universal CCA-3' terminus of tRNA is also a typical characteristic of the molecule. "Why CCA?" is the fundamental unanswered question, but several findings give a comprehensive picture of its origin. Here, the origins and early evolution of tRNA are discussed in terms of various perspectives, including nucleotide ligation, chiral selectivity of amino acids, genetic code evolution, and the organization of the ribosomal peptidyl transferase center (PTC). The proto-tRNA molecules may have evolved not only as adaptors but also as contributors to the composition of the ribosome.

Towards an accurate prediction of the thermal stability of homologous proteins

Protein thermostability has been the focus of growing research interests in the last decades since its understanding and control play important roles in the optimization of a wide series of bioprocesses of academic and industrial importance. The complexity of this issue is rooted in the fact that the mechanisms ensuring thermal resistance are not unique and specific, but rather family- or even protein-dependent. Therefore, and despite the amount of research already accomplished, obtaining fast and precise thermal stability predictions is still a challenge, especially on a large scale. This article deepens the study of protein thermal stability and is focused on the prediction of its best descriptor, the melting temperature Tm. The relations between Tm and a series of factors that are expected to influence the protein stability are analyzed and discussed. Different Tm-prediction methods that utilize these factors, sometimes with additional information about homologous proteins, are introduced, and their individual performances are evaluated. The best methods are based on temperature-dependent statistical potentials, on the environmental temperature of the host organism, on the fraction of charged residues, and on the number of residues. They are combined to build an improved prediction method with significantly increased score. The root mean square deviation between the computed and experimental Tm-values for 45 proteins of known structure from 11 families is about 7°C in cross-validation and decreases to 5°C when 10% outliers are removed. The associated linear correlation coefficients are equal to .91 and .95, respectively.

Improvement in Thermostability of an Achaetomium sp. Strain Xz8 Endopolygalacturonase via the Optimization of Charge-Charge Interactions

Improving enzyme thermostability is of importance for widening the spectrum of application of enzymes. In this study, a structure-based rational design approach was used to improve the thermostability of a highly active, wide-pH-range-adaptable, and stable endopolygalacturonase (PG8fn) from Achaetomium sp. strain Xz8 via the optimization of charge-charge interactions. By using the enzyme thermal stability system (ETSS), two residues--D244 and D299--were inferred to be crucial contributors to thermostability. Single (D244A and D299R) and double (D244A/D299R) mutants were then generated and compared with the wild type. All mutants showed improved thermal properties, in the order D244A < D299R < D244A/D299R. In comparison with PG8fn, D244A/D299R showed the most pronounced shifts in temperature of maximum enzymatic activity (Tmax), temperature at which 50% of the maximal activity of an enzyme is retained (T50), and melting temperature (Tm), of about 10, 17, and 10.2°C upward, respectively, with the half-life (t1/2) extended by 8.4 h at 50°C and 45 min at 55°C. Another distinguishing characteristic of the D244A/D299R mutant was its catalytic activity, which was comparable to that of the wild type (23,000 ± 130 U/mg versus 28,000 ± 293 U/mg); on the other hand, it showed more residual activity (8,400 ± 83 U/mg versus 1,400 ± 57 U/mg) after the feed pelleting process (80°C and 30 min). Molecular dynamics (MD) simulation studies indicated that mutations at sites D244 and D299 lowered the overall root mean square deviation (RMSD) and consequently increased the protein rigidity. This study reveals the importance of charge-charge interactions in protein conformation and provides a viable strategy for enhancing protein stability.

Functional adaptations of the bacterial chaperone trigger factor to extreme environmental temperatures

Trigger factor (TF) is the first molecular chaperone interacting cotranslationally with virtually all nascent polypeptides synthesized by the ribosome in bacteria. Thermal adaptation of chaperone function was investigated in TFs from the Antarctic psychrophile Pseudoalteromonas haloplanktis, the mesophile Escherichia coli and the hyperthermophile Thermotoga maritima. This series covers nearly all temperatures encountered by bacteria. Although structurally homologous, these TFs display strikingly distinct properties that are related to the bacterial environmental temperature. The hyperthermophilic TF strongly binds model proteins during their folding and protects them from heat-induced misfolding and aggregation. It decreases the folding rate and counteracts the fast folding rate imposed by high temperature. It also functions as a carrier of partially folded proteins for delivery to downstream chaperones ensuring final maturation. By contrast, the psychrophilic TF displays weak chaperone activities, showing that these functions are less important in cold conditions because protein folding, misfolding and aggregation are slowed down at low temperature. It efficiently catalyses prolyl isomerization at low temperature as a result of its increased cellular concentration rather than from an improved activity. Some chaperone properties of the mesophilic TF possibly reflect its function as a cold shock protein in E. coli.