Protein thermostability above 100 degreesC: a key role for ionic interactions

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.

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

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

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).

Engineering enhanced thermostability into the Geobacillus pallidus nitrile hydratase

Nitrile hydratases (NHases) are important biocatalysts for the enzymatic conversion of nitriles to industrially-important amides such as acrylamide and nicotinamide. Although thermostability in this enzyme class is generally low, there is not sufficient understanding of its basis for rational enzyme design. The gene expressing the Co-type NHase from the moderate thermophile, Geobacillus pallidus RAPc8 (NRRL B-59396), was subjected to random mutagenesis. Four mutants were selected that were 3 to 15-fold more thermostable than the wild-type NHase, resulting in a 3.4–7.6 ​kJ/mol increase in the activation energy of thermal inactivation at 63 ​°C. High resolution X-ray crystal structures (1.15–1.80 ​Å) were obtained of the wild-type and four mutant enzymes. Mutant 9E, with a resolution of 1.15 ​Å, is the highest resolution crystal structure obtained for a nitrile hydratase to date. Structural comparisons between the wild-type and mutant enzymes illustrated the importance of salt bridges and hydrogen bonds in enhancing NHase thermostability. These additional interactions variously improved thermostability by increased intra- and inter-subunit interactions, preventing cooperative unfolding of α-helices and stabilising loop regions. Some hydrogen bonds were mediated via a water molecule, specifically highlighting the significance of structured water molecules in protein thermostability. Although knowledge of the mutant structures makes it possible to rationalize their behaviour, it would have been challenging to predict in advance that these mutants would be stabilising.

•

Random mutagenesis yields a 15-fold increase in nitrile hydratase thermostability.

•

Salt bridges and hydrogen bonds improves nitrile hydratase thermostability.

•

Water-mediated hydrogen bonds improves protein thermostability.

Biochemical and structural studies of two tetrameric nucleoside 2'-deoxyribosyltransferases from psychrophilic and mesophilic bacteria: Insights into cold-adaptation

Nucleoside 2'-deoxyribosyltransferases (NDTs) catalyze the cleavage of glycosidic bonds of 2'-deoxynucleosides and the following transfer of the 2'-deoxyribose moiety to acceptor nucleobases. Here, we report the crystal structures and biochemical properties of the first tetrameric NDTs: the type I NDT from the mesophilic bacterium Enterococcus faecalis V583 (EfPDT) and the type II NDT from the bacterium Desulfotalea psychrophila (DpNDT), the first psychrophilic NDT. This novel structural and biochemical data permitted an exhaustive comparative analysis aimed to shed light into the basis of the high global stability of the psychrophilic DpNDT, which has a higher melting temperature than EfPDT (58.5 °C versus 54.4 °C) or other mesophilic NDTs. DpNDT possesses a combination of unusual structural motifs not present neither in EfPDT nor any other NDT that most probably contribute to its global stability, in particular, a large aliphatic isoleucine-leucine-valine (ILV) bundle accompanied by a vicinal disulfide bridge and also an intersubunit disulfide bridge, the first described for an NDT. The functional and structural features of DpNDT do not fit the standard features of psychrophilic enzymes, which lead us to consider the implication of (sub)cellular levels together with the protein level in the adaptation of enzymatic activity to low temperatures.

How DNA affects the hyperthermophilic protein Ape10b2 for oligomerization: an investigation using multiple short molecular dynamics simulations

Alba2 is a hyperthermophilic DNA-binding protein, and DNA plays a crucial role in the Alba2 oligomerization process. It is a pity that there is limited research in terms of how DNA affects the conformational change of Alba2 in oligomerization. Herein, we complement the crystal structure of the Ape10b2 (belongs to Alba2)-dsDNA complex (PDB ID: 3U6Y) and employ multiple short molecular dynamics (MSMD) simulations to illuminate the influence of DNA on Ape10b2 at four temperatures (300, 343, 363, and 373 K). Our results indicate that DNA could cause the conformational changes of two important regions (loop1 and loop5), which may be beneficial for protein oligomerization. The results of hydrogen bond analysis show that the increasing number of hydrogen bonds between two monomers of Ape10b2 may also be a favorable factor for oligomerization. In addition, Ape10b2 can stabilize DNA by electrostatic interactions with an increase in temperature, and five residues (Arg40, Arg42, Asn43, Asn45, and Arg46) play a stabilizing role during protein binding to DNA. Our findings could help in understanding the favorable factors leading to protein oligomerization, which contributes to enzyme engineering research from an industrial perspective.

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.

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.

Temperature effect on the SARS-CoV-2: A molecular dynamics study of the spike homotrimeric glycoprotein

Rapid spread of SARS-CoV-2 virus have boosted the need of knowledge about inactivation mechanisms to minimize the impact of COVID-19 pandemic. Recent studies have shown that SARS-CoV-2 virus can be disabled by heating, the exposure time for total inactivation depending on the reached temperature (e.g. more than 45 min at 329 K or less than 5 min at 373 K. In spite of recent crystallographic structures, little is known about the molecular changes induced by the temperature. Here, we unravel the molecular basis of the effect of the temperature over the SARS-CoV-2 spike glycoprotein, which is a homotrimer with three identical monomers, by executing atomistic molecular dynamics (MD) simulations at 298, 310, 324, 338, 358 and 373 K. Furthermore, both the closed down and open up conformational states, which affect the accessibility of receptor binding domain, have been considered. Our results suggest that the spike homotrimer undergoes drastic changes in the topology of the hydrogen bonding interactions and important changes on the secondary structure of the receptor binding domain (RBD), while electrostatic interactions (i.e. salt bridges) are mainly preserved. The proposed inactivation mechanism has important implications for engineering new approaches to fight the SARS-CoV-2 coronavirus, as for example, cleaving or reorganizing the hydrogen bonds through chaotropic agents or nanoparticles with local surface resonant plasmon effect.

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.

Sperm chromatin-condensing protamine enhances SMYD5 thermal stability

Studying thermal stability of proteins not only provides insight into protein structure but also is instrumental in identifying previously unknown interaction partners. We develop a machine learning strategy that combines orthogonal partial least squares regression and stability screening of Silver Bullets Bio library to identify biologically active molecules that enhance protein stability. This strategy proves effective in extracting the stability-enhancing molecules for SMYD5, a histone lysine methyltransferase that regulates chromosome integrity. Protamine, a histone substitute in chromatin condensation during spermatogenesis, is identified as the most influential molecule to enhance SMYD5 thermal stability. We find that the C-terminal poly-glutamic acid tract (poly-E) and a 30-residue insertion in MYND domain (M-insertion), which are unique to SMYD5, regulate the structural stability. However, protamine plays a dominant role in SMYD5 stability, and in the presence of protamine, the poly-E tract or M-insertion loses its ability to affect the stability. The stability-enhancing effect of protamine is SMYD5 specific, and for SMYD2, a closely related homolog, protamine exhibits opposite, destabilizing effects. We find that both SMYD5 and SMYD2 interact with protamine, where SMYD5 interaction is independent of the poly-E tract and M-insertion. Protamine not only helps provide insight into the structure-stability relationships of SMYD5, but also suggests a potential functional link of SMYD5 to spermatogenesis. SMYD5 is a ubiquitously expressed gene with the highest expression in testis, especially in the seminiferous ducts that contain germ cells. Thus, our study opens up avenues that could help delineate major mechanisms underlying chromatin dynamics during spermatogenesis.

X-ray structure of C-phycocyanin from Galdieria phlegrea: Determinants of thermostability and comparison with a C-phycocyanin in the entire phycobilisome

Galdieria phlegrea is a polyextremophilic red alga belonging to Cyanidiophyceae. Galdieria phlegrea C-phycocyanin (GpPC), an abundant light-harvesting pigment with an important role in energy capture and transfer to photosystems, is the C-phycocyanin (C-PC) with the highest thermal stability described so far. GpPC also presents interesting antioxidant and anticancer activities. The X-ray structure of the protein was here solved. GpPC is a [(αβ)₃]₂ hexamer, with the phycocyanobilin chromophore attached to Cys84α, Cys82β and Cys153β. Details of geometry and interaction with solvent of the chromophores are reported. Comparison with the structure of a C-PC in the entire Porphyridium purpureum phycobilisome system reveals that linker polypeptides have a significant effect on the local structure of the chromophores environment. Comparative analyses with the structures of other purified C-PCs, which were carried out including re-refined models of G. sulphuraria C-PC, reveal that GpPC presents a significantly higher number of inter-trimer salt bridges. Notably, the higher number of salt bridges at the (αβ)₃/(αβ)₃ interface is not due to an increased number of charged residues in this region, but to subtle conformational variations of their side chains, which are the result of mutations of close polar and non-polar residues.

Revealing the Structural Contributions to Thermal Adaptation of the TATA-Box Binding Protein: Molecular Dynamics and QSPR Analyses

The TATA-box binding protein (TBP) is an important element of the transcription machinery in archaea and eukaryotic organisms. TBP is expressed in organisms adapted to different temperatures, indicating a robust structure, and experimental studies have shown that the mid-unfolding temperature (Tm) of TBP is directly correlated with the optimal growth temperature (OGT) of the organism. To understand which are the relevant structural requirements for its stability, we present the first structural and dynamic computational study of TBPs, combining molecular dynamics (MD) simulations and a quantitative structure-property relationship (QSPR) over a set of TBPs of organisms adapted to different temperatures. We found that the main structural properties of TBP used to adapt to high temperatures are an increase in the ease of desolvation of charged residues at the surface, an increase in the local resiliency, the presence of Leu clusters in the protein core, and an increase in the loss of hydrophobic packing in the N-terminal subdomain. In view of our results, we consider that TBP is a good model to study thermal adaptation, and our analysis opens the possibility of performing protein engineering on TBPs to study transcription at high or low temperatures.

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

Several organisms, specifically microorganisms survive in a wide range of harsh environments including extreme temperature, pH, and salt concentration. We analyzed systematically a large number of protein sequences with their structures to understand their stability and to discriminate extremophilic proteins from their non-extremophilic orthologs. Our results highlighted that the strategy for the packing of the protein core was influenced by the environmental stresses through substitutive structural events through better ionic interaction. Statistical analysis showed that a significant difference in number and composition of amino acid exist among them. The negative correlation of pairwise sequence alignments and structural alignments indicated that most of the extremophile and non-extremophile proteins didn’t contain any association for maintaining their functional stability. A significant numbers of salt bridges were noticed on the surface of the extremostable proteins. The Ramachandran plot data represented more occurrences of amino acids being present in helix and sheet regions of extremostable proteins. We also found that a significant number of small nonpolar amino acids and moderate number of charged amino acids like Arginine and Aspartic acid represented more nonplanar Omega angles in their peptide bond. Thus, extreme conditions may predispose amino acid composition including geometric variability for molecular adaptation of extremostable proteins against atmospheric variations and associated changes under natural selection pressure. The variation of amino acid composition and structural diversifications in proteins play a major role in evolutionary adaptation to mitigate climate change.

Structural insight into the fungal β-glucosidases and their interactions with organics

Fungal β-glucosidases (BGLs) have unceasingly utilized for industrial applications and recently, they possess a crucial role in bioethanol production. To engineer the BGLs, understanding their structures, intermolecular interactions and molecular docking is requisite, which is carried out in this work based on the glycosyl hydrolase (GH) family. Among 12 BGLs, protein sequence, structure, and conserved sites of GH1 BGLs are evidently diverged to GH3 BGLs. Even biophysical and chemical features of GH1 BGLs are utterly varied from GH3 BGLs, wherein pI, instability index, aliphatic index, surface & buried area, thermostability and thermodynamics are included. On the contrary, aromatic, charged, polar, and hydrophobic residues are significantly higher in GH1 BGLs as compared to that of GH3 BGLs. Moreover, molecular docking of BGLs with 12 substrates and 5 inhibitors revealed that the GH3 BGLs efficiently bound with laminaribose, gentibiose, aryl- and cello-substrates than GH1 BGLs; however, GH3 BGLs are noticeably inhibited by glucose, glucono-δ-lactone, methanetriamine. So, structural insight of BGLs provides an explicit knowledge regarding the catalytic residues, biophysical chemistry and notable binding ligands, which are most important factors for enzyme engineering.

Evaluating the strengths of salt bridges in the CutA1 protein using molecular dynamic simulations: a comparison of different force fields

Ion–ion interactions (salt bridges) between favorable pairs of charged residues are important for the conformational stability of proteins. Molecular dynamic (MD) simulations are useful for elucidating the interactions among charged residues fluctuating in solution. However, the quality of MD results depends strongly on the force fields used. In this study, we compared the strengths of salt bridges among force fields by performing MD simulations using the CutA1 protein (trimer) from the hyperthermophile Pyrococcus horikoshii (PhCutA1), which has an unusually large proportion of charged residues. The force fields Chemistry at HARvard Macromolecular Mechanics (Charmm)27, Assisted Model Building and Energy Refinement (Amber)99sb, Amber14sb, GROningen Molecular Simulation (Gromos)43a1, and Gromos53a6 were used in combination with two different water models, tip3p (for Charmm27, Amber99sb, and Amber14sb) and simple point charge/extended (for Amber99sb, Gromos43a1, and Gromos53a6), yielding a total of six combinations. The RMSDs of all Cα atoms of PhCutA1 were similar among force fields, except for Charmm27, during 400‐ns MD simulations at 300 K; however, the radius of gyration (R_g) was greater for Amber99sb and shorter for Gromos43a1. The average strengths of salt bridges for each positively charged residue did not differ greatly among force fields, but the strengths at specific sites within the structure depended sensitively on the force field used. In the case of the Gromos group, positively charged residues could engage in favorable interactions with many more charged residues than in the other force fields, especially in loop regions; consequently, the apparent strength at each site was lower.

Salt Bridge in Ligand-Protein Complexes-Systematic Theoretical and Statistical Investigations

Although the salt bridge is the strongest among all known noncovalent molecular interactions, no comprehensive studies have been conducted to date to examine its role and significance in drug design. Thus, a systematic study of the salt bridge in biological systems is reported herein, with a broad analysis of publicly available data from Protein Data Bank, DrugBank, ChEMBL, and GPCRdb. The results revealed the distance and angular preferences as well as privileged molecular motifs of salt bridges in ligand-receptor complexes, which could be used to design the strongest interactions. Moreover, using quantum chemical calculations at the MP2 level, the energetic, directionality, and spatial variabilities of salt bridges were investigated using simple model systems mimicking salt bridges in a biological environment. Additionally, natural orbitals for chemical valence (NOCV) combined with the extended-transition-state (ETS) bond-energy decomposition method (ETS-NOCV) were analyzed and indicated a strong covalent contribution to the salt bridge interaction. The present results could be useful for implementation in rational drug design protocols.

From Levinthal's Paradox to the Effects of Cell Environmental Perturbation on Protein Folding

BACKGROUND

The rapidly increasing number of known protein sequences calls for more efficient methods to predict the Three-Dimensional (3D) structures of proteins, thus providing basic knowledge for rational drug design. Understanding the folding mechanism of proteins is valuable for predicting their 3D structures and for designing proteins with new functions and medicinal applications. Levinthal's paradox is that although the astronomical number of conformations possible even for proteins as small as 100 residues cannot be fully sampled, proteins in nature normally fold into the native state within timescales ranging from microseconds to hours. These conflicting results reveal that there are factors in organisms that can assist in protein folding.

METHODS

In this paper, we selected a crowded cell-like environment and temperature, and the top three Posttranslational Modifications (PTMs) as examples to show that Levinthal's paradox does not reflect the folding mechanism of proteins. We then revealed the effects of these factors on protein folding.

RESULTS

The results summarized in this review indicate that a crowded cell-like environment, temperature, and the top three PTMs reshape the Free Energy Landscapes (FELs) of proteins, thereby regulating the folding process. The balance between entropy and enthalpy is the key to understanding the effect of the crowded cell-like environment and PTMs on protein folding. In addition, the stability/flexibility of proteins is regulated by temperature.

CONCLUSION

This paper concludes that the cellular environment could directly intervene in protein folding. The long-term interactions of the cellular environment and sequence evolution may enable proteins to fold efficiently. Therefore, to correctly understand the folding mechanism of proteins, the effect of the cellular environment on protein folding should be considered.

RankProt: A multi criteria-ranking platform to attain protein thermostabilizing mutations and its in vitro applications - Attribute based prediction method on the principles of Analytical Hierarchical Process

Attaining recombinant thermostable proteins is still a challenge for protein engineering. The complexity is the length of time and enormous efforts required to achieve the desired results. Present work proposes a novel and economic strategy of attaining protein thermostability by predicting site-specific mutations at the shortest possible time. The success of the approach can be attributed to Analytical Hierarchical Process and the outcome was a rationalized thermostable mutation(s) prediction tool- RankProt. Briefly the method involved ranking of 17 biophysical protein features as class predictors, derived from 127 pairs of thermostable and mesostable proteins. Among the 17 predictors, ionic interactions and main-chain to main-chain hydrogen bonds were the highest ranked features with eigen value of 0.091. The success of the tool was judged by multi-fold in silico validation tests and it achieved the prediction accuracy of 91% with AUC 0.927. Further, in vitro validation was carried out by predicting thermostabilizing mutations for mesostable Bacillus subtilis lipase and performing the predicted mutations by multi-site directed mutagenesis. The rationalized method was successful to render the lipase thermostable with optimum temperature stability and T_m increase by 20°C and 7°C respectively. Conclusively it can be said that it was the minimum number of mutations in comparison to the number of mutations incorporated to render Bacillus subtilis lipase thermostable, by directed evolution techniques. The present work shows that protein stabilizing mutations can be rationally designed by balancing the biophysical pleiotropy of proteins, in accordance to the selection pressure.

Density discriminates between thermophilic and mesophilic proteins

Despite an intense interest and a remarkable number of studies on the subject, the relationships between thermostability and (primary, secondary and tertiary) structure of proteins are still not fully understood. Here, comparing the protein density - defined by the ratio between the residue number and protein excluded volume - for a set of thermophilic/mesophilic pairs, we provide evidence that this property is connected to the optimal growth temperature. In particular, our results indicate that thermophilic proteins have - in general - a lower density with respect to the mesophilic counterparts, being such a correlation more pronounced for optimal growth temperature differences greater than 40°C. The effect of the protein thermostability changes on the molecular shape is also presented.

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.

Four Inserts within the Catalytic Domain Confer Extra Stability and Activity to Hyperthermostable Pyrolysin from Pyrococcus furiosus

Pyrolysin from the hyperthermophilic archaeon Pyrococcus furiosus is the prototype of the pyrolysin family of the subtilisin-like serine protease superfamily (subtilases). It contains four inserts (IS147, IS29, IS27, and IS8) of unknown function in the catalytic domain. We performed domain deletions and showed that three inserts are either essential (IS147 and IS27) or important (IS8) for efficient maturation of pyrolysin at high temperatures, whereas IS29 is dispensable. The large insert IS147 contains Ca3 and Ca4, two calcium-binding Dx[DN]xDG motifs that are conserved in many pyrolysin-like proteases. Mutagenesis revealed that the Ca3 site contributes to enzyme thermostability and the Ca4 site is necessary for pyrolysin to fold into a maturation-competent conformation. Mature insert-deletion variants were characterized and showed that IS29 and IS8 contribute to enzyme activity and stability, respectively. In the presence of NaCl, pyrolysin undergoes autocleavage at two sites: one within IS29 and the other in IS27 Disrupting the ion pairs in IS27 and IS8 induces autocleavage in the absence of salts. Interestingly, autocleavage products combine noncovalently to form an active, nicked enzyme that is resistant to SDS and urea denaturation. Additionally, a single mutation in IS29 increases resistance to salt-induced autocleavage and further increases enzyme thermostability. Our results suggest that these extra structural elements play a crucial role in adapting pyrolysin to hyperthermal environments.IMPORTANCE Pyrolysin-like proteases belong to the subtilase superfamily and are characterized by large inserts and long C-terminal extensions; however, the role of the inserts in enzyme function is unclear. Our results demonstrate that four inserts in the catalytic domain of hyperthermostable pyrolysin contribute to the folding, maturation, stability, and activity of the enzyme at high temperatures. The modification of extra structural elements in pyrolysin-like proteases is a promising strategy for modulating global structure stability and enzymatic activity of this class of protease.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

Thermodynamics of protein denaturation at temperatures over 100 °C: CutA1 mutant proteins substituted with hydrophobic and charged residues

Although the thermodynamics of protein denaturation at temperatures over 100 °C is essential for the rational design of highly stable proteins, it is not understood well because of the associated technical difficulties. We designed certain hydrophobic mutant proteins of CutA1 from Escherichia coli, which have denaturation temperatures (Td) ranging from 101 to 113 °C and show a reversible heat denaturation. Using a hydrophobic mutant as a template, we successfully designed a hyperthermostable mutant protein (Td = 137 °C) by substituting six residues with charged ones. Thermodynamic analyses of these mutant proteins indicated that the hydrophobic mutants were stabilized by the accumulation of denaturation enthalpy (ΔH) with no entropic gain from hydrophobic solvation around 100 °C, and that the stabilization due to salt bridges resulted from both the increase in ΔH from ion-ion interactions and the entropic effect of the electrostatic solvation over 113 °C. This is the first experimental evidence that has successfully overcome the typical technical difficulties.

Intramolecular hydrogen bonding in the polyextremophilic short-chain dehydrogenase from the archaeon Thermococcus sibiricus and its close structural homologs

The short-chain alcohol dehydrogenase from the archaeon Thermococcus sibiricus (TsAdh319) exhibits adaptation to different kinds of stress: high temperature, high salinity, and the presence of organic solvents and denaturants. Previously a comparison of TsAdh319 with close structural homologs revealed an abnormally large number of charged residues on the surface of TsAdh319 tetramer. We further focused on the analysis of hydrogen bonding of TsAdh319 and its structural homologs from thermophilic and mesophilic organisms as a structural factor of adaptation to extreme environment. The calculation and analysis of the dynamics of hydrogen bonds of different kind were performed. In particular, the intramolecular hydrogen bonds of different kind according to their location and the type of a.a. residues involved in the bond were analyzed. TsAdh319 showed the greatest contribution of charged residues to the formation of surface hydrogen bonds, inner hydrogen bonding, and the bonds between different subunits compared to its structural homologs. Molecular dynamics simulations revealed that, of three enzyme molecules analyzed, TsAdh319 shows the least change in the number of hydrogen bonds of different kinds upon a temperature shift from 27 to 85 °C. The greatest changes were observed for a homologous enzyme from a mesophilic host. Only guanidine hydrochloride being a charged agent was able to deactivate TsAdh319. We suggest that the percentage of charged residues plays a key role in the resistance of TsAdh319 to environmental stress. The analysis shows that salt bridges in TsAdh319 serve as a universal instrument of stabilization under different extreme conditions.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

Structural Basis Underlying the Binding Preference of Human Galectins-1, -3 and -7 for Galβ1-3/4GlcNAc

Galectins represent β-galactoside-binding proteins and are known to bind Galβ1-3/4GlcNAc disaccharides (abbreviated as LN1 and LN2, respectively). Despite high sequence and structural homology shared by the carbohydrate recognition domain (CRD) of all galectin members, how each galectin displays different sugar-binding specificity still remains ambiguous. Herein we provided the first structural evidence of human galectins-1, 3-CRD and 7 in complex with LN1. Galectins-1 and 3 were shown to have higher affinity for LN2 than for LN1, while galectin-7 displayed the reversed specificity. In comparison with the previous LN2-complexed structures, the results indicated that the average glycosidic torsion angle of galectin-bound LN1 (ψ^LN1 ≈ 135°) was significantly differed from that of galectin-bound LN2 (ψ^LN2 ≈ -108°), i.e. the GlcNAc moiety adopted a different orientation to maintain essential interactions. Furthermore, we also identified an Arg-Asp/Glu-Glu-Arg salt-bridge network and the corresponding loop (to position the second Asp/Glu residue) critical for the LN1/2-binding preference.

A novel mechanism of protein thermostability: a unique N-terminal domain confers heat resistance to Fe/Mn-SODs

Superoxide dismutases (SODs), especially thermostable SODs, are widely applied in medical treatments, cosmetics, food, agriculture, and other industries given their excellent antioxidant properties. A novel thermostable cambialistic SOD from Geobacillus thermodenitrificans NG80-2 exhibits maximum activity at 70°C and high thermostability over a broad range of temperatures (20–80°C). Unlike other reported SODs, this enzyme contains an extra repeat-containing N-terminal domain (NTD) of 244 residues adjacent to the conserved functional SODA domain. Deletion of the NTD dramatically decreased its optimum active temperature (OAT) to 30°C and also impaired its thermostability. Conversely, appending the NTD to a mesophilic counterpart from Bacillus subtilis led to a moderately thermophilic enzyme (OAT changed from 30 to 55°C) with improved heat resistance. Temperature-dependant circular dichroism analysis revealed the enhanced conformational stability of SODs fused with this NTD. Furthermore, the NTD also contributes to the stress resistance of host proteins without altering their metal ion specificity or oligomerisation form except for a slight effect on their pH profile. We therefore demonstrate that the NTD confers outstanding thermostability to the host protein. To our knowledge, this is the first discovery of a peptide capable of remarkably improving protein thermostability and provides a novel strategy for bioengineering thermostable SODs.

Protein thermal stability enhancement by designing salt bridges: a combined computational and experimental study

Protein thermal stability is an important factor considered in medical and industrial applications. Many structural characteristics related to protein thermal stability have been elucidated, and increasing salt bridges is considered as one of the most efficient strategies to increase protein thermal stability. However, the accurate simulation of salt bridges remains difficult. In this study, a novel method for salt-bridge design was proposed based on the statistical analysis of 10,556 surface salt bridges on 6,493 X-ray protein structures. These salt bridges were first categorized based on pairing residues, secondary structure locations, and Cα–Cα distances. Pairing preferences generalized from statistical analysis were used to construct a salt-bridge pair index and utilized in a weighted electrostatic attraction model to find the effective pairings for designing salt bridges. The model was also coupled with B-factor, weighted contact number, relative solvent accessibility, and conservation prescreening to determine the residues appropriate for the thermal adaptive design of salt bridges. According to our method, eight putative salt-bridges were designed on a mesophilic β-glucosidase and 24 variants were constructed to verify the predictions. Six putative salt-bridges leaded to the increase of the enzyme thermal stability. A significant increase in melting temperature of 8.8, 4.8, 3.7, 1.3, 1.2, and 0.7°C of the putative salt-bridges N437K–D49, E96R–D28, E96K–D28, S440K–E70, T231K–D388, and Q277E–D282 was detected, respectively. Reversing the polarity of T231K–D388 to T231D–D388K resulted in a further increase in melting temperatures by 3.6°C, which may be caused by the transformation of an intra-subunit electrostatic interaction into an inter-subunit one depending on the local environment. The combination of the thermostable variants (N437K, E96R, T231D and D388K) generated a melting temperature increase of 15.7°C. Thus, this study demonstrated a novel method for the thermal adaptive design of salt bridges through inference of suitable positions and substitutions.

Prediction of the determinants of thermal stability by linear discriminant analysis: the case of the glutamate dehydrogenase protein family

Little is known about the determinants of thermal stability in individual protein families. Most of the knowledge on thermostability comes, in fact, from comparative analyses between large, and heterogeneous, sets of thermo- and mesophilic proteins. Here, we present a multivariate statistical approach aimed to detect signature sequences for thermostability in a single protein family. It was applied to the glutamate dehydrogenase (GDH) family, which is a good model for investigating this peculiar process. The structure of GDH consists of six subunits, each of them organized into two domains. Formation of ion-pair networks on the surface of the protein subunits, or increase in the inter-subunit hydrophobic interactions, have been suggested as important factors for explaining stability at high temperatures. However, identification of the amino acid changes that are involved in this process still remains elusive. Our approach consisted of a linear discriminant analysis on a set of GDH sequences from Archaea and Bacteria (33 thermo- and 36 mesophilic GDHs). It led to detection of 3 amino acid clusters as the putative determinants of thermal stability. They were localized at the subunit interface or in close proximity to the binding site of the NAD(P)(+) coenzyme. Analysis within the clusters led to prediction of 8 critical amino acid sites. This approach could have a wide utility, in the ligth of the notion that each protein family seems to adopt its own strategy for achieving thermostability.

Effects of metal ions on stability and activity of hyperthermophilic pyrolysin and further stabilization of this enzyme by modification of a Ca2+-binding site

Pyrolysin is an extracellular subtilase produced by the marine hyperthermophilic archaeon Pyrococcus furiosus. This enzyme functions at high temperatures in seawater, but little is known about the effects of metal ions on the properties of pyrolysin. Here, we report that the supplementation of Na(+), Ca(2+), or Mg(2+) salts at concentrations similar to those in seawater destabilizes recombinant pyrolysin but leads to an increase in enzyme activity. The destabilizing effect of metal ions on pyrolysin appears to be related to the disturbance of surface electrostatic interactions of the enzyme. In addition, mutational analysis of two predicted high-affinity Ca(2+)-binding sites (Ca1 and Ca2) revealed that the binding of Ca(2+) is important for the stabilization of this enzyme. Interestingly, Asn substitutions at residues Asp818 and Asp820 of the Ca2 site, which is located in the C-terminal extension of pyrolysin, resulted in improvements in both enzyme thermostability and activity without affecting Ca(2+)-binding affinity. These effects were most likely due to the elimination of unfavorable electrostatic repulsion at the Ca2 site. Together, these results suggest that metal ions play important roles in modulating the stability and activity of pyrolysin.

Role of loops connecting secondary structure elements in the stabilization of proteins isolated from thermophilic organisms

It has been recently discovered that the connection of secondary structure elements (ββ-unit, βα- and αβ-units) in proteins follows quite stringent principles regarding the chirality and the orientation of the structural units (Koga et al., Nature 2012;491:222-227). By exploiting these rules, a number of protein scaffolds endowed with a remarkable thermal stability have been designed (Koga et al., Nature 2012;491:222-227). By using structural databases of proteins isolated from either mesophilic or thermophilic organisms, we here investigate the influence of supersecondary associations on the thermal stability of natural proteins. Our results suggest that β-hairpins of proteins from thermophilic organisms are very frequently characterized by shortenings of the loops. Interestingly, this shortening leads to states that display a very strong preference for the most common connectivity of the strands observed in native protein hairpins. The abundance of selective states in these proteins suggests that they may achieve a high stability by adopting a strategy aimed to reduce the possible conformations of the unfolded ensemble. In this scenario, our data indicate that the shortening is effective if it increases the adherence to these rules. We also show that this mechanism may operate in the stabilization of well-known protein folds (thioredoxin and RNase A). These findings suggest that future investigations aimed at defining mechanism of protein stabilization should also consider these effects.

A novel mercuric reductase from the unique deep brine environment of Atlantis II in the Red Sea

A unique combination of physicochemical conditions prevails in the lower convective layer (LCL) of the brine pool at Atlantis II (ATII) Deep in the Red Sea. With a maximum depth of over 2000 m, the pool is characterized by acidic pH (5.3), high temperature (68 °C), salinity (26%), low light levels, anoxia, and high concentrations of heavy metals. We have established a metagenomic dataset derived from the microbial community in the LCL, and here we describe a gene for a novel mercuric reductase, a key component of the bacterial detoxification system for mercuric and organomercurial species. The metagenome-derived gene and an ortholog from an uncultured soil bacterium were synthesized and expressed in Escherichia coli. The properties of their products show that, in contrast to the soil enzyme, the ATII-LCL mercuric reductase is functional in high salt, stable at high temperatures, resistant to high concentrations of Hg(2+), and efficiently detoxifies Hg(2+) in vivo. Interestingly, despite the marked functional differences between the orthologs, their amino acid sequences differ by less than 10%. Site-directed mutagenesis and kinetic analysis of the mutant enzymes, in conjunction with three-dimensional modeling, have identified distinct structural features that contribute to extreme halophilicity, thermostability, and high detoxification capacity, suggesting that these were acquired independently during the evolution of this enzyme. Thus, our work provides fundamental structural insights into a novel protein that has undergone multiple biochemical and biophysical adaptations to promote the survival of microorganisms that reside in the extremely demanding environment of the ATII-LCL.

Engineering proteins for thermostability through rigidifying flexible sites

Engineering proteins for thermostability is an exciting and challenging field since it is critical for broadening the industrial use of recombinant proteins. Thermostability of proteins arises from the simultaneous effect of several forces such as hydrophobic interactions, disulfide bonds, salt bridges and hydrogen bonds. All of these interactions lead to decreased flexibility of polypeptide chain. Structural studies of mesophilic and thermophilic proteins showed that the latter need more rigid structures to compensate for increased thermal fluctuations. Hence flexibility can be an indicator to pinpoint weak spots for enhancing thermostability of enzymes. A strategy has been proven effective in enhancing proteins' thermostability with two steps: predict flexible sites of proteins firstly and then rigidify these sites. We refer to this approach as rigidify flexible sites (RFS) and give an overview of such a method through summarizing the methods to predict flexibility of a protein, the methods to rigidify residues with high flexibility and successful cases regarding enhancing thermostability of proteins using RFS.

Crystal structures of Phanerochaete chrysosporium pyranose 2-oxidase suggest that the N-terminus acts as a propeptide that assists in homotetramer assembly

The flavin-dependent homotetrameric enzyme pyranose 2-oxidase (P2O) is found mostly, but not exclusively, in lignocellulose-degrading fungi where it catalyzes the oxidation of β-d-glucose to the corresponding 2-keto sugar concomitantly with hydrogen peroxide formation during lignin solubilization. Here, we present crystal structures of P2O from the efficient lignocellulolytic basidiomycete Phanerochaete chrysosporium. Structures were determined of wild-type PcP2O from the natural fungal source, and two variants of recombinant full-length PcP2O, both in complex with the slow substrate 3-deoxy-3-fluoro-β-d-glucose. The active sites in PcP2O and P2O from Trametes multicolor (TmP2O) are highly conserved with identical substrate binding. Our structural analysis suggests that the 17 °C higher melting temperature of PcP2O compared to TmP2O is due to an increased number of intersubunit salt bridges. The structure of recombinant PcP2O expressed with its natural N-terminal sequence, including a proposed propeptide segment, reveals that the first five residues of the propeptide intercalate at the interface between A and B subunits to form stabilizing, mainly hydrophobic, interactions. In the structure of mature PcP2O purified from the natural source, the propeptide segment in subunit A has been replaced by a nearby loop in the B subunit. We propose that the propeptide in subunit A stabilizes the A/B interface of essential dimers in the homotetramer and that, upon maturation, it is replaced by the loop in the B subunit to form the mature subunit interface. This would imply that the propeptide segment of PcP2O acts as an intramolecular chaperone for oligomerization at the A/B interface of the essential dimer.

•

Structures of pyranose 2-oxidase from Phanerochaete chrysosporium were determined.

•

The N-terminus may act as a propeptide with a role in homotetramer assembly.

•

A large number of salt bridges between subunits provides thermostability.

•

The substrate is bound in the productive binding mode for oxidation at C2.

Role of key salt bridges in thermostability of G. thermodenitrificans EstGtA2: distinctive patterns within the new bacterial lipolytic enzyme subfamily XIII.2 [corrected]

Bacterial lipolytic enzymes were originally classified into eight different families defined by Arpigny and Jaeger (families I-VIII). Recently, the discovery of new lipolytic enzymes allowed for extending the original classification to fourteen families (I-XIV). We previously reported that G. thermodenitrificans EstGtA2 (access no. AEN92268) belonged to a novel group of bacterial lipolytic enzymes. Here we propose a 15(th) family (family XV) and suggest criteria for the assignation of protein sequences to the N' subfamily. Five selected salt bridges, hallmarks of the N' subfamily (E3/R54, E12/R37, E66/R140, D124/K178 and D205/R220) were disrupted in EstGtA2 using a combinatorial alanine-scanning approach. A set of 14 (R/K→A) mutants was produced, including five single, three double, three triple and three quadruple mutants. Despite a high tolerance to non-conservative mutations for folding, all the alanine substitutions were destabilizing (decreasing T m by 5 to 14°C). A particular combination of four substitutions exceeded this tolerance and prevents the correct folding of EstGtA2, leading to enzyme inactivation. Although other mutants remain active at low temperatures, the accumulation of more than two mutations had a dramatic impact on EstGtA2 activity at high temperatures suggesting an important role of these conserved salt bridge-forming residues in thermostability of lipolytic enzymes from the N' subfamily. We also identified a particular interloop salt bridge in EstGtA2 (D194/H222), located at position i -2 and i -4 residues from the catalytic Asp and His respectively which is conserved in other related bacterial lipolytic enzymes (families IV and XIII) with high tolerance to mutations and charge reversal. We investigated the role of residue identity at position 222 in controlling stability-pH dependence in EstGtA2. The introduction of a His to Arg mutation led to increase thermostability under alkaline pH. Our results suggest primary targets for optimization of EstGtA2 for specific biotechnological purposes.

Sequence analysis and structure prediction of enoyl-CoA hydratase from Avicennia marina: implication of various amino acid residues on substrate-enzyme interactions

Enoyl-CoA hydratase catalyzes the hydration of 2-trans-enoyl-CoA into 3-hydroxyacyl-CoA. The present study focuses on the correlation between the functional and structural aspects of enoyl-CoA hydratase from Avicennia marina. We have used bioinformatics tools to construct and analyze 3D homology models of A. marina enoyl-CoA hydratase (AMECH) bound to different substrates and inhibitors and studied the residues involved in the ligand-enzyme interaction. Structural information obtained from the models was compared with those of the reported crystal structures. We observed that the overall folds were similar; however, AMECH showed few distinct structural changes which include structural variation in the mobile loop, formation and loss of certain interactions between the active site residues and substrates. Some changes were also observed within specific regions of the enzyme. Glu106 is almost completely conserved in sequences of the isomerases/hydratases including AMECH while Glu86 which is the other catalytic residue in most of the isomerases/hydratases is replaced by Gly and shows no interaction with the substrate. Asp114 is located within 4Å distance of the catalytic water which makes it a probable candidate for the second catalytic residue in AMECH. Another prominent feature of AMECH is the presence of structurally distinct mobile loop having a completely different coordination with the hydrophobic binding pocket of acyl portion of the substrate.

Molecular dynamics perspective on the protein thermal stability: a case study using SAICAR synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) Cα, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.