Thermophilic proteins: insight and perspective from in silico experiments

Temperature-dependent iron motion in extremophile rubredoxins - no need for 'corresponding states'

Extremophile organisms are known that can metabolize at temperatures down to - 25 °C (psychrophiles) and up to 122 °C (hyperthermophiles). Understanding viability under extreme conditions is relevant for human health, biotechnological applications, and our search for life elsewhere in the universe. Information about the stability and dynamics of proteins under environmental extremes is an important factor in this regard. Here we compare the dynamics of small Fe-S proteins - rubredoxins - from psychrophilic and hyperthermophilic microorganisms, using three different nuclear techniques as well as molecular dynamics calculations to quantify motion at the Fe site. The theory of 'corresponding states' posits that homologous proteins from different extremophiles have comparable flexibilities at the optimum growth temperatures of their respective organisms. Although 'corresponding states' would predict greater flexibility for rubredoxins that operate at low temperatures, we find that from 4 to 300 K, the dynamics of the Fe sites in these homologous proteins are essentially equivalent.

Decoding the Role of the Global Proteome Dynamics for Cellular Thermal Stability

Molecular mechanisms underlying the thermal response of cells remain elusive. On the basis of the recent result that the short-time diffusive dynamics of the Escherichia coli proteome is an excellent indicator of temperature-dependent bacterial metabolism and death, we used neutron scattering (NS) spectroscopy and molecular dynamics (MD) simulations to investigate the sub-nanosecond proteome mobility in psychro-, meso-, and hyperthermophilic bacteria over a wide temperature range. The magnitude of thermal fluctuations, measured by atomic mean square displacements, is similar among all studied bacteria at their respective thermal cell death. Global roto-translational motions turn out to be the main factor distinguishing the bacterial dynamical properties. We ascribe this behavior to the difference in the average proteome net charge, which becomes less negative for increasing bacterial thermal stability. We propose that the chemical-physical properties of the cytoplasm and the global dynamics of the resulting proteome are fine-tuned by evolution to uphold optimal thermal stability conditions.

Simulation and practical investigation of carbonic anhydrase stability in an industrial solvent system of methyl diethanolamine for carbon dioxide capture

Carbonic anhydrase owing to its potential as an industrial biocatalyst for carbon dioxide sequestration from flue gas has attracted considerable attention in solving global warming problems. A large body of research has been conducted to increase the thermal stability of carbonic anhydrase from different sources against the harsh operational conditions of CO2 capture systems. In contrast to cost-intensive protein engineering methods, solvation with aqueous-organic binary mixtures offers a convenient and economical alternative strategy for retention of protein structure and stability. This study aimed to examine the stabilizing effect of methyl diethanolamine (MDEA) as a component of an aqueous-organic solvent mixture on human carbonic anhydrase II (HCA II) at extreme temperatures. Computational and also spectroscopic examinations were employed for tracking conformational changes and stability evaluation of HCA II in 50:50 (vol %) water: MDEA binary mixture at high temperature. Molecular dynamic (MD) simulation studies predicted the high thermal stability of HCA II in the presence of MDEA. UV absorbance spectra confirmed the thermo-stabilizing effect of the binary solvent mixture on HCA II. While the enzymatic activity of HCA II at 25 °C in the presence of 10, 25, and 50 (vol%) of MDEA was substantially increased, no obvious effect on retention of HCA II activity in the water-MDEA binary solvent mixture at 85 °C was seen. It is shown that the solvation of HCA II in the presence of MDEA could result in the prevention of aggregate formation in high temperatures but not functional stability.Communicated by Ramaswamy H. Sarma.

Biochemical, structural and dynamical characterizations of the lactate dehydrogenase from Selenomonas ruminantium provide information about an intermediate evolutionary step prior to complete allosteric regulation acquisition in the super family of lactate and malate dehydrogenases

In this work, we investigated the lactate dehydrogenase (LDH) from Selenomonas ruminantium (S. rum), an enzyme that differs at key amino acid positions from canonical allosteric LDHs. The wild type (Wt) of this enzyme recognises pyuvate as all LDHs. However, introducing a single point mutation in the active site loop (I85R) allows S. Rum LDH to recognize the oxaloacetate substrate as a typical malate dehydrogenase (MalDH), whilst maintaining homotropic activation as an LDH. We report the tertiary structure of the Wt and I85RLDH mutant. The Wt S. rum enzyme structure binds NADH and malonate, whilst also resembling the typical compact R-active state of canonical LDHs. The structure of the mutant with I85R was solved in the Apo State (without ligand), and shows no large conformational reorganization such as that observed with canonical allosteric LDHs in Apo state. This is due to a local structural feature typical of S. rum LDH that prevents large-scale conformational reorganization. The S. rum LDH was also studied using Molecular Dynamics simulations, probing specific local deformations of the active site that allow the S. rum LDH to sample the T-inactive state. We propose that, with respect to the LDH/MalDH superfamily, the S. rum enzyme possesses a specificstructural and dynamical way to ensure homotropic activation.

In silico characterization of differentially expressed short-read nucleotide sequences identified in dieback stress-induced transcriptomic analysis reveals their role as antimicrobial peptides

We investigated the in silico characterization of short-length nucleotide sequences that were differentially expressed in dieback stress-induced transcriptomic analysis. They displayed homology with C-terminal flanking peptides and defensins-like proteins, revealing their antimicrobial activity. Their predicted fingerprints displayed protein signatures related to antimicrobial peptides. These short-length RGAs have been shown to possess structural motifs such as APLT P-type ATPase, casein kinase II (CK2), protein kinase 3, protein kinase C (PKC), and N-glycosylation site that are the attributes of disease resistance genes. The prediction of arginine and lysine residues in active binding sites in ligand docking analysis prophesied them as antimicrobial peptides due to their strong relation with antimicrobial activity. The in silico structural-functional characterization has predicted their role in resistance against microbial pathogens. Moreover, the predicted antimicrobial peptide regions showed their homology with the signature domain of PR-5-like protein and AMP family Thaumatin.

Structure-Function Relationships of Oxygen Transport Proteins in Marine Invertebrates Enduring Higher Temperatures and Deoxygenation

AbstractPredictions for climate change-to lesser and greater extents-reveal a common scenario in which marine waters are characterized by a deadly trio of stressors: higher temperatures, lower oxygen levels, and acidification. Ectothermic taxa that inhabit coastal waters, such as shellfish, are vulnerable to rapid and prolonged environmental disturbances, such as heatwaves, pollution-induced eutrophication, and dysoxia. Oxygen transport capacity of the hemolymph (blood equivalent) is considered the proximal driver of thermotolerance and respiration in many invertebrates. Moreover, maintaining homeostasis under environmental duress is inextricably linked to the activities of the hemolymph-based oxygen transport or binding proteins. Several protein groups fulfill this role in marine invertebrates: copper-based extracellular hemocyanins, iron-based intracellular hemoglobins and hemerythrins, and giant extracellular hemoglobins. In this brief text, we revisit the distribution and multifunctional properties of oxygen transport proteins, notably hemocyanins, in the context of climate change, and the consequent physiological reprogramming of marine invertebrates.

Life at high temperature observed in vitro upon laser heating of gold nanoparticles

Thermophiles are microorganisms that thrive at high temperature. Studying them can provide valuable information on how life has adapted to extreme conditions. However, high temperature conditions are difficult to achieve on conventional optical microscopes. Some home-made solutions have been proposed, all based on local resistive electric heating, but no simple commercial solution exists. In this article, we introduce the concept of microscale laser heating over the field of view of a microscope to achieve high temperature for the study of thermophiles, while maintaining the user environment in soft conditions. Microscale heating with moderate laser intensities is achieved using a substrate covered with gold nanoparticles, as biocompatible, efficient light absorbers. The influences of possible microscale fluid convection, cell confinement and centrifugal thermophoretic motion are discussed. The method is demonstrated with two species: (i) Geobacillus stearothermophilus, a motile thermophilic bacterium thriving around 65 °C, which we observed to germinate, grow and swim upon microscale heating and (ii) Sulfolobus shibatae, a hyperthermophilic archaeon living at the optimal temperature of 80 °C. This work opens the path toward simple and safe observation of thermophilic microorganisms using current and accessible microscopy tools.

Studying microorganisms at high temperatures is challenging on conventional optical microscopes. Here, the authors introduce the concept of microscale laser heating over the full field of view by using gold nanoparticles as light absorbers, and study thermophile species up to 80 °C.

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

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

Loosely-packed dynamical structures with partially-melted surface being the key for thermophilic argonaute proteins achieving high DNA-cleavage activity

Prokaryotic Argonaute proteins (pAgos) widely participate in hosts to defend against the invasion of nucleic acids. Compared with the CRISPR-Cas system, which requires a specific motif on the target and can only use RNA as guide, pAgos exhibit precise endonuclease activity on any arbitrary target sequence and can use both RNA and DNA as guide, thus rendering great potential for genome editing applications. Hitherto, most in-depth studies on the structure-function relationship of pAgos were conducted on thermophilic ones, functioning at ∼60 to 100°C, whose structures were, however, determined experimentally at much lower temperatures (20-33°C). It remains unclear whether these low-temperature structures can represent the true conformations of the thermophilic pAgos under their physiological conditions. The present work studied three pAgos, PfAgo, TtAgo and CbAgo, whose physiological temperatures differ significantly (95, 75 and 37°C). By conducting thorough experimental and simulation studies, we found that thermophilic pAgos (PfAgo and TtAgo) adopt a loosely-packed structure with a partially-melted surface at the physiological temperatures, largely different from the compact crystalline structures determined at moderate temperatures. In contrast, the mesophilic pAgo (CbAgo) assumes a compact crystalline structure at its optimal function temperature. Such a partially-disrupted structure endows thermophilic pAgos with great flexibility both globally and locally at the catalytic sites, which is crucial for them to achieve high DNA-cleavage activity. To further prove this, we incubated thermophilic pAgos with urea to purposely disrupt their structures, and the resulting cleavage activity was significantly enhanced below the physiological temperature, even at human body temperature. Further testing of many thermophilic Agos present in various thermophilic prokaryotes demonstrated that their structures are generally disrupted under physiological conditions. Therefore, our findings suggest that the highly dynamical structure with a partially-melted surface, distinct from the low-temperature crystalline structure, could be a general strategy assumed by thermophilic pAgos to achieve the high DNA-cleavage activity.

Probing Adaptation of Hydration and Protein Dynamics to Temperature

Protein dynamics is strongly influenced by the surrounding environment and physiological conditions. Here we employ broadband megahertz-to-terahertz spectroscopy to explore the dynamics of water and myoglobin protein on an extended time scale from femto- to nanosecond. The dielectric spectra reveal several relaxations corresponding to the orientational polarization mechanism, including the dynamics of loosely bound, tightly bound, and bulk water, as well as collective vibrational modes of protein in an aqueous environment. The dynamics of loosely bound and bulk water follow non-Arrhenius behavior; however, the dynamics of water molecules in the tightly bound layer obeys the Arrhenius-type relation. Combining molecular simulations and effective-medium approximation, we have determined the number of water molecules in the tightly bound hydration layer and studied the dynamics of protein as a function of temperature. The results provide the important impact of water on the biochemical functions of proteins.

Isolation and thermo-acclimation of thermophilic bacteria in hyperthermophilic fermentation system

Hyperthermophilic microorganisms play a key role in the hyper-thermophilic composting (HTC) technique. However, little information is available about the hyperthermophilic microorganisms prevalent in HTC systems, except for the Calditerricola satsumensis, Calditerricola yamamurae, and Thermaerobacter. To obtain effective hyper-thermophilic microorganisms, a continuous thermo-acclimation of the suitable thermophilic microorganisms was demonstrated in this study. Bacillus thermoamylovorans with high-temperature endurance (70 °C) were newly isolated from sludge composting, and an adequate slow heating rate (2 °C per cycle) was applied to further improve its thermostability. Finally, a strain with a maximum growth temperature of 80 °C was obtained. Moreover, structural and hydrophobic changes in cell proteins, the special amino acid content ratio, and the membrane permeability of the thermophilic bacterium after thermo-acclimation were evaluated for improved thermostability. In addition, the acclimated hyperthermophilic bacterium was further inoculated into the HTC system, and an excellent performance with a maximum operating temperature of 82 °C was observed.

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

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

Structural basis for the hyperthermostability of an archaeal enzyme induced by succinimide formation

Stability of proteins from hyperthermophiles (organisms existing under boiling water conditions) enabled by a reduction of conformational flexibility is realized through various mechanisms. A succinimide (SNN) arising from the post-translational cyclization of the side chains of aspartyl/asparaginyl residues with the backbone amide -NH of the succeeding residue would restrain the torsion angle Ψ and can serve as a new route for hyperthermostability. However, such a succinimide is typically prone to hydrolysis, transforming to either an aspartyl or β-isoaspartyl residue. Here, we present the crystal structure of Methanocaldococcus jannaschii glutamine amidotransferase and, using enhanced sampling molecular dynamics simulations, address the mechanism of its increased thermostability, up to 100°C, imparted by an unexpectedly stable succinimidyl residue at position 109. The stability of SNN109 to hydrolysis is seen to arise from its electrostatic shielding by the side-chain carboxylate group of its succeeding residue Asp110, as well as through n → π^∗ interactions between SNN109 and its preceding residue Glu108, both of which prevent water access to SNN. The stable succinimidyl residue induces the formation of an α-turn structure involving 13-atom hydrogen bonding, which locks the local conformation, reducing protein flexibility. The destabilization of the protein upon replacement of SNN with a Φ-restricted prolyl residue highlights the specificity of the succinimidyl residue in imparting hyperthermostability to the enzyme. The conservation of the succinimide-forming tripeptide sequence (E(N/D)(E/D)) in several archaeal GATases strongly suggests an adaptation of this otherwise detrimental post-translational modification as a harbinger of thermostability.

Biochemical, structural and dynamical studies reveal strong differences in the thermal-dependent allosteric behavior of two extremophilic lactate dehydrogenases

In this work, we combined biochemical and structural investigations with molecular dynamics (MD) simulations to analyze the very different thermal-dependent allosteric behavior of two lactate dehydrogenases (LDH) from thermophilic bacteria. We found that the enzyme from Petrotoga mobilis (P. mob) necessitates an absolute requirement of the allosteric effector (fructose 1, 6-bisphosphate) to ensure functionality. In contrast, even without allosteric effector, the LDH from Thermus thermophilus (T. the) is functional when the temperature is raised. We report the crystal structure of P. mob LDH in the Apo state solved at 1.9 Å resolution. We used this structure and the one from T. the, obtained previously, as a starting point for MD simulations at various temperatures. We found clear differences between the thermal dynamics, which accounts for the behavior of the two enzymes. Our work demonstrates that, within an allosteric enzyme, some areas act as local gatekeepers of signal transmission, allowing the enzyme to populate either the T-inactive or the R-active states with different degrees of stringency.

Deletion of a Peptidylprolyl Isomerase Gene Results in the Inability of Caldicellulosiruptor bescii To Grow on Crystalline Cellulose without Affecting Protein Glycosylation or Growth on Soluble Substrates

Caldicellulosiruptor bescii secretes a large number of complementary multifunctional enzymes with unique activities for biomass deconstruction. The most abundant enzymes in the C. bescii secretome are found in a unique gene cluster containing a glycosyl transferase (GT39) and a putative peptidyl prolyl cis-trans isomerase. Deletion of the glycosyl transferase in this cluster resulted in loss of detectable protein glycosylation in C. bescii, and its activity has been shown to be responsible for the glycosylation of the proline-threonine rich linkers found in many of the multifunctional cellulases. The presence of a putative peptidyl prolyl cis-trans isomerase within this gene cluster suggested that it might also play a role in cellulase modification. Here, we identify this gene as a putative prsA prolyl cis-trans isomerase. Deletion of prsA2 leads to the inability of C. bescii to grow on insoluble substrates such as Avicel, the model cellulose substrate, while exhibiting no differences in phenotype with the wild-type strain on soluble substrates. Finally, we provide evidence that the prsA2 gene is likely needed to increase solubility of multifunctional cellulases and that this unique gene cluster was likely acquired by members of the Caldicellulosiruptor genus with a group of genes to optimize the production and activity of multifunctional cellulases.IMPORTANCE Caldicellulosiruptor has the ability to digest complex plant biomass without pretreatment and have been engineered to convert biomass, a sustainable, carbon neutral substrate, to fuels. Their strategy for deconstructing plant cell walls relies on an interesting class of cellulases consisting of multiple catalytic modules connected by linker regions and carbohydrate binding modules. The best studied of these enzymes, CelA, has a unique deconstruction mechanism. CelA is located in a cluster of genes that likely allows for optimal expression, secretion, and activity. One of the genes in this cluster is a putative isomerase that modifies the CelA protein. In higher eukaryotes, these isomerases are essential for the proper folding of glycoproteins in the endoplasmic reticulum, but little is known about the role of isomerization in cellulase activity. We show that the stability and activity of CelA is dependent on the activity of this isomerase.

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.

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.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes

The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

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

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

Crystal Structure and Conformational Dynamics of Pyrococcus furiosus Prolyl Oligopeptidase

Enzymes in the prolyl oligopeptidase family possess unique structures and substrate specificities that are important for their biological activity and for potential biocatalytic applications. The crystal structures of Pyrococcus furiosus ( Pfu) prolyl oligopeptidase (POP) and the corresponding S477C mutant were determined to 1.9 and 2.2 Å resolution, respectively. The wild type enzyme crystallized in an open conformation, indicating that this state is readily accessible, and it contained bound chloride ions and a prolylproline ligand. These structures were used as starting points for molecular dynamics simulations of Pfu POP conformational dynamics. The simulations showed that large-scale domain opening and closing occurred spontaneously, providing facile substrate access to the active site. Movement of the loop containing the catalytically essential histidine into a conformation similar to those found in structures with fully formed catalytic triads also occurred. This movement was modulated by chloride binding, providing a rationale for experimentally observed activation of POP peptidase catalysis by chloride. Thus, the structures and simulations reported in this study, combined with existing biochemical data, provide a number of insights into POP catalysis.

CAP modifies the structure of a model protein from thermophilic bacteria: mechanisms of CAP-mediated inactivation

Cold atmospheric plasma (CAP) has great potential for sterilization in the food industry, by deactivation of thermophilic bacteria, but the underlying mechanisms are largely unknown. Therefore, we investigate here whether CAP is able to denature/modify protein from thermophilic bacteria. We focus on MTH1880 (MTH) from Methanobacterium thermoautotrophicum as model protein, which we treated with dielectric barrier discharge (DBD) plasma operating in air for 10, 15 and 20 mins. We analysed the structural changes of MTH using circular dichroism, fluorescence and NMR spectroscopy, as well as the thermal and chemical denaturation, upon CAP treatment. Additionally, we performed molecular dynamics (MD) simulations to determine the stability, flexibility and solvent accessible surface area (SASA) of both the native and oxidised protein.

Ion-ion interactions in the denatured state contribute to the stabilization of CutA1 proteins

In order to elucidate features of the denatured state ensembles that exist in equilibrium with the native state under physiological conditions, we performed 1.4-μs molecular dynamics (MD) simulations at 400 K and 450 K using the monomer subunits of three CutA1 mutants from Escherichia coli: an SH-free mutant (Ec0SH) with denaturation temperature (T_d) = 85.6 °C, a hydrophobic mutant (Ec0VV) with T_d = 113.3 °C, and an ionic mutant (Ec0VV_6) with T_d = 136.8 °C. The occupancy of salt bridges by the six substituted charged residues in Ec0VV_6 was 140.1% at 300 K and 89.5% at 450 K, indicating that even in the denatured state, salt bridge occupancy was high, approximately 60% of that at 300 K. From these results, we can infer that proteins from hyperthermophiles with a high ratio of charged residues are stabilized by a decrease in conformational entropy due to ion–ion interactions in the denatured state. The mechanism must be comparable to the stabilization conferred by disulfide bonds within a protein. This suggests that introduction of charged residues, to promote formation of salt bridges in the denatured state, would be a simple way to rationally design stability-enhanced mutants.

Mechanics of Protein Adaptation to High Temperatures

Inspired by Somero's corresponding state principle that relates protein enhanced thermal stability with mechanical rigidity, we deployed state of the art computational techniques (based on atomistic steered molecular dynamics and Hamiltonian-replica exchange simulations) to study the in silico realization of mechanical and thermal unfolding of two homologous Csp proteins that have evolved to thrive in different thermal environments. By complementing recent single-molecule experiments, we unambiguously show that, for these homologues whose structures are very similar, the increased thermal resistance of the thermophilic variant is not associated with an increased mechanical stability. Our approach provides microscopic insights that are otherwise inaccessible to experimental techniques, and explains why the protein weak spots for thermal and mechanical denaturation are distinct.

Critical structural fluctuations of proteins upon thermal unfolding challenge the Lindemann criterion

Internal subnanosecond timescale motions are key for the function of proteins, and are coupled to the surrounding solvent environment. These fast fluctuations guide protein conformational changes, yet their role for protein stability, and for unfolding, remains elusive. Here, in analogy with the Lindemann criterion for the melting of solids, we demonstrate a common scaling of structural fluctuations of lysozyme protein embedded in different environments as the thermal unfolding transition is approached. By combining elastic incoherent neutron scattering and advanced molecular simulations, we show that, although different solvents modify the protein melting temperature, a unique dynamical regime is attained in proximity of thermal unfolding in all solvents that we tested. This solvation shell-independent dynamical regime arises from an equivalent sampling of the energy landscape at the respective melting temperatures. Thus, we propose that a threshold for the conformational entropy provided by structural fluctuations of proteins exists, beyond which thermal unfolding is triggered.

The extraordinary thermal stability of EstA from S. islandicus is independent of post translational modifications

Enzymes from thermophilic and hyper-thermophilic organisms have an intrinsic high stability. Understanding the mechanisms behind their high stability will be important knowledge for the engineering of novel enzymes with high stability. Lysine methylation of proteins is prevalent in Sulfolobus, a genus of hyperthermophilic and acidophilic archaea. Both unspecific and temperature dependent lysine methylations are seen, but the significance of this post-translational modification has not been investigated. Here, we test the effect of eliminating in vivo lysine methylation on the stability of an esterase (EstA). The enzyme was purified from the native host S. islandicus as well as expressed as a recombinant protein in E. coli, a mesophilic host that does not code for any machinery for in vivo lysine methylation. We find that lysine mono methylation indeed has a positive effect on the stability of EstA, but the effect is small. The effect of the lysine methylation on protein stability is secondary to that of protein expression in E. coli, as the E. coli recombinant enzyme is compromised both on stability and activity. We conclude that these differences are not attributed to any covalent difference between the protein expressed in hyperthermophilic versus mesophilic hosts.

Shedding light on the extra thermal stability of thermophilic proteins

An entropic stabilization mechanism has recently gained attention and credibility as the physical ground for the extra thermal stability of globular proteins from thermophilic microorganisms. An empirical result, obtained from the analysis of thermodynamic data for a large set of proteins, strengthens the general reliability of the theoretical approach originally devised to rationalize the occurrence of cold denaturation [Graziano, PCCP 2014, 16, 21755-21767]. It is shown that this theoretical approach can readily account for the entropic stabilization mechanism. On decreasing the conformational entropy gain associated with denaturation, the thermal stability of a model globular protein increases markedly.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

Role of Proteome Physical Chemistry in Cell Behavior

We review how major cell behaviors, such as bacterial growth laws, are derived from the physical chemistry of the cell’s proteins. On one hand, cell actions depend on the individual biological functionalities of their many genes and proteins. On the other hand, the common physics among proteins can be as important as the unique biology that distinguishes them. For example, bacterial growth rates depend strongly on temperature. This dependence can be explained by the folding stabilities across a cell’s proteome. Such modeling explains how thermophilic and mesophilic organisms differ, and how oxidative damage of highly charged proteins can lead to unfolding and aggregation in aging cells. Cells have characteristic time scales. For example, E. coli can duplicate as fast as 2–3 times per hour. These time scales can be explained by protein dynamics (the rates of synthesis and degradation, folding, and diffusional transport). It rationalizes how bacterial growth is slowed down by added salt. In the same way that the behaviors of inanimate materials can be expressed in terms of the statistical distributions of atoms and molecules, some cell behaviors can be expressed in terms of distributions of protein properties, giving insights into the microscopic basis of growth laws in simple cells.

Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins

We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.

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.

Interface matters: the stiffness route to stability of a thermophilic tetrameric malate dehydrogenase

In this work we investigate by computational means the behavior of two orthologous bacterial proteins, a mesophilic and a thermophilic tetrameric malate dehydrogenase (MalDH), at different temperatures. Namely, we quantify how protein mechanical rigidity at different length- and time-scales correlates to protein thermophilicity as commonly believed. In particular by using a clustering analysis strategy to explore the conformational space of the folded proteins, we show that at ambient conditions and at the molecular length-scale the thermophilic variant is indeed more rigid that the mesophilic one. This rigidification is the result of more efficient inter-domain interactions, the strength of which is further quantified via ad hoc free energy calculations. When considered isolated, the thermophilic domain is indeed more flexible than the respective mesophilic one. Upon oligomerization, the induced stiffening of the thermophilic protein propagates from the interface to the active site where the loop, controlling the access to the catalytic pocket, anchors down via an extended network of ion-pairs. On the contrary in the mesophilic tetramer the loop is highly mobile. Simulations at high temperature, could not re-activate the mobility of the loop in the thermophile. This finding opens questions on the similarities of the binding processes for these two homologues at their optimal working temperature and suggests for the thermophilic variant a possible cooperative role of cofactor/substrate.

Vibrational entropy differences between mesophile and thermophile proteins and their use in protein engineering

We recently introduced ENCoM, an elastic network atomic contact model, as the first coarse-grained normal mode analysis method that accounts for the nature of amino acids and can predict the effect of mutations on thermostability based on changes vibrational entropy. In this proof-of-concept article, we use pairs of mesophile and thermophile homolog proteins with identical structures to determine if a measure of vibrational entropy based on normal mode analysis can discriminate thermophile from mesophile proteins. We observe that in around 60% of cases, thermophile proteins are more rigid at equivalent temperatures than their mesophile counterpart and this difference can guide the design of proteins to increase their thermostability through series of mutations. We observe that mutations separating thermophile proteins from their mesophile orthologs contribute independently to a decrease in vibrational entropy and discuss the application and implications of this methodology to protein engineering.

Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains

In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Molecular basis of thermal stability in truncated (2/2) hemoglobins

BACKGROUND

Understanding the molecular mechanism through which proteins are functional at extreme high and low temperatures is one of the key issues in structural biology. To investigate this phenomenon, we have focused on two instructive truncated hemoglobins from Thermobifida fusca (Tf-trHbO) and Mycobacterium tuberculosis (Mt-trHbO); although the two proteins are structurally nearly identical, only the former is stable at high temperatures.

METHODS

We used molecular dynamics simulations at different temperatures as well as thermal melting profile measurements of both wild type proteins and two mutants designed to interchange the amino acid residue, either Pro or Gly, at E3 position.

RESULTS

The results show that the presence of a Pro at the E3 position is able to increase (by 8°) or decrease (by 4°) the melting temperature of Mt-trHbO and Tf-trHbO, respectively. We observed that the ProE3 alters the structure of the CD loop, making it more flexible.

CONCLUSIONS

This gain in flexibility allows the protein to concentrate its fluctuations in this single loop and avoid unfolding. The alternate conformations of the CD loop also favor the formation of more salt-bridge interactions, together augmenting the protein's thermostability.

GENERAL SIGNIFICANCE

These results indicate a clear structural and dynamical role of a key residue for thermal stability in truncated hemoglobins.

Enhancing thermostability and the structural characterization of Microbacterium saccharophilum K-1 β-fructofuranosidase

A β-fructofuranosidase from Microbacterium saccharophilum K-1 (formerly known as Arthrobacter sp. K-1) is useful for producing the sweetener lactosucrose (4(G)-β-D-galactosylsucrose). Thermostability of the β-fructofuranosidase was enhanced by random mutagenesis and saturation mutagenesis. Clones with enhanced thermostability included mutations at residues Thr47, Ser200, Phe447, Phe470, and Pro500. In the highest stability mutant, T47S/S200T/F447P/F470Y/P500S, the half-life at 60 °C was 182 min, 16.5-fold longer than the wild-type enzyme. A comparison of the crystal structures of the full-length wild-type enzyme and three mutants showed that various mechanisms appear to be involved in thermostability enhancement. In particular, the replacement of Phe447 with Val or Pro induced a conformational change in an adjacent residue His477, which results in the formation of a new hydrogen bond in the enzyme. Although the thermostabilization mechanisms of the five residue mutations were explicable on the basis of the crystal structures, it appears to be difficult to predict which amino acid residues should be modified to obtain thermostabilized enzymes.

The molecular basis of conformational instability of the ecdysone receptor DNA binding domain studied by in silico and in vitro experiments

The heterodimer of the ecdysone receptor (EcR) and ultraspiracle (Usp), members of the nuclear receptors superfamily, regulates gene expression associated with molting and metamorphosis in insects. The DNA binding domains (DBDs) of the Usp and EcR play an important role in their DNA-dependent heterodimerization. Analysis of the crystal structure of the UspDBD/EcRDBD heterocomplex from Drosophila melanogaster on the hsp27 gene response element, suggested an appreciable similarity between both DBDs. However, the chemical denaturation experiments showed a categorically lower stability for the EcRDBD in contrast to the UspDBD. The aim of our study was an elucidation of the molecular basis of this intriguing instability. Toward this end, we mapped the EcRDBD amino acid sequence positions which have an impact on the stability of the EcRDBD. The computational protein design and in vitro analyses of the EcRDBD mutants indicate that non-conserved residues within the α-helix 2, forming the EcRDBD hydrophobic core, represent a specific structural element that contributes to instability. In particular, the L58 appears to be a key residue which differentiates the hydrophobic cores of UspDBD and EcRDBD and is the main reason for the low stability of the EcRDBD. Our results might serve as a benchmark for further studies of the intricate nature of the EcR molecule.

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

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

Structure and dynamics of G-protein coupled receptors

G-protein coupled receptors (GPCRs) are seven helical transmembrane proteins that mediate cell-to-cell communication. They also form the largest superfamily of drug targets. Hence detailed studies of the three dimensional structure and dynamics are critical to understanding the functional role of GPCRs in signal transduction pathways, and for drug design. In this chapter we compare the features of the crystal structures of various biogenic amine receptors, such as β1 and β2 adrenergic receptors, dopamine D3 receptor, M2 and M3 muscarinic acetylcholine receptors. This analysis revealed that conserved residues are located facing the inside of the transmembrane domain in these GPCRs improving the efficiency of packing of these structures. The NMR structure of the chemokine receptor CXCR1 without any ligand bound, shows significant dynamics of the transmembrane domain, especially the helical kink angle on the transmembrane helix6. The activation mechanism of the β2-adrenergic receptor has been studied using multiscale computational methods. The results of these studies showed that the receptor without any ligand bound, samples conformations that resemble some of the structural characteristics of the active state of the receptor. Ligand binding stabilizes some of the conformations already sampled by the apo receptor. This was later observed in the NMR study of the dynamics of human β2-adrenergic receptor. The dynamic nature of GPCRs leads to a challenge in obtaining purified receptors for biophysical studies. Deriving thermostable mutants of GPCRs has been a successful strategy to reduce the conformational heterogeneity and stabilize the receptors. This has lead to several crystal structures of GPCRs. However, the cause of how these mutations lead to thermostability is not clear. Computational studies are beginning to shed some insight into the possible structural basis for the thermostability. Molecular Dynamics simulations studying the conformational ensemble of thermostable mutants have shown that the stability could arise from both enthalpic and entropic factors. There are regions of high stress in the wild type GPCR that gets relieved upon mutation conferring thermostability.