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

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

Role of hydrophobic core on the thermal stability of proteins-molecular dynamics simulations on a single point mutant of sso7d

Abstract The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Role of hydrophobic core on the thermal stability of proteins - molecular dynamics simulations on a single point mutant of Sso7d abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 s long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Role of Hydrophobic Core on the Thermal Stability of Proteins—Molecular Dynamics Simulations on a Single Point Mutant of Sso7d

Abstract

The role of salt bridges in chromatin protein Sso7d, from S. solfataricus has previously been shown to be crucial for its unusual high thermal stability. Experimental studies have shown that single site mutation of Sso7d (F31A) leads to a substantial decrease in the thermal stability of this protein due to distortion of the hydrophobic core. In the present study, we have performed a total of 0.2 μs long molecular dynamics (MD) simulations on F31A at room temperature, and at 360 K, close to the melting temperature of the wild type (WT) protein to investigate the role of hydrophobic core on protein stability. Sso7d-WT was shown to be stable at both 300 and 360 K; however, F31A undergoes denaturation at 360 K, consistent with experimental results. The structural and energetic properties obtained using the analysis of MD trajectories indicate that the single mutation results in high flexibility of the protein, and loosening of intramolecular interactions. Correlation between the dynamics of the salt bridges with the structural transitions and the unfolding pathway indicate the importance of both salt bridges and hydrophobic in effecting thermal stability of proteins in general.     

Structural and dynamic effects of alpha-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 and IR spectroscopy in investigation of phase transitions in molecular crystal 4,4'-bis(11-hydroxy-1-undecyloxy)biphenyl

Molecular dynamics (MD) simulations combined with temperature-dependent IR spectroscopic measurements were used to study phase transitions in molecular crystals of the mesogenic diol 4,4'-bis(11-hydroxy-1-undecyloxy)biphenyl. DSC measurements revealed four phase transitions in this molecular crystal at approximately 327.1 K, 389.8 K, 419.1 K and 431.9 K. Analysis of the dynamic trajectories at temperatures of 300 K, 360 K, 400 K and 480 K revealed changes in conformation of the mesogenic diol molecules and consequently changes in crystal packing and crystal structure in the temperature range 300-480 K and enabled us to understand the mechanism of the phase transitions.     

Structural characterization of the functional regions in the archaeal protein Sso7d

Sso7d from the extreme thermophilic crenarchaeon Sulfolobus solfataricus is a multifunctional protein in in vitro assays, whose in vivo role is still puzzling. Crystals of Sso7d in complex with DNA elucidated the protein surface involved in the binding to the nucleic acid, whereas the locations of the Sso7d regions responsible for a chaperone activity in renaturing protein aggregates (i.e., the protein-binding surface and the site of ATPase activity) are still unknown. We identified the regions of Sso7d involved in protein-binding by limited proteolysis experiments associated to advanced mass spectrometric procedures performed on isolated Sso7d and Sso7d in complex with the peptide melittin. By affinity labeling of Sso7d with the ATP analogue 5'-p-fluorosulfonylbenzoyl adenosine and characterization of the labeled tryptic peptides by tandem mass spectrometry, we found that Y7 and K39 are residues involved in ATP binding/hydrolysis. Insights into the positions of the ligands melittin and ATP were achieved by a molecular modeling study; the models obtained were in agreement with most experimental data. A comparison among the complexes of Sso7d with DNA, with melittin, and with ATP showed that the DNA-binding surface and the protein-binding surface overlap, whereas the ATPase site is mostly independent of the binding sites for the nucleic acid and melittin.     

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.     

Electrostatic interactions play an essential role in DNA repair and cold-adaptation of Uracil DNA glycosylase

Life has adapted to most environments on earth, including low and high temperature niches. The increased catalytic efficiency and thermoliability observed for enzymes from organisms living in constantly cold regions when compared to their mesophilic and thermophilic cousins are poorly understood at the molecular level. Uracil DNA glycosylase (UNG) from cod (cUNG) catalyzes removal of uracil from DNA with an increased k_cat and reduced K_m relative to its warm-active human (hUNG) counterpart. Specific issues related to DNA repair and substrate binding/recognition (K_m) are here investigated by continuum electrostatics calculations, MD simulations and free energy calculations. Continuum electrostatic calculations reveal that cUNG has surface potentials that are more complementary to the DNA potential at and around the catalytic site when compared to hUNG, indicating improved substrate binding. Comparative MD simulations combined with free energy calculations using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method show that large opposing energies are involved when forming the enzyme-substrate complexes. Furthermore, the binding free energies obtained reveal that the Michaelis-Menten complex is more stable for cUNG, primarily due to enhanced electrostatic properties, suggesting that energetic fine-tuning of electrostatics can be utilized for enzymatic temperature adaptation. Energy decomposition pinpoints the residual determinants responsible for this adaptation. Figure Electrostatic isosurfaces of cod uracil DNA glycosylase in complex with double stranded DNA     

Partial B-to-A DNA transition upon minor groove binding of protein Sac7d monitored by Raman spectroscopy

Members of the Sso7d/Sac7d protein family and other related proteins are believed to play an important role in DNA packaging and maintenance in archeons. Sso7d/Sac7d are small, abundant, basic, and nonspecific DNA-binding proteins of the hyperthermophilic archeon Sulfolobus. Structures of several complexes of Sso7d/Sac7d with DNA octamers are known. These structures are characterized by sequence unspecific minor groove binding of the proteins and sharp kinking of the double helix. Corresponding Raman vibrational signatures have been identified in this study. A Raman spectroscopic analysis of Sac7d binding to the oligonucleotide decamer d(GAGGCGCCTC)(2) reveals large conformational perturbations in the DNA structure upon complex formation. Perturbed Raman bands are associated with the vibrational modes of the sugar phosphate backbone and frequency shifts of bands assigned to nucleoside vibrations. Large changes in the DNA backbone and partial B- to A-form DNA transitions are indicated that are closely associated with C2'-endo/anti to C3'-endo/anti conversion of the deoxyadenosyl moiety upon Sac7d binding. The major spectral feature of Sac7d binding is kinking of the DNA. Raman markers of minor groove binding do not largely contribute to spectral differences; however, clear indications for minor groove binding come from G-N2 and G-N3 signals that are supported by Trp24 features. Trp24 is the only tryptophan present in Sac7d and binds to guanine N3, as has been demonstrated clearly in X-ray structures of Sac7d-DNA complexes. No changes of the Sac7d secondary structure have been detected upon DNA binding.     

DNA bending,compaction and negative supercoiling by the architectural protein Sso7d of Sulfolobus solfataricus

Members of the Sso7d/Sac7d family are small, abundant, non-specific DNA-binding proteins of the hyperthermophilic Archaea Sulfolobus. Crystal structures of these proteins in complex with oligonucleotides showed that they induce changes in the helical twist and marked DNA bending. On this basis they have been suggested to play a role in organising chromatin structures in these prokaryotes, which lack histones. We report functional in vitro assays to investigate the effects of the observed Sso7d-induced structural modifications on DNA geometry and topology. We show that binding of multiple Sso7d molecules to short DNA fragments induces significant curvature and reduces the stiffness of the complex. Sso7d induces negative supercoiling of DNA molecules of any topology (relaxed, positively or negatively supercoiled) and in physiological conditions of temperature and template topology. Binding of Sso7d induces compaction of positively supercoiled and relaxed DNA molecules, but not of negatively supercoiled ones. Finally, Sso7d inhibits the positive supercoiling activity of the thermophile-specific enzyme reverse gyrase. The proposed biological relevance of these observations is that these proteins might model the behaviour of DNA in constrained chromatin environments.     

Molecular dynamics simulations and MM–PBSA calculations of the lectin from snowdrop (Galanthus nivalis)

Galanthus nivalis agglutinin (GNA), a mannose-specific lectin from snowdrop bulbs, is a member of the monocot mannose-specific lectin family and exhibits antiviral activity toward HIV. In the present study, molecular dynamics (MD) simulations were performed to study the interaction between GNA and its carbohydrate ligand over a specific time span. By analysis of the secondary structures, it was observed that the GNA conformation maintains rather stable along the trajectories and the high fluctuations were only centered on the carbohydrate recognition domains. Our MD simulations also reproduced most of the hydrogen bonds observed in the x-ray crystal structure. Furthermore, the obtained MD trajectories were used to estimate the binding free energy of the complex using the molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) method. It was revealed by the inspection of the binding free energy components that the major contributions to the complex stability arose from electrostatic interactions.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

In vitro DNA binding of the archaeal protein Sso7d induces negative supercoiling at temperatures typical for thermophilic growth.

The topological state of DNA in hyperthermophilic archaea appears to correspond to a linking excess in comparison with DNA in mesophilic organisms. Since DNA binding proteins often contribute to the control of DNA topology by affecting DNA geometry in the presence of DNA topoisomerases, we tested whether the histone-like protein Sso7d from the hyperthermophilic archaeon Sulfolobus solfataricus alters DNA conformation. In ligase-mediated supercoiling assays carried out at 37, 60, 70, 80 and 90 degrees C we found that DNA binding of increasing amounts of Sso7d led to a progressive decrease in plasmid linking number (Lk), producing negative supercoiling. Identical unwinding effects were observed when recombinant non-methylated Sso7d was used. For a given Sso7d concentration the DNA unwinding induced was augmented with increasing temperature. However, after correction for the overwinding effect of high temperature on DNA, plasmids ligated at 60-90 degrees C exhibited similar sigma values at the highest Sso7d concentrations assayed. These results suggest that Sso7d may play a compensatory role in vivo by counteracting the overwinding effect of high temperature on DNA. Additionally, Sso7d unwinding could be involved in the topological changes observed during thermal stress (heat and cold shock), playing an analogous role in crenarchaeal cells to that proposed for HU in bacteria.     

Disruption and formation of surface salt bridges are coupled to DNA binding by the integration host factor: a computational analysis

Revealing the thermodynamic driving force of protein-DNA interactions is crucial to the understanding of factors that dictate the properties and function of protein-DNA complexes. For the binding of DNA to DNA-wrapping proteins, such as the integration host factor (IHF), Record and co-workers proposed that the disruption of a large number of preexisting salt bridges is coupled with the binding process [Holbrook, J. A., et al. (2001) J. Mol. Biol. 310, 379]. To test this proposal, we have conducted explicit solvent MD simulations (multiple ～25-50 ns trajectories for each salt concentration) to examine the behavior of charged residues in IHF, especially concerning their ability to form salt bridges at different salt concentrations. Of the 17 cationic residues noted by Record and co-workers, most are engaged in salt bridge interactions for a significant portion of the trajectories, especially in the absence of salt. This observation suggests that, from a structural point of view, their proposal is plausible. However, the complex behaviors of charged residues observed in the MD simulations also suggest that the unusual thermodynamic characteristics of IHF-DNA binding likely arise from the interplay between complex dynamics of charged residues both in and beyond the DNA binding site. Moreover, a comparison of MD simulations at different salt concentrations suggests that the strong dependence of the IHF-DNA binding enthalpy on salt concentration may not be due to a significant decrease in the number of stable salt bridges in apo IHF at high salt concentrations. In addition to the Hofmeister effects quantified in more recent studies of IHF-DNA binding, we recommend consideration of the variation of the enthalpy change of salt bridge disruption at different salt concentrations. Finally, the simulation study presented here explicitly highlights the fact that the electrostatic properties of DNA-binding proteins can be rather different in the apo and DNA-bound states, which has important implications for the design of robust methods for predicting DNA binding sites in proteins.     

Biochemical characterization of DNA-binding proteins from Pyrobaculum aerophilum and Aeropyrum pernix

Several representatives of the Crenarchaeal branch of the Archaea contain highly abundant, small, positively charged proteins exemplified by the Sso7d protein from Sulfolobus solfataricus. These proteins bind to DNA in a non-sequence-specific manner. Using publicly available genomic sequence information, we identified a second class of small Crenarchaeal DNA-binding proteins represented by the Pyrobaculum aerophilum open reading frame 3192–encoded (Pae3192) protein and its paralogs. We investigated the biochemical properties of the Pae3192 protein and an orthologous protein (Ape1322b) from Aeropyrum pernix in side-by-side experiments with the Sso7d protein. We demonstrate that the recombinant Ape1322b, Pae3192 and Sso7d proteins bind to DNA and that the DNA-protein complexes formed are slightly different for each protein. We show that like Sso7d, Pae3192 constrains negative supercoils in DNA. In addition, we show that all three proteins raise the melting temperature of duplex DNA upon binding. Finally, we present the equilibrium affinity constants and kinetic association constants of each protein for single-stranded and double-stranded DNA.     

Structure prediction and functional analyses of a thermostable lipase obtained from Shewanella putrefaciens

Previous experimental studies on thermostable lipase from Shewanella putrefaciens suggested the maximum activity at higher temperatures, but with little information on its conformational profile. In this study, the three-dimensional structure of lipase was predicted and a 60 ns molecular dynamics (MD) simulation was carried out at temperatures ranging from 300 to 400 K to better understand its thermostable nature at the molecular level. MD simulations were performed in order to predict the optimal activity of thermostable lipase. The results suggested strong conformational temperature dependence. The thermostable lipase maintained its bio-active conformation at 350 K during the 60 ns MD simulations.