Molecular dynamics study of a hyperthermophilic and a mesophilic rubredoxin

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.

Enhanced Thermostability and Enzymatic Activity of Cel6A Variants from Thermobifida fusca by Empirical Domain Engineering (Short Title: Domain Engineering of Cel6A)

Cellulases are a set of lignocellulolytic enzymes, capable of producing eco-friendly low-cost renewable bioethanol. However, low stability and hydrolytic activity limit their wide-scale applicability at the industrial scale. In this work, we report the domain engineering of endoglucanase (Cel6A) of Thermobifida fusca to improve their catalytic activity and thermal stability. Later, enzymatic activity and thermostability of the most efficient variant named as Cel6A.CBC was analyzed by molecular dynamics simulations. This variant demonstrated profound activity against soluble and insoluble cellulosic substrates like filter paper, alkali-treated bagasse, regenerated amorphous cellulose (RAC), and bacterial microcrystalline cellulose. The variant Cel6A.CBC showed the highest catalysis of carboxymethyl cellulose (CMC) and other related insoluble substrates at a pH of 6.0 and a temperature of 60 °C. Furthermore, a sound rationale was observed between experimental findings and molecular modeling of Cel6A.CBC which revealed thermostability of Cel6A.CBC at 26.85, 60.85, and 74.85 °C as well as structural flexibility at 126.85 °C. Therefore, a thermostable derivative of Cel6A engineered in the present work has enhanced biological performance and can be a useful construct for the mass production of bioethanol from plant biomass.

Heat-resistant cytosolic malate dehydrogenases (cMDHs) of thermophilic intertidal snails (genus Echinolittorina): protein underpinnings of tolerance to body temperatures reaching 55°C

Snails of the genus Echinolittorina are among the most heat-tolerant animals; they experience average body temperatures near 41-44°C in summer and withstand temperatures up to at least 55°C. Here, we demonstrate that heat stability of function (indexed by the Michaelis-Menten constant of the cofactor NADH, K_M^NADH) and structure (indexed by rate of denaturation) of cytosolic malate dehydrogenases (cMDHs) of two congeners (E. malaccana and E. radiata) exceeds values previously found for orthologs of this protein from less thermophilic species. The ortholog of E. malaccana is more heat stable than that of E. radiata, in keeping with the congeners' thermal environments. Only two inter-congener differences in amino acid sequence in these 332 residue proteins were identified. In both cases (positions 48 and 114), a glycine in the E. malaccana ortholog is replaced by a serine in the E. radiata protein. To explore the relationship between structure and function and to characterize how amino acid substitutions alter stability of different regions of the enzyme, we used molecular dynamics simulation methods. These computational methods allow determination of thermal effects on fine-scale movements of protein components, for example, by estimating the root mean square deviation in atom position over time and the root mean square fluctuation for individual residues. The minor changes in amino acid sequence favor temperature-adaptive change in flexibility of regions in and around the active sites. Interspecific differences in effects of temperature on fine-scale protein movements are consistent with the differences in thermal effects on binding and rates of heat denaturation.

Structural and dynamical insight into thermally induced functional inactivation of firefly luciferase

Luciferase is the key component of light production in bioluminescence process. Extensive and advantageous application of this enzyme in biotechnology is restricted due to its low thermal stability. Here we report the effect of heating up above Tm on the structure and dynamical properties of luciferase enzyme compared to temperature at 298 K. In this way we demonstrate that the number of hydrogen bonds between N- and C-domain is increased for the free enzyme at 325 K. Increased inter domain hydrogen bonds by three at 325 K suggests that inter domain contact is strengthened. The appearance of simultaneous strong salt bridge and hydrogen bond between K529 and D422 and increased existence probability between R533 and E389 could mechanistically explain stronger contact between N- and C-domain. Mutagenesis studies demonstrated the importance of K529 and D422 experimentally. Also the significant reduction in SASA for experimentally important residues K529, D422 and T343 which are involved in active site region was observed. Principle component analysis (PCA) in our study shows that the dynamical behavior of the enzyme is changed upon heating up which mainly originated from the change of motion modes and associated extent of those motions with respect to 298 K. These findings could explain why heating up of the enzyme or thermal fluctuation of protein conformation reduces luciferase activity in course of time as a possible mechanism of thermal functional inactivation. According to these results we proposed two strategies to improve thermal stability of functional luciferase.

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

A principal component analysis of the dynamics of subdomains and binding sites in human serum albumin

The conformational dynamics of human serum albumin (HSA) was investigated by principal component analysis (PCA) applied to three molecular dynamics trajectories of 200 ns each. The overlap of the essential subspaces spanned by the first 10 principal components (PC) of different trajectories was about 0.3 showing that the PCA based on a trajectory length of 200 ns is not completely convergent for this protein. The contributions of the relative motion of subdomains and of the subdomains (internal) distortion to the first 10 PCs were found to be comparable. Based on the distribution of the first 3 PC, 10 protein conformers are identified showing relative root mean square deviations (RMSD) between 2.3 and 4.6 Å. The main PCs are found to be delocalized over the whole protein structure indicating that the motions of different protein subdomains are coupled. This coupling is considered as being related to the allosteric effects observed upon ligand binding to HSA. On the other hand, the first PC of one of the three trajectories describes a conformational transition of the protein domain I that is close to that experimentally observed upon myristate binding. This is a theoretical support for the older hypothesis stating that changes of the protein onformation favorable to binding can precede the ligand complexation. A detailed all atoms PCA performed on the primary Sites 1 and 2 confirms the multiconformational character of the HSA binding sites as well as the significant coupling of their motions.

A single residue substitution accounts for the significant difference in thermostability between two isoforms of human cytosolic creatine kinase

Creatine kinase (CK) helps maintain homeostasis of intracellular ATP level by catalyzing the reversible phosphotransfer between ATP and phosphocreatine. In humans, there are two cytosolic CK isoforms, the muscle-type (M) and the brain-type (B), which frequently function as homodimers (hMMCK and hBBCK). Interestingly, these isoenzymes exhibit significantly different thermostabilities, despite high similarity in amino acid sequences and tertiary structures. In order to investigate the mechanism of this phenomenon, in this work, we first used domain swapping and site-directed mutagenesis to search for the key residues responsible for the isoenzyme-specific thermostability. Strikingly, the difference in thermostability was found to principally arise from one single residue substitution at position 36 (Pro in hBBCK vs. Leu in hMMCK). We then engaged the molecular dynamics simulations to study the molecular mechanism. The calculations imply that the P36L substitution introduces additional local interactions around residue 36 and thus further stabilizes the dimer interface through a complex interaction network, which rationalizes the observation that hMMCK is more resistant to thermal inactivation than hBBCK. We finally confirmed this molecular explanation through thermal inactivation assays on Asp36 mutants that were proposed to devastate the local interactions and thus the dimer associations in both isoenzymes.

A comparative molecular dynamics study of thermophilic and mesophilic β-fructosidase enzymes

Arabidopsis thaliana cell wall invertase 1 (AtcwINV1) and Thermotoga maritima β-fructosidase (BfrA) are among the best structurally studied members of the glycoside hydrolase family 32. Both enzymes hydrolyze sucrose as the main substrate but differ strongly in their thermal stability. Mesophilic AtcwINV1 and thermophilic BfrA have divergent sequence similarities in the N-terminal five bladed β-propeller catalytic domain (31 %) and the C-terminal β-sandwich domain (15 %) of unknown function. The two enzymes were subjected to 200 ns molecular dynamics simulations at 300 K (27 °C) and 353 K (80 °C). Regular secondary structure regions, but not loops, in AtcwINV1 and BfrA showed no significant fluctuation differences at both temperatures. BfrA was more rigid than AtcwINV1 at 300 K. The simulation at 353 K did not alter the structural stability of BfrA, but did increase the overall flexibility of AtcwINV1 exhibiting the most fluctuating regions in the β-propeller domain. The simulated heat treatment also increased the gyration radius and hydrophobic solvent accessible surface area of the plant enzyme, consistent with the initial steps of an unfolding process. The preservation of the conformational rigidity of BfrA at 353 K is linked to the shorter size of the protein loops. Shortening of BfrA loops appears to be a key mechanism for thermostability.

Adaptation of Extremophilic Proteins with Temperature and Pressure: Evidence from Initiation Factor 6

In this work, we study dynamical properties of an extremophilic protein, Initiation Factor 6 (IF6), produced by the archeabacterium Methanocaldococcus jannascii, which thrives close to deep-sea hydrothermal vents where temperatures reach 80 °C and the pressure is up to 750 bar. Molecular dynamics simulations (MD) and quasi-elastic neutron scattering (QENS) measurements give new insights into the dynamical properties of this protein with respect to its eukaryotic and mesophilic homologue. Results obtained by MD are supported by QENS data and are interpreted within the framework of a fractional Brownian dynamics model for the characterization of protein relaxation dynamics. IF6 from M. jannaschii at high temperature and pressure shares similar flexibility with its eukaryotic homologue from S. cerevisieae under ambient conditions. This work shows for the first time, to our knowledge, that the very common pattern of corresponding states for thermophilic protein adaptation can be extended to thermo-barophilic proteins. A detailed analysis of dynamic properties and of local structural fluctuations reveals a complex pattern for "corresponding" structural flexibilities. In particular, in the case of IF6, the latter seems to be strongly related to the entropic contribution given by an additional, C-terminal, 20 amino-acid tail which is evolutionary conserved in all mesophilic IF6s.

Structure and Dynamics of GeoCyp: A Thermophilic Cyclophilin with a Novel Substrate Binding Mechanism That Functions Efficiently at Low Temperatures

Thermophilic proteins have found extensive use in research and industrial applications because of their high stability and functionality at elevated temperatures while simultaneously providing valuable insight into our understanding of protein folding, stability, dynamics, and function. Cyclophilins, constituting a ubiquitously expressed family of peptidyl-prolyl isomerases with a range of biological functions and disease associations, have been utilized both for conferring stress tolerances and in exploring the link between conformational dynamics and enzymatic function. To date, however, no active thermophilic cyclophilin has been fully biophysically characterized. Here, we determine the structure of a thermophilic cyclophilin (GeoCyp) from Geobacillus kaustophilus, characterize its dynamic motions over several time scales using an array of methodologies that include chemical shift-based methods and relaxation experiments over a range of temperatures, and measure catalytic activity over a range of temperatures to compare its structure, dynamics, and function to those of a mesophilic counterpart, human cyclophilin A (CypA). Unlike those of most thermophile/mesophile pairs, GeoCyp catalysis is not substantially impaired at low temperatures as compared to that of CypA, retaining ~70% of the activity of its mesophilic counterpart. Examination of substrate-bound ensembles reveals a mechanism by which the two cyclophilins may have adapted to their environments through altering dynamic loop motions and a critical residue that acts as a clamp to regulate substrate binding differentially in CypA and GeoCyp. Fast time scale (pico- to nanosecond) dynamics are largely conserved between the two proteins, in accordance with the high degree of structural similarity, although differences do exist in their temperature dependencies. Slower (microsecond) time scale motions are likewise localized to similar regions in the two proteins with some variability in their magnitudes yet do not exhibit significant temperature dependencies in either enzyme.

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.

Structural reorganization renders enhanced metalloprotein stability

The enhanced stability of a mesophilic metalloprotein was assessed using biophysical spectroscopies. Significant local structural interconversions during thermal insult account for a reorganization of the protein scaffold, without disturbing the active metal site. This cushioning mechanism is proposed to be a generic property of metalloproteins contributing to enhanced stability.

Thermophilic proteins: insight and perspective from in silico experiments

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

Flexibility and mobility in mesophilic and thermophilic homologous proteins from molecular dynamics and FoldUnfold method

To function properly protein molecules require both flexibility and rigidity, therefore fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. In this work we used two theoretical approaches to determine flexible regions in four homologous pairs of proteins from thermophilic and mesophilic organisms. Protein pairs chosen in this study were selected to represent four typical folding classes. Our first approach, FoldUnfold, uses amino acid sequence and statistical information on the density of contacts of amino acids in tertiary structures of known globular proteins. The main advantages of such knowledge-based methodology are its computational speed and ability to make predictions in the absence of three-dimensional (3D) structure of a protein. The second approach uses a graph theory-based rigid cluster decomposition termed FIRST, applied together with Molecular Dynamics (MD) simulations of proteins with known structure. While MD simulations are time-consuming, they are the most direct way of studying physical properties of proteins, including their rigidity/flexibility. Flexible regions predicted by both methods in this work were in good agreement with each other. We also showed that high mobility of a site is not necessarily indicative of its high flexibility and vice versa. In our simulations thermophile proteins were less flexible than their mesophilic homologues. Longer flexible loops were found in mesophilic proteins of all classes.

Using molecular dynamics to probe the structural basis for enhanced stability in thermal stable cytochromes P450

High-temperature molecular dynamics (MD) has been used to assess if MD can be employed as a useful tool for probing the structural basis for enhanced stability in thermal stable cytochromes P450. CYP119, the most thermal stable P450 known, unfolds more slowly during 500 K MD simulations than P450s that melt at lower temperatures, P450cam and P450cin. A comparison of the 500 K MD trajectories shows that the Cys ligand loop, a critically important structural feature just under the heme, in both P450cin and P450cam completely unfolds while this region is quite stable in CYP119. In CYP119, this region is stabilized by tight nonpolar interactions involving Tyr26 and Leu308. The corresponding residues in P450cam are Gly and Thr, respectively. The in silico generated Y26A/L308A CYP119 double mutant is substantially less stable than wild-type CYP119, and the Cys ligand loop unfolds in a manner similar to that of P450cam. The MD thus has identified a potential "hot spot" important for stability. As an experimental test of the MD results, the Y26A/L308A double mutant was prepared, and thermal melting curves show that the double mutant exhibits a melting temperature (T(m)) 16 degrees C lower than that of wild-type CYP119. Control mutations that were predicted by MD not to destabilize the protein were also generated, and the experimental melting temperature was not significantly different from that of the wild-type enzyme. Therefore, high-temperature MD is a useful tool in predicting the structural underpinnings of thermal stability in P450s.

Structural study of carboxylesterase from hyperthermophilic bacteria Geobacillus stearothermophilus by molecular dynamics simulation

Carboxylesterases are ubiquitous enzymes with important physiological, industrial and medical applications such as synthesis and hydrolysis of stereo specific compounds, including the metabolic processing of drugs, and antimicrobial agents. Here, we have performed molecular dynamics simulations of carboxylesterase from hyperthermophilic bacterium Geobacillus stearothermophilus (GsEst) for 10ns each at five different temperatures namely at 300K, 343K, 373K, 473K and 500K. Profiles of root mean square fluctuation (RMSF) identify thermostable and thermosensitive regions of GsEst. Unfolding of GsEst initiates at the thermosensitive alpha-helices and proceeds to the thermostable beta-sheets. Five ion-pairs have been identified as critical ion-pairs for thermostability and are maintained stably throughout the higher temperature simulations. A detailed investigation of the active site residues of this enzyme suggests that the geometry of this site is well preserved up to 373K. Furthermore, the hydrogen bonds between Asp188 and His218 of the active site are stably maintained at higher temperatures imparting stability of this site. Radial distribution functions (RDFs) show similar pattern of solvent ordering and water penetration around active site residues up to 373K. Principal component analysis suggests that the motion of the entire protein as well as the active site is similar at 300K, 343K and 373K. Our study may help to identify the factors responsible for thermostability of GsEst that may endeavor to design enzymes with enhanced thermostability.

Hyperfine-shifted (13)C and (15)N NMR signals from Clostridium pasteurianum rubredoxin: extensive assignments and quantum chemical verification

Stable isotope-labeling methods, coupled with novel techniques for detecting fast-relaxing NMR signals, now permit detailed investigations of paramagnetic centers of metalloproteins. We have utilized these advances to carry out comprehensive assignments of the hyperfine-shifted ¹³C and ¹⁵N signals of the rubredoxin from Clostridium pasteurianum (CpRd) in both its oxidized and reduced states. We used residue-specific labeling (by chemical synthesis) and residue-type-selective labeling (by biosynthesis) to assign signals detected by one-dimensional ¹⁵N NMR spectroscopy, to nitrogen atoms near the iron center. We refined and extended these ¹⁵N assignments to the adjacent carbonyl carbons by means of one-dimensional ¹³C[¹⁵N] decoupling difference experiments. We collected paramagnetic-optimized SuperWEFT ¹³C[¹³C] constant time COSY (SW-CT-COSY) data to complete the assignment of ¹³C signals of reduced CpRd. By following these ¹³C signals as the protein was gradually oxidized, we transferred these assignments to carbons in the oxidized state. We have compared these assignments with hyperfine chemical shifts calculated from available X-ray structures of CpRd in its oxidized and reduced forms. The results allow the evaluation of the X-ray structural models as representative of the solution structure of the protein, and they provide a framework for future investigation of the active site of this protein. The methods developed here should be applicable to other proteins that contain a paramagnetic center with high spin and slow electron exchange.

Dynamics in Thermotoga neapolitana adenylate kinase: 15N relaxation and hydrogen-deuterium exchange studies of a hyperthermophilic enzyme highly active at 30 degrees C

Backbone conformational dynamics of Thermotoga neapolitana adenylate kinase in the free form (TNAK) and inhibitor-bound form (TNAK*Ap5A) were investigated at 30 degrees C using (15)N NMR relaxation measurements and NMR monitored hydrogen-deuterium exchange. With kinetic parameters identical to those of Escherichia coli AK (ECAK) at 30 degrees C, TNAK is a unique hyperthermophilic enzyme. These catalytic properties make TNAK an interesting and novel model to study the interplay between protein rigidity, stability, and activity. Comparison of fast time scale dynamics (picosecond to nanosecond) in the open and closed states of TNAK and ECAK at 30 degrees C reveals a uniformly higher rigidity across all domains of TNAK. Within this framework of a rigid TNAK structure, several residues located in the AMP-binding domain and in the core-lid hinge regions display high picosecond to nanosecond time scale flexibility. Together with the recent comparison of ECAK dynamics with those of hyperthermophilic Aquifex aeolicus AK (AAAK), our results provide strong evidence for the role of picosecond to nanosecond time scale fluctuations in both stability and activity. In the slow time scales, TNAK's increased rigidity is not uniform but localized in the AMP-binding and lid domains. The core domain amides of ECAK and TNAK in the open and closed states show comparable protection against exchange. Significantly, the hinges framing the lid domain show similar exchange data in ECAK and TNAK open and closed forms. Our NMR relaxation and hydrogen-deuterium exchange studies therefore suggest that TNAK maintains high activity at 30 degrees C by localizing flexibility to the hinge regions that are key to facilitating conformational changes.

Physical basis of metal-binding specificity in Escherichia coli NikR

In Escherichia coli and other bacteria, nickel uptake is regulated by the transcription factor NikR. Nickel binding at high-affinity sites in E. coli NikR (EcNikR) facilitates EcNikR binding to the nik operon, where it then suppresses transcription of genes encoding the nickel uptake transporter, NikABCDE. A structure of the EcNikR-DNA complex suggests that a second metal-binding site is also present when NikR binds to the nik operon. Moreover, this co-crystal structure raises the question of what metal occupies the second site under physiological conditions: K(+), which is present in the crystal structure, or Ni(2+), which has been proposed to bind to low- as well as high-affinity sites on EcNikR. To determine which ion is preferred at the second metal-binding site and the physical basis for any preference of one ion over another in both the second metal-binding site and the high-affinity sites, we conducted a series of detailed molecular simulations on the EcNikR structure. Simulations that place Ni(2+) at high-affinity sites lead to stable trajectories with realistic ion-ligand distances and geometries, while simulations that place K(+) at these sites lead to conformational changes in the protein that are likely unfavorable for ion binding. By contrast, simulations on the second metal site in the EcNikR-DNA complex lead to stable trajectories with realistic geometries regardless of whether K(+) or Ni(2+) occupies this site. Electrostatic binding free energy calculations, however, suggest that EcNikR binding to DNA is more favorable when the second metal-binding site contains K(+). An analysis of the energetic contributions to the electrostatic binding free energy suggests that, while the interaction between EcNikR and DNA is more favorable when the second site contains Ni(2+), the large desolvation penalty associated with moving Ni(2+) from solution to the relatively buried second site offsets this favorable interaction term. Additional free energy simulations that account for both electrostatic and non-electrostatic effects argue that EcNikR binding to DNA is most favorable when the second site contains a monovalent ion the size of K(+). Taken together, these data suggest that the EcNikR structure is most stable when Ni(2+) occupies high-affinity sites and that EcNikR binding to DNA is more favorable when the second site contains K(+).

Insights into thermal stability of thermophilic nitrile hydratases by molecular dynamics simulation

Thermal stability is of great importance for industrial enzymes. Here we explored the thermal-stable mechanism of thermophilic nitrile hydratases (NHases) utilizing a molecular dynamic simulation. At a nanosecond timescale, profiles of root mean square fluctuation (RMSF) of two thermophilic NHases, 1UGQ and 1V29, under enhancing thermal stress were carried out at 300 K, 320 K, 350 K and 370 K, respectively. Results showed that the region A1 (211-231 aa) and A2 (305-316 aa) in 1UGQ, region B1 (186-192 aa) in 1V29, and most of terminal ends in both enzymes are hyper-sensitive. Salt-bridge analyses revealed that in one hand, salt-bridges contributed to maintaining the rigid structure and stable performance of the thermophilic 1UGQ and 1V29; in the other hand, salt-bridges involved in thermal sensitive regions are relatively weak and prone to be broken at elevated temperature, thereby cannot hold the stable conformation of the spatial neighborhood. In 1V29, region A1 was stabilized by a well-organized hook-hook like cluster with multiple salt-bridge interactions, region A2 was stabilized by two strong salt-bridge interactions of GLU52-ARG332 and GLU334-ARG332. In 1UGQ, the absence of a charged residue decreased its thermal sensitivity of region B1, and the formation of a small beta-sheet containing a stable salt-bridge in C-beta-terminal significantly enhanced its thermal stability. By radius of gyration calculation containing or eliminating the thermal sensitive regions, we quantified the contribution of thermal sensitive regions for thermal sensitivity of 1UGQ and 1V29. Consequently, we presented strategies to improve thermal stability of the industrialized mesophilic NHase by introducing stable salt-bridge interactions into its thermal sensitive regions.

Cavity-creating mutations in Pseudomonas aeruginosa azurin: effects on protein dynamics and stability

Changes in flexibility and structural stability of Pseudomonas aeruginosa azurin in response to cavity-creating mutations were probed by the phosphorescence emission of Trp-48, which was deeply buried in the compact hydrophobic core of the macromolecule, and by measurements of guanidinum hydrochloride unfolding, respectively. Replacement of the bulky side chains Phe-110, Phe-29, and Tyr-108 with the smaller Ala introduced cavities at different distances from the hydrophobic core. The phosphorescence lifetime (tau(0)) of Trp-48, buried inside the protein core, and the acrylamide quenching rate constant (k(q)) were used to monitor local and global flexibility changes induced by the introduction of the cavity. The results of this work demonstrate the following: 1), the effect on core flexibility of the insertion of cavities is not correlated readily to the distance of the cavity from the core; 2), the protein global flexibility results are related to the cavity distance from the packed core of the macromolecule; and 3), the increase in protein flexibility does not correspond necessarily to a comparable destabilizing effect of some mutations.

Dispersion interactions govern the strong thermal stability of a protein

Rubredoxin from the hyperthermophile Pyrococcus furiosus (Pf Rd) is an extremely thermostable protein, which makes it an attractive subject of protein folding and stability studies. A fundamental question arises as to what the reason for such extreme stability is and how it can be elucidated from a complex set of interatomic interactions. We addressed this issue first theoretically through a computational analysis of the hydrophobic core of the protein and its mutants, including the interactions taking place inside the core. Here we show that a single mutation of one of phenylalanine's residues inside the protein's hydrophobic core results in a dramatic decrease in its thermal stability. The calculated unfolding Gibbs energy as well as the stabilization energy differences between a few core residues follows the same trend as the melting temperature of protein variants determined experimentally by microcalorimetry measurements. NMR spectroscopy experiments have shown that the only part of the protein affected by mutation is the reasonably rearranged hydrophobic core. It is hence concluded that stabilization energies, which are dominated by London dispersion, represent the main source of stability of this protein.

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.

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.

Observation of terahertz vibrations in Pyrococcus furiosus rubredoxin via impulsive coherent vibrational spectroscopy and nuclear resonance vibrational spectroscopy – interpretation by molecular mechanics

We have used impulsive coherent vibrational spectroscopy (ICVS) to study the Fe(S-Cys)(4) site in oxidized rubredoxin (Rd) from Pyrococcus furiosus (Pf). In this experiment, a 15 fs visible laser pulse is used to coherently pump the sample to an excited electronic state, and a second <10 fs pulse is used to probe the change in transmission as a function of the time delay. PfRd was observed to relax to the ground state by a single exponential decay with time constants of approximately 255-275 fs. Superimposed on this relaxation are oscillations caused by coherent excitation of vibrational modes in both excited and ground electronic states. Fourier transformation reveals the frequencies of these modes. The strongest ICV mode with 570 nm excitation is the symmetric Fe-S stretching mode near 310 cm(-1), compared to 313 cm(-1) in the low temperature resonance Raman. If the rubredoxin is pumped at 520 nm, a set of strong bands occurs between 20 and 110 cm(-1). Finally, there is a mode at approximately 500 cm(-1) which is similar to features near 508 cm(-1) in blue Cu proteins that have been attributed to excited state vibrations. Normal mode analysis using 488 protein atoms and 558 waters gave calculated spectra that are in good agreement with previous nuclear resonance vibrational spectra (NRVS) results. The lowest frequency normal modes are identified as collective motions of the entire protein or large segments of polypeptide. Motion in these modes may affect the polar environment of the redox site and thus tune the electron transfer functions in rubredoxins.

Cold-active enzymes studied by comparative molecular dynamics simulation

Enzymes from cold-adapted species are significantly more active at low temperatures, even those close to zero Celsius, but the rationale of this adaptation is complex and relatively poorly understood. It is commonly stated that there is a relationship between the flexibility of an enzyme and its catalytic activity at low temperature. This paper gives the results of a study using molecular dynamics simulations performed for five pairs of enzymes, each pair comprising a cold-active enzyme plus its mesophilic or thermophilic counterpart. The enzyme pairs included alpha-amylase, citrate synthase, malate dehydrogenase, alkaline protease and xylanase. Numerous sites with elevated flexibility were observed in all enzymes; however, differences in flexibilities were not striking. Nevertheless, amino acid residues common in both enzymes of a pair (not present in insertions of a structure alignment) are generally more flexible in the cold-active enzymes. The further application of principle component analysis to the protein dynamics revealed that there are differences in the rate and/or extent of opening and closing of the active sites. The results indicate that protein dynamics play an important role in catalytic processes where structural rearrangements, such as those required for active site access by substrate, are involved. They also support the notion that cold adaptation may have evolved by selective changes in regions of enzyme structure rather than in global change to the whole protein.

The effects of cosolutes on protein dynamics: the reversal of denaturant-induced protein fluctuations by trimethylamine N-oxide

The protein stabilizing effects of the small molecule osmolyte, trimethylamine N-oxide, against chemical denaturant was investigated by NMR spin-relaxation measurements and model-free analysis. In the presence of 0.7 M guanidine hydrochloride increased picosecond-nanosecond dynamics are observed in the protein ribonuclease A. These increased fluctuations occur throughout the protein, but the most significant increases in flexibility occur at positions believed to be the first to unfold. Addition of 0.35 M trimethylamine N-oxide to this destabilized form of ribonuclease results in significant rigidification of the protein backbone as assessed by (1)H-(15)N order parameters. Statistically, these order parameters are the same as those measured in native ribonuclease indicating that TMAO reduces the amplitude of backbone fluctuations in a destabilized protein. These data suggest that TMAO restricts the bond vector motions on the protein energy landscape to resemble those motions that occur in the native protein and points to a relation between stability and dynamics in this enzyme.

Molecular Dynamics of Apo-Adenylate Kinase: A Principal Component Analysis

Adenylate kinase from E. coli (AKE) is studied with molecular dynamics. AKE undergoes large-scale motions of its Lid and AMP-binding domains when its open form closes over its substrates, AMP and Mg2+-ATP. The third domain, the Core, is relatively stable during closing. The resulting trajectory is analyzed with a principal component analysis method that decomposes the atom motions into modes ordered by their decreasing contributions to the total protein fluctuation. Simulations at 303 K (normal T) and 500 K (high T) reveal that at both temperatures the first three modes account for 70% of the total fluctuation. The residues that contribute the most to these three modes are concentrated in the Lid and AMP-binding domains. Analysis of the normal T modes indicates that the Lid and AMP-binding domains sample a broad distribution of conformations indicating that AKE is designed to provide its substrates with a large set of conformations. The high T results show that the Lid initially closes toward the Core. Subsequently, the Lid rotates to a new stable conformation that is different from what is observed in the substrate-bound AKE. These results are discussed in the context of experimental data that indicate that adenylate kinases do sample more than one conformational state in solution and that each of these conformational states undergoes substantial fluctuations. A pair of residues is suggested for labeling that would be useful for monitoring distance fluctuations by energy transfer experiments.

Explanation of the stability of thermophilic proteins based on unique micromorphology

Two mesophilic/thermophilic variants of the G-domain of the elongation factor Tu were studied via molecular dynamics simulations. By analyzing the simulation data via the Voronoi space tessellation, we have found that the two proteins have the same macromolecular packing, while the water-exposed surface area is larger for the thermophile. A larger coordination with water is probably due to a peculiar corrugation of the exposed surface of this species. From an enthalpic point of view, the thermophile shows a larger number of intramolecular hydrogen bonds, stronger electrostatic interactions, and a flatter free-energy landscape. Overall, the data suggest that the specific hydration state enhances macromolecular fluctuations but, at the same time, increases thermal stability.

Neutron Scattering Reveals the Dynamic Basis of Protein Adaptation to Extreme Temperature

To explore protein adaptation to extremely high temperatures, two parameters related to macromolecular dynamics, the mean square atomic fluctuation and structural resilience, expressed as a mean force constant, were measured by neutron scattering for hyperthermophilic malate dehydrogenase from Methanococcus jannaschii and a mesophilic homologue, lactate dehydrogenase from Oryctolagus cunniculus (rabbit) muscle. The root mean square fluctuations, defining flexibility, were found to be similar for both enzymes (1.5 A) at their optimal activity temperature. Resilience values, defining structural rigidity, are higher by an order of magnitude for the high temperature-adapted protein (0.15 Newtons/meter for O. cunniculus lactate dehydrogenase and 1.5 Newtons/meter for M. jannaschii malate dehydrogenase). Thermoadaptation appears to have been achieved by evolution through selection of appropriate structural rigidity in order to preserve specific protein structure while allowing the conformational flexibility required for activity.

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.

Identifying and engineering ion pairs in adenylate kinases. Insights from molecular dynamics simulations of thermophilic and mesophilic homologues

Molecular dynamics simulations were performed to study thermal stabilization of proteins via electrostatic interactions of ion pairs. Dynamic motions of four ion pairs previously proposed to be important in thermal stability of adenylate kinase from the thermophile Bacillus stearothermophilus were monitored during the simulation. One of the four ion pairs identified in the crystal structure, Lys180-Asp114, was not maintained in close contact suggesting that the ion pair does not contribute to thermal stability. Among the other three ion pairs, the ion pair Arg116-Glu198 was proposed to be the most important for stability. To predict behaviors of the ion pairs when engineered into a mesophilic homologue to increase stability, in silico mutants of adenylate kinase from the mesophile Bacillus subtilis were generated, and their molecular dynamics simulations were carried out. The ion pairs in the mutant simulations displayed similar behaviors to those in the simulation of the thermophilic protein. To validate the results of the simulations experimentally, the same mutants were produced in vitro and their thermal stabilities were measured using differential scanning calorimetry. In agreement with the simulations, the Lys180-Asp114 did not result in any increase in stability by itself or additive effect with other ion pairs, whereas a mutant with the Arg116-Glu198 exhibited the highest stability among the mutants having one of the four ion pairs. These results provide specific knowledge about stability in adenylate kinases and more generally suggest that molecular dynamics simulations can provide valuable information for identifying and engineering ion pairs.

The Importance of Dynamic Effects on the Enzyme Activity

Onconase (ONC), a member of the RNase A superfamily extracted from oocytes of Rana pipiens, is an effective cancer killer. It is currently used in treatment of various forms of cancer. ONC antitumor properties depend on its ribonucleolytic activity that is low in comparison with other members of the superfamily. The most damaging side effect from Onconase treatment is renal toxicity, which seems to be caused by the unusual stability of the enzyme. Therefore, mutants with reduced thermal stability and/or increased catalytic activity may have significant implications for human cancer chemotherapy. In this context, we have determined the crystal structures of two Onconase mutants (M23L-ONC and C87S,des103-104-ONC) and performed molecular dynamic simulations of ONC and C87S,des103-104-ONC with the aim of explaining on structural grounds the modifications of the activity and thermal stability of the mutants. The results also provide the molecular bases to explain the lower catalytic activity of Onconase compared with RNase A and the unusually high thermal stability of the amphibian enzyme.

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.

Checking the pH-induced conformational transition of prion protein by molecular dynamics simulations: effect of protonation of histidine residues

The role of acidic pH in the conversion of human prion protein to the pathogenic isoform is investigated by means of molecular dynamics simulations, focusing the attention on the effect of protonation of histidine residues on the conformational behavior of human PrPC globular domain. Our simulations reveal a significant loss of alpha-helix content under mildly acidic conditions, due to destructuration of the C-terminal part of HB (thus suggesting a possible involvement of HB into the conformational transition leading to the pathogenic isoform) and a transient lengthening of the native beta-sheet. Protonation of His-187 and His-155 seems to be crucial for the onset of the conformational rearrangement. This finding can be related to the existence of a pathogenic mutation, H187R, which is associated with GSS syndrome. Finally, the relevance of our results for the location of a Cu2+-binding pocket in the C-terminal part of the prion is discussed.

Multiple time scale backbone dynamics of homologous thermophilic and mesophilic ribonuclease HI enzymes

Backbone conformational fluctuations on multiple time scales in a cysteine-free Thermus thermophilus ribonuclease HI mutant (ttRNH(*)) are quantified using (15)N nuclear magnetic spin relaxation. Laboratory-frame relaxation data acquired at 310 K and at static magnetic field strengths of 11.7, 14.1 and 18.8 T are analysed using reduced spectral density mapping and model-free approaches. Chemical exchange line broadening is characterized using Hahn-echo transverse and multiple quantum relaxation data acquired over a temperature range of 290-320 K and at a static magnetic field strength of 14.1 T. Results for ttRNH(*) are compared to previously published data for a mesophilic homologue, Escherichia coli ribonuclease HI (ecRNH). Intramolecular conformational fluctuations on the picosecond-to-nanosecond time scale generally are similar for ttRNH(*) and ecRNH. beta-Strands 3 and 5 and the glycine-rich region are more rigid while the substrate-binding handle region and C-terminal tail are more flexible in ttRNH(*) than in ecRNH. Rigidity in the two beta-strands and the glycine-rich region, located along the periphery of the central beta-sheet, may be associated with the increased thermodynamic stability of the thermophilic enzyme. Chemical exchange line broadening, reflecting microsecond-to-millisecond time scale conformational changes, is more pronounced in ttRNH(*) than in ecRNH, particularly for residues in the handle and surrounding the catalytic site. The temperature dependence of chemical exchange show an increase of approximately 15 kJ/mol in the apparent activation energies for ttRNH(*) residues in the handle compared to ecRNH. Increased activation barriers, coupled with motion between alpha-helices B and C not present in ecRNH, may be associated with the reduced catalytic activity of the thermophilic enzyme at 310 K.

Structural and dynamic effects of α-Helix deletion in Sso7d: Implications for protein thermal stability

Sso7d is a 62-residue protein from the hyperthemophilic archaeon Sulfolobus solfataricus with a denaturation temperature close to 100 degrees C around neutral pH. An engineered form of Sso7d truncated at leucine 54 (L54Delta) is significantly less stable, with a denaturation temperature of 53 degrees C. Molecular dynamics (MD) studies of Sso7d and its truncated form at two different temperatures have been performed. The results of the MD simulations at 300 K indicate that: (1) the flexibility of Sso7d chain at 300 K agrees with that detected from X-ray and NMR structural studies; (2) L54Delta remains stable in the native folded conformation and possesses an overall dynamic behavior similar to that of the parent protein. MD simulations performed at 500 K, 10 ns long, indicate that, while Sso7d is in-silico resistant to high temperature, the truncated variant partially unfolds, revealing the early phases of the thermal unfolding pathway of the protein. Analysis of the trajectories of L54Delta suggests that the unzipping of the N-terminal and C-terminal beta-strands should be the first event of the unfolding pathway, and points out the regions more resistant to thermal unfolding. These findings allow one to understand the role played by specific interactions connecting the two ends of the chain for the high thermal stability of Sso7d, and support recent hypotheses on its folding mechanism emerged from site-directed mutagenesis studies.

Molecular dynamics simulations of metalloproteins

Molecular dynamics simulations are now commonly applied to metalloproteins, despite the challenges introduced by the presence of metal ions. Force field parameters are nowadays available also for these 'exotic' atoms and several biological systems have been successfully studied. Some of the most relevant results and methodological advancements are reviewed.

Molecular dynamics simulations of a K+ channel blocker: Tc1 toxin from Tityus cambridgei

Toxins that block voltage-gated potassium (Kv) channels provide a possible template for improved homology models of the Kv pore. In assessing the interactions of Kv channels and their toxins it is important to determine the dynamic flexibility of the toxins. Multiple 10 ns duration molecular dynamics simulations combined with essential dynamics analysis have been used to explore the flexibility of four different Kv channel-blocking toxins. Three toxins (Tc1, AgTx and ChTx) share a common fold. They also share a common pattern of conformational dynamics, as revealed by essential dynamics analysis of the simulation results. This suggests that some aspects of dynamic behaviour are conserved across a single protein fold class. In each of these three toxins, the residue exhibiting minimum flexibility corresponds to a conserved lysine residue that is suggested to interact with the filter domain of the channel. Thus, comparative simulations reveal functionally important conservation of molecular dynamics as well as protein fold across a family of related toxins.