The effect of nucleic acid geometry on counterion association

Effect of a water molecule on the sugar puckering of uridine, 2'-deoxyuridine, and 2'-O-methyl uridine inserted in duplexes

We used high-level quantum mechanical calculations to determine the pucker (north type or south type) of various compounds: uridine, 2'-deoxyuridine, and 2'-O-methyl uridine. Although the dihedrals of the backbone are set close to their experimental values in double-stranded nucleic acids, calculations using density functional theory show that, in vacuo or in a continuum mimicking the dielectric properties of water, the south puckering conformations of uridine is favored. This contrasts with experimental data: most ribonucleosides inserted into a duplex have the north puckering. We show here that the north puckering is favored when an explicit water molecule is introduced into the calculation. The orientations of the 2' group and of the water molecule have implications for the prevalence of the north puckering. We studied several orientations of the water molecule binding uracil O2 and the 2' group and estimated the energy barriers in the path between the north-to-south conformations. The north puckering is more favored in 2'-OH than in 2'-OCH3 compounds in the presence of the explicit water molecule.

Conformational energetics of stable and metastable states formed by DNA triplet repeat oligonucleotides: implications for triplet expansion diseases

We have embedded the hexameric triplet repeats (CAG)(6) and (CTG)(6) between two (GC)(3) domains to produce two 30-mer hairpins with the sequences d[(GC)(3)(CAG)(6)(GC)(3)] and d[(GC)(3)(CTG)(6)(GC)(3)]. This construct reduces the conformational space available to these repetitive DNA sequences. We find that the (CAG)(6) and (CTG)(6) repeats form stable, ordered, single-stranded structures. These structures are stabilized at 62 degrees C by an average enthalpy per base of 1.38 kcal.mol(-1) for the CAG triplet and 2.87 kcal.mol(-1) for the CTG triplet, while being entropically destabilized by 3.50 cal.K(-1).mol(-1) for the CAG triplet and 7.6 cal.K(-1).mol(-1) for the CTG triplet. Remarkably, these values correspond, respectively, to 1/3 (for CAG) and 2/3 (for CTG) of the enthalpy and entropy per base values associated with Watson-Crick base pairs. We show that the presence of the loop structure kinetically inhibits duplex formation from the two complementary 30-mer hairpins, even though the duplex is the thermodynamically more stable state. Duplex formation, however, does occur at elevated temperatures. We propose that this thermally induced formation of a more stable duplex results from thermal disruption of the single-stranded order, thereby allowing the complementary domains to associate (perhaps via "kissing hairpins"). Our melting profiles show that, once duplex formation has occurred, the hairpin intermediate state cannot be reformed, consistent with our interpretation of kinetically trapped hairpin structures. The duplex formed by the two complementary oligonucleotides does not have any unusual optical or thermodynamic properties. By contrast, the very stable structures formed by the individual single-stranded triplet repeat sequences are thermally and thermodynamically unusual. We discuss this stable, triplet repeat, single-stranded structure and its interconversion with duplex in terms of triplet expansion diseases.

DAPI binding to the DNA minor groove: a continuum solvent analysis

A continuum solvent model based on the generalized Born (GB) or finite-difference Poisson-Boltzmann (FDPB) approaches has been employed to compare the binding of 4'-6-diamidine-2-phenyl indole (DAPI) to the minor groove of various DNA sequences. Qualitative agreement between the results of GB and FDPB approaches as well as between calculated and experimentally observed trends regarding the sequence specificity of DAPI binding to B-DNA was obtained. Calculated binding energies were decomposed into various contributions to solvation and DNA-ligand interaction. DNA conformational adaptation was found to make a favorable contribution to the calculated total interaction energy but did not change the DAPI binding affinity ranking of different DNA sequences. The calculations indicate that closed complex formation is mainly driven by nonpolar contributions and was found to be disfavored electrostatically due to a desolvation penalty that outbalances the attractive Coulomb interaction. The calculated penalty was larger for DAPI binding to GC-rich sequences compared with AT-rich target sequences and generally larger for the FDPB vs the GB continuum model. A radial interaction profile for DAPI at different distances from the DNA minor groove revealed an electrostatic energy minimum a few Angstroms farther away from the closed binding geometry. The calculated electrostatic interaction up to this distance is attractive and it may stabilize a nonspecific binding arrangement.

Sodium and chlorine ions as part of the DNA solvation shell

The distribution of sodium and chlorine ions around DNA is presented from two molecular dynamics simulations of the DNA fragment d(C(5)T(5)). (A(5)G(5)) in explicit solvent with 0.8 M additional NaCl salt. One simulation was carried out for 10 ns with the CHARMM force field that keeps the DNA structure close to A-DNA, the other for 12 ns with the AMBER force field that preferentially stabilizes B-DNA conformations (, Biophys. J. 75:134-149). From radial distributions of sodium and chlorine ions a primary ion shell is defined. The ion counts and residence times of ions within this shell are compared between conformations and with experiment. Ordered sodium ion sites were found in minor and major grooves around both A and B-DNA conformations. Changes in the surrounding hydration structure are analyzed and implications for the stabilization of A-DNA and B-DNA conformations are discussed.

Modeling high-resolution hydration patterns in correlation with DNA sequence and conformation

Hydration around the DNA fragment d(C5T5).(A5G5) is presented from two molecular dynamics simulations of 10 and 12 ns total simulation time. The DNA has been simulated as a flexible molecule with both the CHARMM and AMBER force fields in explicit solvent including counterions and 0.8 M additional NaCl salt. From the previous analysis of the DNA structure B-DNA conformations were found with the AMBER force-field and A-DNA conformations with CHARMM parameters. High-resolution hydration patterns are compared between the two conformations and between C.G and T.A base-pairs from the homopolymeric parts of the simulated sequence. Crystallographic results from a statistical analysis of hydration sites around DNA crystal structures compare very well with the simulation results. Differences between the crystal sites and our data are explained by variations in conformation, sequence, and limitations in the resolution of water sites by crystal diffraction. Hydration layers are defined from radial distribution functions and compared with experimental results. Excellent agreement is found when the measured experimental quantities are compared with the equivalent distribution of water molecules in the first hydration shell. The number of water molecules bound to DNA was found smaller around T.A base-pairs and around A-DNA as compared to B-DNA. This is partially offset by a larger number of water molecules in hydrophobic contact with DNA around T.A base-pairs and around A-DNA. The numbers of water molecules in minor and major grooves have been correlated with helical roll, twist, and inclination angles. The data more fully explain the observed B-->A transition at low humidity.