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.     

Dynamics and unfolding pathways of a hyperthermophilic and a mesophilic rubredoxin.

Molecular dynamics simulations in solution are performed for a rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus (RdPf) and one from the mesophilic organism Desulfovibrio vulgaris (RdDv). The two proteins are simulated at four temperatures: 300 K, 373 K, 473 K (two sets), and 500 K; the various simulations extended from 200 ps to 1,020 ps. At room temperature, the two proteins are stable, remain close to the crystal structure, and exhibit similar dynamic behavior; the RMS residue fluctuations are slightly smaller in the hyperthermophilic protein. An analysis of the average energy contributions in the two proteins is made; the results suggest that the intraprotein energy stabilizes RdPf relative to RdDv. At 373 K, the mesophilic protein unfolds rapidly (it begins to unfold at 300 ps), whereas the hyperthermophilic does not unfold over the simulation of 600 ps. This is in accord with the expected stability of the two proteins. At 473 K, where both proteins are expected to be unstable, unfolding behavior is observed within 200 ps and the mesophilic protein unfolds faster than the hyperthermophilic one. At 500 K, both proteins unfold; the hyperthermophilic protein does so faster than the mesophilic protein. The unfolding behavior for the two proteins is found to be very similar. Although the exact order of events differs from one trajectory to another, both proteins unfold first by opening of the loop region to expose the hydrophobic core. This is followed by unzipping of the beta-sheet. The results obtained in the simulation are discussed in terms of the factors involved in flexibility and thermostability.     

A comparative molecular dynamics study of thermophilic and mesophilic ribonuclease HI enzymes

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.     

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.     

Dynamic arrangement of ion pairs and individual contributions to the thermal stability of the cofactor-binding domain of glutamate dehydrogenase from Thermotoga maritima

The dynamics of a hyperthermophilic protein fragment in a water environment, as studied by performing molecular dynamics (MD) simulations at various temperatures, is compared to the dynamical behavior of a homologous mesophilic protein simulated under identical conditions. The effects on the stability of the spatial arrangement and mobility of the charged residues in solution were quantified by calculating free energy changes upon salt bridge formation in these proteins. Electrostatic free energy terms derived from a thermodynamic cycle were obtained by solving the linearized Poisson-Boltzmann equation for a series of protein conformations generated by MD simulations and placed subsequently in a continuum solvent medium. Our results show that the ion pairs are electrostatically stabilizing in most of the cases, but their individual contributions vary significantly. The greater contribution of the charged residues to the stability of the hyperthermophilic protein as compared with the mesophilic counterpart was evidenced only by the calculations that included conformations sampled at 343 and 373 K. The "dynamic" structure of the hyperthermophilic protein fragment simulated at elevated temperatures reveals an optimum placement of the ionizable residues within the protein structure as well as the role of their cooperative interactions in promoting thermal stability. The thermodynamic properties such as electrostatic free energy differences, configurational entropies, and specific heat capacities calculated in the dynamic context of the protein structure provided new insight into the mechanism of protein thermostabilization.     

Effects of environment on the structure of Pyrococcus furiosus rubredoxin: a molecular dynamics study

Molecular dynamics simulations based on a 0.95-A resolution crystal structure of Pyrococcus furiosus have been performed to elucidate the effects of the environment on the structure of rubredoxin, and proteins in general. Three 1-ns simulations are reported here: two crystalline state simulations at 123 and 300 K, and a solution state simulation at 300 K. These simulations show that temperature has a greater impact on the protein structure than the close molecular contacts of the crystal matrix in rubredoxin, although both have an effect on its dynamic properties. These results indicate that differences between NMR solution structures and X-ray crystal structures will be relatively minor if they are done at similar temperatures. In addition, the crystal simulations appears to mimic previous crystallographic experiments on the effects of cryo-temperature on temperature factors, and might provide a useful tool in the structural analysis of protein structures solved at cryo-temperatures.     

Enhanced thermal stability achieved without increased conformational rigidity at physiological temperatures: spatial propagation of differential flexibility in rubredoxin hybrids

The extreme thermal stability of proteins from hyperthermophilic organisms is widely believed to arise from an increased conformational rigidity in the native state. In apparent contrast to this paradigm, both Pyrococcus furiosus (Pf) rubredoxin, the most thermostable protein characterized to date, and its Clostridium pasteurianum (Cp) mesophile homolog undergo a transient conformational opening of their multi-turn segments, which is more favorable in hyperthermophile proteins below room temperature. Substitution of the hyperthermophile multi-turn sequence into the mesophile protein sequence yields a hybrid, (14-33(Pf)) Cp, that exhibits a 12 degrees increase in its reversible thermal unfolding transition midpoint. Nuclear magnetic resonance (NMR) magnetization transfer-based hydrogen exchange was used to monitor backbone conformational dynamics in the subsecond time regime. Despite the substantially increased thermostability, flexibility throughout the entire main chain of the more thermostable hybrid is equal to or greater than that of the wild type mesophile rubredoxin near its normal growth temperature. In comparison to the identical core residues of the (14-33(Pf)) Cp rubredoxin hybrid, six spatially clustered residues in the parental mesophile protein exhibit a substantially larger temperature dependence of exchange. The exchange behavior of these six residues closely matches that observed in the multi-turn segment, consistent with a more extensive conformational process. These six core residues exhibit a much weaker temperature dependence of exchange in the (14-33(Pf)) Cp hybrid, similar to that observed for the multi-turn segment in its parental Pf rubredoxin. These results suggest that differential temperature dependence of flexibility can underlie variations in thermostability observed for mesophile versus hyperthermophile homologs.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson–Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.     

Crystal structure of rubredoxin from Pyrococcus furiosus at 0.95 Å resolution, and the structures of N-terminal methionine and formylmethionine variants of Pf Rd. Contributions of N-terminal interactions to thermostability

The high-resolution crystal structure of the small iron-sulfur protein rubredoxin (Rd) from the hyperthermophilic archeon Pyrococcus furiosus (Pf) is reported in this paper, together with those of its methionine ([_0M]Pf Rd) and formylmethionine (f[_0M]Pf Rd) variants. These studies were conducted to assess the consequences of the presence or absence of a salt bridge between the amino terminal nitrogen of Ala1 and the side chain of Glu14 to the structure and stability of this rubredoxin. The structure of wild-type Pf Rd was solved to a resolution of 0.95?Å and refined by full-matrix least-squares techniques to a crystallographic agreement factor of 12.8% [F>2σ(F) data, 25?617 reflections], while those of the [_0M]Pf and f[_0M]Pf Rd variants were solved at slightly lower resolutions (1.1?Å, R=11.5%, 17?213 reflections; 1.2?Å, R=13.7%, 12?478 reflections, respectively). The quality of the data was such that about half of the hydrogen atoms of the protein were clearly visible. All three structures were ultimately refined using the program SHELXL-93 with anisotropic atomic displacement parameters for all non-hydrogen protein atoms, and calculated hydrogen positions included but not refined. In this paper we also report thermostability data for all three forms of Pf Rd, and show that they follow the sequence wild-type >[_0M]Pf>formyl[_0M]Pf. Comparison of the three Pf Rd structures in the N-terminal region show that the structures of wild-type Pf Rd and f[_0M]Pf are rather similar, while that of [_0M]Pf Rd shows a number of additional hydrogen bonds involving the extra methionine group. While the salt bridge between the Ala1 amino group and the Glu14 carboxylate is not the primary determinant of the thermostability of Pf Rd, alterations to the amino terminus do have a moderate influence on the thermostability of this protein.     

Molecular dynamics simulations of the thermal stability of tteRBP and ecRBP

Molecular dynamics simulations were performed for investigating the thermal stability of the extremely thermophilic Thermoanaerobacter tengcongensis ribose binding protein (tteRBP) and the mesophilic homologous Escherichia coli ribose binding protein (ecRBP). The simulations for the two proteins were carried out under the room temperature (300?K) and the optimal activity temperature (tteRBP 375?K and ecRBP 329?K), respectively. The comparative analyses of the trajectories show that the two proteins have stable overall structures at the two temperatures; further analyses indicate that they both have strong side-chain interactions and different backbone flexibilities at the different temperatures. The tteRBP 375?K and ecRBP 329?K have stronger internal motion and higher flexibility than tteRBP 300?K and ecRBP 300?K, respectively, it is noted that the flexibility of tteRBP is much higher than that of ecRBP at the two temperatures. Therefore, tteRBP 375?K can adapt to high temperature due to its higher flexibility of backbone. Combining with the researches by Cuneo et al., it is concluded that the side-chain interactions and flexibility of backbone are both the key factors to maintain thermal stability of the two proteins.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:22     

Reduced temperature dependence of collective conformational opening in a hyperthermophile rubredoxin.

Spatially localized differences in the conformational dynamics of the rubredoxins from the hyperthermophile Pyrococcus furiosus (Pf) and the mesophile Clostridium pasteurianum (Cp) are monitored via amide exchange measurements. As shown previously for the hyperthermophile protein, nearly all backbone amides of the Cp rubredoxin exhibit EX(2) hydrogen exchange kinetics with conformational opening rates of >1 s(-)(1). Significantly slower amide exchange is observed for Pf rubredoxin in the region surrounding the metal site and the proximal end of the three-stranded beta-sheet, while for the rest of the structure, the exchange rates at 23 degrees C are similar for both proteins. For the multiple-turn region comprising residues 14-32 in both rubredoxins, the uniformity of both the exchange rate constants and the values of the activation energy at the slowly exchanging sites is consistent with a model of solvent exposure via a subglobal cooperative conformational opening. In contrast to the common expectation of increased rigidity in the hyperthermophile proteins, below room temperature Pf rubredoxin exhibits a larger apparent flexibility in this multiple-turn region. The smaller enthalpy for the conformational opening process of this region in Pf rubredoxin reflects the much weaker temperature dependence of the underlying conformational equilibrium in the hyperthermophile protein.     

A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

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.     

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.     

Adaptation to high temperatures through macromolecular dynamics by neutron scattering

Work on the relationship between hyperthermophile protein dynamics, stability and activity is reviewed. Neutron spectroscopy has been applied to measure and compare the macromolecular dynamics of various hyperthermophilic and mesophilic proteins, under different conditions. First, molecular dynamics have been analyzed for the hyperthermophile malate dehydrogenase from Methanococcus jannaschii and a mesophilic homologue, the lactate dehydrogenase from Oryctolagus cunniculus (rabbit) muscle. The neutron scattering approach has provided independent measurements of the global flexibility and structural resilience of each protein, and it has been demonstrated that macromolecular dynamics represents one of the molecular mechanisms of thermoadaptation. The resilience was found to be higher for the hyperthermophilic protein, thus ensuring similar flexibilities in both enzymes at their optimal activity temperature. Second, the neutron method has been developed to quantify the average macromolecular flexibility and resilience within the natural crowded environment of the cell, and mean macromolecular motions have been measured in vivo in psychrophile, mesophile, thermophile and hyperthermophile bacteria. The macromolecular resilience in bacteria was found to increase with adaptation to high temperatures, whereas flexibility was maintained within narrow limits, independent of physiological temperature for all cells in their active state. Third, macromolecular motions have been measured in free and immobilized dihydrofolate reductase from Escherichia coli. The immobilized mesophilic enzyme has increased stability and decreased activity, so that its properties are changed to resemble those of a thermophilic enzyme. Quasi-elastic neutron scattering measurements have also been performed to probe the protein motions. Compared to the free enzyme, the average height of the activation free energy barrier to local motions was found to be increased by 0.54 kcal.mol(-1) in the immobilized dihydrofolate reductase, a value that is of the same order as expected from the theoretical rate equation.     

Insights into the thermal stabilization and conformational transitions of DNA by hyperthermophile protein Sso7d: molecular dynamics simulations and MM-PBSA analysis

In the assembly of DNA-protein complex, the DNA kinking plays an important role in nucleoprotein structures and gene regulation. Molecular dynamics (MD) simulations were performed on specific protein-DNA complexes in this study to investigate the stability and structural transitions of DNA depending on temperature. Furthermore, we introduced the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) approach to analyze the interactions between DNA and protein in hyperthermophile. Focused on two specific Sso7d-DNA complexes (PDB codes: 1BNZ and 1BF4), we performed MD simulations at four temperatures (300, 360, 420, and 480?K) and MM-PBSA at 300 and 360?K to illustrate detailed information on the changes of DNA. Our results show that Sso7d stabilizes DNA duplex over a certain temperature range and DNA molecules undergo B-like to A-like form transitions in the binary complex with the temperature increasing, which are consistent with the experimental data. Our work will contribute to a better understanding of protein-DNA interaction.     

Molecular dynamics effects on protein electrostatics

Electrostatic calculations have been carried out on a number of structural conformers of tuna cytochrome c. Conformers were generated using molecular dynamics simulations with a range of solvent simulating, macroscopic dielectric formalisms, and one solvent model that explicitly included solvent water molecules. Structures generated using the lowest dielectric models were relatively tight, with side chains collapsed on the surface, while those from the higher dielectric models had more internal and external fluidity, with surface side chains exploring a fuller range of conformational space. The average structure generated with the explicitly solvated model corresponded most closely with the crystal structure. Individual pK values, overall titration curves, and electrostatic potential surfaces were calculated for average structures and structures along each simulation. Differences between structural conformers within each simulation give rise to substantial changes in calculated local electrostatic interactions, resulting in pK value fluctuations for individual sites in the protein that vary by 0.3-2.0 pK units from the calculated time average. These variations are due to the thermal side chain reorientations that produce fluctuations in charge site separations. Properties like overall titration curves and pH dependent stability are not as sensitive to side chain fluctuations within a simulation, but there are substantial effects between simulations due to marked differences in average side chain behavior. These findings underscore the importance of proper dielectric formalism in molecular dynamics simulations when used to generate alternate solution structures from a crystal structure, and suggest that conformers significantly removed from the average structure have altered electrostatic properties that may prove important in episodic protein properties such as catalysis.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

Protein flexibility: multiple molecular dynamics simulations of insulin chain B

Multiple molecular dynamics simulations totaling more than 100 ns were performed on chain B of insulin in explicit solvent at 300 K and 400 K. Despite some individual variations, a comparison of the protein dynamics of each simulation showed similar trends and most structures were consistent with NMR experimental values, even at the elevated temperature. The importance of packing interactions in determining the conformational transitions of the protein was observed, sometimes resulting in conformations induced by localized hydrophobic interactions. The high temperature simulation generated a more diverse range of structures with similar elements of secondary structure and populated conformations to the simulations at room temperature. A broad sampling of the conformational space of insulin chain B illustrated a wide range of conformational states with many transitions at room temperature in addition to the conformational states observed experimentally. The T-state conformation associated with insulin activity was consistently present and a possible mechanism of behavior was suggested.