Molecular dynamics simulations of the hyperthermophilic protein sac7d from Sulfolobus acidocaldarius: contribution of salt bridges to thermostability

Alteration of Chain-Length Selectivity and Thermostability of Rhizopus oryzae Lipase via Virtual Saturation Mutagenesis Coupled with Disulfide Bond Design

Rhizopus oryzae lipase (ROL) is one of the most important enzymes used in the food, biofuel, and pharmaceutical industries. However, the highly demanding conditions of industrial processes can reduce its stability and activity. To seek a feasible method to improve both the catalytic activity and the thermostability of this lipase, first, the structure of ROL was divided into catalytic and noncatalytic regions by identifying critical amino acids in the crevice-like binding pocket. Second, a mutant screening library aimed at improvement of ROL catalytic performance by virtual saturation mutagenesis of residues in the catalytic region was constructed based on Rosetta's Cartesian_ddg protocol. A double mutant, E265V/S267W (with an E-to-V change at residue 265 and an S-to-W change at residue 267), with markedly improved catalytic activity toward diverse chain-length fatty acid esters was identified. Then, computational design of disulfide bonds was conducted for the noncatalytic amino acids of E265V/S267W, and two potential disulfide bonds, S61C-S115C and E190C-E238C, were identified as candidates. Experimental data validated that the variant E265V/S267W/S61C-S115C/E190C-E238C had superior stability, with an increase of 8.5°C in the melting temperature and a half-life of 31.7 min at 60°C, 4.2-fold longer than that of the wild-type enzyme. Moreover, the variant improved the lipase activity toward five 4-nitrophenyl esters by 1.5 to 3.8 times, exhibiting a potential to modify the catalytic efficiency. IMPORTANCE Rhizopus oryzae lipase (ROL) is very attractive in biotechnology and industry as a safe and environmentally friendly biocatalyst. Functional expression of ROL in Escherichia coli facilitates effective high-throughput screening for positive variants. This work highlights a method to improve both selectivity and thermostability based on a combination of virtual saturation mutagenesis in the substrate pocket and disulfide bond prediction in the noncatalytic region. Using the method, ROL thermostability and activity to diverse 4-nitrophenyl esters could be substantially improved. The strategy of rational introduction of multiple mutations in different functional domains of the enzyme is a great prospect in the modification of biocatalysts.

Modulation of Flexible Loops in Catalytic Cavities Reveals the Thermal Activation Mechanism of a Glycogen-Debranching Enzyme

Some thermophilic enzymes not only exhibit high thermostability at high temperatures but also have an activation effect by thermal incubation. However, the correlations between temperature-induced structural modulation and thermal activation are still unclear. In this study, we selected a thermophilic glycogen-debranching enzyme from Saccharolobus solfataricus STB09 (SsGDE), which was a promising starch-debranching enzyme with a thermal activation property at temperatures ranging from 50 to 70 °C, to explore the thermal activation mechanism. Molecular dynamics simulations were performed for SsGDE at 30, 50, or 70 °C to reveal the temperature dependence of structure modulation and catalytic function. The results revealed that four loops (loop1 313-337, loop2 399-418, loop3 481-513, and loop4 540-574) in SsGDE were reshaped, which made the catalytic cavity more open. The internal residues, including the catalytic triad Asp3631, Glu399, and Asp471, could be exposed, due to the structural modulation, to exert catalytic functions. We proposed that the thermal activation effect of SsGDE was closely associated with the temperature-induced modulation of the catalytic cavity, which paved the way for further engineering enzymes to achieve higher catalytic performance and stability.

Overview of strategies for developing high thermostability industrial enzymes: Discovery, mechanism, modification and challenges

Biocatalysts such as enzymes are environmentally friendly and have substrate specificity, which are preferred in the production of various industrial products. However, the strict reaction conditions in industry including high temperature, organic solvents, strong acids and bases and other harsh environments often destabilize enzymes, and thus substantially compromise their catalytic functions, and greatly restrict their applications in food, pharmaceutical, textile, bio-refining and feed industries. Therefore, developing industrial enzymes with high thermostability becomes very important in industry as thermozymes have more advantages under high temperature. Discovering new thermostable enzymes using genome sequencing, metagenomics and sample isolation from extreme environments, or performing molecular modification of the existing enzymes with poor thermostability using emerging protein engineering technology have become an effective means of obtaining thermozymes. Based on the thermozymes as biocatalytic chips in industry, this review systematically analyzes the ways to discover thermostable enzymes from extreme environment, clarifies various interaction forces that will affect thermal stability of enzymes, and proposes different strategies to improve enzymes' thermostability. Furthermore, latest development in the thermal stability modification of industrial enzymes through rational design strategies is comprehensively introduced from structure-activity relationship point of view. Challenges and future research perspectives are put forward as well.

Electric Field-Mediated Fibronectin-Hydroxyapatite Interaction: A Molecular Insight

In experimental research-driven biomaterials science, the influence of different material properties (elastic stiffness, surface energy, etc.) and, to a relatively lesser extent, biophysical stimulation (electric/magnetic) on cell-material interactions has been extensively investigated. Despite the central importance of protein adsorption on cell-material interactions, the quantitative analysis to probe into the role of physicochemical factors in protein adsorption remains largely unexplored in biomaterials science. In recent studies, the critical role of electric field stimulation toward the modulation of cell functionality in implantable biomaterials has been experimentally demonstrated. Given this background, we investigated the influence of external electric field stimulation (upto 1.00 V/nm) on fibronectin (FN) adsorption on a hydroxyapatite (HA) (001) surface at 300 K using the all-atom molecular dynamics (MD) simulation method. FN adsorption was found to be governed by attractive electrostatic interactions, which changed with the electric field strength. Nonmonotonous changes in the structural integrity of FN were recorded with the change in the field strength and direction. This can be attributed to the spatial rearrangement of the positions of local charges and the global structural changes of proteins. The dipole moment vectors of FN, water, and HA quantitatively exhibited a similar pattern of orienting themselves parallel to the field direction, with field strength-dependent increase in their magnitudes. No significant change has been recorded for the radial distribution function of water surrounding FN. Field-dependent variation in the salt bridge nets and the number of hydrogen bonds between FN and HA were also examined. One of the important results in the context of cell-material interaction is that the RGD (Arg-Gly-Asp) sequence of FN was exposed to the solvent side when the field was applied along an outward direction perpendicular to the HA (001) surface. In summary, the present study provides molecular insights into the influence of electric field stimulation on phenomenological interactions involved in FN adsorption on the HA surface.

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

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

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.

Peptide Bond Isomerization in High-Temperature Simulations

Force fields for molecular simulation are generally optimized to model macromolecules such as proteins at ambient temperature and pressure. Nevertheless, elevated temperatures are frequently used to enhance conformational sampling, either during system setup or as a component of an advanced sampling technique such as temperature replica exchange. Because macromolecular force fields are now put upon to simulate temperatures and time scales that greatly exceed their original design specifications, it is appropriate to re-evaluate whether these force fields are up to the task. Here, we quantify the rates of peptide bond isomerization in high-temperature simulations of three octameric peptides and a small fast-folding protein. We show that peptide octamers with and without proline residues undergo cis/trans isomerization every 1-5 ns at 800 K with three classical atomistic force fields (AMBER99SB-ILDN, CHARMM22/CMAP, and OPLS-AA/L). On the low microsecond time scale, these force fields permit isomerization of nonprolyl peptide bonds at temperatures ≥500 K, and the CHARMM22/CMAP force field permits isomerization of prolyl peptide bonds ≥400 K. Moreover, the OPLS-AA/L force field allows chiral inversion about the Cα atom at 800 K. Finally, we show that temperature replica exchange permits cis peptide bonds developed at 540 K to subsequently migrate back to the 300 K ensemble, where cis peptide bonds are present in 2 ± 1% of the population of Trp-cage TC5b, including up to 4% of its folded state. Further work is required to assess the accuracy of cis/trans isomerization in the current generation of protein force fields.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

Thermodynamics of protein denaturation at temperatures over 100 °C: CutA1 mutant proteins substituted with hydrophobic and charged residues

Although the thermodynamics of protein denaturation at temperatures over 100 °C is essential for the rational design of highly stable proteins, it is not understood well because of the associated technical difficulties. We designed certain hydrophobic mutant proteins of CutA1 from Escherichia coli, which have denaturation temperatures (Td) ranging from 101 to 113 °C and show a reversible heat denaturation. Using a hydrophobic mutant as a template, we successfully designed a hyperthermostable mutant protein (Td = 137 °C) by substituting six residues with charged ones. Thermodynamic analyses of these mutant proteins indicated that the hydrophobic mutants were stabilized by the accumulation of denaturation enthalpy (ΔH) with no entropic gain from hydrophobic solvation around 100 °C, and that the stabilization due to salt bridges resulted from both the increase in ΔH from ion-ion interactions and the entropic effect of the electrostatic solvation over 113 °C. This is the first experimental evidence that has successfully overcome the typical technical difficulties.

Physicochemical Properties of Ion Pairs of Biological Macromolecules

Ion pairs (also known as salt bridges) of electrostatically interacting cationic and anionic moieties are important for proteins and nucleic acids to perform their function. Although numerous three-dimensional structures show ion pairs at functionally important sites of biological macromolecules and their complexes, the physicochemical properties of the ion pairs are not well understood. Crystal structures typically show a single state for each ion pair. However, recent studies have revealed the dynamic nature of the ion pairs of the biological macromolecules. Biomolecular ion pairs undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. This dynamic behavior is reasonable in light of the fundamental concepts that were established for small ions over the last century. In this review, we introduce the physicochemical concepts relevant to the ion pairs and provide an overview of the recent advancement in biophysical research on the ion pairs of biological macromolecules.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

Insights into the unfolding pathway and identification of thermally sensitive regions of phytase from Aspergillus niger by molecular dynamics simulations

Thermal stability is of great importance in the application of commercial phytases. Phytase A (PhyA) is a monomeric protein comprising twelve α-helices and ten β-sheets. Comparative molecular dynamics (MD) simulations (at 310, 350, 400, and 500 K) revealed that the thermal stability of PhyA from Aspergillus niger (A. niger) is associated with its conformational rigidity. The most thermally sensitive regions were identified as loops 8 (residues 83-106), 10 (161-174), 14 (224-230), 17 (306-331), and 24 (442-444), which are present on the surface of the protein. It was observed that solvent-exposed loops denature before or show higher flexibility than buried residues. We observed that PhyA begins to unfold at loops 8 and 14, which further extends to loop 24 at the C-terminus. The intense movement of loop 8 causes the helix H2 and beta-sheet B3 to fluctuate at high temperature. The high flexibility of the H2, H10, and H12 helices at high temperature resulted in complete denaturation. The high mobility of loop 14 easily transfers to the adjacent helices H7, H8, and H9, which fluctuate and partially unfold at high temperature (500 K). It was also observed that the salt bridges Asp110-Lys149, Asp205-Lys277, Asp335-Arg136, Asp416-Arg420, and Glu387-Arg400 are important influences on the structural stability but not the thermostability, as the lengths of these salt bridges did not increase with rising temperature. The salt bridges Glu125-Arg163, Asp299-Arg136, Asp266-Arg219, Asp339-Lys278, Asp335-Arg136, and Asp424-Arg428 are all important for thermostability, as the lengths of these bridges increased dramatically with increasing temperature. Here, for the first time, we have computationally identified the thermolabile regions of PhyA, and this information could be used to engineer novel thermostable phytases. Numerous homologous phytases of fungal as well as bacterial origin are known, and these homologs show high sequence similarity. Our findings could prove useful in attempts to increase the thermostability of homologous phytases via protein engineering.

The role of proline substitutions within flexible regions on thermostability of luciferase

Improving the stability of firefly luciferase has been a critical issue for its wider industrial applications. Studies about hyperthermophile proteins show that flexibility could be an effective indicator to find out weak spots to engineering thermostability of proteins. However, the relationship among flexibility, activity and stability in most of proteins is unclear. Proline is the most rigid residue and can be introduced to rigidify flexible regions to enhance thermostability of proteins. We firstly apply three different methods, molecular dynamics (MD) simulation, B-FITTER and framework rigidity optimized dynamics algorithm (FRODA) to determine the flexible regions of Photinus pyralis luciferase: Fragment 197-207; Fragment 471-481 and Fragment 487-495. Then, introduction of proline is used to rigidify these flexible regions. Two mutants D476P and H489P within most flexible regions are finally designed. In the results, H489P mutant shows improved thermostability while maintaining its catalytic efficiency compared to that of wild type luciferase. Flexibility analysis confirms that the overall rigidity and local rigidity of H489P mutant are greatly strengthened. D476P mutant shows decreased thermosatbility and the reason for this is elucidated at the molecular level. S307P mutation is randomly chosen outside the flexible regions as a control. Thermostability analysis shows that S307P mutation has decreased kinetic stability and enhanced thermodynamic stability.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

GENERAL SIGNIFICANCE

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

Engineering proteins for thermostability through rigidifying flexible sites

Engineering proteins for thermostability is an exciting and challenging field since it is critical for broadening the industrial use of recombinant proteins. Thermostability of proteins arises from the simultaneous effect of several forces such as hydrophobic interactions, disulfide bonds, salt bridges and hydrogen bonds. All of these interactions lead to decreased flexibility of polypeptide chain. Structural studies of mesophilic and thermophilic proteins showed that the latter need more rigid structures to compensate for increased thermal fluctuations. Hence flexibility can be an indicator to pinpoint weak spots for enhancing thermostability of enzymes. A strategy has been proven effective in enhancing proteins' thermostability with two steps: predict flexible sites of proteins firstly and then rigidify these sites. We refer to this approach as rigidify flexible sites (RFS) and give an overview of such a method through summarizing the methods to predict flexibility of a protein, the methods to rigidify residues with high flexibility and successful cases regarding enhancing thermostability of proteins using RFS.

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.

Testing of the GROMOS Force-Field Parameter Set 54A8: Structural Properties of Electrolyte Solutions, Lipid Bilayers, and Proteins

The GROMOS 54A8 force field [Reif et al. J. Chem. Theory Comput.2012, 8, 3705–3723] is the first of its kind to contain nonbonded parameters for charged amino acid side chains that are derived in a rigorously thermodynamic fashion, namely a calibration against single-ion hydration free energies. Considering charged moieties in solution, the most decisive signature of the GROMOS 54A8 force field in comparison to its predecessor 54A7 can probably be found in the thermodynamic equilibrium between salt-bridged ion pair formation and hydration. Possible shifts in this equilibrium might crucially affect the properties of electrolyte solutions or/and the stability of (bio)molecules. It is therefore important to investigate the consequences of the altered description of charged oligoatomic species in the GROMOS 54A8 force field. The present study focuses on examining the ability of the GROMOS 54A8 force field to accurately model the structural properties of electrolyte solutions, lipid bilayers, and proteins. It is found that (i) aqueous electrolytes involving oligoatomic species (sodium acetate, methylammonium chloride, guanidinium chloride) reproduce experimental salt activity derivatives for concentrations up to 1.0 m (1.0-molal) very well, and good agreement between simulated and experimental data is also reached for sodium acetate and methylammonium chloride at 2.0 m concentration, while not even qualitative agreement is found for sodium chloride throughout the whole range of examined concentrations, indicating a failure of the GROMOS 54A7 and 54A8 force-field parameter sets to correctly account for the balance between ion–ion and ion–water binding propensities of sodium and chloride ions; (ii) the GROMOS 54A8 force field reproduces the liquid crystalline-like phase of a hydrated DPPC bilayer at a pressure of 1 bar and a temperature of 323 K, the area per lipid being in agreement with experimental data, whereas other structural properties (volume per lipid, bilayer thickness) appear underestimated; (iii) the secondary structure of a range of different proteins simulated with the GROMOS 54A8 force field at pH 7 is maintained and compatible with experimental NMR data, while, as also observed for the GROMOS 54A7 force field, α-helices are slightly overstabilized with respect to 3₁₀-helices; (iv) with the GROMOS 54A8 force field, the side chains of arginine, lysine, aspartate, and glutamate residues appear slightly more hydrated and present a slight excess of oppositely-charged solution components in their vicinity, whereas salt-bridge formation properties between charged residues at the protein surface, as assessed by probability distributions of interionic distances, are largely equivalent in the GROMOS 54A7 and 54A8 force-field parameter sets.

Crowding induces differences in the diffusion of thermophilic and mesophilic proteins: a new look at neutron scattering results

The dynamical basis underlying the increased thermal stability of thermophilic proteins remains uncertain. Here, we challenge the new paradigm established by neutron scattering experiments in solution, in which the adaptation of thermophilic proteins to high temperatures lies in the lower sensitivity of their flexibility to temperature changes. By means of a combination of molecular dynamics and Brownian dynamics simulations, we report a reinterpretation of those experiments and show evidence that under crowding conditions, such as in vivo, thermophilic and homolog mesophilic proteins have diffusional properties with different thermal behavior.

Molecular dynamics of mesophilic-like mutants of a cold-adapted enzyme: insights into distal effects induced by the mutations

Networks and clusters of intramolecular interactions, as well as their "communication" across the three-dimensional architecture have a prominent role in determining protein stability and function. Special attention has been dedicated to their role in thermal adaptation. In the present contribution, seven previously experimentally characterized mutants of a cold-adapted α-amylase, featuring mesophilic-like behavior, have been investigated by multiple molecular dynamics simulations, essential dynamics and analyses of correlated motions and electrostatic interactions. Our data elucidate the molecular mechanisms underlying the ability of single and multiple mutations to globally modulate dynamic properties of the cold-adapted α-amylase, including both local and complex unpredictable distal effects. Our investigation also shows, in agreement with the experimental data, that the conversion of the cold-adapted enzyme in a warm-adapted variant cannot be completely achieved by the introduction of few mutations, also providing the rationale behind these effects. Moreover, pivotal residues, which are likely to mediate the effects induced by the mutations, have been identified from our analyses, as well as a group of suitable candidates for protein engineering. In fact, a subset of residues here identified (as an isoleucine, or networks of mesophilic-like salt bridges in the proximity of the catalytic site) should be considered, in experimental studies, to get a more efficient modification of the features of the cold-adapted enzyme.

Molecular simulations on the thermal stabilization of DNA by hyperthermophilic chromatin protein Sac7d, and associated conformational transitions

Sac7d belongs to a family of chromosomal proteins, which are crucial for thermal stabilization of DNA at higher growth temperatures. It is capable of binding DNA nonspecifically, and is responsible for the increase in the melting temperature of DNA in the bound form up to 85 °C. Molecular dynamics (MD) simulations were performed at different temperatures on two protein-DNA complexes of Sac7d. Various structural and energetic parameters were calculated to examine the DNA stability and to investigate the conformational changes in DNA and the protein-DNA interactions. Room temperature simulations indicated very good agreement with the experimental structures. The protein structure is nearly unchanged at both 300 and 360 K, and only up to five base pairs of the DNA are stabilized by Sac7d at 360 K. However, the MD simulations on DNA alone systems show that they lose their helical structures at 360 K further supporting the role of Sac7d in stabilizing the oligomers. At higher temperatures (420 and 480 K), DNA undergoes denaturation in the presence and the absence of the protein. The DNA molecules were found to undergo B- to A-form transitions consistent with experimental studies, and the extent of these transitions are examined in detail. The extent of sampling B- and A-form regions was found to show temperature and sequence dependence. Multiple MD simulations yielded similar results validating the proposed model. Interaction energy calculations corresponding to protein-DNA binding indicates major contribution due to DNA backbone, explaining the nonspecific interactions of Sac7d.

Energetic decomposition with the generalized-born and Poisson-Boltzmann solvent models: lessons from association of G-protein components

Continuum electrostatic models have been shown to be powerful tools in providing insight into the energetics of biomolecular processes. While the Poisson-Boltzmann (PB) equation provides a theoretically rigorous approach to computing electrostatic free energies of solution in such a model, computational cost makes its use for large ensembles of states impractical. The generalized-Born (GB) approximation provides a much faster alternative, although with a weaker theoretical framework. While much attention has been given to how GB recapitulates PB energetics for the overall stability of a biomolecule or the affinity of a complex, little attention has been given to how the contributions of individual functional groups are captured by the two methods. Accurately capturing these individual electrostatic components is essential both for the development of a mechanistic understanding of biomolecular processes and for the design of variant sequences and structures with desired properties. Here, we present a detailed comparison of the group-wise decomposition of both PB and GB electrostatic free energies of binding, using association of various components of the heterotrimeric-G-protein complex as a model. We find that, while net binding free energies are strongly correlated in the two models, the correlations of individual group contributions are highly variable; in some cases, strong correlation is seen, while in others, there is essentially none. Structurally, the GB model seems to capture the magnitude of direct, short-range electrostatic interactions quite well but performs more poorly with moderate-range "action-at-a-distance" interactions--GB has a tendency to overestimate solvent screening over moderate distances, and to underestimate the costs of desolvating charged groups somewhat removed from the binding interface. Despite this, however, GB does seem to be quite effective as a predictor of those groups that will be computed to be most significant in a PB-based model.

A combined solid-state NMR and MD characterization of the stability and dynamics of the HET-s(218-289) prion in its amyloid conformation

The three-dimensional structure of amyloid fibrils of the prion-forming part of the HET-s protein [HET-s(218-289)], as determined by solid-state NMR, contains rigid and remarkably well-ordered parts, as witnessed by the narrow solid-state NMR line widths for this system. On the other hand, high-resolution magic-angle-spinning (HRMAS) NMR results have shown that HET-s(218-289) amyloid fibrils contain highly flexible parts as well. Here, we further explore this unexpected behaviour using solid-state NMR and molecular dynamics (MD). The NMR data provide new information on order and dynamics in the rigid and flexible parts of HET-s(218-289), respectively. The MD study addresses whether or not small multimers, in an amyloid conformation, are stable on the 10 ns timescale of the MD run and provides insight into the dynamic parameters on the nanosecond timescale. The atom-positional, root-mean-squared fluctuations (RMSFs) and order parameters S(2) obtained are in agreement with the NMR data. A flexible loop and the N terminus exhibit dynamics on the ps-ns timescale, whereas the hydrophobic core of HET-s(218-289) is rigid. The high degree of order in the core region of HET-s(218-289) amyloids, as observed in the MD simulations, is in agreement with the narrow, solid-state, NMR lines. Finally, we employed MD to predict the behaviour of the salt-bridge network in HET-s(218-289), which cannot be obtained easily by experiment. Simulations at different temperatures indicated that the network is highly dynamic and that it contributes to the thermostability of the HET-s(218-289) amyloids.

Exploring the conserved water site and hydration of a coiled-coil trimerisation motif: a MD simulation study

The solvent structure and dynamics around ccbeta-p, a 17-residue peptide that forms a parallel three-stranded alpha-helical coiled coil in solution, was analysed through 10 ns explicit solvent molecular dynamics (MD) simulations at 278 and 330 K. Comparison with two corresponding simulations of the monomeric form of ccbeta-p was used to investigate the changes of hydration upon coiled-coil formation. Pronounced peaks in the solvent density distribution between residues Arg8 and Glu13 of neighbouring helices show the presence of water bridges between the helices of the ccbeta-p trimer; this is in agreement with the water sites observed in X-ray crystallography experiments. Interestingly, this water site is structurally conserved in many three-stranded coiled coils and, together with the Arg and Glu residues, forms part of a motif that determines three-stranded coiled-coil formation. Our findings show that little direct correlation exists between the solvent density distribution and the temporal ordering of water around the trimeric coiled coil. The MD-calculated effective residence times of up to 40 ps show rapid exchange of surface water molecules with the bulk phase, and indicate that the solvent distribution around biomolecules requires interpretation in terms of continuous density distributions rather than in terms of discrete molecules of water. Together, our study contributes to understanding the principles of three-stranded coiled-coil formation.

Single-mutation-induced stability loss in protein lysozyme

Recent NMR experiments have revealed that a single residue mutation W62G on protein hen's-egg white lysozyme can cause a dramatic loss of long-range interactions and protein stability; however, the molecular mechanism for this surprising phenomenon is not completely clear. In this mini-review, we have summarized some of our recent work on the molecular mechanism with large-scale molecular modelling, and also utilized a new wavelet method to analyse the local structural clusters present in both the wild-type and mutant folding trajectories. These extensive MD (Molecular Dynamics) simulations (10+ μs) were performed in 8 M urea, mimicking the experimental condition. Detailed analyses revealed that the Trp62 residue is the key to a co-operative long-range interaction within the wild-type protein: it acts as a bridge between neighbouring basic residues, mainly arginine residues, through π-type hydrogen bonds or π-cation interactions to form an Arg-Trp-Arg ‘sandwich-like’ local structure. The local cluster near Trp62 further extends its interaction to other clusters, such as the one near Trp111, through Arg112, which is involved in such an Arg-Trp-Arg bridging structure, thus achieving the long-range interactions for the wild-type. On the other hand, the mutant does not have this bridging effect and forms much less local clusters or contacts, and therefore results in a much less stable structure. Overall, these findings not only support the general conclusions of the experiment, but also provide a detailed but somewhat different molecular picture of the disruption of the long-range interactions.

Carboxyl pK(a) values, ion pairs, hydrogen bonding, and the pH-dependence of folding the hyperthermophile proteins Sac7d and Sso7d

Sac7d and Sso7d are homologous, hyperthermophile proteins with a high density of charged surface residues and potential ion pairs. To determine the relative importance of specific amino acid side-chains in defining the stability and function of these Archaeal chromatin proteins, pK(a) values were measured for the acidic residues in both proteins using (13)C NMR chemical shifts. The stability of Sso7d enabled titrations to pH 1 under low-salt conditions. Two aspartate residues in Sso7d (D16 and D35) and a single glutamate residue (G54) showed significantly perturbed pK(a) values in low salt, indicating that the observed pH-dependence of stability was primarily due to these three residues. The pH-dependence of backbone amide NMR resonances demonstrated that perturbation of all three pK(a) values was primarily the result of side-chain to backbone amide hydrogen bonds. Few of the significantly perturbed acidic pK(a) values in Sac7d and Sso7d could be attributed to primarily ion pair or electrostatic interactions. A smaller perturbation of E48 (E47 in Sac7d) was ascribed to an ion pair interaction that may be important in defining the DNA binding surface. The small number (three) of significantly altered pK(a) values was in good agreement with a linkage analysis of the temperature, pH, and salt-dependence of folding. The linkage of the ionization of two or more side-chains to protein folding led to apparent cooperativity in the pH-dependence of folding, although each group titrated independently with a Hill coefficient near unity. These results demonstrate that the acid pH-dependence of protein stability in these hyperthermophile proteins is due to independent titration of acidic residues with pK(a) values perturbed primarily by hydrogen bonding of the side-chain to the backbone. This work demonstrates the need for caution in using structural data alone to argue the importance of ion pairs in stabilizing hyperthermophile proteins.

Configurational entropy elucidates the role of salt-bridge networks in protein thermostability

Detailed knowledge of how networks of surface salt bridges contribute to protein thermal stability is essential not only to understand protein structure and function but also to design thermostable proteins for industrial applications. Experimental studies investigating thermodynamic stability through measurements of free energy associated with mutational alterations in proteins provide only macroscopic evidence regarding the structure of salt-bridge networks and assessment of their contribution to protein stability. Using explicit-solvent molecular dynamics simulations to provide insight on the atomic scale, we investigate here the structural stability, defined in terms of root-mean-square fluctuations, of a short polypeptide designed to fold into a stable trimeric coiled coil with a well-packed hydrophobic core and an optimal number of intra- and interhelical surface salt bridges. We find that the increase of configurational entropy of the backbone and side-chain atoms and decreased pair correlations of these with increased temperature are consistent with nearly constant atom-positional root-mean-square fluctuations, increased salt-bridge occupancies, and stronger electrostatic interactions in the coiled coil. Thus, our study of the coiled coil suggests a mechanism in which well-designed salt-bridge networks could accommodate stochastically the disorder of increased thermal motion to produce 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.

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.

Destruction of long-range interactions by a single mutation in lysozyme

We propose a mechanism, based on a > or =10-micros molecular dynamics simulation, for the surprising misfolding of hen egg-white lysozyme caused by a single mutation (W62G). Our simulations of the wild-type and mutant lysozymes in 8 M urea solution at biological temperature (with both pH 2 and 7) reveal that the mutant structure is much less stable than that of the wild type, with the mutant showing larger fluctuations and less native-like contacts. Analysis of local contacts reveals that the Trp-62 residue is the key to a cooperative long-range interaction within the wild type, where it acts like a bridge between two neighboring basic residues. Thus, a native-like cluster or nucleation site can form near these residues in the wild type but not in the mutant. The time evolution of the secondary structure also exhibits a quicker loss of the beta-sheets in the mutant than in the wild type, whereas some of the alpha-helices persist during the entire simulation in both the wild type and the mutant in 8 M urea (even though the tertiary structures are basically all gone). These findings, while supporting the general conclusions of a recent experimental study by Dobson and coworkers [Klein-Seetharam J, Oikama M, Grimshaw SB, Wirmer J, Duchardt E, Ueda T, Imoto T, Smith LJ, Dobson CM, Schwalbe H (2002) Science 295:1719-1722], provide a detailed but different molecular picture of the misfolding mechanism.

Thermal Denaturing of Mutant Lysozyme with Both the OPLSAA and the CHARMM Force Fields

Biomolecular simulations enabled by massively parallel supercomputers such as BlueGene/L promise to bridge the gap between the currently accessible simulation time scale and the experimental time scale for many important protein folding processes. In this study, molecular dynamics simulations were carried out for both the wild-type and the mutant hen lysozyme (TRP62GLY) to study the single mutation effect on lysozyme stability and misfolding. Our thermal denaturing simulations at 400-500 K with both the OPLSAA and the CHARMM force fields show that the mutant structure is indeed much less stable than the wild-type, which is consistent with the recent urea denaturing experiment (Dobson et al. Science 2002, 295, 1719-1722; Nature 2003, 424, 783-788). Detailed results also reveal that the single mutation TRP62GLY first induces the loss of native contacts in the beta-domain region of the lysozyme protein at high temperatures, and then the unfolding process spreads into the alpha-domain region through Helix C. Even though the OPLSAA force field in general shows a more stable protein structure than does the CHARMM force field at high temperatures, the two force fields examined here display qualitatively similar results for the misfolding process, indicating that the thermal denaturing of the single mutation is robust and reproducible with various modern force fields.

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.