Thermodynamic stability and folding of proteins from hyperthermophilic organisms

Cellular and epigenetic perspective of protein stability and its implications in the biological system

Protein stability is a fundamental prerequisite in both experimental and therapeutic applications. Current advancements in high throughput experimental techniques and functional ontology approaches have elucidated that impairment in the structure and stability of proteins is intricately associated with the cause and cure of several diseases. Therefore, it is paramount to deeply understand the physical and molecular confounding factors governing the stability of proteins. In this review article, we comprehensively investigated the evolution of protein stability, examining its emergence over time, its relationship with organizational aspects and the experimental methods used to understand it. Furthermore, we have also emphasized the role of Epigenetics and its interplay with post-translational modifications (PTMs) in regulating the stability of proteins.

Genomic estimates of mutation and substitution rates contradict the evolutionary speed hypothesis of the latitudinal diversity gradient

The latitudinal diversity gradient (LDG) refers to a decrease in biodiversity from the equator to the poles. The evolutionary speed hypothesis, backed by the metabolic theory of ecology, asserts that nucleotide mutation and substitution rates per site per year are higher and thereby speciation rates are higher at higher temperatures, generating the LDG. However, prior empirical investigations of the relationship between the temperature and mutation or substitution rate were based on a few genes and the results were mixed. We here revisit this relationship using genomic data. No significant correlation between the temperature and mutation rate is found in 13 prokaryotes or in 107 eukaryotes. An analysis of 234 diverse trios of bacterial taxa indicates that the synonymous substitution rate is not significantly associated with the growth temperature. The same data, however, reveal a significant negative association between the nonsynonymous substitution rate and temperature, which is explainable by a larger fraction of detrimental nonsynonymous mutations at higher temperatures due to a stronger demand for protein stability. We conclude that the evolutionary speed hypothesis of the LDG is unsupported by genomic data and advise that future mechanistic studies of the LDG should focus on other hypotheses.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Structural and functional adaptation in extremophilic microbial α-amylases

Maintaining stable native conformation of a protein under a given ecological condition is the prerequisite for survival of organisms. Extremophilic bacteria and archaea have evolved to adapt under extreme conditions of temperature, pH, salt, and pressure. Molecular adaptations of proteins under these conditions are essential for their survival. These organisms have the capability to maintain stable, native conformations of proteins under extreme conditions. The enzymes produced by the extremophiles are also known as extremozyme, which are used in several industries. Stability and functionality of extremozymes under varying temperature, pH, and solvent conditions are the most desirable requirement of industry. α-Amylase is one of the most important enzymes used in food, pharmaceutical, textile, and detergent industries. This enzyme is produced by diverse microorganisms including various extremophiles. Therefore, understanding its stability is important from fundamental as well as an applied point of view. Each class of extremophiles has a distinctive set of dominant non-covalent interactions which are important for their stability. Static information obtained by comparative analysis of amino acid sequence and atomic resolution structure provides information on the prevalence of particular amino acids or a group of non-covalent interactions. Protein folding studies give the information about thermodynamic and kinetic stability in order to understand dynamic aspect of molecular adaptations. In this review, we have summarized information on amino acid sequence, structure, stability, and adaptability of α-amylases from different classes of extremophiles.

Differential effect of polyol and sugar osmolytes on the refolding of homologous alpha amylases: A comparative study

Polyol and sugar osmolytes are known to enhance the stability of proteins, however, their role in assisting protein folding is not well understood. We asked whether these osmolytes have the same effect during refolding of a pair of thermophilic and mesophilic proteins. Herein, we have chosen α-amylases from Bacillus licheniformis (BLA) and Bacillus amyloliquefaciens (BAA) as thermophilic like and mesophilic counterparts respectively, having similar structures but differing thermostability. The effect of a series of polyols with varying number of -OH groups from 2 to 6 (Ethylene glycol, glycerol, erythritol, xylitol and sorbitol) and sugars (trehalose and sucrose) has been studied on the refolding of BLA and BAA. Our study demonstrates that glycerol, sorbitol and trehalose are the efficient cosolvents for BAA refolding, while comparatively less effective for BLA. Urea induced destabilization of BLA and BAA is differently compensated by polyol and sugar osmolytes during refolding. This suggests that the early species formed during BLA and BAA refolding are differently susceptible towards urea, indicating differential nature of their refolding pathways. Addition of trehalose at different times during refolding showed that the presence of trehalose is essential at the early stages of refolding. It is one of the first systematic study wherein the comparative effect of polyol and sugar assisted refolding of thermophilic and mesophilic protein has been carried out. The study highlights the differential effect of protein-osmolyte interactions during refolding of thermophilic and mesophilic proteins which may have implications in protein formulations, refolding and inhibition of aggregation.

Thermal processing implications on the digestibility of meat, fish and seafood proteins

Thermal processing is an inevitable part of the processing and preparation of meat and meat products for human consumption. However, thermal processing techniques, both commercial and domestic, induce modifications in muscle proteins which can have implications for their digestibility. The nutritive value of muscle proteins is closely related to their digestibility in the gastrointestinal tract and is determined by the end products that it presents in the assimilable form (amino acids and small peptides) for the absorption. The present review examines how different thermal processing techniques, such as sous-vide, microwave, stewing, roasting, boiling, frying, grilling, and steam cooking, affect the digestibility of muscle proteins in the gastrointestinal tract. By altering the functional and structural properties of muscle proteins, thermal processing has the potential to influence the digestibility negatively or positively, depending on the processing conditions. Thermal processes such as sous-vide can induce favourable changes, such as partial unfolding or exposure of cleavage sites, in muscle proteins and improve their digestibility whereas processes such as stewing and roasting can induce unfavourable changes, such as protein aggregation, severe oxidation, cross linking or increased disulfide (S-S) content and decrease the susceptibility of proteins during gastrointestinal digestion. The review examines how the underlying mechanisms of different processing conditions can be translated into higher or lower protein digestibility in detail. This review expands the current understanding of muscle protein digestion and generates knowledge that will be indispensable for optimizing the digestibility of thermally processed muscle foods for maximum nutritional benefits and optimal meal planning.

The Molecular Determinants of Thermoadaptation: Methanococcales as a Case Study

Previous reports have shown that environmental temperature impacts proteome evolution in Bacteria and Archaea. However, it is unknown whether thermoadaptation mainly occurs via the sequential accumulation of substitutions, massive horizontal gene transfers, or both. Measuring the real contribution of amino acid substitution to thermoadaptation is challenging, because of confounding environmental and genetic factors (e.g., pH, salinity, genomic G + C content) that also affect proteome evolution. Here, using Methanococcales, a major archaeal lineage, as a study model, we show that optimal growth temperature is the major factor affecting variations in amino acid frequencies of proteomes. By combining phylogenomic and ancestral sequence reconstruction approaches, we disclose a sequential substitutional scheme in which lysine plays a central role by fine tuning the pool of arginine, serine, threonine, glutamine, and asparagine, whose frequencies are strongly correlated with optimal growth temperature. Finally, we show that colonization to new thermal niches is not associated with high amounts of horizontal gene transfers. Altogether, although the acquisition of a few key proteins through horizontal gene transfer may have favored thermoadaptation in Methanococcales, our findings support sequential amino acid substitutions as the main factor driving thermoadaptation.

Structural and thermodynamic analysis of factors governing the stability and thermal folding/unfolding of SazCA

Molecular basis of protein stability at different temperatures is a fundamental problem in protein science that is substantially far from being accurately and quantitatively solved as it requires an explicit knowledge of the temperature dependence of folding free energy of amino acid residues. In the present study, we attempted to gain insights into the thermodynamic stability of SazCA and its implications on protein folding/unfolding. We report molecular dynamics simulations of water solvated SazCA in a temperature range of 293-393 K to study the relationship between the thermostability and flexibility. Our structural analysis shows that the protein maintains the highest structural stability at 353 K and the protein conformations are highly flexible at temperatures above 353 K. Larger exposure of hydrophobic surface residues to the solvent medium for conformations beyond 353 K were identified from H-bond analysis. Higher number of secondary structure contents exhibited by SazCA at 353 K corroborated the conformations at 353 K to exhibit the highest thermal stability. The analysis of thermodynamics of protein stability revealed that the conformations that denature at higher melting temperatures tend to have greater maximum thermal stability. Our analysis shows that 353 K conformations have the highest melting temperature, which was found to be close to the experimental optimum temperature. The enhanced protein stability at 353 K due the least value of heat capacity at unfolding suggested an increase in folding. Comparative Gibbs free energy analysis and funnel shaped energy landscape confirmed a transition in folding/unfolding pathway of SazCA at 353 K.

A conserved π-helix plays a key role in thermoadaptation of catalysis in the glycoside hydrolase family 4

Here, we characterize the role of a π-helix in the molecular mechanisms underlying thermoadaptation in the glycoside hydrolase family 4 (GH4). The interspersed π-helix present in a subgroup is evolutionarily related to a conserved α-helix in other orthologs by a single residue insertion/deletion event. The insertional residue, Phe407, in a hyperthermophilic α-glucuronidase, makes specific interactions across the inter-subunit interface. In order to establish the sequence-structure-stability implications of the π-helix, the wild-type and the deletion variant (Δ407) were characterized. The variant showed a significant lowering of melting temperature and optimum temperature for the highest activity. Crystal structures of the proteins show a transformation of the π-helix to a continuous α-helix in the variant, identical to that in orthologs lacking this insertion. Thermodynamic parameters were determined from stability curves representing the temperature dependence of unfolding free energy. Though the proteins display maximum stabilities at similar temperatures, a higher melting temperature in the wild-type is achieved by a combination of higher enthalpy and lower heat capacity of unfolding. Comparisons of the structural changes, and the activity and thermodynamic profiles allow us to infer that specific non-covalent interactions, and the existence of residual structure in the unfolded state, are crucial determinants of its thermostability. These features permit the enzyme to balance the preservation of structure at a higher temperature with the thermodynamic stability required for optimum catalysis.

Survival of the cheapest: how proteome cost minimization drives evolution

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

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

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

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.

Remote oxidative modifications induced by oxygen free radicals modify T/R allosteric equilibrium of a hyperthermophilic lactate dehydrogenase

L-Lactate dehydrogenase (LDH) is a model protein allowing to shed light on the fundamental molecular mechanisms that drive the acquisition, evolution and regulation of enzyme properties. In this study, we test the hypothesis of a link between thermal stability of LDHs and their capacity against unfolding induced by reactive oxygen species (ROS) generated by γ-rays irradiation. By using circular dichroism spectroscopy, we analysed that high thermal stability of a thermophilic LDH favours strong resistance against ROS-induced unfolding, in contrast to its psychrophilic and mesophilic counterparts that are less resistant. We suggest that a protein's phenotype linking strong thermal stability and resistance against ROS damages would have been a selective evolutionary advantage. We also find that the enzymatic activity of the thermophilic LDH that is strongly resistant against ROS-unfolding is very sensitive to inactivation by irradiation. To address this counter-intuitive observation, we combined mass spectrometry analyses and enzymatic activity measurements. We demonstrate that the dramatic change on LDH activity was linked to remote chemical modifications away from the active site, that change the equilibrium between low-affinity tense (T-inactive) and high-affinity relaxed (R-active) forms. We found the T-inactive thermophilic enzyme obtained after irradiation can recover its LDH activity by addition of the allosteric effector 1, 6 fructose bis phosphate. We analyse our data within the general framework of allosteric regulation, which requires that an enzyme in solution populates a large diversity of dynamically-interchanging conformations. Our work demonstrates that the radiation-induced inactivation of an enzyme is controlled by its dynamical properties.

Computer-Assisted Recombination (CompassR) Teaches us How to Recombine Beneficial Substitutions from Directed Evolution Campaigns

A main remaining challenge in protein engineering is how to recombine beneficial substitutions. Systematic recombination studies show that poorly performing variants are usually obtained after recombination of 3 to 4 beneficial substitutions. This limits researchers in exploiting nature's potential in generating better enzymes. The Computer-assisted Recombination (CompassR) strategy provides a selection guide for beneficial substitutions that can be recombined to gradually improve enzyme performance by analysis of the relative free energy of folding (ΔΔGfold ). The performance of CompassR was evaluated by analysis of 84 recombinants located on 13 positions of Bacillus subtilis lipase A. The finally obtained variant F17S/V54K/D64N/D91E had a 2.7-fold improved specific activity in 18.3 % (v/v) 1-butyl-3-methylimidazolium chloride ([BMIM][Cl]). In essence, the deducted CompassR rule allows recombination of beneficial substitutions in an iterative manner and empowers researchers to generate better enzymes in a time-efficient manner.

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.

Thermophilic Proteins as Versatile Scaffolds for Protein Engineering

Literature from the past two decades has outlined the existence of a trade-off between protein stability and function. This trade-off creates a unique challenge for protein engineers who seek to introduce new functionality to proteins. These engineers must carefully balance the mutation-mediated creation and/or optimization of function with the destabilizing effect of those mutations. Subsequent research has shown that protein stability is positively correlated with "evolvability" or the ability to support mutations which bestow new functionality on the protein. Since the ultimate goal of protein engineering is to create and/or optimize a protein's function, highly stable proteins are preferred as potential scaffolds for protein engineering. This review focuses on the application potential for thermophilic proteins as scaffolds for protein engineering. The relatively high inherent thermostability of these proteins grants them a great deal of mutational robustness, making them promising scaffolds for various protein engineering applications. Comparative studies on the evolvability of thermophilic and mesophilic proteins have strongly supported the argument that thermophilic proteins are more evolvable than mesophilic proteins. These findings indicate that thermophilic proteins may represent the scaffold of choice for protein engineering in the future.

Review: Engineering of thermostable enzymes for industrial applications

The catalytic properties of some selected enzymes have long been exploited to carry out efficient and cost-effective bioconversions in a multitude of research and industrial sectors, such as food, health, cosmetics, agriculture, chemistry, energy, and others. Nonetheless, for several applications, naturally occurring enzymes are not considered to be viable options owing to their limited stability in the required working conditions. Over the years, the quest for novel enzymes with actual potential for biotechnological applications has involved various complementary approaches such as mining enzyme variants from organisms living in extreme conditions (extremophiles), mimicking evolution in the laboratory to develop more stable enzyme variants, and more recently, using rational, computer-assisted enzyme engineering strategies. In this review, we provide an overview of the most relevant enzymes that are used for industrial applications and we discuss the strategies that are adopted to enhance enzyme stability and/or activity, along with some of the most relevant achievements. In all living species, many different enzymes catalyze fundamental chemical reactions with high substrate specificity and rate enhancements. Besides specificity, enzymes also possess many other favorable properties, such as, for instance, cost-effectiveness, good stability under mild pH and temperature conditions, generally low toxicity levels, and ease of termination of activity. As efficient natural biocatalysts, enzymes provide great opportunities to carry out important chemical reactions in several research and industrial settings, ranging from food to pharmaceutical, cosmetic, agricultural, and other crucial economic sectors.

Thermal stabilization of the deglycating enzyme Amadoriase I by rational design

Amadoriases are a class of FAD-dependent enzymes that are found in fungi, yeast and bacteria and that are able to hydrolyze glycated amino acids, cleaving the sugar moiety from the amino acidic portion. So far, engineered Amadoriases have mostly found practical application in the measurement of the concentration of glycated albumin in blood samples. However, these engineered forms of Amadoriases show relatively low absolute activity and stability levels, which affect their conditions of use. Therefore, enzyme stabilization is desirable prior to function-altering molecular engineering. In this work, we describe a rational design strategy based on a computational screening method to evaluate a library of potentially stabilizing disulfide bonds. Our approach allowed the identification of two thermostable Amadoriase I mutants (SS03 and SS17) featuring a significantly higher T₅₀ (55.3 °C and 60.6 °C, respectively) compared to the wild-type enzyme (52.4 °C). Moreover, SS17 shows clear hyperstabilization, with residual activity up to 95 °C, whereas the wild-type enzyme is fully inactive at 55 °C. Our computational screening method can therefore be considered as a promising approach to expedite the design of thermostable enzymes.

Cell-wide analysis of protein thermal unfolding reveals determinants of thermostability

Temperature-induced cell death is thought to be due to protein denaturation, but the determinants of thermal sensitivity of proteomes remain largely uncharacterized. We developed a structural proteomic strategy to measure protein thermostability on a proteome-wide scale and with domain-level resolution. We applied it to Escherichia coli, Saccharomyces cerevisiae, Thermus thermophilus, and human cells, yielding thermostability data for more than 8000 proteins. Our results (i) indicate that temperature-induced cellular collapse is due to the loss of a subset of proteins with key functions, (ii) shed light on the evolutionary conservation of protein and domain stability, and (iii) suggest that natively disordered proteins in a cell are less prevalent than predicted and (iv) that highly expressed proteins are stable because they are designed to tolerate translational errors that would lead to the accumulation of toxic misfolded species.

Protein folding at extreme temperatures: Current issues

The range of temperatures compatible with life is currently estimated from -25°C, as exemplified by metabolically active bacteria between sea ice crystals, and up to 122°C in hydrothermal vents as exemplified by the archaeon Methanopyrus kandleri. In the context of protein folding, as soon as a polypeptide emerges from the ribosome, it is exposed to the effects of environmental temperatures. Recent investigations have shown that the rate of protein folding is not adapted to extreme temperatures and should be very fast at high temperature and low in cold environments. This lack of adaptation is driven by kinetic constraints on protein stability. To counteract the deleterious effects of fast protein folding in hyperthermophiles, chaperones such as the Trigger Factor hold and slow down the rate of folding intermediates. Prolyl isomerization, a rate-limiting step in the folding of many proteins, is strongly temperature-dependent and impairs folding of psychrophilic proteins in the cold. This is compensated by reduction of the proline content in cold-adapted proteins, by an increased number of prolyl isomerases encoded in the genome of psychrophilic microorganisms and by overexpression of prolyl isomerases under low temperature cultivation. After folding, the native state is reached and although extremophilic proteins share the same fold, dramatic differences in stability have been recorded by differential scanning calorimetry.

An optimized strategy to measure protein stability highlights differences between cold and hot unfolded states

Macromolecular crowding ought to stabilize folded forms of proteins, through an excluded volume effect. This explanation has been questioned and observed effects attributed to weak interactions with other cell components. Here we show conclusively that protein stability is affected by volume exclusion and that the effect is more pronounced when the crowder's size is closer to that of the protein under study. Accurate evaluation of the volume exclusion effect is made possible by the choice of yeast frataxin, a protein that undergoes cold denaturation above zero degrees, because the unfolded form at low temperature is more expanded than the corresponding one at high temperature. To achieve optimum sensitivity to changes in stability we introduce an empirical parameter derived from the stability curve. The large effect of PEG 20 on cold denaturation can be explained by a change in water activity, according to Privalov's interpretation of cold denaturation.

Derivative of Extremophilic 50S Ribosomal Protein L35Ae as an Alternative Protein Scaffold

Small antibody mimetics, or alternative binding proteins (ABPs), extend and complement antibody functionality with numerous applications in research, diagnostics and therapeutics. Given the superiority of ABPs, the last two decades have witnessed development of dozens of alternative protein scaffolds (APSs) for the design of ABPs. Proteins from extremophiles with their high structural stability are especially favorable for APS design. Here, a 10X mutant of the 50S ribosomal protein L35Ae from hyperthermophilic archaea Pyrococcus horikoshii has been probed as an APS. A phage display library of L35Ae 10X was generated by randomization of its three CDR-like loop regions (repertoire size of 2×108). Two L35Ae 10X variants specific to a model target, the hen egg-white lysozyme (HEL), were isolated from the resulting library using phage display. The affinity of these variants (L4 and L7) to HEL ranges from 0.10 μM to 1.6 μM, according to surface plasmon resonance data. While L4 has 1-2 orders of magnitude lower affinity to HEL homologue, bovine α-lactalbumin (BLA), L7 is equally specific to HEL and BLA. The reference L35Ae 10X is non-specific to both HEL and BLA. L4 and L7 are more resistant to denaturation by guanidine hydrochloride compared to the reference L35Ae 10X (mid-transition concentration is higher by 0.1-0.5 M). Chemical crosslinking experiments reveal an increased propensity of L4 and L7 to multimerization. Overall, the CDR-like loop regions of L35Ae 10X represent a proper interface for generation of functional ABPs. Hence, L35Ae is shown to extend the growing family of protein scaffolds dedicated to the design of novel binding proteins.

Probing the Folding-Unfolding Transition of a Thermophilic Protein, MTH1880

The folding mechanism of typical proteins has been studied widely, while our understanding of the origin of the high stability of thermophilic proteins is still elusive. Of particular interest is how an atypical thermophilic protein with a novel fold maintains its structure and stability under extreme conditions. Folding-unfolding transitions of MTH1880, a thermophilic protein from Methanobacterium thermoautotrophicum, induced by heat, urea, and GdnHCl, were investigated using spectroscopic techniques including circular dichorism, fluorescence, NMR combined with molecular dynamics (MD) simulations. Our results suggest that MTH1880 undergoes a two-state N to D transition and it is extremely stable against temperature and denaturants. The reversibility of refolding was confirmed by spectroscopic methods and size exclusion chromatography. We found that the hyper-stability of the thermophilic MTH1880 protein originates from an extensive network of both electrostatic and hydrophobic interactions coordinated by the central β-sheet. Spectroscopic measurements, in combination with computational simulations, have helped to clarify the thermodynamic and structural basis for hyper-stability of the novel thermophilic protein MTH1880.

Polarity Alteration of a Calcium Site Induces a Hydrophobic Interaction Network and Enhances Cel9A Endoglucanase Thermostability

Structural calcium sites control protein thermostability and activity by stabilizing native folds and changing local conformations. Alicyclobacillus acidocaldarius survives in thermal-acidic conditions and produces an endoglucanase Cel9A (AaCel9A) which contains a calcium-binding site (Ser465 to Val470) near the catalytic cleft. By superimposing the Ca(2+)-free and Ca(2+)-bounded conformations of the calcium site, we found that Ca(2+) induces hydrophobic interactions between the calcium site and its nearby region by driving a conformational change. The hydrophobic interactions at the high-B-factor region could be enhanced further by replacing the surrounding polar residues with hydrophobic residues to affect enzyme thermostability and activity. Therefore, the calcium-binding residue Asp468 (whose side chain directly ligates Ca(2+)), Asp469, and Asp471 of AaCel9A were separately replaced by alanine and valine. Mutants D468A and D468V showed increased activity compared with those of the wild type with 0 mM or 10 mM Ca(2+) added, whereas the Asp469 or Asp471 substitution resulted in decreased activity. The D468A crystal structure revealed that mutation D468A triggered a conformational change similar to that induced by Ca(2+) in the wild type and developed a hydrophobic interaction network between the calcium site and the neighboring hydrophobic region (Ala113 to Ala117). Mutations D468V and D468A increased 4.5°C and 5.9°C, respectively, in melting temperature, and enzyme half-life at 75°C increased approximately 13 times. Structural comparisons between AaCel9A and other endoglucanases of the GH9 family suggested that the stability of the regions corresponding to the AaCel9A calcium site plays an important role in GH9 endoglucanase catalysis at high temperature.

An acidified thermostabilizing mini-peptide derived from the carboxyl extension of the larger isoform of the plant Rubisco activase

Thermostable fusion peptide partners are valuable in engineering thermostability in proteins. We evaluated the Arabidopsis counterpart (AtRAce) and an acidified derivative (mRAce) of the conserved carboxyl extension (RAce) of plant Rubisco activase (RCA) for their thermostabilizing properties in Escherichia coli and Saccharomyces cerevisiae using a protein fusion strategy. We used AtRAce and mRAce as fusion tails for the thermolabile protein RCA2 from Arabidopsis thaliana and Nicotiana tabacum. The homologous fusion of AtRAce with Arabidopsis RCA2 and the heterologous fusion of AtRAce with tobacco RCA2 increased the thermostability of both proteins. The acidified derivative mRAce conferred greater thermostability upon both proteins as compared with AtRAce. Moreover, mRAce enhanced the thermostability of other two thermolabile proteins from Jatropha curcas: the cytosolic ascorbate peroxidase 1 (JcAPX1) and the TATA-box binding protein isoform 1 (JcTBP1). We further report - for the first time - that JcTBP1 mediates heat tolerance in vivo in yeast. Thus, our study identifies a C-terminal acidic mini-peptide - the acidified derivative mRAce - with potential uses in improving the thermostability of heat-labile proteins and their associated heat tolerance in host organisms.

Extremophilic 50S Ribosomal RNA-Binding Protein L35Ae as a Basis for Engineering of an Alternative Protein Scaffold

Due to their remarkably high structural stability, proteins from extremophiles are particularly useful in numerous biological applications. Their utility as alternative protein scaffolds could be especially valuable in small antibody mimetic engineering. These artificial binding proteins occupy a specific niche between antibodies and low molecular weight substances, paving the way for development of innovative approaches in therapeutics, diagnostics, and reagent use. Here, the 50S ribosomal RNA-binding protein L35Ae from the extremophilic archaea Pyrococcus horikoshii has been probed for its potential to serve as a backbone in alternative scaffold engineering. The recombinant wild type L35Ae has a native-like secondary structure, extreme thermal stability (mid-transition temperature of 90°C) and a moderate resistance to the denaturation by guanidine hydrochloride (half-transition at 2.6 M). Chemical crosslinking and dynamic light scattering data revealed that the wild type L35Ae protein has a propensity for multimerization and aggregation correlating with its non-specific binding to a model cell surface of HEK293 cells, as evidenced by flow cytometry. To suppress these negative features, a 10-amino acid mutant (called L35Ae 10X) was designed, which lacks the interaction with HEK293 cells, is less susceptible to aggregation, and maintains native-like secondary structure and thermal stability. However, L35Ae 10X also shows lowered resistance to guanidine hydrochloride (half-transition at 2.0M) and is more prone to oligomerization. This investigation of an extremophile protein’s scaffolding potential demonstrates that lowered resistance to charged chemical denaturants and increased propensity to multimerization may limit the utility of extremophile proteins as alternative scaffolds.

Thermal stabilization of dihydrofolate reductase using monte carlo unfolding simulations and its functional consequences

Design of proteins with desired thermal properties is important for scientific and biotechnological applications. Here we developed a theoretical approach to predict the effect of mutations on protein stability from non-equilibrium unfolding simulations. We establish a relative measure based on apparent simulated melting temperatures that is independent of simulation length and, under certain assumptions, proportional to equilibrium stability, and we justify this theoretical development with extensive simulations and experimental data. Using our new method based on all-atom Monte-Carlo unfolding simulations, we carried out a saturating mutagenesis of Dihydrofolate Reductase (DHFR), a key target of antibiotics and chemotherapeutic drugs. The method predicted more than 500 stabilizing mutations, several of which were selected for detailed computational and experimental analysis. We find a highly significant correlation of r = 0.65–0.68 between predicted and experimentally determined melting temperatures and unfolding denaturant concentrations for WT DHFR and 42 mutants. The correlation between energy of the native state and experimental denaturation temperature was much weaker, indicating the important role of entropy in protein stability. The most stabilizing point mutation was D27F, which is located in the active site of the protein, rendering it inactive. However for the rest of mutations outside of the active site we observed a weak yet statistically significant positive correlation between thermal stability and catalytic activity indicating the lack of a stability-activity tradeoff for DHFR. By combining stabilizing mutations predicted by our method, we created a highly stable catalytically active E. coli DHFR mutant with measured denaturation temperature 7.2°C higher than WT. Prediction results for DHFR and several other proteins indicate that computational approaches based on unfolding simulations are useful as a general technique to discover stabilizing mutations.

All-atom molecular simulations have provided valuable insight into the workings of molecular machines and the folding and unfolding of proteins. However, commonly employed molecular dynamics simulations suffer from a limitation in accessible time scale, making it difficult to model large-scale unfolding events in a realistic amount of simulation time without employing unrealistically high temperatures. Here, we describe a rapid all-atom Monte Carlo simulation approach to simulate unfolding of the essential bacterial enzyme Dihydrofolate Reductase (DHFR) and all possible single point-mutants. We use these simulations to predict which mutants will be more thermodynamically stable (i.e., reside more often in the native folded state vs. the unfolded state) than the wild-type protein, and we confirm our predictions experimentally, creating several highly stable and catalytically active mutants. Thermally stable active engineered proteins can be used as a starting point in directed evolution experiments to evolve new functions on the background of this additional “reservoir of stability.” The stabilized enzyme may be able to accumulate a greater number of destabilizing yet functionally important mutations before unfolding, protease digestion, and aggregation abolish its activity.

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.

Formation of high-order oligomers by a hyperthemostable Fe-superoxide dismutase (tcSOD)

Hyperthermostable proteins are highly resistant to various extreme conditions. Many factors have been proposed to contribute to their ultrahigh structural stability. Some thermostable proteins have larger oligomeric size when compared to their mesophilic homologues. The formation of compact oligomers can minimize the solvent accessible surface area and increase the changes of Gibbs free energy for unfolding. Similar to mesophilic proteins, hyperthermostable proteins also face the problem of unproductive aggregation. In this research, we investigated the role of high-order oligomerization in the fight against aggregation by a hyperthermostable superoxide dismutase identified from Tengchong, China (tcSOD). Besides the predominant tetramers, tcSOD could also form active high-order oligomers containing at least eight subunits. The dynamic equilibrium between tetramers and high-order oligomers was not significantly affected by pH, salt concentration or moderate temperature. The secondary and tertiary structures of tcSOD remained unchanged during heating, while cross-linking experiments showed that there were conformational changes or structural fluctuations at high temperatures. Mutational analysis indicated that the last helix at the C-terminus was involved in the formation of high-order oligomers, probably via domain swapping. Based on these results, we proposed that the reversible conversion between the active tetramers and high-order oligomers might provide a buffering system for tcSOD to fight against the irreversible protein aggregation pathway. The formation of active high-order oligomers not only increases the energy barrier between the native state and unfolded/aggregated state, but also provides the enzyme the ability to reproduce the predominant oligomers from the active high-order oligomers.

Protein thermodynamics can be predicted directly from biological growth rates

Life on Earth is capable of growing from temperatures well below freezing to above the boiling point of water, with some organisms preferring cooler and others hotter conditions. The growth rate of each organism ultimately depends on its intracellular chemical reactions. Here we show that a thermodynamic model based on a single, rate-limiting, enzyme-catalysed reaction accurately describes population growth rates in 230 diverse strains of unicellular and multicellular organisms. Collectively these represent all three domains of life, ranging from psychrophilic to hyperthermophilic, and including the highest temperature so far observed for growth (122°C). The results provide credible estimates of thermodynamic properties of proteins and obtain, purely from organism intrinsic growth rate data, relationships between parameters previously identified experimentally, thus bridging a gap between biochemistry and whole organism biology. We find that growth rates of both unicellular and multicellular life forms can be described by the same temperature dependence model. The model results provide strong support for a single highly-conserved reaction present in the last universal common ancestor (LUCA). This is remarkable in that it means that the growth rate dependence on temperature of unicellular and multicellular life forms that evolved over geological time spans can be explained by the same model.

Circular dichroism and fluorescence spectroscopy of cysteinyl-tRNA synthetase from Halobacterium salinarum ssp. NRC-1 demonstrates that group I cations are particularly effective in providing structure and stability to this halophilic protein

Proteins from extremophiles have the ability to fold and remain stable in their extreme environment. Here, we investigate the presence of this effect in the cysteinyl-tRNA synthetase from Halobacterium salinarum ssp. NRC-1 (NRC-1), which was used as a model halophilic protein. The effects of salt on the structure and stability of NRC-1 and of E. coli CysRS were investigated through far-UV circular dichroism (CD) spectroscopy, fluorescence spectroscopy, and thermal denaturation melts. The CD of NRC-1 CysRS was examined in different group I and group II chloride salts to examine the effects of the metal ions. Potassium was observed to have the strongest effect on NRC-1 CysRS structure, with the other group I salts having reduced strength. The group II salts had little effect on the protein. This suggests that the halophilic adaptations in this protein are mediated by potassium. CD and fluorescence spectra showed structural changes taking place in NRC-1 CysRS over the concentration range of 0-3 M KCl, while the structure of E. coli CysRS was relatively unaffected. Salt was also shown to increase the thermal stability of NRC-1 CysRS since the melt temperature of the CysRS from NRC-1 was increased in the presence of high salt, whereas the E. coli enzyme showed a decrease. By characterizing these interactions, this study not only explains the stability of halophilic proteins in extremes of salt, but also helps us to understand why and how group I salts stabilize proteins in general.

An efficient strategy for small-scale screening and production of archaeal membrane transport proteins in Escherichia coli

Background

Membrane proteins play a key role in many fundamental cellular processes such as transport of nutrients, sensing of environmental signals and energy transduction, and account for over 50% of all known drug targets. Despite their importance, structural and functional characterisation of membrane proteins still remains a challenge, partially due to the difficulties in recombinant expression and purification. Therefore the need for development of efficient methods for heterologous production is essential.

Methodology/Principal Findings

Fifteen integral membrane transport proteins from Archaea were selected as test targets, chosen to represent two superfamilies widespread in all organisms known as the Major Facilitator Superfamily (MFS) and the 5-Helix Inverted Repeat Transporter superfamily (5HIRT). These proteins typically have eleven to twelve predicted transmembrane helices and are putative transporters for sugar, metabolite, nucleobase, vitamin or neurotransmitter. They include a wide range of examples from the following families: Metabolite-H⁺-symporter; Sugar Porter; Nucleobase-Cation-Symporter-1; Nucleobase-Cation-Symporter-2; and neurotransmitter-sodium-symporter. Overproduction of transporters was evaluated with three vectors (pTTQ18, pET52b, pWarf) and two Escherichia coli strains (BL21 Star and C43 (DE3)). Thirteen transporter genes were successfully expressed; only two did not express in any of the tested vector-strain combinations. Initial trials showed that seven transporters could be purified and six of these yielded quantities of ≥ 0.4 mg per litre suitable for functional and structural studies. Size-exclusion chromatography confirmed that two purified transporters were almost homogeneous while four others were shown to be non-aggregating, indicating that they are ready for up-scale production and crystallisation trials.

Conclusions/Significance

Here, we describe an efficient strategy for heterologous production of membrane transport proteins in E. coli. Small-volume cultures (10 mL) produced sufficient amount of proteins to assess their purity and aggregation state. The methods described in this work are simple to implement and can be easily applied to many more membrane proteins.

Dynamics measured by neutron scattering correlates with the organization of bioenergetics complexes in natural membranes from hyperthermophile and mesophile bacteria

Various models on membrane structure and organization of proteins and complexes in natural membranes emerged during the last years. However, the lack of systematic dynamical studies to complement structural investigations hindered the establishment of a more complete picture of these systems. Elastic incoherent neutron scattering gives access to the dynamics on a molecular level and was applied to natural membranes extracted from the hyperthermophile Aquifex aeolicus and the mesophile Wolinella succinogenes bacteria. The results permitted to extract a hierarchy of dynamic flexibility and atomic resilience within the samples, which correlated with the organization of proteins in bioenergetics complexes and the functionality of the membranes.

Atypical features of Thermus thermophilus succinate:quinone reductase

The Thermus thermophilus succinate:quinone reductase (SQR), serving as the respiratory complex II, has been homologously produced under the control of a constitutive promoter and subsequently purified. The detailed biochemical characterization of the resulting wild type (wt-rcII) and His-tagged (rcII-His₈-SdhB and rcII-SdhB-His₆) complex II variants showed the same properties as the native enzyme with respect to the subunit composition, redox cofactor content and sensitivity to the inhibitors malonate, oxaloacetate, 3-nitropropionic acid and nonyl-4-hydroxyquinoline-N-oxide (NQNO). The position of the His-tag determined whether the enzyme retained its native trimeric conformation or whether it was present in a monomeric form. Only the trimer exhibited positive cooperativity at high temperatures. The EPR signal of the [2Fe-2S] cluster was sensitive to the presence of substrate and showed an increased rhombicity in the presence of succinate in the native and in all recombinant forms of the enzyme. The detailed analysis of the shape of this signal as a function of pH, substrate concentration and in the presence of various inhibitors and quinones is presented, leading to a model for the molecular mechanism that underlies the influence of succinate on the rhombicity of the EPR signal of the proximal iron-sulfur cluster.

Psychrophilic enzymes: from folding to function and biotechnology

Psychrophiles thriving permanently at near-zero temperatures synthesize cold-active enzymes to sustain their cell cycle. Genome sequences, proteomic, and transcriptomic studies suggest various adaptive features to maintain adequate translation and proper protein folding under cold conditions. Most psychrophilic enzymes optimize a high activity at low temperature at the expense of substrate affinity, therefore reducing the free energy barrier of the transition state. Furthermore, a weak temperature dependence of activity ensures moderate reduction of the catalytic activity in the cold. In these naturally evolved enzymes, the optimization to low temperature activity is reached via destabilization of the structures bearing the active site or by destabilization of the whole molecule. This involves a reduction in the number and strength of all types of weak interactions or the disappearance of stability factors, resulting in improved dynamics of active site residues in the cold. These enzymes are already used in many biotechnological applications requiring high activity at mild temperatures or fast heat-inactivation rate. Several open questions in the field are also highlighted.

Thermal stabilization of an endoglucanase by cyclization

An intein-driven protein splicing approach allowed for the covalent linkage between the N- and C-termini of a polypeptide chain to create circular variants of the endo-β-1,3-1,4-glucanase, LicA, from Bacillus licheniformis. Two circular variants, LicA-C1 and LicA-C2, which have connecting loops of 20 and 14 amino acids, respectively, showed catalytic activities that are approximately two and three times higher, respectively, compared to that of the linear LicA (LicA-L1). The thermal stability of the circular variants was significantly increased compared to the linear form. Whereas the linear glucanase lost half of its activity after 3 min at 65 °C, the two circular variants have 6-fold (LicA-C1) and 16-fold (LicA-C2) increased half-life time of inactivation. In agreement with this, fluorescence spectroscopy and differential scanning calorimetry studies revealed that circular enzymes undergo structural changes at higher temperatures compared to that of the linear form. The effect of calcium on the conformational stability and function of the circular LicAs was also investigated, and we observed that the presence of calcium ions results in increased thermal stability. The impact of the length of the designed loops on thermal stability of the circular proteins is discussed, and it is suggested that cyclization may be an efficient strategy for the increased stability of proteins.

Role of disulfide bonds in conformational stability and folding of 5'-deoxy-5'-methylthioadenosine phosphorylase II from the hyperthermophilic archaeon Sulfolobus solfataricus

Sulfolobus solfataricus 5'-deoxy-5'-melthylthioadenosine phosphorylase II (SsMTAPII), is a hyperthermophilic hexameric protein with two intrasubunit disulfide bonds (C138-C205 and C200-C262) and a CXC motif (C259-C261). To get information on the role played by these covalent links in stability and folding, the conformational stability of SsMTAPII and C262S and C259S/C261S mutants was studied by thermal and guanidinium chloride (GdmCl)-induced unfolding and analyzed by fluorescence spectroscopy, circular dichroism, and SDS-PAGE. No thermal unfolding transition of SsMTAPII can be obtained under nonreducing conditions, while in the presence of the reducing agent Tris-(2-carboxyethyl) phosphine (TCEP), a Tm of 100°C can be measured demonstrating the involvement of disulfide bridges in enzyme thermostability. Different from the wild-type, C262S and C259S/C261S show complete thermal denaturation curves with sigmoidal transitions centered at 102°C and 99°C respectively. Under reducing conditions these values decrease by 4°C and 8°C respectively, highlighting the important role exerted by the CXC disulfide on enzyme thermostability. The contribution of disulfide bonds to the conformational stability of SsMTAPII was further assessed by GdmCl-induced unfolding experiments carried out under reducing and nonreducing conditions. Thermal unfolding was found to be reversible if the protein was heated in the presence of TCEP up to 90°C but irreversible above this temperature because of aggregation. In analogy, only chemical unfolding carried out in the presence of reducing agents resulted in a reversible process suggesting that disulfide bonds play a role in enzyme denaturation. Thermal and chemical unfolding of SsMTAPII occur with dissociation of the native hexameric state into denatured monomers, as indicated by SDS-PAGE.