Molecular dynamics simulation of hydrated d(CGGGTACCCG)4 as a four-way DNA Holliday junction and comparison with the crystallographic structure

Molecular dynamics (MD) simulation of decamer sequence (CGGGTACCCG)(4) as a four-way Holliday junction is reported here for 15.0ns at three different temperatures 100, 200 and 300K, respectively, using AMBER force field. Particle mesh Ewald method has been utilised to deal long-range interaction potentials. After MD simulation, various parameters of the junction model including backbone and helical parameters have been worked out and the dynamical pathway is discussed. Structural analysis and geometrical calculations were carried out through X3DNA. The computational results obtained are compared with the previously reported crystallographic outcomes. The width and depth of the major and minor grooves of the duplex of the four arms of the DNA junction have been calculated. The variations in the C1-C1 distances between the two complementary strands are discussed in detail. A close observation of the results reveals that the conformation of the average simulated structure at low temperature is of B' form and the structural integrity of the DNA junction having a twofold sequence symmetry is temperature dependent. It also seems that besides the other parameters (i.e. presence of ions, solvents, etc.), temperature may be playing a key role in preserving the structural integrity of the DNA junction.