Abstract
In this study, we report what we believe to be the first multiscale simulation of the dynamic relaxation of DNA supercoils by human topoisomerase IB (topo IB). We leverage our previous molecular dynamics calculations of the free energy landscape describing the interaction between a short DNA fragment and topo IB. Herein, this landscape is used to prescribe boundary conditions for a computational, elastodynamic continuum rod model of a long length of supercoiled DNA. The rod model, which accounts for the nonlinear bending, twisting, and electrostatic interaction of the (negatively charged) DNA backbone, is extended to include the hydrodynamic drag induced by the surrounding physiological buffer. Simulations for a 200-bp-long DNA supercoil in complex with topo IB reveal a relaxation timescale of ∼0.1-1.0 μs. The relaxation follows a sequence of cascading reductions in the supercoil linking number (Lk), twist (Tw), and writhe (Wr) that follow companion cascading reductions in the supercoil elastic and electrostatic energies. The novel (to our knowledge) multiscale modeling method may enable simulations of the entire experimental setup that measures DNA supercoiling and relaxation via single molecule magnetic trapping.