A Comparative Molecular Dynamics Study of Thermophilic and Mesophilic Ribonuclease HI Enzymes

Abstract

We studied a pair of homologous thermophilic and mesophilic ribonuclease HI enzymes by molecular dynamics simulations. Each protein was subjected to three 5 ns simulations in explicit water at both 310 K and 340 K. The thermophilic enzyme showed larger overall positional fluctuations at both temperatures, while only the mesophilic enzyme at the higher temperature showed significant instability. When the temperature is changed, the relative flexibility of different local segments on the two proteins changed differently. Principal component analysis showed that the simulations of the two proteins explored largely overlapping regions in the conformational space. However, at 340 K, the collective structure variations of the thermophilic protein are different from those of the mesophilic protein. Our results, although not in accordance with the view that hyperthermostability of proteins may originate from their conformational rigidity, are consistent with several recent experimental and simulation studies which showed that thermophilic proteins may be conformationally more flexible than their mesophilic counterparts. The decorrelation between conformational rigidity and hyperthermostability may be attributed to the temperature dependence and long range nature of electrostatic interactions that play more important roles in the structural stability of thermophilic proteins.     

Relation between stability, dynamics and enzyme activity in 3-phosphoglycerate kinases from yeast and Thermus thermophilus.

3-Phosphoglycerate kinases from yeast and the extreme thermophilic bacterium Thermus thermophilus HB8 have been used as models for investigating the relationship between stability, dynamics and activity. It was found that while at a given temperature the thermophilic protein is more stable, its conformational dynamics as measured by the ability of acrylamide to quench the fluorescence of a buried tryptophan as well as its specific activity, are both lower than for the mesophilic protein. As the temperature is increased, the thermodynamic stability of the thermophilic protein approaches that of the mesophilic protein at its working temperature. Its conformational dynamics and specific activity however were both shown to increase, until at the physiologically operational temperature, they become similar to those of the mesophilic enzyme at its operational temperature. These results confirm the proposal that a direct relationship and balance holds between thermodynamic stability, dynamics and specific activity in globular proteins. They demonstrate also the constraining effect of increased stability upon conformational dynamics and enzyme activity.     

Role of loop dynamics in thermal stability of mesophilic and thermophilic adenylosuccinate synthetase: a molecular dynamics and normal mode analysis study

Enzymes from thermophiles are poorly active at temperatures at which their mesophilic homologs exhibit high activity and attain corresponding active states at high temperatures. In this study, comparative molecular dynamics (MD) simulations, supplemented by normal mode analysis, have been performed on an enzyme Adenylosuccinate synthetase (AdSS) from E. coli (mesophilic) and P. horikoshii (thermophilic) systems to understand the effects of loop dynamics on thermal stability of AdSS. In mesophilic AdSS, both ligand binding and catalysis are facilitated through the coordinated movement of five loops on the protein. The simulation results suggest that thermophilic P. horikoshii preserves structure and catalytic function at high temperatures by using the movement of only a subset of loops (two out of five) for ligand binding and catalysis unlike its mesophilic counterpart in E. coli. The pre-arrangement of the catalytic residues in P. horikoshii is well-preserved and salt bridges remain stable at high temperature (363K). The simulations suggest a general mechanism (including pre-arrangement of catalytic residues, increased polar residue content, stable salt bridges, increased rigidity, and fewer loop movements) used by thermophilic enzymes to preserve structure and be catalytically active at elevated temperatures.     

Comparative molecular dynamics simulation studies for determining factors contributing to the thermostability of chemotaxis protein “CheY”

Comparative molecular dynamics simulations of chemotaxis protein “CheY” from thermophilic origin Thermotoga maritima and its mesophilic counterpart Salmonella enterica have been performed for 10?ns each at 300 and 350?K, and 20?ns each at 400 and 450?K. The trajectories were analyzed in terms of different factors like root-mean-square deviation, root-mean-square fluctuation, radius of gyration, solvent accessible surface area, H-bonds, salt bridge content, and protein–solvent interactions which indicate distinct differences between the two of them. The two proteins also follow dissimilar unfolding pathways. The overall flexibility calculated by the trace of the diagonalized covariance matrix displays similar flexibility of both the proteins near their optimum growth temperatures. However, at higher temperatures mesophilic protein shows increased overall flexibility than its thermophilic counterpart. Principal component analysis also indicates that the essential subspaces explored by the simulations of two proteins at different temperatures are nonoverlapping and they show significantly different directions of motion. However, there are significant overlaps within the trajectories and similar direction of motions are observed for both proteins at 300?K. Overall, the mesophilic protein leads to increased conformational sampling of the phase space than its thermophilic counterpart. This is the first ever study of thermostability of CheY protein homologs by using protein dynamism as a main impact. Our study might be used as a model for studying the molecular basis of thermostability of two homologous proteins from two organisms living at different temperatures with less visible differences.     

Solution structure and backbone dynamics of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin mutant: a molecular analysis of its reduced thermal stability

No general strategy for thermostability has been yet established, because the extra stability of thermophiles appears to be the sum of different cumulative stabilizing interactions. In addition, the increase of conformational rigidity observed in many thermophilic proteins, which in some cases disappears when mesophilic and thermophilic proteins are compared at their respective physiological temperatures, suggests that evolutionary adaptation tends to maintain corresponding states with respect to conformational flexibility. In this study, we accomplished a structural analysis of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin (BacTrx) mutant, which has reduced heat resistance with respect to the thermostable wild-type. Furthermore, we have also achieved a detailed study, carried out at 25, 45, and 65 degrees C, of the backbone dynamics of both the BacTrx and its K18G/R82E mutant. Our findings clearly indicate that the insertion of the two mutations causes a loss of energetically favorable long-range interactions and renders the secondary structure elements of the double mutants more similar to those of the mesophilic Escherichia coli thioredoxin. Moreover, protein dynamics analysis shows that at room temperature the BacTrx, as well as the double mutant, are globally as rigid as the mesophilic thioredoxins; differently, at 65 degrees C, which is in the optimal growth temperature range of A. acidocaldarius, the wild-type retains its rigidity while the double mutant is characterized by a large increase of the amplitude of the internal motions. Finally, our research interestingly shows that fast motions on the pico- to nanosecond time scale are not detrimental to protein stability and provide an entropic stabilization of the native state. This study further confirms that protein thermostability is reached through diverse stabilizing interactions, which have the key role to maintain the structural folding stable and functional at the working temperature.     

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.     

A rigidifying salt-bridge favors the activity of thermophilic enzyme at high temperatures at the expense of low-temperature activity

Background

Thermophilic enzymes are often less active than their mesophilic homologues at low temperatures. One hypothesis to explain this observation is that the extra stabilizing interactions increase the rigidity of thermophilic enzymes and hence reduce their activity. Here we employed a thermophilic acylphosphatase from Pyrococcus horikoshii and its homologous mesophilic acylphosphatase from human as a model to study how local rigidity of an active-site residue affects the enzymatic activity.

Methods and Findings

Acylphosphatases have a unique structural feature that its conserved active-site arginine residue forms a salt-bridge with the C-terminal carboxyl group only in thermophilic acylphosphatases, but not in mesophilic acylphosphatases. We perturbed the local rigidity of this active-site residue by removing the salt-bridge in the thermophilic acylphosphatase and by introducing the salt-bridge in the mesophilic homologue. The mutagenesis design was confirmed by x-ray crystallography. Removing the salt-bridge in the thermophilic enzyme lowered the activation energy that decreased the activation enthalpy and entropy. Conversely, the introduction of the salt-bridge to the mesophilic homologue increased the activation energy and resulted in increases in both activation enthalpy and entropy. Revealed by molecular dynamics simulations, the unrestrained arginine residue can populate more rotamer conformations, and the loss of this conformational freedom upon the formation of transition state justified the observed reduction in activation entropy.

Conclusions

Our results support the conclusion that restricting the active-site flexibility entropically favors the enzymatic activity at high temperatures. However, the accompanying enthalpy-entropy compensation leads to a stronger temperature-dependency of the enzymatic activity, which explains the less active nature of the thermophilic enzymes at low temperatures.     

Investigations of the thermostability of rubredoxin models using molecular dynamics simulations.

The affects of differences in amino acid sequence on the temperature stability of the three-dimensional structure of the small beta-sheet protein, rubredoxin (Rd), was revealed when a set of homology models was subjected to molecular dynamics simulations at relatively high temperatures. Models of Rd from the hyperthermophile, Pyrococcus furiosus (Pf), an organism that grows optimally at 100 degrees C, were compared to three mesophilic Rds of known X-ray crystal structure. Simulations covering the limits of known Rd thermostabilities were carried out at temperatures of 300 K, 343 K, 373 K, and 413 K. They suggest that Rd stability is correlated with structural dynamics. Because the dynamic behavior of three Pf Rd models was consistently different from the dynamic behavior of the three mesophilic Rd structures, detailed analysis of the temperature-dependent dynamic behavior was carried out. The major differences between the models of the protein from the hyperthermophile and the others were: (1) an obvious temperature-dependent transition in the mesophilic structures not seen with the Pf Rd models, (2) consistent AMBER energy for the Pf Rd due to differences in nonbonded interaction terms, (3) less variation in the average conformations for the Pf Rd models with temperature, and (4) the presence of more extensive secondary structure for the Pf Rd models. These unsolvated dynamics simulations support a simple, general hypothesis to explain the hyperthermostability of Pf Rd. Its structure simplifies the conformational space to give a single minimum accessible over an extreme range of temperatures, whereas the mesophilic proteins sample a more complex conformational space with two or more minima over the same temperature range.     

Model-free analysis of a thermophilic Fe(7)S(8) protein compared with a mesophilic Fe(4)S(4) protein

15N T(1), T(2) and (1)H-(15)N NOE were measured for the thermophilic Fe(7)S(8) protein from Bacillus schlegelii and for the Fe(4)S(4) HiPIP protein from Chromatium vinosum, which is a mesophilic protein. The investigation was performed at 276, 300, and 330 K at 11.7 T for the former, whereas only the 298 K data at 14.1 T for the latter were acquired. The data were analyzed with the model-free protocol after correcting the measured parameters for the effect of paramagnetism, because both proteins are paramagnetic. Both thermophilic and mesophilic proteins are quite rigid, with an average value of the generalized order parameter S2at room temperature of 0.92 and 0.94 for Fe(7)S(8) and Fe(4)S(4) proteins, respectively. The analyzed nitrogens for the Fe(7)S(8) protein showed a significant decrease in S2with increasing temperature, and at the highest temperature >70% of the residues had an internal correlation time. This research shows that subnanosecond rigidity is not related to thermostability and provides an estimate of the effect of increasing temperature on this time scale.     

Exploiting the Link between Protein Rigidity and Thermostability for Data‐Driven Protein Engineering

Understanding and exploiting the relationship between microscopic structure and macroscopic stability is important for developing strategies to improve protein stability at high temperatures. The thermostability of proteins has been repeatedly linked to an enhanced structural rigidity of the folded native state. In the current study, the rigidity of protein structures from mesophilic and thermophilic organisms along a thermal unfolding trajectory is directly probed. In order to perform this, protein structures were modeled as constraint networks, and the rigidity in these networks was quantified using the Floppy Inclusion and Rigid Substructure Topography (FIRST) method. During the thermal unfolding, a phase transition was observed that defines the rigidity percolation threshold and corresponds to the folded‐unfolded transition in protein folding. Using concepts from percolation theory and network science, a higher phase transition temperature was observed for ca. two‐thirds of the proteins from thermophilic organisms compared to their mesophilic counterparts, when applied to a data set of 20 pairs of homologues. From both the analysis of the microstructure of the constraint networks and monitoring the macroscopic behavior during the thermal unfolding, direct evidence was found for the “corresponding states” concept, which states that mesophilic and thermophilic enzymes are in corresponding states of similar flexibility at their respective optimal temperature. Finally, the current approach facilitated the identification of structural features from which a destabilization of the structure originates upon thermal unfolding. These predictions show a good agreement with the experimental data. Therefore, the information might be exploited in data‐driven protein engineering by pointing to residues that should be varied to obtain a protein with higher thermostability.     

Effect of molecular confinement on internal enzyme dynamics: frequency domain fluorometry and molecular dynamics simulation studies

The tryptophanyl emission decay of the mesophilic beta-galactosidase from Aspergillus oryzae free in buffer and entrapped in agarose gel is investigated as a function of temperature and compared to that of the hyperthermophilic enzyme from Sulfolobus solfataricus. Both enzymes are tetrameric proteins with a large number of tryptophanyl residues, so the fluorescence emission can provide information on the conformational dynamics of the overall protein structure rather than that of the local environment. The tryptophanyl emission decays are best fitted by bimodal Lorentzian distributions. The long-lived component is ascribed to close, deeply buried tryptophanyl residues with reduced mobility; the short-lived one arises from tryptophanyl residues located in more flexible external regions of each subunit, some of which are involved in forming the catalytic site. The center of both lifetime distribution components at each temperature increases when going from the free in solution mesophilic enzyme to the gel-entrapped and hyperthermophilic enzyme, thus indicating that confinement of the mesophilic enzyme in the agarose gel limits the freedom of the polypeptide chain. A more complex dependence is observed for the distribution widths. Computer modeling techniques are used to recognize that the catalytic sites are similar for the mesophilic and hyperthermophilic beta-galactosidases. The effect due to gel entrapment is considered in dynamic simulations by imposing harmonic restraints to solvent-exposed atoms of the protein with the exclusion of those around the active site. The temperature dependence of the tryptophanyl fluorescence emission decay and the dynamic simulation confirm that more rigid structures, as in the case of the immobilized and/or hyperthermophilic enzyme, require higher temperatures to achieve the requisite conformational dynamics for an effective catalytic action and strongly suggest a link between conformational rigidity and enhanced thermal stability.     

Stability and rigidity/flexibility—Two sides of the same coin?

Protein molecules require both flexibility and rigidity for functioning. The fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. We have determined flexible regions for four homologous pairs from thermophilic and mesophilic organisms by two methods: the fast FoldUnfold which uses amino acid sequence and the time consuming MDFirst which uses three-dimensional structures. We demonstrate that both methods allow determining flexible regions in protein structure. For three of the four thermophile–mesophile pairs of proteins, FoldUnfold predicts practically the same flexible regions which have been found by the MD/First method. As expected, molecular dynamics simulations show that thermophilic proteins are more rigid in comparison to their mesophilic homologues. Analysis of rigid clusters and their decomposition provides new insights into protein stability. It has been found that the local networks of salt bridges and hydrogen bonds in thermophiles render their structure more stable with respect to fluctuations of individual contacts. Such network includes salt bridge triads Agr-Glu-Lys and Arg-Glu-Arg, or salt bridges (such as Arg-Glu) connected with hydrogen bonds. This ionic network connects alpha helices and rigidifies the structure. Mesophiles can be characterized by stand alone salt bridges and hydrogen bonds or small ionic clusters. Such difference in the network of salt bridges results in different flexibility of homologous proteins. Combining both approaches allows characterizing structural features in atomic detail that determine the rigidity/flexibility of a protein structure. This article is a part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

A simple strategy for the purification of large thermophilic proteins overexpressed in mesophilic microorganisms: application to multimeric enzymes from Thermus sp. strain T2 expressed in Escherichia coli

The heating of protein preparations of mesophilic organism (e.g., E. coli) produces the obliteration of all soluble multimeric proteins from this organism. In this way, if a multimeric enzyme from a thermophilic microorganism is expressed in these mesophilic hosts, the only large protein remaining soluble in the preparation after heating is the thermophilic enzyme. These large proteins may be then selectively adsorbed on lowly activated anionic exchangers, enabling their full purification in just these two simple steps. This strategy has been applied to the purification of an alpha-galactosidase and a beta-galactosidase from Thermus sp. strain T2, both expressed in E. coli, achieving the almost full purification of both enzymes in only these two simple steps. This very simple strategy seems to be of general applicability to the purification of any thermophilic multimeric enzyme expressed in a mesophilic host.     

A Comparison of Structural and Evolutionary Attributes of Escherichia coli and Thermus thermophilus Small Ribosomal Subunits: Signatures of Thermal Adaptation

Here we compare the structural and evolutionary attributes of Thermus thermophilus and Escherichia coli small ribosomal subunits (SSU). Our results indicate that with few exceptions, thermophilic 16S ribosomal RNA (16S rRNA) is densely packed compared to that of mesophilic at most of the analogous spatial regions. In addition, we have located species-specific cavity clusters (SSCCs) in both species. E. coli SSCCs are numerous and larger compared to T. thermophilus SSCCs, which again indicates densely packed thermophilic 16S rRNA. Thermophilic ribosomal proteins (r-proteins) have longer disordered regions than their mesophilic homologs and they experience larger disorder-to-order transitions during SSU-assembly. This is reflected in the predicted higher conformational changes of thermophilic r-proteins compared to their mesophilic homologs during SSU-assembly. This high conformational change of thermophilic r-proteins may help them to associate with the 16S ribosomal RNA with high complementary interfaces, larger interface areas, and denser molecular contacts, compared to those of mesophilic. Thus, thermophilic protein-rRNA interfaces are tightly associated with 16S rRNA than their mesophilic homologs. Densely packed 16S rRNA interior and tight protein-rRNA binding of T. thermophilus (compared to those of E. coli) are likely the signatures of its thermal adaptation. We have found a linear correlation between the free energy of protein-RNA interface formation, interface size, and square of conformational changes, which is followed in both prokaryotic and eukaryotic SSU. Disorder is associated with high protein-RNA interface polarity. We have found an evolutionary tendency to maintain high polarity (thereby disorder) at protein-rRNA interfaces, than that at rest of the protein structures. However, some proteins exhibit exceptions to this general trend.     

Structural features of cold-adapted dimeric GH2 β-D-galactosidase from Arthrobacter sp. 32cB

Crystal structures of cold-adapted β-d-galactosidase (EC 3.2.1.23) from the Antarctic bacterium Arthrobacter sp. 32cB (ArthβDG) have been determined in an unliganded form resulting from diffraction experiments conducted at 100 K (at resolution 1.8 Å) and at room temperature (at resolution 3.0 Å). A detailed comparison of those two structures of the same enzyme was performed in order to estimate differences in their molecular flexibility and rigidity and to study structural rationalization for the cold-adaptation of the investigated enzyme. Furthermore, a comparative analysis with structures of homologous enzymes from psychrophilic, mesophilic, and thermophilic sources has been discussed to elucidate the relationship between structure and cold-adaptation in a wider context. The performed studies confirm that the structure of cold-adapted ArthβDG maintains balance between molecular stability and structural flexibility, which can be observed independently on the temperature of conducted X-ray diffraction experiments. Obtained information about proper protein function under given conditions provide a guideline for rational engineering of proteins in terms of their temperature optimum and thermal stability.     

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

Protein rigidity and thermophilic adaptation

We probe the hypothesis of corresponding states, according to which homologues from mesophilic and thermophilic organisms are in corresponding states of similar rigidity and flexibility at their respective optimal temperatures. For this, the local distribution of flexible and rigid regions in 19 pairs of homologous proteins from meso- and thermophilic organisms is analyzed and related to activity characteristics of the enzymes by constraint network analysis (CNA). Two pairs of enzymes are considered in more detail: 3-isopropylmalate dehydrogenase and thermolysin-like protease. By comparing microscopic stability features of homologues with the help of stability maps, introduced for the first time, we show that adaptive mutations in enzymes from thermophilic organisms maintain the balance between overall rigidity, important for thermostability, and local flexibility, important for activity, at the appropriate working temperature. Thermophilic adaptation in general leads to an increase of structural rigidity but conserves the distribution of functionally important flexible regions between homologues. This finding provides direct evidence for the hypothesis of corresponding states. CNA thereby implicitly captures and unifies many different mechanisms that contribute to increased thermostability and to activity at high temperatures. This allows to qualitatively relate changes in the flexibility of active site regions, induced either by a temperature change or by the introduction of mutations, to experimentally observed losses of the enzyme function. As for applications, the results demonstrate that exploiting the principle of corresponding states not only allows for successful thermostability optimization but also for guiding experiments in order to improve enzyme activity in protein engineering.     

Molecular dynamics study of a hyperthermophilic and a mesophilic rubredoxin

In recent years, increased interest in the origin of protein thermal stability has gained attention both for its possible role in understanding the forces governing the folding of a protein and for the design of new highly stable engineered biocatalysts. To study the origin of thermostability, we have performed molecular dynamics simulations of two rubredoxins, from the mesophile Clostridium pasteurianum and from the hyperthermophile Pyrococcus furiosus. The simulations were carried out at two temperatures, 300 and 373 K, for each molecule. The length of the simulations was within the range of 6-7.2 ns. The rubredoxin from the hyperthermophilic organism was more flexible than its mesophilic counterpart at both temperatures; however, the overall flexibility of both molecules at their optimal growth temperature was the same, despite 59% sequence homology. The conformational space sampled by both molecules was larger at 300 K than at 373 K. The essential dynamics analysis showed that the principal overall motions of the two molecules are significantly different. On the contrary, each molecule showed similar directions of motion at both temperatures.     

Protein dynamics and stability: the distribution of atomic fluctuations in thermophilic and mesophilic dihydrofolate reductase derived using elastic incoherent neutron scattering

The temperature dependence of the dynamics of mesophilic and thermophilic dihydrofolate reductase is examined using elastic incoherent neutron scattering. It is demonstrated that the distribution of atomic displacement amplitudes can be derived from the elastic scattering data by assuming a (Weibull) functional form that resembles distributions seen in molecular dynamics simulations. The thermophilic enzyme has a significantly broader distribution than its mesophilic counterpart. Furthermore, although the rate of increase with temperature of the atomic mean-square displacements extracted from the dynamic structure factor is found to be comparable for both enzymes, the amplitudes are found to be slightly larger for the thermophilic enzyme. Therefore, these results imply that the thermophilic enzyme is the more flexible of the two.