Influence of the dynamic positions of cations on the structure of the DNA minor groove: sequence-dependent effects

Binding of Phosphate Species to Ca2+ and Mg2+ in Aqueous Solution

Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H2PO4-, HPO42-, PO43- coordinating with Ca2+ and Mg2+. First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca2+ and Mg2+ generally show similar trends at longer metal-phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation-anion interactions.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

Changes in Microenvironment Modulate the B- to A-DNA Transition

B- to A-DNA transition is known to be sensitive to the macroscopic properties of the solution, such as salt and ethanol concentrations. Microenvironmental effects on DNA conformational transition have been broadly studied. Providing an intuitive picture of how DNA responds to environmental changes is, however, still needed. Analyzing the chemical equilibrium of B-to-A DNA transition at critical concentrations, employing explicit-solvent simulations, is envisioned to help understand such microenvironmental effects. In the present study, free-energy calculations characterizing the B- to A-DNA transition and the distribution of cations were carried out in solvents with different ethanol concentrations. With the addition of ethanol, the most stable structure of DNA changes from the B- to A-form, in agreement with previous experimental observation. In 60% ethanol, a chemical equilibrium is found, showing reversible transition between B- and A-DNA. Analysis of the microenvironment around DNA suggests that with the increase of ethanol concentration, the cations exhibit a significant tendency to move toward the backbone, and mobility of water molecules around the major groove and backbone decreases gradually, leading eventually to a B-to-A transition. The present results provide a free-energy view of DNA microenvironment and of the role of cation motion in the conformational transition.

Molecular Dynamics Simulations of Dna and Its Complexes

This article describes how classical molecular simulation methods are being used to gain a molecular-level understanding of the interaction mechanisms responsible for DNA–ligand recognition, and that govern the response of DNA to ligand binding. Case studies using a variety of different ligands—including small pharmaceutical drugs, proteins and lipids—are used to illustrate the power of modern molecular dynamics simulation methods for understanding how we may control the function and structure of DNA.

New Insights into the Interactions of a DNA Oligonucleotide with mPEGylated-PAMAM by Circular Dichroism and Solution NMR

Dendrimers are well-defined, highly branched, synthetic three-dimensional molecules with a large number of reactive end groups. PAMAM dendrimers form stable complexes with DNA chemistries and constitute an important class of nonviral, cationic vectors in gene delivery. The aim of this study is to examine the interactions of a 12 bp DNA oligonucletide with PAMAM-G2 and mPEG- b-PAMAM-G3 having eight surface amine groups under physiological conditions, using constant DNA concentration but varying dendrimer concentration. 1D ³¹P NMR, 2D NOESY, and CD spectroscopic methods were employed to investigate the interactions between the dendrimer and the DNA. The CD experiments carried out with a constant DNA concentration of 25 μM and dendrimer concentrations from 0 to 100 μM indicated minimal change to the chirality of the DNA for both types of dendrimers. While the PAMAM-G2 dendrimer caused aggregation of the majority of the DNA, the 2D NMR data of the DNA with an mPEG- b-PAMAM-G3 dendrimer indicated general broadening of the 1D ³¹P peaks from the DNA phosphates, a small number of ¹H chemical shift perturbations (CSPs), and reduction of specific ¹H-¹H NOE intensities. These data suggest there is minimal structural alteration of the DNA in the complex and indicate preferential binding of the dendrimer to the central AATT region of the DNA sequence. The results herein are the first such results demonstrating a soluble DNA complex with the mPEG- b-PAMAM-G3 dendrimer analyzed by multidimensional NMR.

Interaction of Na+, K+, Mg2+ and Ca2+ counter cations with RNA

Data on the location of alkaline and alkaline earth ions at RNA from crystallography, spectroscopy and computational modeling are reviewed.

Exploring ion induced folding of a single-stranded DNA oligomer from molecular simulation studies

One crucial issue in DNA hydration is the effect of salts on its conformational features. This has relevance in biology as cations present in the cellular environment shield the negative charges on the DNA backbone, thereby reducing the repulsive force between them. By screening the negative charges along the backbone, cations stabilize the folded structure of DNA. To study the effect of the added salt on single-stranded DNA (ss-DNA) conformations, we have performed room temperature molecular dynamics simulations of an aqueous solution containing the ss-DNA dodecamer with the 5'-CGCGAATTCGCG-3' sequence in the presence of 0.2, 0.5, and 0.8 M NaCl. Our calculations reveal that in the presence of the salt, the DNA molecule forms more collapsed coil-like conformations due to the screening of negative charges along the backbone. Additionally, we demonstrated that the formation of an octahedral inner-sphere complex by the strongly bound ion plays an important role in the stabilization of such folded conformation of DNA. Importantly, it is found that ion-DNA interactions can also explain the formation of non-sequential base stackings with longer lifetimes. Such non-sequential base stackings further stabilize the collapsed coil-like folded form of the DNA oligomer.

Recent findings on the mechanisms involved in tenofovir resistance

Since its approval for clinical use in 2001, tenofovir (TFV) has become one of the most frequently prescribed nucleotide analogues used in combination with other antiretroviral agents against HIV-1 infection. Although reverse transcriptase inhibitors (RTIs) including TFV have been shown to be highly potent with reasonable safety profiles in the clinic, drug resistance hinders the effectiveness of current therapies and even causes treatment failure. Therefore, understanding the resistance mechanisms of RT and exploring the potential antiviral synergy between the different RTIs in combination therapies against the resistance mechanisms would greatly improve the long-term efficacy of existing and future regimens. We have studied the pyrophosphorolytic removal of TFV, a major resistance mechanism that RT utilizes, from two different viral sequences and observed interesting outcomes associated with the sequence context. Furthermore, addition of efavirenz, a non-nucleoside RTI, inhibits this removal process confirming the synergistic antiviral effects. This article highlights our recently published work on the viral sequence context contributing to the study of anti-HIV drug resistance in conjunction with the benefits of combining various RTIs that may have been neglected previously.

Free energy of binding of cationic metal complexes to AuNPs through electron-transfer processes

AuNP effects on the oxidation reaction of [Ru(NH3)5pz](2+) with S2O8(2-) were studied. Experimental results show the effects produced by gold nanoparticles of different sizes, which are then discussed by using an approach based on the two-state model. Changes in the observed reactivity are explained by a change in the degree of association of the reactants to the nanoclusters, which depends on the negative surface potential of the gold nanoparticles. TEM and zeta potential measurements show that this potential determines the binding strength of one of the reactants ([Ru(NH3)5pz](2+)) to the citrate gold surface. The number of binding sites on a citrate nanoparticle receptor has also been determined. The experimental results confirm that an electron transfer reaction can be used as a probe for the determination of the free energy of binding of cationic metal complexes to gold nanoparticles.

Roles of the amino group of purine bases in the thermodynamic stability of DNA base pairing

The energetic aspects of hydrogen-bonded base-pair interactions are important for the design of functional nucleotide analogs and for practical applications of oligonucleotides. The present study investigated the contribution of the 2-amino group of DNA purine bases to the thermodynamic stability of oligonucleotide duplexes under different salt and solvent conditions, using 2'-deoxyriboinosine (I) and 2'-deoxyribo-2,6-diaminopurine (D) as non-canonical nucleotides. The stability of DNA duplexes was changed by substitution of a single base pair in the following order: G•C > D•T ≈ I•C > A•T > G•T > I•T. The apparent stabilization energy due to the presence of the 2-amino group of G and D varied depending on the salt concentration, and decreased in the water-ethanol mixed solvent. The effects of salt concentration on the thermodynamics of DNA duplexes were found to be partially sequence-dependent, and the 2-amino group of the purine bases might have an influence on the binding of ions to DNA through the formation of a stable base-paired structure. Our results also showed that physiological salt conditions were energetically favorable for complementary base recognition, and conversely, low salt concentration media and ethanol-containing solvents were effective for low stringency oligonucleotide hybridization, in the context of conditions employed in this study.

Balancing the interactions of ions, water, and DNA in the Drude polarizable force field

Recently we presented a first-generation all-atom Drude polarizable force field for DNA based on the classical Drude oscillator model, focusing on optimization of key dihedral angles followed by extensive validation of the force field parameters. Presently, we describe the procedure for balancing the electrostatic interactions between ions, water, and DNA as required for development of the Drude force field for DNA. The proper balance of these interactions is shown to impact DNA stability and subtler conformational properties, including the conformational equilibrium between the BI and BII states, and the A and B forms of DNA. The parametrization efforts were simultaneously guided by gas-phase quantum mechanics (QM) data on small model compounds and condensed-phase experimental data on the hydration and osmotic properties of biologically relevant ions and their solutions, as well as theoretical predictions for ionic distribution around DNA oligomer. In addition, fine-tuning of the internal base parameters was performed to obtain the final DNA model. Notably, the Drude model is shown to more accurately reproduce counterion condensation theory predictions of DNA charge neutralization by the condensed ions as compared to the CHARMM36 additive DNA force field, indicating an improved physical description of the forces dictating the ionic solvation of DNA due to the explicit treatment of electronic polarizability. In combination with the polarizable DNA force field, the availability of Drude polarizable parameters for proteins, lipids, and carbohydrates will allow for simulation studies of heterogeneous biological systems.

Ab initio molecular dynamics of Na⁺ and Mg²⁺ countercations at the backbone of RNA in water solution

The interactions between sodium or magnesium ions and phosphate groups of the RNA backbone represented as dinucleotide fragments in water solution have been studied using ab initio Born-Oppenheimer molecular dynamics. All systems have been simulated at 300 and 320 K. Sodium ions have mobility higher than that of the magnesium ions and readily change their position with respect to the phosphate groups, from directly bonded to completely solvated state, with a rough estimate of the lifetime of bonded Na(+) of about 20-30 ps. The coordination number of the sodium ions frequently changes in irregular intervals ranging from several femtoseconds to about 10 ps with the most frequently encountered coordination number five, followed by six. The magnesium ion is stable both as directly bonded to an oxygen atom from the phosphate group and completely solvated by water. In both states the Mg(2+) ion has exactly six oxygen atoms in the first coordination shell; moreover, during the whole simulation of more than 100 ps no exchange of ligand in the first coordination shells has been observed. Solvation of the terminal phosphate oxygen atoms by water molecules forming hydrogen bonds in different locations of the ions is also discussed. The stability of the system containing sodium ions essentially does not depend on the position of the ions with respect to the phosphate groups.

The ABCs of molecular dynamics simulations on B-DNA, circa 2012

This article provides a retrospective on the ABC initiative in the area of all-atom molecular dynamics (MD) simulations including explicit solvent on all tetranucleotide steps of duplex B-form DNA duplex, ca. 2012. The ABC consortium has completed two phases of simulations, the most current being a set of 50-100 trajectories based on the AMBER ff99 force field together with the parmbsc0 modification. Some general perspectives on the field of MD on DNA and sequence effects on DNA structure are provided, followed by an overview our MD results, including a detailed comparison of the ff99/parmbsc0 results with crystal and NMR structures available for d(CGCGAATTCGCG). Some projects inspired by or related to the ABC initiative and database are also reviewed, including methods for the trajectory analyses, informatics of dealing with the large database of results, compressions of trajectories for efficacy of distribution, DNA solvation by water and ions, parameterization of coarse-grained models with applications and gene finding and genome annotation.

Physiologically relevant concentrations of NaCl and KCl increase DNA photocleavage by an N-substituted 9-aminomethylanthracene dye

This paper describes the synthesis of a new 9-aminomethylanthracene dye N-substituted with a pyridinylpolyamine side chain (4). The effects of NaCl and KCl on anthracene/DNA interactions were then studied, with the goal of simulating the conditions of high ionic strength that a DNA photosensitizer might encounter in the cell nucleus (~150 mM of NaCl and 260 mM of KCl). As exemplified by methylene blue (5), the expected effect of increasing ionic strength is to decrease DNA binding and photocleavage yields. In contrast, the addition of 150 mM of NaCl in combination with 260 mM of KCl to photocleavage reactions containing micromolar concentrations of 4 triggers the conversion of supercoiled, nicked, and linear forms of pUC19 plasmid into a highly degraded band of DNA fragments (350 nm hν, pH 7.0). Circular dichroism spectra point to a correlation between salt-induced unwinding of the DNA helix and the increase in DNA photocleavage yields. The results of circular dichroism, UV-vis absorption, fluorescence emission, thermal denaturation, and photocleavage inhibition experiments suggest that the combination of salts causes a change in the DNA binding mode of 4 from intercalation to an external interaction. This in turn leads to an increase in the anthracene-sensitized production of DNA-damaging reactive oxygen species.

Induced topological changes in DNA complexes: influence of DNA sequences and small molecule structures

Heterocyclic diamidines are compounds with antiparasitic properties that target the minor groove of kinetoplast DNA. The mechanism of action of these compounds is unknown, but topological changes to DNA structures are likely to be involved. In this study, we have developed a polyacrylamide gel electrophoresis-based screening method to determine topological effects of heterocyclic diamidines on four minor groove target sequences: AAAAA, TTTAA, AAATT and ATATA. The AAAAA and AAATT sequences have the largest intrinsic bend, whereas the TTTAA and ATATA sequences are relatively straight. The changes caused by binding of the compounds are sequence dependent, but generally the topological effects on AAAAA and AAATT are similar as are the effects on TTTAA and ATATA. A total of 13 compounds with a variety of structural differences were evaluated for topological changes to DNA. All compounds decrease the mobility of the ATATA sequence that is consistent with decreased minor groove width and bending of the relatively straight DNA into the minor groove. Similar, but generally smaller, effects are seen with TTTAA. The intrinsically bent AAAAA and AAATT sequences, which have more narrow minor grooves, have smaller mobility changes on binding that are consistent with increased or decreased bending depending on compound structure.

On the molecular basis of uracil recognition in DNA: comparative study of T-A versus U-A structure, dynamics and open base pair kinetics

Uracil (U) can be found in DNA as a mismatch paired either to adenine (A) or to guanine (G). Removal of U from DNA is performed by a class of enzymes known as uracil–DNA–glycosylases (UDG). Recent studies suggest that recognition of U–A and U–G mismatches by UDG takes place via an extra-helical mechanism. In this work, we use molecular dynamics simulations to analyze the structure, dynamics and open base pair kinetics of U–A base pairs relative to their natural T–A counterpart in 12 dodecamers. Our results show that the presence of U does not alter the local conformation of B-DNA. Breathing dynamics and base pair closing kinetics are only weakly dependent on the presence of U versus T, with open T–A and U–A pairs lifetimes in the nanosecond timescale. Additionally, we observed spontaneous base flipping in U–A pairs. We analyze the structure and dynamics for this event and compare the results to available crystallographic data of open base pair conformations. Our results are in agreement with both structural and kinetic data derived from NMR imino proton exchange measurements, providing the first detailed description at the molecular level of elusive events such as spontaneous base pair opening and flipping in mismatched U–A sequences in DNA. Based on these results, we propose that base pair flipping can occur spontaneously at room temperature via a 3-step mechanism with an open base pair intermediate. Implications for the molecular basis of U recognition by UDG are discussed.

Factors affecting the yields of C1' and C5' oxidation products in radiation-damaged DNA: the indirect effect

This study reports the effects of denaturation and deoxygenation on radiation-induced formation of 2-deoxyribonolactone (2-dL) and 5'-aldehyde (5'-Ald) lesions in highly polymerized DNA. The radiation-chemical yields of 2-dL were determined through quantification of its dephosphorylation product 5-methylenefuranone (5MF). The formation of 5'-Ald was monitored qualitatively through the release of furfural (Fur) under the same conditions. The yields of 2-dL were found to be 7.3 ± 0.3 nmol J(-1), or about 18% of the yield of free base release measured in the same samples. Denaturation increased the efficiency of 2-dL formation approximately twofold while deoxygenation resulted in a fourfold decrease. The release of Fur is about twofold lower than that of 5MF in aerated native DNA samples and is further reduced by denaturation of the DNA. Unlike 5MF, the formation of Fur requires the presence of molecular oxygen, which is consistent with peroxyl radical-mediated oxidation of C5' radicals into 5'-Ald. In contrast, the existence of an oxygen-independent pathway of 2-dL formation suggests that C1' sugar radicals can also be oxidized by radiation-produced oxidizing intermediates such as electron-loss centers on guanines.

Minor groove orientation of the KWKK peptide tethered via the N-terminal amine to the acrolein-derived 1,N2-gamma-hydroxypropanodeoxyguanosine lesion with a trimethylene linkage

DNA−protein conjugates are potentially repaired via proteolytic digestion to DNA−peptide conjugates. The latter have been modeled with the amino-terminal lysine of the peptide KWKK conjugated via a trimethylene linkage to the N²-dG amine positioned in 5′-d(GCTAGCXAGTCC)-3′·5′-d(GGACTCGCTAGC)-3′ (X = N²-dG−trimethylene link−KWKK). This linkage is a surrogate for the reversible linkage formed by the γ-OH-1,N²-propanodeoxyguanosine (γ-OH-PdG) adduct. This conjugated KWKK stabilizes the DNA. Amino acids K²⁶, W²⁷, K²⁸, and K²⁹ are in the minor groove. The W²⁷ indolyl group does not intercalate into the DNA. The G⁷N² amine and the K²⁶ N-terminal amine nitrogens are in the trans configuration with respect to the C_α or C_γ of the trimethylene tether, respectively. The structure of this DNA−KWKK conjugate is discussed in the context of its biological processing.

DNA curvature and flexibility in vitro and in vivo

It has been more than 50 years since the elucidation of the structure of double-helical DNA. Despite active research and progress in DNA biology and biochemistry, much remains to be learned in the field of DNA biophysics. Predicting the sequence-dependent curvature and flexibility of DNA is difficult. Applicability of the conventional worm-like chain polymer model of DNA has been challenged. The fundamental forces responsible for the remarkable resistance of DNA to bending and twisting remain controversial. The apparent 'softening' of DNA measured in vivo in the presence of kinking proteins and superhelical strain is incompletely understood. New methods and insights are being applied to these problems. This review places current work on DNA biophysics in historical context and illustrates the ongoing interplay between theory and experiment in this exciting field.

Monovalent cation binding in the minor groove of DNA A-tracts

The binding of five different monovalent cations to DNA oligomers containing A-tracts, runs of four or more contiguous adenine residues, has been assessed by capillary electrophoresis, using the Replacement Ion method. In this method, a nonbinding cation in the background electrolyte is gradually replaced by a binding cation, keeping the ionic strength of the solution constant. Monovalent cation binding reduces the effective charge of an A-tract-containing oligomer, decreasing its free solution mobility. The cations bind in the A-tract minor groove, because the binding site can be blocked by the minor groove binding drug netropsin. Li(+), NH(4)(+), and Tris(+) ions bind to A-tracts with similar affinities; the binding of Na(+) ions is weaker, and K(+) ion binding is highly variable. Each A-tract appears to bind one monovalent cation upon saturation of the binding site(s). For a given cation, the apparent dissociation constants depend on A-tract sequence and orientation, but not on the phasing of the A-tracts with respect to the helix repeat. Differences in the cooperativity of binding of the various cations to A-tracts with different sequences suggest that monovalent cation binding may be coupled with a conformational transition leading to the formation of the characteristic narrow minor groove A-tract structure.

Molecular dynamics of a kappaB DNA element: base flipping via cross-strand intercalative stacking in a microsecond-scale simulation

The sequence-dependent structural variability and conformational dynamics of DNA play pivotal roles in many biological milieus, such as in the site-specific binding of transcription factors to target regulatory elements. To better understand DNA structure, function, and dynamics in general, and protein···DNA recognition in the ‘κB’ family of genetic regulatory elements in particular, we performed molecular dynamics simulations of a 20-bp DNA encompassing a cognate κB site recognized by the proto-oncogenic ‘c-Rel’ subfamily of NF-κB transcription factors. Simulations of the κB DNA in explicit water were extended to microsecond duration, providing a broad, atomically detailed glimpse into the structural and dynamical behavior of double helical DNA over many timescales. Of particular note, novel (and structurally plausible) conformations of DNA developed only at the long times sampled in this simulation—including a peculiar state arising at ≈0.7 μs and characterized by cross-strand intercalative stacking of nucleotides within a longitudinally sheared base pair, followed (at ≈1 μs) by spontaneous base flipping of a neighboring thymine within the A-rich duplex. Results and predictions from the microsecond-scale simulation include implications for a dynamical NF-κB recognition motif, and are amenable to testing and further exploration via specific experimental approaches that are suggested herein.

A solid-state 23Na NMR study of monovalent cation binding to double-stranded DNA at low relative humidity

We report a solid-state (23)Na NMR study of monovalent cation (Li(+), Na(+), K(+), Rb(+), Cs(+) and NH(4) (+)) binding to double-stranded calf thymus DNA (CT DNA) at low relative humidity, ca 0-10%. Results from (23)Na--(31)P rotational echo double resonance (REDOR) NMR experiments firmly establish that, at low relative humidity, monovalent cations are directly bound to the phosphate group of CT DNA and are partially dehydrated. On the basis of solid-state (23)Na NMR titration experiments, we obtain quantitative thermodynamic parameters concerning the cation-binding affinity for the phosphate group of CT DNA. The free energy difference (DeltaG degrees ) between M(+) and Na(+) ions is as follows: Li(+) (-1.0 kcal mol(-1)), K(+) (7.2 kcal mol(-1)), NH(4) (+) (1.0 kcal mol(-1)), Rb(+) (4.5 kcal mol(-1)) and Cs(+) (1.5 kcal mol(-1)). These results suggest that, at low relative humidity, the binding affinity of monovalent cations for the phosphate group of CT DNA follows the order: Li(+) > Na(+) > NH(4) (+) > Cs(+) > Rb(+) > K(+). This sequence is drastically different from that observed for CT DNA in solution. This discrepancy is attributed to the different modes of cation binding in dry and wet states of DNA. In the wet state of DNA, cations are fully hydrated. Our results suggest that the free energy balance between direct cation-phosphate contact and dehydration interactions is important. The reported experimental results on relative ion-binding affinity for the DNA backbone may be used for testing theoretical treatment of cation-phosphate interactions in DNA.

Modifying the N-terminus of polyamides: PyImPyIm has improved sequence specificity over f-ImPyIm

Seven N-terminus modified derivatives of a previously published minor-groove binding polyamide (f-ImPyIm, 1) were synthesized and the biochemical and biophysical chemistry evaluated. These compounds were synthesized with the aim of attaining a higher level of sequence selectivity over f-ImPyIm (1), a previously published strong minor-groove binder. Two compounds possessing a furan or a benzofuran moiety at the N-terminus showed a footprint of 0.5microM at the cognate ACGCGT site (determined by DNase I footprinting); however, the specificity of these compounds was not improved. In contrast, PyImPyIm (4) produced a footprint of 0.5microM but showed a superior specificity using the same technique. When evaluated by thermal melting experiments and circular dichroism using ACGCGT and the non-cognate AAATTT sequence, all compounds were shown to bind in the minor-groove of DNA and stabilize the cognate sequence much better than the non-cognate (except for the non-amido-compound that did not bind either sequence, as expected). PyImPyIm (4) was interesting as the DeltaT(m) for this compound was only 4 degrees C but the footprint was very selective. No binding was observed for this compound with a third DNA (non-cognate, ACCGGT). ITC studies on compound 4 showed exothermic binding with ACGCGT and no heat change was observed for titrating the compound to the other two DNA sequences. The heat capacity (DeltaC(p)) of the PIPI/ACGCGT complex calculated from the hydrophobic interactions and SASA calculations was comparable to the experimental value obtained from ITC (-146calmol(-1)K(-1)). SPR results provided confirmation of the sequence specificity of PyImPyIm (4), with a K(eq) value determined to be 7.1x10(6) M(-1) for the cognate sequence and no observable binding to AAATTT and ACCGGT. Molecular dynamic simulations affirmed that PyImPyIm (4) binds as a dimer in an overlapped conformation, and it fits snugly in the minor-groove of the ACGCGT oligonucleotide. PyImPyIm (4) is an especially interesting molecule, because although the binding affinity is slightly reduced, the specificity with respect to f-ImPyIm (1) is significantly improved.

Antiparasitic compounds that target DNA

Designed, synthetic heterocyclic diamidines have excellent activity against eukaryotic parasites that cause diseases such as sleeping sickness and leishmania and adversely affect millions of people each year. The most active compounds bind specifically and strongly in the DNA minor groove at AT sequences. The compounds enter parasite cells rapidly and appear first in the kinetoplast that contains the mitochondrial DNA of the parasite. With time the compounds are also generally seen in the cell nucleus but are not significantly observed in the cytoplasm. The kinetoplast decays over time and disappears from the mitochondria of treated cells. At this point the compounds begin to be observed in other regions of the cell, such as the acidocalcisomes. The cells typically die in 24-48h after treatment. Active compounds appear to selectively target extended AT sequences and induce changes in kinetoplast DNA minicircles that cause a synergistic destruction of the catenated kinetoplast DNA network and cell death.

Sequence and length dependent thermodynamic differences in heterocyclic diamidine interactions at AT base pairs in the DNA minor groove

With the goal of developing a better understanding of the antiparasitic biological action of DB75, we have evaluated its interaction with duplex alternating and nonalternating sequence AT polymers and oligomers. These DNAs provide an important pair of sequences in a detailed thermodynamic analysis of variations in interaction of DB75 with AT sites. The results for DB75 binding to the alternating and nonalternating AT sequences are quite different at the fundamental thermodynamic level. Although the Gibbs energies are similar, the enthalpies for DB75 binding with poly(dA).poly(dT) and poly(dA-dT).poly(dA-dT) are +3.1 and -4.5 kcal/mol, respectively, while the binding entropies are 41.7 and 15.2 cal/mol.K, respectively. The underlying thermodynamics of binding to AT sites in the minor groove plays a key role in the recognition process. It was also observed that DB75 binding with poly(dA).poly(dT) can induce T.A.T triplet formation and the compound binds strongly to the dT.dA.dT triplex.

A review of the role of the sequence-dependent electrostatic landscape in DNA alkylation patterns

Alkylating agents, including environmental and endogenous carcinogens and DNA targeting antineoplastic agents, that adduct DNA via intermediates with significant cationic charge show a sequence selectively in their covalent bonding to nucleobases. The resulting patterns of alkylation eventually contribute to the agent-dependent distributions and types of mutations. The origin of the regioselective modification of DNA by electrophiles has been attributed to steric and/or electronic factors, but attempts to mechanistically model and predict alkylation patterns have had limited success. In this review, we present data consistent with the role of the intrinsic sequence-dependent electrostatic landscape (SDEL) in DNA that modulates the equilibrium binding of cations and the bonding of reactive charged alkylating agents to atoms that line the floor of the major groove of DNA.

Conformations and Dynamics of the Phosphodiester Backbone of a DNA Fragment That Bears a Strong Topoisomerase II Cleavage Site

The dynamics of the DNA phosphodiester backbone conformations have been studied for a strong topoisomerase II cleavage site (site 22) using molecular dynamics simulations in explicit water and in the presence of sodium ions. We investigated the backbone motions and more particularly the BI/BII transitions involving the epsilon and zeta angles. The consensus cleavage site is adjacent to the phosphate which shows the most important phosphodiester backbone flexibility in the sequence. We infer that these latter properties could be responsible for the preferential cleavage at this site possibly through the perturbation of the cleavage/ligation activities of the topoisomerase II. More generally, the steps pur-pur and pyr-pur are those presenting the highest BII contents. Relations are observed between the backbone phosphodiester BI/BII transitions and the flexibility of the deoxyribose sugar and the helical parameters such as roll. The roll is sequence dependent when the related phosphate is in the BI form, whereas this appears not to be true when it is in the BII form. The BI/BII transitions are associated with water migration, and new relations are observed with counterions. Indeed, it is observed that a strong coupling exists between the BII form and the presence of sodium ions near the adjacent sugar deoxyribose. The presence of sodium ions in the O4' surroundings or their binding could assist the BI to BII transition by furnishing energy. The implications of these new findings and, namely, their importance in the context of the sequence-dependent behavior of BI/BII transitions will be investigated in future studies.

The DNA structure responds differently to physiological concentrations of K(+) or Na(+)

The influence of monovalent cations on DNA conformation and readout is an open question. This NMR study of DNA with either Na(+) or K(+) at physiological concentrations shows that the nature of the cation affects the (31)P chemical shifts (deltaP) and the sequential distances H2'(i)-H6/8(i+1), H2"(i)-H6/8(i+1), and H6/8(i)-H6/8(i+1). The deltaP and distance variations ascertain that the nature of the cation affects the DNA overall structure, i.e. both the conformational equilibria between the backbone BI (epsilon-zeta <0 degrees ) and BII (epsilon-zeta >0 degrees ) states and the helical parameters, via their strong mechanical coupling. These results reveal that Na(+) and K(+) interactions with DNA are different and sequence-dependent. These ions modulate the overall intrinsic properties of DNA, and possibly its packaging and readout.

Physical and Structural Basis for the Strong Interactions of the -ImPy- Central Pairing Motif in the Polyamide f-ImPyIm

The polyamide f-ImPyIm has a higher affinity for its cognate DNA than either the parent analogue, distamycin A (10-fold), or the structural isomer, f-PyImIm (250-fold), has for its respective cognate DNA sequence. These findings have led to the formulation of a two-letter polyamide "language" in which the -ImPy- central pairings associate more strongly with Watson-Crick DNA than -PyPy-, -PyIm-, and -ImIm-. Herein, we further characterize f-ImPyIm and f-PyImIm, and we report thermodynamic and structural differences between -ImPy- (f-ImPyIm) and -PyIm- (f-PyImIm) central pairings. DNase I footprinting studies confirmed that f-ImPyIm is a stronger binder than distamycin A and f-PyImIm and that f-ImPyIm preferentially binds CGCG over multiple competing sequences. The difference in the binding of f-ImPyIm and f-PyImIm to their cognate sequences was supported by the Na(+)-dependent nature of DNA melting studies, in which significantly higher Na(+) concentrations were needed to match the ability of f-ImPyIm to stabilize CGCG with that of f-PyImIm stabilizing CCGG. The selectivity of f-ImPyIm beyond the four-base CGCG recognition site was tested by circular dichroism and isothermal titration microcalorimetry, which shows that f-ImPyIm has marginal selectivity for (A.T)CGCG(A.T) over (G.C)CGCG(G.C). In addition, changes adjacent to this 6 bp binding site do not affect f-ImPyIm affinity. Calorimetric studies revealed that binding of f-ImPyIm, f-PyImIm, and distamycin A to their respective hairpin cognate sequences is exothermic; however, changes in enthalpy, entropy, and heat capacity (DeltaC(p)) contribute differently to formation of the 2:1 complexes for each triamide. Experimental and theoretical determinations of DeltaC(p) for binding of f-ImPyIm to CGCG were in good agreement (-142 and -177 cal mol(-)(1) K(-)(1), respectively). (1)H NMR of f-ImPyIm and f-PyImIm complexed with their respective cognate DNAs confirmed positively cooperative formation of distinct 2:1 complexes. The NMR results also showed that these triamides bind in the DNA minor groove and that the oligonucleotide retains the B-form conformation. Using minimal distance restraints from the NMR experiments, molecular modeling and dynamics were used to illustrate the structural complementarity between f-ImPyIm and CGCG. Collectively, the NMR and ITC experiments show that formation of the 2:1 f-ImPyIm-CGCG complex achieves a structure more ordered and more thermodynamically favored than the structure of the 2:1 f-PyImIm-CCGG complex.

Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules

Many interesting dynamic properties of biological molecules cannot be simulated directly using molecular dynamics because of nanosecond time scale limitations. These systems are trapped in potential energy minima with high free energy barriers for large numbers of computational steps. The dynamic evolution of many molecular systems occurs through a series of rare events as the system moves from one potential energy basin to another. Therefore, we have proposed a robust bias potential function that can be used in an efficient accelerated molecular dynamics approach to simulate the transition of high energy barriers without any advance knowledge of the location of either the potential energy wells or saddle points. In this method, the potential energy landscape is altered by adding a bias potential to the true potential such that the escape rates from potential wells are enhanced, which accelerates and extends the time scale in molecular dynamics simulations. Our definition of the bias potential echoes the underlying shape of the potential energy landscape on the modified surface, thus allowing for the potential energy minima to be well defined, and hence properly sampled during the simulation. We have shown that our approach, which can be extended to biomolecules, samples the conformational space more efficiently than normal molecular dynamics simulations, and converges to the correct canonical distribution.

Computer modeling demonstrates that electrostatic attraction of nucleosomal DNA is mediated by histone tails

We conducted molecular dynamics computer simulations of charged histone tail-DNA interactions in systems mimicking nucleosome core particles (NCP) . In a coarse-grained model, the NCP is modeled as a negatively charged spherical particle with flexible polycationic histone tails attached to it in a dielectric continuum with explicit mobile counterions and added salt. The size, charge, and distribution of the tails relative to the core were built to mimick real NCP. In this way, we incorporate attractive ion-ion correlation effects due to fluctuations in the ion cloud and the attractive entropic and energetic tail-bridging effects. In agreement with experimental data, increase of monovalent salt content from salt-free to physiological concentration leads to the formation of NCP aggregates; likewise, in the presence of MgCl2, the NCPs form condensed systems via histone-tail bridging and accumulation of counterions. More detailed mechanisms of the histone tail-DNA interactions and dynamics have been obtained from all-atom molecular dynamics simulations (including water), comprising three DNA 22-mers and 14 short fragments of the H4 histone tail (amino acids 5-12) carrying three positive charges on lysine+ interacting with DNA. We found correlation of the DNA-DNA distance with the presence and association of the histone tail between the DNA molecules.

Similarities and differences in interaction of K+ and Na+ with condensed ordered DNA. A molecular dynamics computer simulation study

Four 20 ns molecular dynamics simulations have been performed with two counterions, K⁺ or Na⁺, at two water contents, 15 or 20 H₂O per nucleotide. A hexagonal simulation cell comprised of three identical DNA decamers [d(5′-ATGCAGTCAG) × d(5′-TGACTGCATC)] with periodic boundary condition along the DNA helix was used. The simulation setup mimics the DNA state in oriented DNA fibers or in crystals of DNA oligomers. Variation of counterion nature and water content do not alter averaged DNA structure. K⁺ and Na⁺ binding to DNA are different. K⁺ binds to the electronegative sites of DNA bases in the major and the minor grooves, while Na⁺ interacts preferentially with the phosphate groups. Increase of water causes a shift of both K⁺ and Na⁺ from the first hydration shell of O1P/O2P and of the DNA bases in the minor groove with lesser influence for the cation binding to the bases in the major groove. Mobility of both water and cations in the K–DNA systems is faster than in the Na–DNA systems: Na⁺ organizes and immobilizes water structure around itself and near DNA while for K⁺ water is less organized and more dynamic.

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps

Molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide basepair steps are reported. The objective is to obtain the calculated dynamical structure for at least two copies of each case, use the results to examine issues with regard to convergence and dynamical stability of MD on DNA, and determine the significance of sequence context effects on all unique dinucleotide steps. This information is essential to understand sequence effects on DNA structure and has implications on diverse problems in the structural biology of DNA. Calculations were carried out on the 136 cases embedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All simulations were carried out using a well-defined state-of-the-art MD protocol, the AMBER suite of programs, and the parm94 force field. In a previous article (Beveridge et al. 2004. Biophysical Journal. 87:3799-3813), the research design, details of the simulation protocol, and informatics issues were described. Preliminary results from 15 ns MD trajectories were presented for the d(CpG) step in all 10 unique sequence contexts. The results indicated the sequence context effects to be small for this step, but revealed that MD on DNA at this length of trajectory is subject to surprisingly persistent cooperative transitions of the sugar-phosphate backbone torsion angles alpha and gamma. In this article, we report detailed analysis of the entire trajectory database and occurrence of various conformational substates and its impact on studies of context effects. The analysis reveals a possible direct correspondence between the sequence-dependent dynamical tendencies of DNA structure and the tendency to undergo transitions that "trap" them in nonstandard conformational substates. The difference in mean of the observed basepair step helicoidal parameter distribution with different flanking sequence sometimes differs by as much as one standard deviation, indicating that the extent of sequence effects could be significant. The observations reveal that the impact of a flexible dinucleotide such as CpG could extend beyond the immediate basepair neighbors. The results in general provide new insight into MD on DNA and the sequence-dependent dynamical structural characteristics of DNA.

Nuclear aggregates of polyamines are supramolecular structures that play a crucial role in genomic DNA protection and conformation

In a previous study we showed that natural polyamines interact in the nuclear environment with phosphate groups to form molecular aggregates [nuclear aggregates of polyamines (NAPs)] with estimated molecular mass values of 8000, 4800 and 1000 Da. NAPs were found to interact with genomic DNA, influence its conformation and interfere with the action of nucleases. In the present work, we demonstrated that NAPs protect naked genomic DNA from DNase I, whereas natural polyamines (spermine, spermidine and putrescine) fail to do so. In the context of DNA protection, NAPs induced noticeable changes in DNA conformation, which were revealed by temperature-dependent modifications of DNA electrophoretic properties. In addition, we presented, for NAPs, a structural model of polyamine aggregation into macropolycyclic compounds. We believe that NAPs are the sole biological forms by which polyamines efficiently protect genomic DNA against DNase I, while maintaining its dynamic structure.

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds specific sites deep inside the groove it can reside there for very long periods of time. According to our simulation the presence of Na(+) inside the groove does not have a dramatic influence in the structure or dynamics of the duplex DNA. The ability of current MD simulations to obtain equilibrated pictures of the counterion atmosphere around DNA is critically discussed.

Structural selectivity of aromatic diamidines

Competition dialysis was used to study the interactions of 13 substituted aromatic diamidine compounds with 13 nucleic acid structures and sequences. The results show a striking selectivity of these compounds for the triplex structure poly dA:(poly dT)2, a novel aspect of their interaction with nucleic acids not previously described. The triplex selectivity of selected compounds was confirmed by thermal denaturation studies. Triplex selectivity was found to be modulated by the location of amidine substituents on the core phenyl-furan-phenyl ring scaffold. Molecular models were constructed to rationalize the triplex selectivity of DB359, the most selective compound in the series. Its triplex selectivity was found to arise from optimal ring stacking on base triplets, along with proper positioning of its amidine substituents to occupy the minor and the major-minor grooves of the triplex. New insights into the molecular recognition of nucleic acid structures emerged from these studies, adding to the list of available design principles for selectively targeting DNA and RNA.

Energetic diversity of DNA minor-groove recognition by small molecules displayed through some model ligand-DNA systems

Energetics of interactions occurring in the model ligand-DNA systems constituted from distamycin A (DST), netropsin (NET) and the oligomeric duplexes d(GCAAGTTGCGATATACG)d(CGTATATCGCAACTTGC)=D#1 and d(GCAAGTTGCGAAAAACG)d(CGTTTTTCGCAACTTGC)=D#2 was studied by spectropolarimetry, UV-absorption spectroscopy and isothermal titration calorimetry. Model analysis of the measured signals was applied to describe individual and competitive binding in terms of populations of various species in the solution. Our results reveal several unprecedented ligand-DNA binding features. DST binds to the neighboring 5'-AAGTT-3' and 5'-ATATA-3' sites of D#1 statistically in a 2:1 binding mode. By contrast, its association to D#2 appears to be a 2:1 binding event only at the DST/D#2 molar ratios between 0 and 2 while its further binding to D#2 may be considered as a step-by-step binding to the unoccupied 5'-AAAAA-3' sites resulting first in DST3D#2 and finally in DST4D#2 complex formation. Competition between DST and NET binding shows that for the most part DST displaces NET from its complexes with D#1 and D#2. In contrast to the obligatory 1:1 binding of DST to the ligand-free 5'-AAAAA-3' sites observed at DST/5'-AAAAA-3' <1 the displacement of NET bound to the 5'-AAAAA-3' sites by added DST occurs even at the smallest additions of DST in a 2:1 manner. NET can also displace DST molecules but only those bound monomerically to the 5'-AAAAA-3' sites of DST3D#2. Actually, only half of these molecules can be displaced due to the simultaneous rebinding of the displaced DST to the unreacted 5'-AAAAA-3' sites in DST3D#2. Binding of DST and NET to D#1 and D#2 is an enthalpy driven process accompanied by large unfavorable (DST), small (NET) or large favorable (NET binding to 5'-AAAAA-3') entropy contributions and negative deltaCP degrees that are reasonably close to deltaCP degrees predicted from the calculated changes in solvent-accessible surface areas that accompany complex formation. Although various modes of DST and NET binding within D#1 and D#2 are characterized by significant energetic differences they seem to be governed by the same driving forces; the hydrophobic transfer of ligand from the solution into the duplex binding site and the accompanying specific non-covalent ligand-DNA and ligand-ligand interactions occurring within the DNA minor groove.

Ion motions in molecular dynamics simulations on DNA

Counterions play a significant role in DNA structure and function, and molecular dynamics (MD) simulations offer the prospect of detailed description of the dynamical structure of ions at the molecular level. However, the motions of mobile counterions are notably slow to converge in MD on DNA. Obtaining accurate and reliable MD simulations requires knowing just how much sampling is required for convergence of each of the properties of interest. To address this issue, MD on a d(CGCGAATTCGCG) duplex in a dilute aqueous solution of water and 22 Na+ counterions was performed until convergence was achieved. The calculated first shell ion occupancies and DNA-Na+ radial distribution functions were computed as a function of time to assess convergence, and compared with relaxation times of the DNA internal parameters shift, slide, rise, tilt, roll, and twist. The sequence dependence of fractional occupancies of ions in the major and minor grooves of the DNA is examined, and the possibility of correlation between ion proximity and DNA minor groove widths is investigated.