Nucleic Acids Research · 2021 · 139 citations · 78 references
Natural SciencesBiophysicsMolecular BiologyComputational BiophysicsComputational ChemistryMolecular BiophysicsBase StackingComputational Molecular BiophysicsDna ChemistryDna StructureMolecular DynamicsNucleic Acid DynamicsBase Stacking ModelsBiomolecular SimulationDna Duplex StabilityDna Unwinding MechanismIntercalation Effects On Base StackingMolecular Dynamics Of ComplexesMolecular Dynamics Of Nucleic AcidsMolecular Dynamics Simulations Of Nucleic Acids
Visualization of double stranded DNA in gels with the binding of the fluorescent dye ethidium bromide has been a basic experimental technique in any molecular biology laboratory for >40 years. The interaction between ethidium and double stranded DNA has been observed to be an intercalation between base pairs with strong experimental evidence. This presents a unique opportunity for computational chemistry and biomolecular simulation techniques to benchmark and assess their models in order to see if the theory can reproduce experiments and ultimately provide new insights. We present molecular dynamics simulations of the interaction of ethidium with two different double stranded DNA models. The first model system is the classic sequence d(CGCGAATTCGCG)2 also known as the Drew-Dickerson dodecamer. We found that the ethidium ligand binds mainly stacked on, or intercalated between, the terminal base pairs of the DNA with little to no interaction with the inner base pairs. As the intercalation at the terminal CpG steps is relatively rapid, the resultant DNA unwinding, rigidification, and increased stability of the internal base pair steps inhibits further intercalation. In order to reduce these interactions and to provide a larger groove space, a second 18-mer DNA duplex system with the sequence d(GCATGAACGAACGAACGC) was tested. We computed molecular dynamics simulations for 20 independent replicas with this sequence, each with ∼27 μs of sampling time. Results show several spontaneous intercalation and base-pair eversion events that are consistent with experimental observations. The present work suggests that extended MD simulations with modern DNA force fields and optimized simulation codes are allowing the ability to reproduce unbiased intercalation events that we were not able to previously reach due to limits in computing power and the lack of extensively tested force fields and analysis tools.
78
VMD: Visual molecular dynamics
William Humphrey, Andrew Dalke, Klaus Schulten · Journal of Molecular Graphics · 1996 · 66.6K citations
A smooth particle mesh Ewald method
Ulrich Essmann, L. Perera, Max L. Berkowitz et al. · The Journal of Chemical Physics · 1995 · 22.9K citations
Development and testing of a general amber force field
Junmei Wang, Romain M. Wolf, James W. Caldwell et al. · Journal of Computational Chemistry · 2004 · 19.6K citations