Stabilizing salt-bridge enhances protein thermostability by reducing the heat capacity change of unfolding

PEGylation versus glycosylation: effect on the thermodynamics and thermostability of crisantaspase

Thermostability is an important and desired feature of therapeutic proteins and is critical for the success or failure of protein drugs development. It can be increased by PEGylation-binding of poly(ethylene glycol) moieties-or glycosylation-post-translational modification to add glycans. Here, the thermostability and thermodynamic parameters of native, PEGylated, and glycosylated versions of the antileukemic enzyme crisantaspase were investigated. First-order kinetics was found to describe the irreversible deactivation process. Activation energy of the enzyme-catalyzed reaction (E*) was estimated for native, PEGylated, and glycosylated enzyme (10.2, 14.8, and 18.8 kJ mol-1 respectively). Half-life decreased progressively with increasing temperature, and longer half-life was observed for PEG-crisantaspase (87.74 min) at 50 °C compared to the native form (9.79 min). The activation energy of denaturation of PEG-crisantaspase (307.1 kJ mol-1) was higher than for crisantaspase (218.1 kJ mol-1) and Glyco-crisantaspase (120.0 kJ mol-1), which means that more energy is required to overcome the energy barrier of the unfolding process. According to our results, PEG-crisantaspase is more thermostable than its native form, while Glyco-crisantaspase is more thermosensitive.

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.

Stability of collagen heterotrimer with same charge pattern and different charged residue identities

Salt bridges are important factors in maintaining the stability of proteins, and their contribution to protein folding has received much attention. Although the interaction energies, or stabilizing contributions, of individual salt bridges have been measured in various proteins, a systematic assessment of various types of salt bridges in a relatively uniform environment is still a valuable analysis. Here, we used a collagen heterotrimer as a host-guest platform to construct 48 heterotrimers with the same charge pattern. A variety of salt bridges were formed between the oppositely charged residues Lys, Arg, Asp, and Glu. The melting temperature (Tm) of the heterotrimers was measured with circular dichroism. The atomic structures of 10 salt bridges were shown in three x-ray crystals of heterotrimer. Molecular dynamics simulation based on the crystal structures indicated that strong, intermediate, and weak salt bridges have distinctive N-O distances. A linear regression model was used to predict the stability of heterotrimers with high accuracy (R2 = 0.93). We developed an online database to help readers understand how a salt bridge stabilizes collagen. This work will help us better understand the stabilizing mechanism of salt bridges in collagen folding and provide a new strategy to design collagen heterotrimers.

Comparative Modeling and Analysis of Extremophilic D-Ala-D-Ala Carboxypeptidases

Understanding the molecular adaptations of organisms to extreme environments requires a comparative analysis of protein structure, function, and dynamics across species found in different environmental conditions. Computational studies can be particularly useful in this pursuit, allowing exploratory studies of large numbers of proteins under different thermal and chemical conditions that would be infeasible to carry out experimentally. Here, we perform such a study of the MEROPS family S11, S12, and S13 proteases from psychophilic, mesophilic, and thermophilic bacteria. Using a combination of protein structure prediction, atomistic molecular dynamics, and trajectory analysis, we examine both conserved features and trends across thermal groups. Our findings suggest a number of hypotheses for experimental investigation.

Computational analysis of the effect of Gly100Ala mutation on the thermostability of SazCA

Maintaining the protein stability upon mutation is a challenging task in protein engineering. In the present computational study, we induced a single point Gly100Ala mutation in SazCA and examined the factors governing the stability and flexibility of the mutated form, and compared it to that of the wildtype using molecular dynamics simulations. We observed higher structural stability and lesser residual mobility in the mutated SazCA. Improved H-bonding due to Gly100Ala was observed. Ala100 was responsible for the increased helical contents in the mutated SazCA while Gly100 compromised the secondary structure contents in the wildtype. A strong network of salt bridges and high local ordering of the solvent molecules at the protein surface contributed to the enhanced stability of the mutated protein. Our simulations conclusively highlight Gly100Ala mutation as a step towards designing a more robust and thermostable SazCA.Communicated by Ramaswamy H. Sarma.

Thermo-L-Asparaginases: From the Role in the Viability of Thermophiles and Hyperthermophiles at High Temperatures to a Molecular Understanding of Their Thermoactivity and Thermostability

L-asparaginase (L-ASNase) is a vital enzyme with a broad range of applications in medicine, food industry, and diagnostics. Among various organisms expressing L-ASNases, thermophiles and hyperthermophiles produce enzymes with superior performances-stable and heat resistant thermo-ASNases. This review is an attempt to take a broader view on the thermo-ASNases. Here we discuss the position of thermo-ASNases in the large family of L-ASNases, their role in the heat-tolerance cellular system of thermophiles and hyperthermophiles, and molecular aspects of their thermoactivity and thermostability. Different types of thermo-ASNases exhibit specific L-asparaginase activity and additional secondary activities. All products of these enzymatic reactions are associated with diverse metabolic pathways and are important for mitigating heat stress. Thermo-ASNases are quite distinct from typical mesophilic L-ASNases based on structural properties, kinetic and activity profiles. Here we attempt to summarize the current understanding of the molecular mechanisms of thermo-ASNases' thermoactivity and thermostability, from amino acid composition to structural-functional relationships. Research of these enzymes has fundamental and biotechnological significance. Thermo-ASNases and their improved variants, cloned and expressed in mesophilic hosts, can form a large pool of enzymes with valuable characteristics for biotechnological application.

Understanding the Role of Intramolecular Ion-Pair Interactions in Conformational Stability Using an Ab Initio Thermodynamic Cycle

Intramolecular ion-pair interactions yield shape and functionality to many molecules. With proper orientation, these interactions overcome steric factors and are responsible for the compact structures of several peptides. In this study, we present a thermodynamic cycle based on isoelectronic and alchemical mutation to estimate the intramolecular ion-pair interaction energy. We determine these energies for 26 benchmark molecules with common ion-pair combinations and compare them with results obtained using intramolecular symmetry-adapted perturbation theory. For systems with long linkers, the ion-pair energies evaluated using both approaches deviate by less than 2.5% in the vacuum phase. The thermodynamic cycle based on density functional theory facilitates calculations of salt-bridge interactions in model tripeptides with continuum/microsolvation modeling and four large peptides: 1EJG (crambin), 1BDK (bradykinin), 1L2Y (a mini-protein with a tryptophan cage), and 1SCO (a toxin from the scorpion venom).

Identification and Mutation Analysis of Nonconserved Residues on the TIM-Barrel Surface of GH5_5 Cellulases for Catalytic Efficiency and Stability Improvement

Exploring the potential functions of nonconserved residues on the outer side of α-helices and systematically optimizing them are pivotal for their application in protein engineering. Based on the evolutionary structural conservation analysis of GH5_5 cellulases, a practical molecular improvement strategy was developed. Highly variable sites on the outer side of the α-helices of the GH5_5 cellulase from Aspergillus niger (AnCel5A) were screened, and 14 out of the 34 highly variable sites were confirmed to exert a positive effect on the activity. After the modular combination of the positive mutations, the catalytic efficiency of the mutants was further improved. By using CMC-Na as the substrate, the catalytic efficiency and specific activity of variant AnCel5A_N193A/T300P/D307P were approximately 2.0-fold that of AnCel5A (227 ± 21 versus 451 ± 43 ml/s/mg and 1,726 ± 19 versus 3,472 ± 42 U/mg, respectively). The half-life (t_1/2) of variant AnCel5A_N193A/T300P/D307P at 75°C was 2.36 times that of AnCel5A. The role of these sites was successfully validated in other GH5_5 cellulases. Computational analyses revealed that the flexibility of the loop 6-loop 7-loop 8 region was responsible for the increased catalytic performance. This work not only illustrated the important role of rapidly evolving positions on the outer side of the α-helices of GH5_5 cellulases but also revealed new insights into engineering the proteins that nature left as clues for us to find. IMPORTANCE A comprehensive understanding of the residues on the α-helices of the GH5_5 cellulases is important for catalytic efficiency and stability improvement. The main objective of this study was to use the evolutionary conservation and plasticity of the TIM-barrel fold to probe the relationship between nonconserved residues on the outer side of the α-helices and the catalytic efficiency of GH5_5 cellulases by conducting structure-guided protein engineering. By using a four-step nonconserved residue screening strategy, the functional role of nonconserved residues on the outer side of the α-helices was effectively identified, and a variant with superior performance and capability was constructed. Hence, this study proved the effectiveness of this strategy in engineering GH5_5 cellulases and provided a potential competitor for industrial applications. Furthermore, this study sheds new light on engineering TIM-barrel proteins.

A Statistical Analysis of the Sequence and Structure of Thermophilic and Non-Thermophilic Proteins

Thermophilic proteins have various practical applications in theoretical research and in industry. In recent years, the demand for thermophilic proteins on an industrial scale has been increasing; therefore, the engineering of thermophilic proteins has become a hot direction in the field of protein engineering. However, the exact mechanism of thermostability of proteins is not yet known, for engineering thermophilic proteins knowing the basis of thermostability is necessary. In order to understand the basis of the thermostability in proteins, we have made a statistical analysis of the sequences, secondary structures, hydrogen bonds, salt bridges, DHA (Donor-Hydrogen-Accepter) angles, and bond lengths of ten pairs of thermophilic proteins and their non-thermophilic orthologous. Our findings suggest that polar amino acids contribute to thermostability in proteins by forming hydrogen bonds and salt bridges which provide resistance against protein denaturation. Short bond length and a wider DHA angle provide greater bond stability in thermophilic proteins. Moreover, the increased frequency of aromatic amino acids in thermophilic proteins contributes to thermal stability by forming more aromatic interactions. Additionally, the coil, helix, and loop in the secondary structure also contribute to thermostability.

A comprehensive in silico exploration of the impacts of missense variants on two different conformations of human pirin protein

BACKGROUND

Pirin, a member of the cupin superfamily, is an iron-binding non-heme protein. It acts as a coregulator of several transcription factors, especially the members of NFκB transcription factor family. Based on the redox state of its iron cofactor, it can assume two different conformations and thereby act as a redox sensor inside the nucleus. Previous studies suggested that pirin may be associated with cancer, inflammatory diseases as well as COVID-19 severities. Hence, it is important to explore the pathogenicity of its missense variants. In this study, we used a number of in silico tools to investigate the effects of missense variants of pirin on its structure, stability, metal cofactor binding affinity and interactions with partner proteins. In addition, we used protein dynamics simulation to elucidate the effects of selected variants on its dynamics. Furthermore, we calculated the frequencies of haplotypes containing pirin missense variants across five major super-populations (African, Admixed American, East Asian, European and South Asian).

RESULTS

Among a total of 153 missense variants of pirin, 45 were uniformly predicted to be pathogenic. Of these, seven variants can be considered for further experimental studies. Variants R59P and L116P were predicted to significantly destabilize and damage pirin structure, substantially reduce its affinity to its binding partners and alter pirin residue fluctuation profile via changing the flexibility of several key residues. Additionally, variants R59Q, F78V, G98D, V151D and L220P were found to impact pirin structure and function in multiple ways. As no haplotype was identified to be harboring more than one missense variant, further interrogation of the individual effects of these seven missense variants is highly recommended.

CONCLUSIONS

Pirin is involved in the transcriptional regulation of several genes and can play an important role in inflammatory responses. The variants predicted to be pathogenic in this study may thus contribute to a better understanding of the underlying molecular mechanisms of various inflammatory diseases. Future studies should be focused on clarifying if any of these variants can be used as disease biomarkers.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s42269-022-00917-7.

Thermostabilization mechanisms in thermophilic versus mesophilic three-helix bundle proteins

The engineered three‐helix bundle, UVF, is thermostabilized entropically due to heightened, native‐state dynamics. However, it is unclear whether this thermostabilization strategy is observed in natural proteins from thermophiles. We performed all‐atom, explicit solvent molecular dynamics simulations of two three‐helix bundles from thermophilic H. butylicus (2lvsN and 2lvsC) and compared their dynamics to a mesophilic three‐helix bundle, the Engrailed homeodomain (EnHD). Like UVF, 2lvsC had heightened native dynamics, which it maintained without unfolding at 100°C. Shortening and rigidification of loops in 2lvsN and 2lvsC and increased surface hydrogen bonds in 2lvsN were observed, as is common in thermophilic proteins. A buried disulfide and salt bridge in 2lvsN and 2lvsC, respectively, provided some stabilization, and addition of a homologous disulfide bond in EnHD slowed unfolding. The transferability and commonality of stabilization strategies among members of the three‐helix bundle fold suggest that these strategies may be general and deployable in designing thermostable proteins.

New Variants of the Cytochrome P450 2R1 (CYP2R1) Gene in Individuals with Severe Vitamin D-Activating Enzyme 25(OH)D Deficiency

BACKGROUND

Vitamin D is a fat-soluble cholesterol derivative found in two forms, vitamin D2, and vitamin D3. Cytochrome P450 2R1 (CYP2R1) encoded by the CYP2R1 gene is the major hydroxylase that activates vitamin D by catalyzing the formation of 25-hydroxyvitamin D (25(OH)D).

METHODS

We collected 89 (100%) subjects, 46 of which (51.69%) had a documented severe deficiency of 25(OH)D (<10 ng/mL) and 43 (48.31%) in the control group with documented optimum levels of 25(OH)D (>30 ng/mL). We performed Sanger sequencing of three selected fragments of the CYP2R1 gene (Ch11: 14878000-14878499; Ch11: 14880058-14880883 and Ch11: 14885321-14886113) that affect the binding of substrates to this enzyme and analyzed the possible involvement of genetic variation in these regions with an increased risk of vitamin D deficiency in healthy Polish individuals.

RESULTS

Two substitutions were found within the three fragments. Bioinformatic analysis suggested that one of these (NC_000011.10: g.14878291G>A) may influence the structure and function of CYP2R1.

CONCLUSIONS

Variant NC_000011.10: g.14878291G>A may have a perturbing effect on heme binding in the active site of CYP2R1 and on the function of 25-hydroxylase and probably affects the concentration of 25(OH)D in vivo. We intend to perform functional verification in a larger patient population to confirm and extend these results.

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.

Surface salt bridges contribute to the extreme thermal stability of an FN3-like domain from a thermophilic bacterium

This study uses differential scanning calorimetry, X-ray crystallography, and molecular dynamics simulations to investigate the structural basis for the high thermal stability (melting temperature 97.5°C) of a FN3-like protein domain from thermophilic bacteria Thermoanaerobacter tengcongensis (FN3tt). FN3tt adopts a typical FN3 fold with a three-stranded beta sheet packing against a four-stranded beta sheet. We identified three solvent exposed arginine residues (R23, R25, and R72), which stabilize the protein through salt bridge interactions with glutamic acid residues on adjacent strands. Alanine mutation of the three arginine residues reduced melting temperature by up to 22°C. Crystal structures of the wild type (WT) and a thermally destabilized (∆Tm -19.7°C) triple mutant (R23L/R25T/R72I) were found to be nearly identical, suggesting that the destabilization is due to interactions of the arginine residues. Molecular dynamics simulations showed that the salt bridge interactions in the WT were stable and provided a dynamical explanation for the cooperativity observed between R23 and R25 based on calorimetry measurements. In addition, folding free energy changes computed using free energy perturbation molecular dynamics simulations showed high correlation with melting temperature changes. This work is another example of surface salt bridges contributing to the enhanced thermal stability of thermophilic proteins. The molecular dynamics simulation methods employed in this study may be broadly useful for in silico surface charge engineering of proteins.

Intrinsic basis of thermostability of prolyl oligopeptidase from Pyrococcus furiosus

Salt-bridges play a key role in the thermostability of proteins adapted in stress environments whose intrinsic basis remains to be understood. We find that the higher hydrophilicity of PfP than that of HuP is due to the charged but not the polar residues. The primary role of these residues is to enhance the salt-bridges and their ME. Unlike HuP, PfP has made many changes in its intrinsic property to strengthen the salt-bridge. First, the desolvation energy is reduced by directing the salt-bridge towards the surface. Second, it has made bridge-energy more favorable by recruiting energetically advantageous partners with high helix-propensity among the six possible salt-bridge pairs. Third, ME-residues that perform intricate interactions have increased their energy contribution by making major changes in their binary properties. The use of salt-bridge partners as ME-residues, and ME-residues' overlapping usage, predominant in helices, and energetically favorable substitution are some of the favorable features of PfP compared to HuP. These changes in PfP reduce the unfavorable, increase the favorable ME-energy. Thus, the per salt-bridge stability of PfP is greater than that of HuP. Further, unfavorable target ME-residues can be identified whose mutation can increase the stability of salt-bridge. The study applies to other similar systems.

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Compost Samples from Different Temperature Zones as a Model to Study Co-occurrence of Thermophilic and Psychrophilic Bacterial Population: a Metagenomics Approach

In this study, using a metagenomic approach, we explore the bacterial diversity of compost sites categorized based on their ambient temperatures. The two sites were Reckong Peo in the lower Himalayas and Tambaram in the southern region of the country, namely, CPR and CT. Following assembly of the raw reads from shotgun metagenomics, similarity hits were generated using NCBI BLAST + and SILVA database. A total of 1463 and 1483 species were annotated from CPR and CT. A species-level annotation was performed using a python-based literature search pipeline revealing their growth characteristics. Thermophiles Thermomonospora curvata and Thermus scotoductus were among the prominent species in CT. CPR too was seen abundant with Acidothermus cellulolyticus and Moorella thermoacetica, constituting 10% of the population. Nearly 3% of the identified species in the site CPR were psychrophilic. Although found higher in CPR, psychrophilic species were identified in CT too. Flavobacterium and Psychrobacter spp. were present in both sites without any significant changes in their relative distribution contrary to the thermophilic species abundance (z = - 4.3). Akin to the sequenced samples, database-derived metagenomes also showed similar distribution of thermophiles and psychrophiles. Identifying such peculiar prevalence of extremophiles can be central to understanding extended growth temperatures.

Parallel molecular mechanisms for enzyme temperature adaptation

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

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.

Expression and characterization of a thermotolerant and pH-stable hyaluronate lyase from Thermasporomyces composti DSM22891

Hyaluronate lyases have received extensive attention due to their applications in medical science, drug and biochemical engineering. However, few thermotolerant and pH-stable hyaluronate lyases have been found. In this study, hyaluronate lyase TcHly8B from Thermasporomyces composti DSM22891 was expressed in Escherichia coli BL21(DE3), purified, and characterized. Phylogenetic analysis revealed that TcHly8B belonged to a new subfamily in PL8. The molecular mass of recombinant TcHly8B determined by SDS-PAGE was approximately 86 kDa. The optimal temperature of TcHly8B was 70 °C, which was higher than that of previously reported hyaluronate lyases. TcHly8B was very stable at temperatures from 0 to 60 °C. The optimal pH of TcHly8B was 6.6. It could retain more than 80% of its original enzyme activity after incubation for 12 h in the pH range of 3.0-10.6. TcHly8B degraded hyaluronic acid into unsaturated disaccharides as the end products. The amino acid sequence and structure analysis of TcHly8B demonstrated that the amino acid composition and salt bridges might contribute to the thermostability of TcHly8B. Overall, this study provides an excellent example for the discovery of thermotolerant hyaluronate lyases and can be applied to the industrialized production and basic research of hyaluronate oligosaccharides.

Probing conformational hotspots for the recognition and intervention of protein complexes by lysine reactivity profiling

Probing the conformational and functional hotspot sites within aqueous native protein complexes is still a challenging task. Herein, a mass spectrometry (MS)-based two-step isotope labeling-lysine reactivity profiling (TILLRP) strategy is developed to quantify the reactivities of lysine residues and probe the molecular details of protein–protein interactions as well as evaluate the conformational interventions by small-molecule active compounds. The hotspot lysine sites that are crucial to the SARS-CoV-2 S1–ACE2 combination could be successfully probed, such as S1 Lys⁴¹⁷ and Lys⁴⁴⁴. Significant alteration of the reactivities of lysine residues at the interaction interface of S1-RBD Lys³⁸⁶–Lys⁴⁶² was observed during the formation of complexes, which might be utilized as indicators for investigating the S1-ACE2 dynamic recognition and intervention at the molecular level in high throughput.

A mass spectrometry-based two-step isotope labeling-lysine reactivity profiling strategy is developed to probe the molecular details of protein–protein interactions and evaluate the conformational interventions by small-molecule active compounds.

Enzymatic characterization, molecular dynamics simulation, and application of a novel Bacillus licheniformis laccase

A new laccase gene from newly isolated Bacillus licheniformis TCCC 111219 was actively expressed in Escherichia coli. This recombinant laccase (rLAC) exhibited a high stability towards a wide pH range and high temperatures. 170% of the initial activity was detected at pH 10.0 after 10-d incubation, and 60% of the initial activity was even kept after 2-h incubation at 70 °C. It indicated that only single type of extreme environment, such as strong alkaline environment (300 K, pH 12) or high temperature (370 K, pH 7), did not show obvious impact on the structural stability of rLAC during molecular dynamics simulation process. But the four loop regions of rLAC where the active site is situated were seriously destroyed when strong alkaline and high temperature environment existed simultaneously (370 K, pH 12) because of the damage of hydrogen bonds and salt bridges. Moreover, this thermo- and alkaline-stable enzyme could efficiently decolorize the structurally differing azo, triphenylmethane, and anthraquinone dyes with appropriate mediator at pH 3.0, 7.0, and 9.0 at 60 °C. These rare characteristics suggested its high potential in industrial applications to decolorize textile dyeing effluent.

Alpha-Carbonic Anhydrases from Hydrothermal Vent Sources as Potential Carbon Dioxide Sequestration Agents: In Silico Sequence, Structure and Dynamics Analyses

With the increase in CO2 emissions worldwide and its dire effects, there is a need to reduce CO2 concentrations in the atmosphere. Alpha-carbonic anhydrases (α-CAs) have been identified as suitable sequestration agents. This study reports the sequence and structural analysis of 15 α-CAs from bacteria, originating from hydrothermal vent systems. Structural analysis of the multimers enabled the identification of hotspot and interface residues. Molecular dynamics simulations of the homo-multimers were performed at 300 K, 363 K, 393 K and 423 K to unearth potentially thermostable α-CAs. Average betweenness centrality (BC) calculations confirmed the relevance of some hotspot and interface residues. The key residues responsible for dimer thermostability were identified by comparing fluctuating interfaces with stable ones, and were part of conserved motifs. Crucial long-lived hydrogen bond networks were observed around residues with high BC values. Dynamic cross correlation fortified the relevance of oligomerization of these proteins, thus the importance of simulating them in their multimeric forms. A consensus of the simulation analyses used in this study suggested high thermostability for the α-CA from Nitratiruptor tergarcus. Overall, our novel findings enhance the potential of biotechnology applications through the discovery of alternative thermostable CO2 sequestration agents and their potential protein design.

On the Thermal Stability of O6-Methylguanine-DNA Methyltransferase from Archaeon Pyrococcus kodakaraensis by Molecular Dynamics Simulations

We have employed molecular dynamics simulations to analyze the thermal stability of the O⁶-methylguanine-DNA methyltransferase (MGMT) protein, both hyperthermophilic archaeon Pyrococcus kodakaraensis (Pk-MGMT) and its mesophilic homologue pair, obtained from enterobacterium Escherichia coli (AdaC). This theoretical study was done at three different temperatures: 302, 371, and 450 K. The molecular dynamics has been performed in explicit aqueous solvent during a period of time of 95 ns, including periodic boundary conditions and constant pressure. The same procedure has been used for both proteins, and each simulation has been carried out by triplicate. Hence, we performed 18 simulations. In this way, we have done different analyses to explore the factors that may affect the thermal stability of Pk-MGMT. The structural behavior was analyzed using indicators such as root-mean-square deviation, radius of gyration, solvent-accessible surface area, hydrogen bonds, native contacts, secondary structure, and salt bridge formation. The results showed that when the temperature increases, the global atomic fluctuations increase too, which suggests that both proteins lose thermal stability, but as expected, this fact is highlighted in AdaC. Moreover, the contacts of the native state in AdaC are considerably lower than those found in Pk-MGMT at 450 K. Also, the structural studies showed that conserved and nonconserved salt bridges kept close contacts with the Pk-MGMT protein at high temperatures. These interaction types act as molecular staples and are mainly responsible to provide thermostability to the hyperthermophilic protein.

Lysosomal protein thermal stability does not correlate with cellular half-life: global observations and a case study of tripeptidyl-peptidase 1

Late-infantile neuronal ceroid lipofuscinosis (LINCL) is a neurodegenerative lysosomal storage disorder caused by mutations in the gene encoding the protease tripeptidyl-peptidase 1 (TPP1). Progression of LINCL can be slowed or halted by enzyme replacement therapy, where recombinant human TPP1 is administered to patients. In this study, we utilized protein engineering techniques to increase the stability of recombinant TPP1 with the rationale that this may lengthen its lysosomal half-life, potentially increasing the potency of the therapeutic protein. Utilizing multiple structure-based methods that have been shown to increase the stability of other proteins, we have generated and evaluated over 70 TPP1 variants. The most effective mutation, R465G, increased the melting temperature of TPP1 from 55.6°C to 64.4°C and increased its enzymatic half-life at 60°C from 5.4 min to 21.9 min. However, the intracellular half-life of R465G and all other variants tested in cultured LINCL patient-derived lymphoblasts was similar to that of WT TPP1. These results provide structure/function insights into TPP1 and indicate that improving in vitro thermal stability alone is insufficient to generate TPP1 variants with improved physiological stability. This conclusion is supported by a proteome-wide analysis that indicates that lysosomal proteins have higher melting temperatures but also higher turnover rates than proteins of other organelles. These results have implications for similar efforts where protein engineering approaches, which are frequently evaluated in vitro, may be considered for improving the physiological properties of proteins, particularly those that function in the lysosomal environment.

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.

Disulfide Bond Engineering of an Endoglucanase from Penicillium verruculosum to Improve Its Thermostability

Endoglucanases (EGLs) are important components of multienzyme cocktails used in the production of a wide variety of fine and bulk chemicals from lignocellulosic feedstocks. However, a low thermostability and the loss of catalytic performance of EGLs at industrially required temperatures limit their commercial applications. A structure-based disulfide bond (DSB) engineering was carried out in order to improve the thermostability of EGLII from Penicillium verruculosum. Based on in silico prediction, two improved enzyme variants, S127C-A165C (DSB2) and Y171C-L201C (DSB3), were obtained. Both engineered enzymes displayed a 15⁻21% increase in specific activity against carboxymethylcellulose and β-glucan compared to the wild-type EGLII (EGLII-wt). After incubation at 70 °C for 2 h, they retained 52⁻58% of their activity, while EGLII-wt retained only 38% of its activity. At 80 °C, the enzyme-engineered forms retained 15⁻22% of their activity after 2 h, whereas EGLII-wt was completely inactivated after the same incubation time. Molecular dynamics simulations revealed that the introduced DSB rigidified a global structure of DSB2 and DSB3 variants, thus enhancing their thermostability. In conclusion, this work provides an insight into DSB protein engineering as a potential rational design strategy that might be applicable for improving the stability of other enzymes for industrial applications.

Biochemical and thermodynamic analyses of energy conversion in extremophiles

A variety of extreme environments, characterized by extreme values of various physicochemical parameters (temperature, pressure, salinity, pH, and so on), are found on Earth. Organisms that favorably live in such extreme environments are called extremophiles. All living organisms, including extremophiles, must acquire energy to maintain cellular homeostasis, including extremophiles. For energy conversion in harsh environments, thermodynamically useful reactions and stable biomolecules are essential. In this review, I briefly summarize recent studies of extreme environments and extremophiles living in these environments and describe energy conversion processes in various extremophiles based on my previous research. Furthermore, I discuss the correlation between the biological system of electrotrophy, a third biological energy acquisition system, and the mechanism underlying microbiologically influenced corrosion. These insights into energy conversion in extremophiles may improve our understanding of the "limits of life". Abbreviations: PPi: pyrophosphate; PPase: pyrophosphatase; ITC: isothermal titration microcalorimetry; SVNTase: Shewanella violacea 5'-nucleotidase; SANTase: Shewanella amazonensis 5'-nucleotidase.

Key Role of the Local Hydrophobicity in the East Patch of Plastocyanins on Their Thermal Stability and Redox Properties

Understanding the molecular basis of the thermal stability and functionality of redox proteins has important practical applications. Here, we show a distinct thermal dependence of the spectroscopic and electrochemical properties of two plastocyanins from the thermophilic cyanobacterium Phormidium laminosum and their mesophilic counterpart from Synechocystis sp. PCC 6803, despite the similarity of their molecular structures. To explore the origin of these differences, we have mimicked the local hydrophobicity in the east patch of the thermophilic protein by replacing a valine of the mesophilic plastocyanin by isoleucine. Interestingly, the resulting mutant approaches the thermal stability, redox thermodynamics, and dynamic coupling of the flexible site motions of the thermophilic protein, indicating the existence of a close connection between the hydrophobic packing of the east patch region of plastocyanin and the functional control and stability of the oxidized and reduced forms of the protein.

The Effects of One Amino Acid Substitutions at the C-Terminal Region of Thermostable L2 Lipase by Computational and Experimental Approach

The substitutions of the amino acid at the predetermined critical point at the C-terminal of L2 lipase may increase its thermostability and enzymatic activity, or even otherwise speed up the unfolding of the protein structure. The C-terminal of most proteins is often flexible and disordered. However, some protein functions are directly related to flexibility and play significant role in enzyme reaction. The critical point for mutation of L2 lipase structure was predicted at the position 385 of the L2 sequence, and the best three mutants were determined based on I-Mutant2.0 software. The best three mutants were S385E, S385I and S385V. The effects of the substitution of the amino acids at the critical point were analysed with molecular dynamics simulation by using Yet Another Scientific Artificial Reality Application software. The predicted mutant L2 lipases were found to have lower root mean square deviation value as compared to L2 lipase. It was indicated that all the three mutants had higher compactness in the structure, consequently enhanced the stability. Root mean square fluctuation analysis showed that the flexibility of L2 lipase was reduced by mutations. Purified S385E lipase had an optimum temperature of 80 °C in Tris-HCl pH 8. The highest enzymatic activity of purified S385E lipase was obtained at 80 °C temperature in Tris-HCl pH 8, while for L2 lipase it was at 70 °C in Glycine-NaOH pH 9. The thermal stability of S385V lipase was enhanced as compared to other protein since that the melting point (T _m) value was at 85.96 °C. S385I lipase was more thermostable compared to recombinant L2 lipase and other mutants at temperature 60 °C within 16 h preincubation.

Archaeal DNA on the histone merry-go-round

How did the nucleosome, the fundamental building block of all eukaryotic chromatin, evolve? This central question has been impossible to address because the four core histones that make up the protein core of the nucleosome are so highly conserved in all eukaryotes. With the discovery of small, minimalist histone-like proteins in most known archaea, the likely origin of histones was identified. We recently determined the structure of an archaeal histone-DNA complex, revealing that archaeal DNA topology and protein-DNA interactions are astonishingly similar compared to the eukaryotic nucleosome. This was surprising since most archaeal histones form homodimers which consist only of the minimal histone fold and are devoid of histone tails and extensions. Unlike eukaryotic H2A-H2B and H3-H4 heterodimers that assemble into octameric particles wrapping ~ 150 bp DNA, archaeal histones form polymers around which DNA coils in a quasi-continuous superhelix. At any given point, this superhelix has the same geometry as nucleosomal DNA. This suggests that the architectural role of histones (i.e. the ability to bend DNA into a nucleosomal superhelix) was established before archaea and eukaryotes diverged, while the ability to form discrete particles, together with signaling functions of eukaryotic chromatin (i.e. epigenetic modifications) were secondary additions.

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

Discovering how organisms and their proteins adapt to extreme conditions is a complicated process. Every condition has its own set of adaptations that make it uniquely stable in its environment. The purpose of our review is to discuss what is known in the extremophilic community about protein adaptations. To simplify our mission, we broke the extremophiles into three broad categories: thermophiles, halophiles and psychrophiles. While there are crossover organisms- organisms that exist in two or more extremes, like heat plus acid or cold plus pressure, most of them have a primary adaptation that is within one of these categories which tends to be the most easily identifiable one. While the generally known adaptations are still accepted, like thermophilic proteins have increased ionic interactions and a hardier hydrophobic core, halophilic proteins have a large increase in acidic amino acids and amino acid/peptide insertions and psychrophiles have a much more open structure and reduced ionic interactions, some new information has come to light. Thermophilic stability can be improved by increased subunit-subunit or subunit-cofactor interactions. Halophilic proteins have reversible folding when in the presence of salt. Psychrophilic proteins have an increase in cavities that not only decrease the formation of ice, but also increase flexibility under low temperature conditions. In a proof of concept experiment, we applied what is currently known about adaptations to a well characterized protein, malate dehydrogenase (MDH). While this protein has been profiled in the literature, we are applying our adaptation predictions to its sequence and structure to see if the described adaptations apply. Our analysis demonstrates that thermophilic and halophilic adaptations fit the corresponding MDHs very well. However, because the number of psychrophiles MDH sequences and structures is low, our analysis on psychrophiles is inconclusive and needs more information. By discussing known extremophilic adaptations and applying them to a random, conserved protein, we have found that general adaptations are conserved and can be predicted in proposed extremophilic proteins. The present field of extremophile adaptations is discovering more and more ways organisms and their proteins have adapted. The more that is learned about protein adaptation, the closer we get to custom proteins, designed to fit any extreme and solve some of the world's most pressing environmental problems.

Dynamic structure mediates halophilic adaptation of a DNA polymerase from the deep-sea brines of the Red Sea

The deep-sea brines of the Red Sea are remote and unexplored environments characterized by high temperatures, anoxic water, and elevated concentrations of salt and heavy metals. This environment provides a rare system to study the interplay between halophilic and thermophilic adaptation in biologic macromolecules. The present article reports the first DNA polymerase with halophilic and thermophilic features. Biochemical and structural analysis by Raman and circular dichroism spectroscopy showed that the charge distribution on the protein’s surface mediates the structural balance between stability for thermal adaptation and flexibility for counteracting the salt-induced rigid and nonfunctional hydrophobic packing. Salt bridge interactions via increased negative and positive charges contribute to structural stability. Salt tolerance, conversely, is mediated by a dynamic structure that becomes more fixed and functional with increasing salt concentration. We propose that repulsive forces among excess negative charges, in addition to a high percentage of negatively charged random coils, mediate this structural dynamism. This knowledge enabled us to engineer a halophilic version of Thermococcus kodakarensis DNA polymerase.—Takahashi, M., Takahashi, E., Joudeh, L. I., Marini, M., Das, G., Elshenawy, M. M., Akal, A., Sakashita, K., Alam, I., Tehseen, M., Sobhy, M. A., Stingl, U., Merzaban, J. S., Di Fabrizio, E., Hamdan, S. M. Dynamic structure mediates halophilic adaptation of a DNA polymerase from the deep-sea brines of the Red Sea.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Biochemical and structural insights into a thermostable cellobiohydrolase from Myceliophthora thermophila

Cellobiohydrolases hydrolyze cellulose, a linear polymer with glucose monomers linked exclusively by β-1,4 glycosidic linkages. The widespread hydrogen bonding network tethers individual cellulose polymers forming crystalline cellulose, which prevent the access of hydrolytic enzymes and water molecules. The most abundant enzyme secreted by Myceliophthora thermophila M77 in response to the presence of biomass is the cellobiohydrolase MtCel7A, which is composed by a GH7-catalytic domain (CD), a linker, and a CBM1-type carbohydrate-binding module. GH7 cellobiohydrolases have been studied before, and structural models have been proposed. However, currently available GH7 crystal structures only define separate catalytic domains and/or cellulose-binding modules and do not include the full-length structures that are involved in shaping the catalytic mode of operation. In this study, we determined the 3D structure of catalytic domain using X-ray crystallography and retrieved the full-length enzyme envelope via small-angle X-ray scattering (SAXS) technique. The SAXS data reveal a tadpole-like molecular shape with a rigid linker connecting the CD and CBM. Our biochemical studies show that MtCel7A has higher catalytic efficiency and thermostability as well as lower processivity when compared to the well-studied TrCel7A from Trichoderma reesei. Based on a comparison of the crystallographic structures of CDs and their molecular dynamic simulations, we demonstrate that MtCel7A has considerably higher flexibility than TrCel7A. In particular, loops that cover the active site are more flexible and undergo higher conformational fluctuations, which might account for decreased processivity and enhanced enzymatic efficiency. Our statistical coupling analysis suggests co-evolution of amino acid clusters comprising the catalytic site of MtCel7A, which correlate with the steps in the catalytic cycle of the enzyme.

DATABASE

The atomic coordinates and structural factors of MtCel7A have been deposited in the Protein Data Bank with accession number 5W11.

Effects of Trimethylamine-N-oxide on the Conformation of Peptides and its Implications for Proteins

To provide insights into the stabilizing mechanisms of trimethylamine-N-oxide (TMAO) on protein structures, we perform all-atom molecular dynamics simulations of peptides and the Trp-cage miniprotein. The effects of TMAO on the backbone and charged residues of peptides are found to stabilize compact conformations, whereas effects of TMAO on nonpolar residues lead to peptide swelling. This suggests competing mechanisms of TMAO on proteins, which accounts for hydrophobic swelling, backbone collapse, and stabilization of charge-charge interactions. These mechanisms are observed in Trp cage.

Determinants of thermostability in the cytochrome P450 fold

Cytochromes P450 are found throughout the biosphere in a wide range of environments, serving a multitude of physiological functions. The ubiquity of the P450 fold suggests that it has been co-opted by evolution many times, and likely presents a useful compromise between structural stability and conformational flexibility. The diversity of substrates metabolized and reactions catalyzed by P450s makes them attractive starting materials for use as biocatalysts of commercially useful reactions. However, process conditions impose different requirements on enzymes to those in which they have evolved naturally. Most natural environments are relatively mild, and therefore most P450s have not been selected in Nature for the ability to withstand temperatures above ~40°C, yet industrial processes frequently require extended incubations at much higher temperatures. Thus, there has been considerable interest and effort invested in finding or engineering thermostable P450 systems. Numerous P450s have now been identified in thermophilic organisms and analysis of their structures provides information as to mechanisms by which the P450 fold can be stabilized. In addition, protein engineering, particularly by directed or artificial evolution, has revealed mutations that serve to stabilize particular mesophilic enzymes of interest. Here we review the current understanding of thermostability as it applies to the P450 fold, gleaned from the analysis of P450s characterized from thermophilic organisms and the parallel engineering of mesophilic forms for greater thermostability. We then present a perspective on how this information might be used to design stable P450 enzymes for industrial application. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Complex network analysis of thermostable mutants of Bacillus subtilis Lipase A

Three-dimensional structures of proteins that regulate their functions can be modelled using complex network based approaches for understanding the structure-function relationship. The six mutants of the protein Lipase A from Bacillus subtilis, harbouring 2 to 12 mutations, retain their function at higher temperatures with negligible variation in their overall three-dimensional crystallographic structures. This enhanced thermostability of the mutants questions the structure-function paradigm. In this paper, a coarse-grained complex network approach is used to elucidate the structural basis of enhanced thermostability in the mutant proteins, by uncovering small but significant local changes distributed throughout the structure, rendering stability to the mutants at higher temperatures. Community structure analysis of the six mutant protein networks uncovers the specific reorganisations among the nodes/residues that occur, in absence of overall structural variations, which induce enhanced rigidity underlying the increased thermostability. This study offers a novel and significant application of complex network analysis that proposes to be useful in the understanding and designing of thermostable proteins.

Binding Properties of Split tRNA to the C-terminal Domain of Methionyl-tRNA Synthetase of Nanoarchaeum equitans

The C-terminal domain of methionyl-tRNA synthetase (MetRS-C) from Nanoarchaeum equitans is homologous to a tRNA-binding protein consisting of 111 amino acids (Trbp111) from Aquifex aeolicus. The crystal structure of MetRS-C showed that it existed as a homodimer, and that each monomer possessed an oligonucleotide/oligosaccharide-binding fold (OB-fold). Analysis using a quartz crystal microbalance indicated that MetRS-C freshly isolated from N. equitans was bound to tRNA. However, binding of the split 3'-half tRNA species was stronger than that of the 5'-half species. The T-loop and the 3'-end regions of the split 3'-half tRNA were found to be responsible for the binding. The minimum structure for binding to MetRS-C might be a minihelix-like stem-loop with single-stranded 3'-terminus. After successive duplications of such a small hairpin structure with the assistance of a Trbp-like structure, the interaction of the T-loop region of the 3'-half with a Trbp-like structure could have been evolutionarily replaced by RNA-RNA interactions, along with many combinational tertiary interactions, to form the modern tRNA structure.

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.