Denaturation of Drew-Dickerson DNA in a high salt concentration medium: Molecular dynamics simulations

How global DNA unwinding causes non-uniform stress distribution and melting of DNA

DNA unwinding is an important process that controls binding of proteins, gene expression and melting of double-stranded DNA. In a series of all-atom MD simulations on two DNA molecules containing a transcription start TATA-box sequence we demonstrate that application of a global restraint on the DNA twisting dramatically changes the coupling between helical parameters and the distribution of deformation energy along the sequence. Whereas only short range nearest-neighbor coupling is observed in the relaxed case, long-range coupling is induced in the globally restrained case. With increased overall unwinding the elastic deformation energy is strongly non-uniformly distributed resulting ultimately in a local melting transition of only the TATA box segment during the simulations. The deformation energy tends to be stored more in cytidine/guanine rich regions associated with a change in conformational substate distribution. Upon TATA box melting the deformation energy is largely absorbed by the melting bubble with the rest of the sequences relaxing back to near B-form. The simulations allow us to characterize the structural changes and the propagation of the elastic energy but also to calculate the associated free energy change upon DNA unwinding up to DNA melting. Finally, we design an Ising model for predicting the local melting transition based on empirical parameters. The direct comparison with the atomistic MD simulations indicates a remarkably good agreement for the predicted necessary torsional stress to induce a melting transition, for the position and length of the melted region and for the calculated associated free energy change between both approaches.

Asymmetric base-pair opening drives helicase unwinding dynamics

The opening of a Watson-Crick double helix is required for crucial cellular processes, including replication, repair, and transcription. It has long been assumed that RNA or DNA base pairs are broken by the concerted symmetric movement of complementary nucleobases. By analyzing thousands of base-pair opening and closing events from molecular simulations, here, we uncover a systematic stepwise process driven by the asymmetric flipping-out probability of paired nucleobases. We demonstrate experimentally that such asymmetry strongly biases the unwinding efficiency of DNA helicases toward substrates that bear highly dynamic nucleobases, such as pyrimidines, on the displaced strand. Duplex substrates with identical thermodynamic stability are thus shown to be more easily unwound from one side than the other, in a quantifiable and predictable manner. Our results indicate a possible layer of gene regulation coded in the direction-dependent unwindability of the double helix.

Effect of gold nanoparticle on stability of the DNA molecule: A study of molecular dynamics simulation

An understanding of the mechanism of DNA interactions with gold nanoparticles is useful in today medicine applications. We have performed a molecular dynamics simulation on a B-DNA duplex (CCTCAGGCCTCC) in the vicinity of a gold nanoparticle with a truncated octahedron structure composed of 201 gold atoms (diameter ∼1.8 nm) to investigate gold nanoparticle (GNP) effects on the stability of DNA. During simulation, the nanoparticle is closed to DNA and phosphate groups direct the particles into the major grooves of the DNA molecule. Because of peeling and untwisting states that are occur at end of DNA, the nucleotide base lies flat on the surface of GNP. The configuration entropy is estimated using the covariance matrix of atom-positional fluctuations for different bases. The results show that when a gold nanoparticle has interaction with DNA, entropy increases. The results of conformational energy and the hydrogen bond numbers for DNA indicated that DNA becomes unstable in the vicinity of a gold nanoparticle. The radial distribution function was calculated for water hydrogen-phosphate oxygen pairs. Almost for all nucleotide, the presence of a nanoparticle around DNA caused water molecules to be released from the DNA duplex and cations were close to the DNA.

Evaluation of the Gibbs Free Energy Changes and Melting Temperatures of DNA/DNA Duplexes Using Hybridization Enthalpy Calculated by Molecular Dynamics Simulation

A molecular dynamics simulation approach was applied for the prediction of the thermal stability of oligonucleotide duplexes. It was shown that the enthalpy of the DNA/DNA complex formation could be calculated using this approach. We have studied the influence of various simulation parameters on the secondary structure and the hybridization enthalpy value of Dickerson-Drew dodecamer. The optimal simulation parameters for the most reliable prediction of the enthalpy values were determined. The thermodynamic parameters (enthalpy and entropy changes) of a duplex formation were obtained experimentally for 305 oligonucleotides of various lengths and GC-content. The resulting database was studied with molecular dynamics (MD) simulation using the optimized simulation parameters. Gibbs free energy changes and the melting temperatures were evaluated using the experimental correlation between enthalpy and entropy changes of the duplex formation and the enthalpy values calculated by the MD simulation. The average errors in the predictions of enthalpy, the Gibbs free energy change, and the melting temperature of oligonucleotide complexes were 11%, 10%, and 4.4 °C, respectively. We have shown that the molecular dynamics simulation gives a possibility to calculate the thermal stability of native DNA/DNA complexes a priori with an unexpectedly high accuracy.

Investigation of changes in the arrangement of water molecules and salt ions surrounding different atoms of the DNA molecule during the melting process: a molecular dynamics simulation study

A molecular dynamics simulation was performed on a B-DNA duplex (CGCGAATTGCGC) at different temperatures. The DNA was immerged in a saltwater medium with 1 mol/L NaCl concentration. The arrangements of water molecules and cations around the different atoms of DNA on the melting pathway were investigated. Almost for all atoms of the DNA by double helix → single-stranded transition, the water molecules released from the DNA duplex and cations were close to single-stranded DNA, but this behavior was not clearly seen at melting temperatures. Therefore, release of water molecules and cations approaching the DNA by the increase of temperature does not have any effect on the sharpness of the transition curve. Most of the water molecules and cations were found to be around the negatively charged phosphate oxygen atoms. The number of water molecules released from the first shell hydration upon melting in the minor groove was higher than in the major groove, and intrusion of cations into the minor groove after melting was higher than into the major groove. The hydrations of imino protons were different from each other and were dependent on DNA bases.

Temperature effect on poly(dA).poly(dT): molecular dynamics simulation studies of polymeric and oligomeric constructs

Understanding unwinding and melting of double helical DNA is very important to characterize role of DNA in replication, transcription, translation etc. Sequence dependent melting thermodynamics is used extensively for detecting promoter regions but melting studies are generally done for short oligonucleotides. This study reports several molecular dynamics (MD) simulations of homopolymeric poly(dA).poly(dT) as regular oligonucleotide fragments as well as its corresponding polymeric constructs with water and charge-neutralizing counterions at different temperatures ranging from 300 to 400 K. We have eliminated the end-effect or terminal peeling propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The dynamic properties such as basepairing and stacking geometry, groove width, backbone conformational parameters, bending, distribution of counter ions and number of hydrogen bonds of oligomeric and polymeric constructs of poly(dA).poly(dT) have been analyzed. The oligomer shows terminal fraying or peeling effect at temperatures above 340 K. The polymer shows partial melting at elevated temperatures although complete denaturations of basepairs do not take place. The analysis of cross strand hydrogen bonds shows that the number of N-H···O hydrogen bonds increases with increase in temperature while C-H···O hydrogen bond frequencies decrease with temperature. Restructuring of counterions in the minor groove with temperature appear as initiation of melting in duplex structures.

Effect of temperature on DNA double helix: An insight from molecular dynamics simulation

The three-dimensional structure of DNA contains various sequence-dependent structural information, which control many cellular processes in life, such as replication, transcription, DNA repair, etc. For the above functions, DNA double helices need to unwind or melt locally, which is different from terminal melting, as often seen in molecular dynamics (MD) simulations or even in many DNA crystal structures. We have carried out detailed MD simulations of DNA double helices of regular oligonucleotide fragments as well as in polymeric constructs with water and charge-neutralizing counter-ions at several different temperatures. We wanted to eliminate the end-effect or terminal melting propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The polymeric construct is expected to allow us to see local melting at elevated temperatures. Comparative structural analysis of oligonucleotides and its corresponding virtual polymer at various temperatures ranging from 300 K to 400 K is discussed. The general behaviour, such as volume expansion coefficients of both the simulations show high similarity, indicating polymeric construct, does not give many artificial constraints. Local melting of a polymer, even at elevated temperature, may need a high nucleation energy that was not available in the short (7 ns) simulations. We expected to observe such nucleation followed by cooperative melting of the polymers in longer MD runs. Such simulations of different polymeric sequences would facilitate us to predict probable melting origins in a polymeric DNA.