Transition between B-DNA and Z-DNA: free energy landscape for the B-Z junction propagation

Canonical, right-handed B-DNA can be transformed into noncanonical, left-handed Z-DNA in vitro at high salt concentrations or in vivo under physiological conditions. The molecular mechanism of this drastic conformational transition is still unknown despite numerous studies. Inspired by the crystal structure of a B-Z junction and the previous zipper model, we show here, with the aid of molecular dynamics simulations, that a stepwise propagation of a B-Z junction is a highly probable pathway for the B-Z transition. In this paper, the movement of a B-Z junction by a two-base-pair step in a double-strand nonamer, [d(GpCpGpCpGpCpGpCpG)](2), is considered. Targeted molecular dynamics simulations and umbrella sampling for this transition resulted in a transition pathway with a free energy barrier of 13 kcal/mol. This barrier is much more favorable than those obtained from previous atomistic simulations that lead to concerted transitions of the whole strands. The free energy difference between B-DNA and Z-DNA evaluated from our simulation is 0.9 kcal/mol per dinucleotide unit, which is consistent with previous experiments. The current computation thus strongly supports the proposal that the B-Z transition involves a relatively fast extension of B-DNA or Z-DNA by sequential propagation of B-Z junctions once nucleation of junctions is established.