Abstract:DNA minicircles are closed double-stranded DNA (dsDNA) fragments that have been demonstrated to be an important experimental tool to understand supercoiled, or stressed, DNA mechanics, such as nucleosome positioning and DNA-protein interactions. Specific minicircles can be simulated using Molecular Dynamics (MD) simulation. However, the enormous sequence space makes it unfeasible to exhaustively explore the sequence-dependent mechanics of DNA minicircles using either experiment or MD. For linear fragments, the cgNA+ model, a computationally efficient sequence-dependent coarse-grained model using enhanced Curves+ internal coordinates (rigid base plus rigid phosphate) of double-stranded nucleic acids (dsNAs), predicts highly accurate nonlocal sequence-dependent equilibrium distributions for an arbitrary sequence when compared with MD simulations. This article addresses the problem of modeling sequence-dependent topologically closed and, therefore, stressed fragments of dsDNA. We introduce cgNA+min, a computational approach within the cgNA+ framework, which extends the cgNA+ model applicability to compute the sequence-dependent energy minimising configurations of covalently closed dsDNA minicircles of various lengths and linking numbers (Lk). The main idea is to derive the appropriate chain rule to express the cgNA+ energy in absolute coordinates involving quaternions where the closure condition is simple to handle. We also present a semi-analytic method for efficiently computing sequence-dependent initial minicircles having arbitrary Lk and length. For different classes and lengths of sequences, we demonstrate that the dsDNA minicircle energies computed using cgNA+min agree well with the energies approximated from experimentally measured J-factor values. Finally, we present the minicircle shape, energy, and multiplicity of Lk for more than 120K random DNA sequences of different lengths.

Unexpected Dissociation Energetics Of The Na+ Counterion From Gc Motifs In Dna Hole-Migration

Metastable Hydrogen-bonds Featuring Negative Dissociation Energies in Protein-bound DNA in Hole Migration

Wrapping/Unwrapping Transition of Double-Stranded Dna in Dna-Nanosphere Complexes Induced by Multivalent Anions

Insights Into Guanine Radical Cation Deprotonation Using the Quantum Mechanics and Quantum Mechanics/Molecular Mechanics (ABEEM) Methods

Positively and Negatively Hydrated Counterions in Molecular Dynamics Simulations of DNA Double Helix

DNA minor-groove binder Hoechst 33258 destabilizes base-pairing adjacent to its binding site

Molecular dynamics simulations reveal the balance of forces governing the formation of a guanine tetrad—a common structural unit of G-quadruplex DNA

Multiple-charge transfer and trapping in DNA dimers

Negative Electron Transfer Collision-Induced Dissociation of G-Quadruplexes: Uncovering the Guanine Radical Anion Loss Pathway

Communication: Asymmetrical Cation Movements Through G-quadruplex DNA.

Quantitative Effects of the Loop Region on Topology, Thermodynamics, and Cation Binding of DNA G-quadruplexes

Attractive Forces between Cation Condensed DNA Double Helices

Unexpected large charge transfer rate mediated by adenine in twisted DNA structure

Structuring of counterions around dna double helix: a molecular dynamics study

Predicting Site-Binding Modes of Ions and Water to Nucleic Acids Using Molecular Solvation Theory

Excess-electron Capture and Energy Transfer to Bulk Water for Aqueous DNA Nucleotide

cgNA+min: computation of sequence-dependent dsDNA energy-minimising minicircle configurations

Slow G-Quadruplex Conformation Rearrangement and Accessibility Change Induced by Potassium in Human Telomeric Single-Stranded DNA

Counterion atmosphere around DNA double helix: trapping of counterions at the nanoscale

Simulation Study on Dynamics of A- to B-form Transition in Aqueous DNA Solution: Effect of Alkali Metal Counterions

Multivalent Counterions Induced Attraction Between DNA Polyelectrolytes