Molecular dynamics simulation of homo-DNA: The role of crystal packing in duplex conformation

Abstract

The (4'→6')-linked DNA homolog 2',3'-dideoxy-β-D-glucopyranosyl nucleic acid (dideoxy-glucose nucleic acid or homo-DNA) exhibits stable self-pairing of the Watson–Crick and reverse-Hoogsteen types, but does not cross-pair with DNA. Molecular modeling and NMR solution studies of homo-DNA duplexes pointed to a conformation that was nearly devoid of a twist and a stacking distance in excess of 4.5 Å. By contrast, the crystal structure of the homo-DNA octamer dd(CGAATTCG) revealed a right-handed duplex with average values for helical twist and rise of ca. 15° and 3.8 Å, respectively. Other key features of the structure were strongly inclined base-pair and backbone axes in the duplex with concomitant base-pair slide and cross-strand stacking, and the formation of a dimer across a crystallographic dyad with inter-duplex base swapping. To investigate the conformational flexibility of the homo-DNA duplex and a potential influence of lattice interactions on its geometry, we used molecular dynamics (MD) simulations of the crystallographically observed dimer of duplexes and an isolated duplex in the solution state. The dimer of duplexes showed limited conformational flexibility, and key parameters such as helical rise, twist, and base-pair slide exhibited only minor fluctuations. The single duplex was clearly more flexible by comparison and underwent partial unwinding, albeit without significant lengthening. Thus, base stacking was preserved in the isolated duplex and two adenosines extruded from the stack in the dimer of duplexes were reinserted into the duplex and pair with Ts in a Hoogsteen mode. Our results confirmed that efficient stacking in homo-DNA seen in the crystal structure of a dimer of duplexes was maintained in the separate duplex. Therefore, lattice interactions did not account for the different geometries of the homo-DNA duplex in the crystal and earlier models that resembled inclined ladders with large base-pair separations that precluded efficient stacking.