Do ultrastable proteins from hyperthermophiles have high or low conformational rigidity?

Efficient Base-Catalyzed Kemp Elimination in an Engineered Ancestral Enzyme

The routine generation of enzymes with completely new active sites is a major unsolved problem in protein engineering. Advances in this field have thus far been modest, perhaps due, at least in part, to the widespread use of modern natural proteins as scaffolds for de novo engineering. Most modern proteins are highly evolved and specialized and, consequently, difficult to repurpose for completely new functionalities. Conceivably, resurrected ancestral proteins with the biophysical properties that promote evolvability, such as high stability and conformational diversity, could provide better scaffolds for de novo enzyme generation. Kemp elimination, a non-natural reaction that provides a simple model of proton abstraction from carbon, has been extensively used as a benchmark in de novo enzyme engineering. Here, we present an engineered ancestral β-lactamase with a new active site that is capable of efficiently catalyzing Kemp elimination. The engineering of our Kemp eliminase involved minimalist design based on a single function-generating mutation, inclusion of an extra polypeptide segment at a position close to the de novo active site, and sharply focused, low-throughput library screening. Nevertheless, its catalytic parameters (k_cat/K_M~2·10⁵ M⁻¹ s⁻¹, k_cat~635 s⁻¹) compare favorably with the average modern natural enzyme and match the best proton-abstraction de novo Kemp eliminases that are reported in the literature. The general implications of our results for de novo enzyme engineering are discussed.

Improved thermostability of D-allulose 3-epimerase from Clostridium bolteae ATCC BAA-613 by proline residue substitution

d-allulose, a rare sugar that is scarce in nature, exerts several beneficial effects and has commercial potential. d-allulose 3-epimerase (DAEase) plays a vital role in catalyzing the isomerization from d-fructose to d-allulose. However, the industrial application of DAEase for d-allulose production is hindered by its poor long-term thermostability. In the present research, we introduced a proline residue (i) to restrict its spatial conformation and (ii) to reduce the entropy of the unfolded state of DAEase. The t_1/2 value of the double-site Clostridium bolteae DAEase mutant Cb-51P/89P was prolonged to 58 min at 55 °C, a 2.32-fold increase compared with wild-type DAEase. The manipulation did not cause obvious changes in the enzymatic properties, including optimum pH, optimal temperature, optimum metal ion, and enzymatic activity. As the accumulation of multiple small effects through proline substitution could dramatically improve the thermostability of the mutant protein, our method to improve the thermostability while roughly retaining the original enzymatic properties is promising.

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

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

Photosynthetic adaptation to polar life: Energy balance, photoprotection and genetic redundancy

The persistent low temperature that characterize polar habitats combined with the requirement for light for all photoautotrophs creates a conundrum. The absorption of too much light at low temperature can cause an energy imbalance that decreases photosynthetic performance that has a negative impact on growth and can affect long-term survival. The goal of this review is to survey the mechanism(s) by which polar photoautotrophs maintain cellular energy balance, that is, photostasis to overcome the potential for cellular energy imbalance in their low temperature environments. Photopsychrophiles are photosynthetic organisms that are obligately adapted to low temperature (0⁰- 15 °C) but usually die at higher temperatures (≥20 °C). In contrast, photopsychrotolerant species can usually tolerate and survive a broad range of temperatures (5⁰- 40 °C). First, we summarize the basic concepts of excess excitation energy, energy balance, photoprotection and photostasis and their importance to survival in polar habitats. Second, we compare the photoprotective mechanisms that underlie photostasis and survival in aquatic cyanobacteria and green algae as well as terrestrial Antarctic and Arctic plants. We show that polar photopsychrophilic and photopsychrotolerant organisms attain energy balance at low temperature either through a regulated reduction in the efficiency of light absorption or through enhanced capacity to consume photosynthetic electrons by the induction of O₂ as an alternative electron acceptor. Finally, we compare the published genomes of three photopsychrophilic and one photopsychrotolerant alga with five mesophilic green algae including the model green alga, Chlamydomonas reinhardtii. We relate our genomic analyses to photoprotective mechanisms that contribute to the potential attainment of photostasis. Finally, we discuss how the observed genomic redundancy in photopsychrophilic genomes may confer energy balance, photoprotection and resilience to their harsh polar environment. Primary production in aquatic, Antarctic and Arctic environments is dependent on diverse algal and cyanobacterial communities. Although mosses and lichens dominate the Antarctic terrestrial landscape, only two extant angiosperms exist in the Antarctic. The identification of a single 'molecular key' to unravel adaptation of photopsychrophily and photopsychrotolerance remains elusive. Since these photoautotrophs represent excellent biomarkers to assess the impact of global warming on polar ecosystems, increased study of these polar photoautotrophs remains essential.

Current Status of Mining, Modification, and Application of Cellulases in Bioactive Substance Extraction

Cellulases have been used to extract bioactive ingredients from medical plants; however, the poor enzymatic properties of current cellulases significantly limit their application. Two strategies are expected to address this concern: (1) new cellulase gene mining strategies have been promoted, optimized, and integrated, thanks to the improvement of gene sequencing, genomic data, and algorithm optimization, and (2) known cellulases are being modified, thanks to the development of protein engineering, crystal structure data, and computing power. Here, we focus on mining strategies and provide a systemic overview of two approaches based on sequencing and function. Strategies based on protein structure modification, such as introducing disulfide bonds, proline, salt bridges, N-glycosylation modification, and truncation of loop structures, have already been summarized. This review discusses four aspects of cellulase-assisted extraction. Initially, cellulase alone was used to extract bioactive substances, and later, mixed enzyme systems were developed. Physical methods such as ultrasound, microwave, and high hydrostatic pressure have assisted in improving extraction efficiency. Cellulase changes the structure of biomolecules during the extraction process to convert them into effective ingredients with better activity and bioavailability. The combination of cellulase with other enzymes and physical technologies is a promising strategy for future extraction applications.

An effective computational-screening strategy for simultaneously improving both catalytic activity and thermostability of α-l-rhamnosidase

Catalytic efficiency and thermostability are the two most important characteristics of enzymes. However, it is always tough to improve both catalytic efficiency and thermostability of enzymes simultaneously. In the present study, a computational strategy with double-screening steps was proposed to simultaneously improve both catalysis efficiency and thermostability of enzymes; and a fungal α-l-rhamnosidase was used to validate the strategy. As the result, by molecular docking and sequence alignment analysis within the binding pocket, seven mutant candidates were predicted with better catalytic efficiency. By energy variety analysis, A355N, S356Y, and D525N among the seven mutant candidates were predicted with better thermostability. The expression and characterization results showed the mutant D525N had significant improvements in both enzyme activity and thermostability. Molecular dynamics simulations indicated that the mutations located within the 5 Å range of the catalytic domain, which could improve root mean squared deviation, electrostatic, Van der Waal interaction, and polar salvation values, and formed water bridge between the substrate and the enzyme. The study indicated that the computational strategy based on the binding energy, conservation degree and mutation energy analyses was effective to develop enzymes with better catalysis and thermostability, providing practical approach for developing industrial enzymes.

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.

Protein thermal stability

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

Entropic contribution to enhanced thermal stability in the thermostable P450 CYP119

The enhanced thermostability of thermophilic proteins with respect to their mesophilic counterparts is often attributed to the enthalpy effect, arising from strong interactions between protein residues. Intuitively, these strong interresidue interactions will rigidify the biomolecules. However, the present work utilizing neutron scattering and solution NMR spectroscopy measurements demonstrates a contrary example that the thermophilic cytochrome P450, CYP119, is much more flexible than its mesophilic counterpart, CYP101A1, something which is not apparent just from structural comparison of the two proteins. A mechanism to explain this apparent contradiction is that higher flexibility in the folded state of CYP119 increases its conformational entropy and thereby reduces the entropy gain during denaturation, which will increase the free energy needed for unfolding and thus stabilize the protein. This scenario is supported by thermodynamic data on the temperature dependence of unfolding free energy, which shows a significant entropic contribution to the thermostability of CYP119 and lends an added dimension to enhanced stability, previously attributed only to presence of aromatic stacking interactions and salt bridge networks. Our experimental data also support the notion that highly thermophilic P450s such as CYP119 may use a mechanism that partitions flexibility differently from mesophilic P450s between ligand binding and thermal stability.

Protein stability and dynamics influenced by ligands in extremophilic complexes - a molecular dynamics investigation

In this study, we explore the structural and dynamic adaptations of the Tryptophan synthase α-subunit in a ligand bound state in psychrophilic, mesophilic and hyperthermophilic organisms at different temperatures by MD simulations. We quantify the global and local fluctuations in the 40 ns time scale by analyzing the root mean square deviation/fluctuations. The distinct behavior of the active site and loop 6 is observed with the elevation of temperature. Protein stability relies more on electrostatic interactions, and these interactions might be responsible for the stability of varying temperature evolved proteins. The paper also focuses on the effect of temperature on protein dynamics and stability governed by the distinct behavior of the ligand associated with its retention, binding and dissociation over the course of time. The integration of principle component analysis and a free energy landscape was useful in identifying the conformational space accessible to ligand bound homologues and how the presence of the ligand alters the conformational and dynamic properties of the protein.

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.

On the indirect relationship between protein dynamics and enzyme activity

The behaviors of simple thermal systems have been well studied in physical chemistry and the principles obtained from such studies have been applied to complex thermal systems, such as proteins and enzymes. But the simple application of such principles is questionable and may lead to mistakes under some circumstances. In enzymology, the transition state theory of chemical reactions has been accepted as a fundamental theory, but the role of protein dynamics in enzyme catalysis is controversial in the context of transition state theory. By studying behaviors of complex thermal systems, we have revised the Arrhenius equation and transition state theory and our model is validated in enzymology. Formally speaking, the revised Arrhenius equation is apparently similar to a conventional Arrhenius equation, but the physical meanings of its parameters differ from that of traditional forms in principle. Within this model, the role of protein dynamics in enzyme catalysis is well defined and quantified.

Pnserpin: A Novel Serine Protease Inhibitor from Extremophile Pyrobaculum neutrophilum

Serine protease inhibitors (serpins) are native inhibitors of serine proteases, constituting a large protein family with members spread over eukaryotes and prokaryotes. However, only very few prokaryotic serpins, especially from extremophiles, have been characterized to date. In this study, Pnserpin, a putative serine protease inhibitor from the thermophile Pyrobaculum neutrophilum, was overexpressed in Escherichia coli for purification and characterization. It irreversibly inhibits chymotrypsin-, trypsin-, elastase-, and subtilisin-like proteases in a temperature range from 20 to 100 °C in a concentration-dependent manner. The stoichiometry of inhibition (SI) of Pnserpin for proteases decreases as the temperature increases, indicating that the inhibitory activity of Pnserpin increases with the temperature. SDS-PAGE (sodium dodecyl sulfate polyacrylamide gel electrophoresis) showed that Pnserpin inhibits proteases by forming a SDS-resistant covalent complex. Homology modeling and molecular dynamic simulations predicted that Pnserpin can form a stable common serpin fold. Results of the present work will help in understanding the structural and functional characteristics of thermophilic serpin and will broaden the current knowledge about serpins from extremophiles.

Peculiarities of thermal denaturation of OmpF porin from Yersinia ruckeri

Irreversible denaturation of membrane proteins in detergent solutions is similar to unfolding of water-soluble multidomain proteins and represents a complex, multistage process.

Tuning protein mechanics through an ionic cluster graft from an extremophilic protein

Proteins from extremophilic organisms provide excellent model systems to determine the role of non-covalent interactions in defining protein stability and dynamics as well as being attractive targets for the development of robust biomaterials. Hyperthermophilic proteins have a prevalence of salt bridges, relative to their mesophilic homologues, which are thought to be important for enhanced thermal stability. However, the impact of salt bridges on the mechanical properties of proteins is far from understood. Here, a combination of protein engineering, biophysical characterisation, single molecule force spectroscopy (SMFS) and molecular dynamics (MD) simulations directly investigates the role of salt bridges in the mechanical stability of two cold shock proteins; BsCSP from the mesophilic organism Bacillus subtilis and TmCSP from the hyperthermophilic organism Thermotoga maritima. Single molecule force spectroscopy shows that at ambient temperatures TmCSP is mechanically stronger yet, counter-intuitively, its native state can withstand greater deformation before unfolding (i.e. it is mechanically soft) compared with BsCSP. MD simulations were used to identify the location and quantify the population of salt bridges, and reveal that TmCSP contains a larger number of highly occupied salt bridges than BsCSP. To test the hypothesis that salt-bridges endow these mechanical properties on the hyperthermophilic CSP, a charged triple mutant (CTM) variant of BsCSP was generated by grafting an ionic cluster from TmCSP into the BsCSP scaffold. As expected CTM is thermodynamically more stable and mechanically softer than BsCSP. We show that a grafted ionic cluster can increase the mechanical softness of a protein and speculate that it could provide a mechanical recovery mechanism and that it may be a design feature applicable to other proteins.

Toward rational thermostabilization of Aspergillus oryzae cutinase: Insights into catalytic and structural stability

Cutinases are powerful hydrolases that can cleave ester bonds of polyesters such as poly(ethylene terephthalate) (PET), opening up new options for enzymatic routes for polymer recycling and surface modification reactions. Cutinase from Aspergillus oryzae (AoC) is promising owing to the presence of an extended groove near the catalytic triad which is important for the orientation of polymeric chains. However, the catalytic efficiency of AoC on rigid polymers like PET is limited by its low thermostability; as it is essential to work at or over the glass transition temperature (Tg) of PET, that is, 70 °C. Consequently, in this study we worked toward the thermostabilization of AoC. Use of Rosetta computational protein design software in conjunction with rational design led to a 6 °C improvement in the thermal unfolding temperature (Tm) and a 10-fold increase in the half-life of the enzyme activity at 60 °C. Surprisingly, thermostabilization did not improve the rate or temperature optimum of enzyme activity. Three notable findings are presented as steps toward designing more thermophilic cutinase: (a) surface salt bridge optimization produced enthalpic stabilization, (b) mutations to proline reduced the entropy loss upon folding, and (c) the lack of a correlative increase in the temperature optimum of catalytic activity with thermodynamic stability suggests that the active site is locally denatured at a temperature below the Tm of the global structure.

Improving the Thermostability and Activity of a Thermophilic Subtilase by Incorporating Structural Elements of Its Psychrophilic Counterpart

The incorporation of the structural elements of thermostable enzymes into their less stable counterparts is generally used to improve enzyme thermostability. However, the process of engineering enzymes with both high thermostability and high activity remains an important challenge. Here, we report that the thermostability and activity of a thermophilic subtilase were simultaneously improved by incorporating structural elements of a psychrophilic subtilase. There were 64 variable regions/residues (VRs) in the alignment of the thermophilic WF146 protease, mesophilic sphericase, and psychrophilic S41. The WF146 protease was subjected to systematic mutagenesis, in which each of its VRs was replaced with those from S41 and sphericase. After successive rounds of combination and screening, we constructed the variant PBL5X with eight amino acid residues from S41. The half-life of PBL5X at 85°C (57.1 min) was approximately 9-fold longer than that of the wild-type (WT) WF146 protease (6.3 min). The substitutions also led to an increase in the apparent thermal denaturation midpoint temperature (Tm) of the enzyme by 5.5°C, as determined by differential scanning calorimetry. Compared to the WT, PBL5X exhibited high caseinolytic activity (25 to 95°C) and high values of Km and kcat (25 to 80°C). Our study may provide a rational basis for developing highly stable and active enzymes, which are highly desired in industrial applications.

Structural Rigidity and Protein Thermostability in Variants of Lipase A from Bacillus subtilis

Understanding the origin of thermostability is of fundamental importance in protein biochemistry. Opposing views on increased or decreased structural rigidity of the folded state have been put forward in this context. They have been related to differences in the temporal resolution of experiments and computations that probe atomic mobility. Here, we find a significant (p = 0.004) and fair (R² = 0.46) correlation between the structural rigidity of a well-characterized set of 16 mutants of lipase A from Bacillus subtilis (BsLipA) and their thermodynamic thermostability. We apply the rigidity theory-based Constraint Network Analysis (CNA) approach, analyzing directly and in a time-independent manner the statics of the BsLipA mutants. We carefully validate the CNA results on macroscopic and microscopic experimental observables and probe for their sensitivity with respect to input structures. Furthermore, we introduce a robust, local stability measure for predicting thermodynamic thermostability. Our results complement work that showed for pairs of homologous proteins that raising the structural stability is the most common way to obtain a higher thermostability. Furthermore, they demonstrate that related series of mutants with only a small number of mutations can be successfully analyzed by CNA, which suggests that CNA can be applied prospectively in rational protein design aimed at higher thermodynamic thermostability.

Discriminating between stabilizing and destabilizing protein design mutations via recombination and simulation

Accuracy of current computational protein design (CPD) methods is limited by inherent approximations in energy potentials and sampling. These limitations are often used to qualitatively explain design failures; however, relatively few studies provide specific examples or quantitative details that can be used to improve future CPD methods. Expanding the design method to include a library of sequences provides data that is well suited for discriminating between stabilizing and destabilizing design elements. Using thermophilic endoglucanase E1 from Acidothermus cellulolyticus as a model enzyme, we computationally designed a sequence with 60 mutations. The design sequence was rationally divided into structural blocks and recombined with the wild-type sequence. Resulting chimeras were assessed for activity and thermostability. Surprisingly, unlike previous chimera libraries, regression analysis based on one- and two-body effects was not sufficient for predicting chimera stability. Analysis of molecular dynamics simulations proved helpful in distinguishing stabilizing and destabilizing mutations. Reverting to the wild-type amino acid at destabilized sites partially regained design stability, and introducing predicted stabilizing mutations in wild-type E1 significantly enhanced thermostability. The ability to isolate stabilizing and destabilizing elements in computational design offers an opportunity to interpret previous design failures and improve future CPD methods.

Understanding the thermostability and activity of Bacillus subtilis lipase mutants: insights from molecular dynamics simulations

Improving the thermostability of industrial enzymes is an important protein engineering challenge. Point mutations, induced to increase thermostability, affect the structure and dynamics of the target protein in several ways and thus can also affect its activity. There appears to be no general rules for improving the thermostabilty of enzymes without adversely affecting their enzymatic activity. We report MD simulations, of wild type Bacillus subtilis lipase (WT) and its six progressively thermostable mutants (2M, 3M, 4M, 6M, 9M, and 12M), performed at different temperatures, to address this issue. Less thermostable mutants (LTMs), 2M to 6M, show WT-like dynamics at all simulation temperatures. However, the two more thermostable mutants (MTMs) show the required flexibility at appropriate temperature ranges and maintain conformational stability at high temperature. They show a deep and rugged free-energy landscape, confining them within a near-native conformational space by conserving noncovalent interactions, and thus protecting them from possible aggregation. In contrast, the LTMs having marginally higher thermostabilities than WT show greater probabilities of accessing non-native conformations, which, due to aggregation, have reduced possibilities of reverting to their respective native states under refolding conditions. Our analysis indicates the possibility of nonadditive effects of point mutations on the conformational stability of LTMs.

Are coarse-grained models apt to detect protein thermal stability? The case of OPEP force field

We present the first investigation of the kinetic and thermodynamic stability of two homologous thermophilic and mesophilic proteins based on the coarse-grained model OPEP. The object of our investigation is a pair of G-domains of relatively large size, 200 amino acids each, with an experimental stability gap of about 40 K. The OPEP force field is able to maintain stable the fold of these relatively large proteins within the hundrend-nanosecond time scale without including external constraints. This makes possible to characterize the conformational landscape of the folded protein as well as to explore the unfolding. In agreement with all-atom simulations used as a reference, we show that the conformational landscape of the thermophilic protein is characterized by a larger number of substates with slower dynamics on the network of states and more resilient to temperature increase. Moreover, we verify the stability gap between the two proteins using replica-exchange simulations and estimate a difference between the melting temperatures of about 23 K, in fair agreement with experiment. The detailed investigation of the unfolding thermodynamics, allows to gain insight into the mechanism underlying the enhanced stability of the thermophile relating it to a smaller heat capacity of unfolding.

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.

Supercritical carbon dioxide and its potential as a life-sustaining solvent in a planetary environment

Supercritical fluids have different properties compared to regular fluids and could play a role as life-sustaining solvents on other worlds. Even on Earth, some bacterial species have been shown to be tolerant to supercritical fluids. The special properties of supercritical fluids, which include various types of selectivities (e.g., stereo-, regio-, and chemo-selectivity) have recently been recognized in biotechnology and used to catalyze reactions that do not occur in water. One suitable example is enzymes when they are exposed to supercritical fluids such as supercritical carbon dioxide: enzymes become even more stable, because they are conformationally rigid in the dehydrated state. Furthermore, enzymes in anhydrous organic solvents exhibit a "molecular memory", i.e., the capacity to "remember" a conformational or pH state from being exposed to a previous solvent. Planetary environments with supercritical fluids, particularly supercritical carbon dioxide, exist, even on Earth (below the ocean floor), on Venus, and likely on Super-Earth type exoplanets. These planetary environments may present a possible habitat for exotic life.

How conformational flexibility stabilizes the hyperthermophilic elongation factor G-domain

Proteins from thermophilic organisms are stable and functional well above ambient temperature. Understanding the molecular mechanism underlying such a resistance is of crucial interest for many technological applications. For some time, thermal stability has been assumed to correlate with high mechanical rigidity of the protein matrix. In this work we address this common belief by carefully studying a pair of homologous G-domain proteins, with their melting temperatures differing by 40 K. To probe the thermal-stability content of the two proteins we use extensive simulations covering the microsecond time range and employ several different indicators to assess the salient features of the conformational landscape and the role of internal fluctuations at ambient condition. At the atomistic level, while the magnitude of fluctuations is comparable, the distribution of flexible and rigid stretches of amino-acids is more regular in the thermophilic protein causing a cage-like correlation of amplitudes along the sequence. This caging effect is suggested to favor stability at high T by confining the mechanical excitations. Moreover, it is found that the thermophilic protein, when folded, visits a higher number of conformational substates than the mesophilic homologue. The entropy associated with the occupation of the different substates and the thermal resilience of the protein intrinsic compressibility provide a qualitative insight on the thermal stability of the thermophilic protein as compared to its mesophilic homologue. Our findings potentially open the route to new strategies in the design of thermostable proteins.

A comparative study of ribosomal proteins: linkage between amino acid distribution and ribosomal assembly

BACKGROUND

Assembly of the ribosome from its protein and RNA constituents must occur quickly and efficiently in order to synthesize the proteins necessary for all cellular activity. Since the early 1960's, certain characteristics of possible assembly pathways have been elucidated, yet the mechanisms that govern the precise recognition events remain unclear.We utilize a comparative analysis to investigate the amino acid composition of ribosomal proteins (r-proteins) with respect to their role in the assembly process. We compared small subunit (30S) r-protein sequences to those of other housekeeping proteins from 560 bacterial species and searched for correlations between r-protein amino acid content and factors such as assembly binding order, environmental growth temperature, protein size, and contact with ribosomal RNA (rRNA) in the 30S complex.

RESULTS

We find r-proteins have a significantly high percent of positive residues, which are highly represented at rRNA contact sites. An inverse correlation between the percent of positive residues and r-protein size was identified and is mainly due to the content of Lysine residues, rather than Arginine. Nearly all r-proteins carry a net positive charge, but no statistical correlation between the net charge and the binding order was detected. Thermophilic (high-temperature) r-proteins contain increased Arginine, Isoleucine, and Tyrosine, and decreased Serine and Threonine compared to mesophilic (lower-temperature), reflecting a known distinction between thermophiles and mesophiles, possibly to account for protein thermostability. However, this difference in amino acid content does not extend to rRNA contact sites, as the proportions of thermophilic and mesophilic contact residues are not significantly different.

CONCLUSIONS

Given the significantly higher level of positively charged residues in r-proteins and at contact sites, we conclude that ribosome assembly relies heavily on an electrostatic component of interaction. However, the binding order of r-proteins in assembly does not appear to depend on these electrostatics interactions. Additionally, because thermophiles and mesophiles exhibit significantly different amino acid compositions in their sequences but not in the identities of contact sites, we conclude that this electrostatic component of interaction is insensitive to temperature and is not the determining factor differentiating the temperature sensitivity of ribosome assembly.

Hyaluronic acid nanoparticles titrate the viscoelastic properties of viscosupplements

Hyaluronic acid (HA) is a glycosaminoglycan with diverse biomedical applications including viscosupplementation of synovial fluid for the treatment osteoarthritis. Current HA viscosupplements such as Synvisc, Orthovisc, and Hyalgan have shown positive effects of reducing pain and improving joint function. The therapeutic efficacy, however, is highly transient, and these viscous fluids suffer from poor injectability. HA nanoparticles were found to modify the rheological properties of a model of the HA viscosupplement Orthovisc. Nanoparticles were successfully synthesized from 17 and 1500 kDa HA. Nanoparticle suspensions of HA were studied at different concentrations and in blends with the model viscosupplement. Nanoparticles made from 1500 kDa HA reduced the viscosupplement viscosity and elasticity to a much greater degree than nanoparticles made from 17 kDa HA. The difference in the nanoparticle effect on viscoelasticity suggested that nanoparticles made from 17 kDa HA may have dangling surface polymers that facilitated interactions with HA in solution. This hypothesis was supported by the greater compressibility of 17 kDa nanoparticles as determined by ultrasonic vibrational spectroscopy. Rheological investigations showed that the viscoelasticity of viscosupplements could be discretely titrated by modulating the concentration and type of HA nanoparticle additive (hard sphere or hairy). Thus, the injectability of viscosupplements may be enhanced while maintaining high elasticity.

An in silico method for designing thermostable variant of a dimeric mesophilic protein based on its 3D structure

Designing proteins with enhanced thermostability has been a major interest of protein engineering because of its potential industrial applications. Here, we have presented a computational method for designing dimeric thermostable protein based on rational mutations on a mesophilic protein. Experimental and structural data indicate that the surface stability of a protein is a major factor controlling denaturation of a protein and ion-pairs are most efficient in enhancing the stability of the surfaces of the monomers and the interface between them. Our mutation based strategy is to first identify several polar or charged residues on the protein surface, interacting weakly with the rest of the protein and then replacing the side-chains of suitable neighboring residues to increase the interaction between these two residues. In stabilizing the interface, mutation is done in the interface for forming an ion pairs between the monomers. Application of this design strategy to a homo-dimeric protein and a hetero-dimeric protein as examples has produced excellent results. In both the cases the designed mutated proteins including the individual monomers and the interfaces were found to be considerably more stable than the respective mesophilic proteins as judged by self-energies and residue-wise interaction patterns. This method is easily applicable to any multi-meric proteins.

Comparative analysis to identify determinants of changing life style in Thermosynechococcus elongatus BP-1, a thermophilic cyanobacterium

A comparative genomics analysis among all forty whole genome sequences available for cyanobacteria (3 thermophiles- Thermosynechococcus elongatus BP-1, Synechococcus sp. JA-2-3B'a (2-13), Synechococcus sp. JA-3-3Ab and 37 mesophiles) was performed to identify genomic and proteomic factors responsible for the behaviour of T. elongatus BP-1, a thermophilic unicellular cyanobacterium with optimum growth temperature [OGT] of 55°C. Majority of genomic and proteomic characteristics for this cyanobacterium indicated contrasting features indicating its mesophilic behaviour while the role of mutational biasness and selection pressure is thought to be responsible for high OGT. Contradictory results were obtained for T. elongatus for synonymous codon usage, CvP-bias and amino acid composition with respect to thermophilic behaviour. Calculated J2 index is lowest among all cyanobacterial genomes. Except for proline and termination codons, T. elongatus showed synonymous codon usage pattern which is expected for mesophiles. Results indicated that among cyanobacterial genomes, majority of genomic and proteomic determinants put T. elongatus very close to mesophiles and the whole genome of this organism represents continuous gain of mesophilic rather than thermophilic behavior.

Hyperstability and substrate promiscuity in laboratory resurrections of Precambrian β-lactamases

We report a sequence reconstruction analysis targeting several Precambrian nodes in the evolution of class-A β-lactamases and the preparation and experimental characterization of their encoded proteins. Despite extensive sequence differences with the modern enzymes (~100 amino acid differences), the proteins resurrected in the laboratory properly fold into the canonical lactamase structure. The encoded proteins from 2-3 billion years (Gyr)-old β-lactamase sequences undergo cooperative two-state thermal denaturation and display very large denaturation temperature enhancements (~35 °C) relative to modern β-lactamases. They degrade different antibiotics in vitro with catalytic efficiencies comparable to that of an average modern enzyme. This enhanced substrate promiscuity is not accompanied by significant changes in the active-site region as seen in static X-ray structures, suggesting a plausible role for dynamics in the evolution of function in these proteins. Laboratory resurrections of 2-3 Gyr-old β-lactamases also endowed modern microorganisms with significant levels of resistance toward a variety of antibiotics, opening up the possibility of performing laboratory replays of the molecular tape of lactamase evolution. Overall, these results support the notions that Precambrian life was thermophilic and that proteins can evolve from substrate-promiscuous generalists into specialists during the course of natural evolution. They also highlight the biotechnological potential of laboratory resurrection of Precambrian proteins, as both high stability and enhanced promiscuity (likely contributors to high evolvability) are advantageous features in protein scaffolds for molecular design and laboratory evolution.

The effect of protein composition on hydration dynamics

Water dynamics at the surface of two homologous proteins with different thermal resistances is found to be unaffected by the different underlying amino-acid compositions, and when proteins are folded it responds similarly to temperature variations. Upon unfolding the water dynamics slowdown with respect to bulk decreases by a factor of two. Our findings are explained by the dominant topological perturbation induced by the protein on the water hydrogen bond dynamics.

Without salt, the 'thermophilic' protein Mth10b is just mesophilic

Most proteins from thermophiles or hyperthermophiles are intrinsically thermostable. However, though Methanobacterium thermoautotrophicum ΔH is a thermophilic archaeon with an optimal growth temperature of 65 °C, Mth10b, an atypical member the Sac10b protein family from M. thermoautotrophicum ΔH, seems not intrinsically thermostable. In this work, to clarify the molecular mechanism of Mth10b remaining stable under its physiological conditions, the thermodynamic properties of Mth10b were studied through equilibrium unfolding experiments performed at pH 7.0 monitored by circular dichroism (CD) spectra in detail. Our work demonstrated that Mth10b is not intrinsically thermostable and that due to the masking effect upon the large numbers of destabilizing electrostatic repulsions resulting from the extremely uneven distribution of charged residues over the surface of Mth10b, salt can contribute to the thermostability of Mth10b greatly. Considering that the intracellular salt concentration is high to 0.7 M, we concluded that salt is the key extrinsic factor to Mth10b remaining stable under its physiological conditions. In other word, without salt, 'thermophilic' protein Mth10b is just a mesophilic one.

Excipients differentially influence the conformational stability and pretransition dynamics of two IgG1 monoclonal antibodies

Since immunoglobulins are conformationally dynamic molecules in solution, we studied the effect of stabilizing and destabilizing excipients on the conformational stability and dynamics of two IgG1 monoclonal antibodies (mAbs; mAb-A and mAb-B) using a variety of biophysical approaches. Even though the two mAbs are of the same IgG1 subtype, the unfolding patterns, aggregation behavior, and pretransition dynamics of these two antibodies were strikingly different in response to external perturbations such as pH, temperature, and presence of excipients. Sucrose and arginine were identified as stabilizers and destabilizers, respectively, on the basis of their influence on conformational stability for both the IgG1 mAbs. The two excipients, however, had distinct effective concentrations and different effects on the conformational stability and pretransition dynamics of the two mAbs as measured by a combination of differential scanning calorimetry, high-resolution ultrasonic spectroscopy, and red-edge excitation shift fluorescence studies. Stabilizing concentrations of sucrose were found to decrease the internal motions of mAb-B, whereas arginine marginally increased its adiabatic compressibility in the pretransition region. Both sucrose and arginine did not influence the pretransition dynamics of mAb-A. The potential reasons for such differences in excipient effects between two IgG1 mAbs are discussed.

Functional validation of hydrophobic adaptation to physiological temperature in the small heat shock protein αA-crystallin

Small heat shock proteins (sHsps) maintain cellular homeostasis by preventing stress and disease-induced protein aggregation. While it is known that hydrophobicity impacts the ability of sHsps to bind aggregation-prone denaturing proteins, the complex quaternary structure of globular sHsps has made understanding the significance of specific changes in hydrophobicity difficult. Here we used recombinant protein of the lenticular sHsp α A-crystallin from six teleost fishes environmentally adapted to temperatures ranging from -2°C to 40°C to identify correlations between physiological temperature, protein stability and chaperone-like activity. Using sequence and structural modeling analysis we identified specific amino acid differences between the warm adapted zebrafish and cold adapted Antarctic toothfish that could contribute to these correlations and validated the functional consequences of three specific hydrophobicity-altering amino acid substitutions in αA-crystallin. Site directed mutagenesis of three residues in the zebrafish (V62T, C143S, T147V) confirmed that each impacts either protein stability or chaperone-like activity or both, with the V62T substitution having the greatest impact. Our results indicate a role for changing hydrophobicity in the thermal adaptation of α A-crystallin and suggest ways to produce sHsp variants with altered chaperone-like activity. These data also demonstrate that a comparative approach can provide new information about sHsp function and evolution.

Temperature-induced transitions in disordered proteins probed by NMR spectroscopy

Intrinsically disordered proteins are abundant in nature and perform many important physiological functions. Multidimensional NMR spectroscopy has been crucial for the understanding of the conformational properties of disordered proteins and is increasingly used to probe their conformational ensembles. Compared to folded proteins, disordered proteins are more malleable and more easily perturbed by environmental factors. Accordingly, the experimental conditions and especially the temperature modify the structural and functional properties of disordered proteins. NMR spectroscopy allows analysis of temperature-induced structural changes at residue resolution using secondary chemical shift analysis, paramagnetic relaxation enhancement, and residual dipolar couplings. This chapter discusses practical aspects of NMR studies of temperature-induced structural changes in disordered proteins.

Relationship between stability and flexibility in the most flexible region of Photinus pyralis luciferase

Firefly luciferase is a protein with a large N-terminal and a small C-terminal domain. B-factor analysis shows that its C-terminal is much more flexible than its N-terminal. Studies on hyperthermophile proteins have been shown that the increased thermal stability of hyperthermophile proteins is due to their enhanced conformational rigidity and the relationship between flexibility, stability and function in most of proteins is on debate. Two mutations (D474K and D476N) in the most flexible region of firefly luciferase were designed. Thermostability analysis shows that D476N mutation doesn't have any significant effect but D474K mutation destabilized protein. On the other hand, flexibility analysis using dynamic quenching and limited proteolysis demonstrates that D474K mutation became much more flexible than wild type although D476N doesn't have any significant difference. Intrinsic and ANS fluorescence studies demonstrate that D476N mutation is brought about by structural changes without significant effect on thermostability and flexibility. Molecular modeling reveals that disruption of a salt bridge between D(474) and K(445) accompanying with some H-bond deletion may be involved in destabilization of D474K mutant.

Thermophilic proteins: insight and perspective from in silico experiments

Proteins from thermophilic and hyperthermophilic organisms are stable and function at high temperatures (50-100 °C). The importance of understanding the microscopic mechanisms underlying this thermal resistance is twofold: it is key for acquiring general clues on how proteins maintain their fold stable and for targeting those medical and industrial applications that aim at designing enzymes that can work under harsh conditions. In this tutorial review we first provide the general background of protein thermostability by specifically focusing on the structural and thermodynamic peculiarities; next, we discuss how computational studies based on Molecular Dynamics simulations can broaden and refine our knowledge on such special class of proteins.

How do thermophilic proteins and proteomes withstand high temperature?

We attempt to understand the origin of enhanced stability in thermophilic proteins by analyzing thermodynamic data for 116 proteins, the largest data set achieved to date. We compute changes in entropy and enthalpy at the convergence temperature where different driving forces are maximally decoupled, in contrast to the majority of previous studies that were performed at the melting temperature. We find, on average, that the gain in enthalpy upon folding is lower in thermophiles than in mesophiles, whereas the loss in entropy upon folding is higher in mesophiles than in thermophiles. This implies that entropic stabilization may be responsible for the high melting temperature, and hints at residual structure or compactness of the denatured state in thermophiles. We find a similar trend by analyzing a homologous set of proteins classified based only on the optimum growth temperature of the organisms from which they were extracted. We find that the folding free energy at the temperature of maximal stability is significantly more favorable in thermophiles than in mesophiles, whereas the maximal stability temperature itself is similar between these two classes. Furthermore, we extend the thermodynamic analysis to model the entire proteome. The results explain the high optimal growth temperature in thermophilic organisms and are in excellent quantitative agreement with full thermal growth rate data obtained in a dozen thermophilic and mesophilic organisms.

Temperature dependence of fast carbonyl backbone dynamics in chicken villin headpiece subdomain

Temperature-dependence of protein dynamics can provide information on details of the free energy landscape by probing the characteristics of the potential responsible for the fluctuations. We have investigated the temperature-dependence of picosecond to nanosecond backbone dynamics at carbonyl carbon sites in chicken villin headpiece subdomain protein using a combination of three NMR relaxation rates: (13)C' longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, (13)C'/(13)C'-(13)C(α) CSA/dipolar and (13)C'/(13)C'-(15)N CSA/dipolar. Order parameters have been extracted using the Lipari-Szabo model-free approach assuming a separation of the time scales of internal and molecular motions in the 2-16°C temperature range. There is a gradual deviation from this assumption from lower to higher temperatures, such that above 16°C the separation of the time scales is inconsistent with the experimental data and, thus, the Lipari-Szabo formalism can not be applied. While there are variations among the residues, on the average the order parameters indicate a markedly steeper temperature dependence at backbone carbonyl carbons compared to that probed at amide nitrogens in an earlier study. This strongly advocates for probing sites other than amide nitrogen for accurate characterization of the potential and other thermodynamics characteristics of protein backbone.