Mechanisms of DNA Hybridization: Transition Path Analysis of a Simulation-Informed Markov Model

TL;DR

This study combines simulations and Markov models to analyze DNA hybridization mechanisms, revealing a dominant zipper-like pathway and emphasizing the complexity of transition states beyond simple base pairing counts.

Contribution

It introduces a multi-step computational framework integrating Brownian dynamics, Markov State Models, and Transition Path Theory to elucidate DNA hybridization pathways.

Findings

Identifies a dominant zipper-like hybridization pathway

Shows the number of base pairs alone doesn't define the transition state

Provides quantitative thermodynamic and kinetic insights

Abstract

Complementary DNA strands in solution reliably hybridize to form stable duplexes. We study the kinetics of the hybridization process and the mechanisms by which two initially isolated strands come together to form a stable double helix. We adopt a multi-step computational approach. First, we perform a large number of Brownian dynamics simulations of the hybridization process using the coarse-grained oxDNA2 model. Second, we use these simulations to construct a Markov State Model of DNA dynamics that uses a state decomposition based on the inter-strand hydrogen bonding pattern. Third, we take advantage of Transition Path Theory to obtain quantitative information about the thermodynamic and dynamic properties of the hybridization process. We find that while there is a large ensemble of possible hybridization pathways there is a single dominant mechanism in which an initial base pair forms close to either end of the nascent double helix, and the remaining bases pair sequentially in a zipper-like fashion. We also show that the number of formed base pairs by itself is insufficient to describe the transition state of the hybridization process.