Search

DNA hybridization plays a central role in biology and, increasingly, in materials science. Yet, there is no precedent for examining the pathways by which specific single-stranded DNA sequences interact to assemble into a double helix. A detailed model of DNA is adopted in this work to examine such pathways and to determine the role of sequence, if any, on DNA hybridization. Transition path sampling simulations reveal that DNA rehybridization is prompted by a distinct nucleation event involving molecular sites with approximately four bases pairing with partners slightly offset from those involved in ideal duplexation. Nucleation is promoted in regions with repetitive base pair sequence motifs, which yield multiple possibilities for finding complementary base partners. Repetitive sequences follow a nonspecific pathway to renaturation consistent with a molecular "slithering" mechanism, whereas random sequences favor a restrictive pathway involving the formation of key base pairs before renaturation fully ensues.

A mesoscale model of DNA is presented (3SPN.1), extending the scheme previously developed by our group. Each nucleotide is mapped onto three interaction sites. Solvent is accounted for implicitly through a medium-effective dielectric constant and electrostatic interactions are treated at the level of Debye-Hückel theory. The force field includes a weak, solvent-induced attraction, which helps mediate the renaturation of DNA. Model parameterization is accomplished through replica exchange molecular dynamics simulations of short oligonucleotide sequences over a range of composition and chain length. The model describes the melting temperature of DNA as a function of composition as well as ionic strength, and is consistent with heat capacity profiles from experiments. The dependence of persistence length on ionic strength is also captured by the force field. The proposed model is used to examine the renaturation of DNA. It is found that a typical renaturation event occurs through a nucleation step, whereby an interplay between repulsive electrostatic interactions and colloidal-like attractions allows the system to undergo a series of rearrangements before complete molecular reassociation occurs.

The conformation and diffusion of a single DNA molecule confined between two parallel plates are examined using both single molecule experiments and Brownian dynamics simulations accounting for hydrodynamic interactions. The degree of chain stretching and the diffusivity are characterized as a function of the chain confinement and the channel geometry. Good agreement is found between the simulations, experiments, and scaling theory predictions.