Localization of Divalent Metal Ions in the Minor Groove of DNA A-Tracts

Effective Interactions between Double-Stranded DNA Molecules in Aqueous Electrolyte Solutions: Effects of Molecular Architecture and Counterion Valency

A computational investigation of the effects of molecular topology, namely, linear and circular, as well as counterion valency, on the ensuing pairwise effective interactions between DNA molecules in an unlinked state is presented. Umbrella sampling simulations have been performed through the introduction of bias potential along a reaction coordinate defined as the distance between the centers-of-mass of pairs of DNA molecules, and effective pair interaction potentials have been computed by employing the weighted histogram analysis method. An interesting comparison can be drawn between the different DNA topologies studied here, especially with regard to the contrasting effects of divalent counterions on the effective pair potentials: while DNA-DNA repulsion in short center-of-mass distances decreases significantly in the presence of divalent counterion-ions (as compared to monovalent ions) for linear DNA, the opposite effect occurs for the DNA minicircles. This can be attributed to the fact that linear DNA fragments can easily adopt relative orientations that minimize electrostatic and steric repulsions by rotating relative to one another and by exhibiting more pronounced bending due to the presence of free ends.

Insights into DNA solvation found in protein-DNA structures

The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.

Diffusion NMR-based comparison of electrostatic influences of DNA on various monovalent cations

Counterions are important constituents for the structure and function of nucleic acids. Using 7Li and 133Cs nuclear magnetic resonance (NMR) spectroscopy, we investigated how ionic radii affect the behavior of counterions around DNA through diffusion measurements of Li+ and Cs+ ions around a 15-bp DNA duplex. Together with our previous data on 23Na+ and 15NH4+ ions around the same DNA under the same conditions, we were able to compare the dynamics of four different monovalent ions around DNA. From the apparent diffusion coefficients at varied concentrations of DNA, we determined the diffusion coefficients of these cations inside and outside the ion atmosphere around DNA (Db and Df, respectively). We also analyzed ionic competition with K+ ions for the ion atmosphere and assessed the relative affinities of these cations for DNA. Interestingly, all cations (i.e., Li+, Na+, NH4+, and Cs+) analyzed by diffusion NMR spectroscopy exhibited nearly identical Db/Df ratios despite the differences in their ionic radii, relative affinities, and diffusion coefficients. These results, along with the theoretical relationship between diffusion and entropy, suggest that the entropy change due to the release of counterions from the ion atmosphere around DNA is also similar regardless of the monovalent ion types. These findings and the experimental diffusion data on the monovalent ions are useful for examination of computational models for electrostatic interactions or ion solvation.

Dynamics of Cations around DNA and Protein as Revealed by 23Na Diffusion NMR Spectroscopy

Counterions are vital for the structure and function of biomolecules. However, the behavior of counterions remains elusive due to the difficulty in characterizing mobile ions. Here, we demonstrate that the dynamics of cations around biological macromolecules can be revealed by ²³Na diffusion nuclear magnetic resonance (NMR) spectroscopy. NMR probe hardware capable of generating strong magnetic field gradients enables ²³Na NMR-based diffusion measurements for Na⁺ ions in solutions of biological macromolecules and their complexes. The dynamic properties of Na⁺ ions interacting with the macromolecules can be investigated using apparent ²³Na diffusion coefficients measured under various conditions. Our diffusion data clearly show that Na⁺ ions retain high mobility within the ion atmosphere around DNA. The ²³Na diffusion NMR method also permits direct observation of the release of Na⁺ ions from nucleic acids upon protein-nucleic acid association. The entropy change due to the ion release can be estimated from the diffusion data.

Insight into the thermal stability of DNA in hydrated ionic liquids from multi-wavelength UV resonance Raman experiments

The utility of ionic liquids (ILs) as alternative solvents for stabilizing and preserving the native structure of DNA over the long term may be envisaged for biotechnological and biomedical applications in the near future. The delicate balance between the stabilizing and destabilizing effects of IL-mediated interactions with the structure of DNA is complex and is still not well understood. This work reports a fundamental study dealing with the effect exerted by cations and anions in imidazolium-based ILs on the thermal structural stability of large nucleic acid molecules. Multi-wavelength UV resonance Raman spectroscopy is used for selectively detecting heat-induced structural transitions of DNA localized on specific base tracts. Our study reveals the establishment of preferential interactions between the imidazolium cations of ILs and the guanine bases in the DNA groove that lead to more effective stacking between the guanine bases even at high temperatures. Interestingly, we observe that this trend for ILs sharing the same chloride anion is further enhanced as the alkyl chain on the imidazolium cation gets shorter. The results from the present investigation lead to a more comprehensive view of the IL-mediated interactions with A-T and G-C base pairs during thermal unfolding.

Experimental approaches for investigating ion atmospheres around nucleic acids and proteins

Ionic interactions are crucial to biological functions of DNA, RNA, and proteins. Experimental research on how ions behave around biological macromolecules has lagged behind corresponding theoretical and computational research. In the 21st century, quantitative experimental approaches for investigating ionic interactions of biomolecules have become available and greatly facilitated examinations of theoretical electrostatic models. These approaches utilize anomalous small-angle X-ray scattering, atomic emission spectroscopy, mass spectrometry, or nuclear magnetic resonance (NMR) spectroscopy. We provide an overview on the experimental methodologies that can quantify and characterize ions within the ion atmospheres around nucleic acids, proteins, and their complexes.

Dynamics of Ionic Interactions at Protein-Nucleic Acid Interfaces

Molecular association of proteins with nucleic acids is required for many biological processes essential to life. Electrostatic interactions via ion pairs (salt bridges) of nucleic acid phosphates and protein side chains are crucial for proteins to bind to DNA or RNA. Counterions around the macromolecules are also key constituents for the thermodynamics of protein–nucleic acid association. Until recently, there had been only a limited amount of experiment-based information about how ions and ionic moieties behave in biological macromolecular processes. In the past decade, there has been significant progress in quantitative experimental research on ionic interactions with nucleic acids and their complexes with proteins. The highly negatively charged surfaces of DNA and RNA electrostatically attract and condense cations, creating a zone called the ion atmosphere. Recent experimental studies were able to examine and validate theoretical models on ions and their mobility and interactions with macromolecules. The ionic interactions are highly dynamic. The counterions rapidly diffuse within the ion atmosphere. Some of the ions are released from the ion atmosphere when proteins bind to nucleic acids, balancing the charge via intermolecular ion pairs of positively charged side chains and negatively charged backbone phosphates. Previously, the release of counterions had been implicated indirectly by the salt-concentration dependence of the association constant.

Recently, direct detection of counterion release by NMR spectroscopy has become possible and enabled more accurate and quantitative analysis of the counterion release and its entropic impact on the thermodynamics of protein–nucleic acid association. Recent studies also revealed the dynamic nature of ion pairs of protein side chains and nucleic acid phosphates. These ion pairs undergo transitions between two major states. In one of the major states, the cation and the anion are in direct contact and form hydrogen bonds. In the other major state, the cation and the anion are separated by water. Transitions between these states rapidly occur on a picosecond to nanosecond time scale. When proteins interact with nucleic acids, interfacial arginine (Arg) and lysine (Lys) side chains exhibit considerably different behaviors. Arg side chains show a higher propensity to form rigid contacts with nucleotide bases, whereas Lys side chains tend to be more mobile at the molecular interfaces. The dynamic ionic interactions may facilitate adaptive molecular recognition and play both thermodynamic and kinetic roles in protein–nucleic acid interactions.

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.

A review of heavy metal cation binding to deoxyribonucleic acids for the creation of chemical sensors

Various human activities lead to the pollution of ground, drinking, and wastewater with toxic metals. It is well known that metal ions preferentially bind to DNA phosphate backbones or DNA nucleobases, or both. Foreman et al. (Environ Toxicol Chem 30(8):1810-1818, 2011) reported the use of a DNA-dye based assay suitable for use as a toxicity test for potable environmental water. They compared the results of this test with the responses of live-organism bioassays. The DNA-based demonstrated that the loss of SYBR Green I fluorescence dye bound to calf thymus DNA was proportional to the toxicity of the water sample. However, this report raised questions about the mechanism that formed the basis of this quasi-quantitatively test. In this review, we identify the unique and preferred DNA-binding sites of individual metals. We show how highly sensitive and selective DNA-based sensors can be designed that contain multiple binding sites for 21 heavy metal cations that bind to DNA and change its structure, consistent with the release of the DNA-bound dye.

NMR Structural Studies on a DNA Four-Way Junction: Stacking Preference and Localization of the Metal-ion Binding Site

Abstract The stacking preference of a DNA Four-Way junction (4H), with a novel junction sequence, has been determined in the presence of magnesium ions as well as in the presence of cobalt(III)hexammine ions by means of NMR spectroscopy. In both cases this 4H has a strong preference (>80%) to fold in an A/D-stacked conformer. NOESY spectra showed intermolecular NOE contacts between 4H protons and the ammine protons of the cobalt(III)hexammine complex. These contacts define the metal-ion binding site, located in the vicinity of the junction. The position is similar to the observed site in a recent X-ray structure of a RNA/DNA hybrid 4H and consistent with the position deduced from an uranyl ion photoprobing study on 4Hs with different sequences.

Ion distributions around left- and right-handed DNA and RNA duplexes: a comparative study

The ion atmosphere around nucleic acids is an integral part of their solvated structure. However, detailed aspects of the ionic distribution are difficult to probe experimentally, and comparative studies for different structures of the same sequence are almost non-existent. Here, we have used large-scale molecular dynamics simulations to perform a comparative study of the ion distribution around (5'-CGCGCGCGCGCG-3')2 dodecamers in solution in B-DNA, A-RNA, Z-DNA and Z-RNA forms. The CG sequence is very sensitive to ionic strength and it allows the comparison with the rare but important left-handed forms. The ions investigated include Na(+), K(+) and Mg(2 +), with various concentrations of their chloride salts. Our results quantitatively describe the characteristics of the ionic distributions for different structures at varying ionic strengths, tracing these differences to nucleic acid structure and ion type. Several binding pockets with rather long ion residence times are described, both for the monovalent ions and for the hexahydrated Mg[(H2O)6](2+) ion. The conformations of these binding pockets include direct binding through desolvated ion bridges in the GpC steps in B-DNA and A-RNA; direct binding to backbone oxygens; binding of Mg[(H2O)6](2+) to distant phosphates, resulting in acute bending of A-RNA; tight 'ion traps' in Z-RNA between C-O2 and the C-O2' atoms in GpC steps; and others.

DNA conformational flexibility study using phosphate backbone neutralization model

Due to the critical role of DNA in the processes of the cell cycle, the structural and physicochemical properties of DNA have long been of concern. In the present work, the effect of interplay between the DNA duplex and metal ions in solution on the DNA structure and conformational flexibility is studied by comparing the structure and dynamic conformational behavior of a duplex in a normal form and its “null isomer” using molecular dynamics methods. It was found that the phosphate neutralization changes the cation atmosphere around the DNA duplex greatly, increases the major groove width, decreases the minor groove width, and reduces the global bending direction preference. We also noted that the probability of BI phosphate linkages increases significantly because of the charge reduction in the backbone phosphate groups. More importantly, we found that the electrostatic effect on the DNA conformational flexibility is dependent on the sequence; that is, the phosphate backbone neutralization induces the global dynamic bending to be less flexible for GC-rich sequences but more flexible for AT-rich sequences.

Interactions of Mn2+ with a non-self-complementary Z-type DNA duplex

Crystal structures of the hexanucleotide d(CACGCG)·d(CGCGTG) were determined in two crystal lattices when different concentrations of the counterion Mn2+ were used in crystallization. The availability of Mn2+ during the crystallization process appears to play an important role in inducing different crystal packings that lead to crystals belonging to the two space groups P2(1) and P6(5). Analysis of the molecular interactions of Mn2+ with the Z-form duplexes shows direct coordination to the purine residues G and A.

Structure of the tetradecanucleotide d(CCCCGGTACCGGGG)2 as an A-DNA duplex

The crystal structure of the tetradecanucleotide sequence d(CCCCGGTACCGGGG)(2) has been determined at 2.5 Å resolution in the tetragonal space group P4(1). This sequence was designed with the expectation of a four-way junction. However, the sequence crystallized as an A-DNA duplex and represents more than one full turn of the A-helix. The crystallographic asymmetric unit consists of one tetradecanucleotide duplex. The structural parameters of the A-type DNA duplex structure and the crystal-packing arrangement are described. One Mn(2+) ion was identified with direct coordination to the N7 position of G(13) and a water molecule at the major-groove side of the C(2)·G(13) base pair.

Infrared multiple photon dissociation action spectroscopy and theoretical studies of diethyl phosphate complexes: effects of protonation and sodium cationization on structure

The gas-phase structures of deprotonated, protonated, and sodium-cationized complexes of diethyl phosphate (DEP) including [DEP - H](-), [DEP + H](+), [DEP + Na](+), and [DEP - H + 2Na](+) are examined via infrared multiple photon dissociation (IRMPD) action spectroscopy using tunable IR radiation generated by a free electron laser, a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS) with an electrospray ionization (ESI) source, and theoretical electronic structure calculations. Measured IRMPD spectra are compared to linear IR spectra calculated at the B3LYP/6-31G(d,p) level of theory to identify the structures accessed in the experimental studies. For comparison, theoretical studies of neutral complexes are also performed. These experiments and calculations suggest that specific geometric changes occur upon the binding of protons and/or sodium cations, including changes correlating to nucleic acid backbone geometry, specifically P-O bond lengths and ∠OPO bond angles. Information from these observations may be used to gain insight into the structures of more complex systems, such as nucleotides and solvated nucleic acids.

Unexpected dissociation energetics of the Na(+) counterion from GC motifs in DNA hole-migration

We present here a theoretical investigation of the electronic and energetic properties of Na(+)GC, a DNA motif bound to a sodium ion (Na(+)) at the N(7) and O(6) sites of guanine (G), and its hole-trapped derivative [Na(+)GC](+) using density functional theory calculations. Normally, Na(+)GC has positive dissociation energies along various dissociation channels. However, hole-trapping of the Na(+)GC motif can lead to an unusual energetic phenomenon. Hole-trapping can reduce not only the dissociation barrier by destabilizing the Na(+)GC motif to a metastable state, but also the dissociation energy of the Na(+)N(7)/O(6) bond with an unexpected change from a positive to a negative value (61.51 versus-16.18 kcal mol(-1)). This unexpected negative dissociation energy phenomenon implies that this motif can store energy (∼16 kcal mol(-1)) in the Na(+)N(7)/O(6) bond zone due to hole-trapping. The topological properties of electron densities and the Laplacian values at the bond critical points indicate that this energetic phenomenon mainly originates from additional electrostatic repulsions between two moieties linked via a high-energy bond (Na(+)N(7)/O(6)). Proton transfer from G induced by hole-trapping can expand the negative dissociation energy zone to both Na(+)N(7)/O(6) and Watson-Crick (WC) H-bond zones. Similar phenomena can be observed for the Na(+) binding at the minor groove. Solvation of the hole-trapped Na(+)GC motif can change the negative dissociation energies by varying degrees, depending on the solvent-binding sites and the polarity of the solvents.

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.

Intensities of DNA ion-phosphate modes in the low-frequency Raman spectra

The Raman intensities of counterion vibrations with respect to the phosphate groups of the double-helix backbone (ion-phosphate modes) in the low-frequency spectra (< 200 cm(-1)) of B -DNA with different alkali metal counterions have been calculated using the model for DNA conformational vibrations and the valence-optic approach. The results have showed that the spectra of DNA with heavy counterions (Rb(+) and Cs(+)) differ from the spectra of DNA with light counterions (Na(+) and K(+)). The calculated spectra of DNA with heavy counterions are characterized by intensive modes of ion-phosphate vibrations that form one united band near 115 cm(-1). Ion-phosphate modes in the spectra of DNA with light counterions are characterized by higher frequencies (near 180 cm(-1)) and much lower intensity. Our calculations explain why the modes of ion-phosphate vibrations are observed in Cs-DNA spectra rather than in Na-DNA. The determined sensitivity of the intensities of DNA low-frequency spectra to the counterion type proves the existence of the ion-phosphate modes.

The unique structure of A-tracts and intrinsic DNA bending

Short runs of adenines are a ubiquitous DNA element in regulatory regions of many organisms. When runs of 4–6 adenine base pairs (‘A-tracts’) are repeated with the helical periodicity, they give rise to global curvature of the DNA double helix, which can be macroscopically characterized by anomalously slow migration on polyacrylamide gels. The molecular structure of these DNA tracts is unusual and distinct from that of canonical B-DNA. We review here our current knowledge about the molecular details of A-tract structure and its interaction with sequences flanking them of either side and with the environment. Various molecular models were proposed to describe A-tract structure and how it causes global deflection of the DNA helical axis. We review old and recent findings that enable us to amalgamate the various findings to one model that conforms to the experimental data. Sequences containing phased repeats of A-tracts have from the very beginning been synonymous with global intrinsic DNA bending. In this review, we show that very often it is the unique structure of A-tracts that is at the basis of their widespread occurrence in regulatory regions of many organisms. Thus, the biological importance of A-tracts may often be residing in their distinct structure rather than in the global curvature that they induce on sequences containing them.

Remarkable metal counterion effect on the internucleotide J-couplings and chemical shifts of the N-H...N hydrogen bonds in the W-C base pairs

The effects of metal ion binding on the (2h) J(NN)-coupling and delta( (1)H)/Deltadelta( (15)N) chemical shifts of N-H...N H-bond units in internucleotide base pairs were explored by a combination of density functional theory calculations and molecular dynamics (MD) simulations. Results indicate that the NMR parameters vary considerably upon cation binding to the natural GC or AT base pairs, and thus can be used to identify the status of the base pairs, if cation-perturbed. The basic trend is that cation perturbation causes (2h) J(NN) to increase, Deltadelta( (15)N) to decrease, and delta( (1)H) to shift upfield for GC, and in the opposite directions for AT. The magnitudes of variation are closely related to the Lewis acidity of the metal ions. For both base pair series (M(z+)GC and M(z+)AT), these NMR parameters are linearly correlated among themselves. Their values depend strongly on the energy gaps (Delta(ELP-->sigma*)) and the second-order interaction energies ( E(2)) between the donor N lone pair (LP(N)) and the acceptor sigma* N-H localized NBO orbitals. In addition, the (2h) J NN changes are also sensitive to the amount of sigma charge transfer from LP(N) to sigma*(N-H) NBOs or from the purine to the pyrimidine moieties. The different trends are a consequence of the different H-bond patterns combined with the polarization effect of the metal ions in the cationized M(z+)AT series, M(z+) <-- A --> T, and the cationized GC series, M(z+) <-- G <-- C. The predicted cation-induced systematic trends of (2h) J(NN) and delta( (15)N, (1)H) in N-H...N H-bond units may provide a new approach to the determination of H-bond structure and strength in Watson-Crick base pairs, and provide an alternative probe of the heterogeneity of DNA sequences.

Molecular modeling of nucleic acid structure: setup and analysis

The last in a set of units by these authors, this unit addresses some important remaining questions about molecular modeling of nucleic acids. It describes how to choose an appropriate molecular mechanics force field; how to set up and equilibrate the system for accurate simulation of a nucleic acid in an explicit solvent by molecular dynamics or Monte Carlo simulation; and how to analyze molecular dynamics trajectories.

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.

The stability and characteristics of a DNA Holliday junction switch

A Holliday junction (HJ) consists of four DNA double helices, with a branch point discontinuity at the intersection of the component strands. At low ionic strength, the HJ adopts an open conformation, with four widely spaced arms, primarily due to strong electrostatic repulsion between the phosphate groups on the backbones. At high ionic strength, screening of this repulsion induces a switch to a more compact (closed) junction conformation. Fluorescent labelling with dyes placed on the HJ arms allows this conformational switch to be detected optically using fluorescence resonance energy transfer (FRET), producing a sensitive fluorescent output of the switch state. This paper presents a systematic and quantitative survey of the switch characteristics of such a labelled HJ. A short HJ (arm length 8 bp) is shown to be prone to dissociation at low switching ion concentration, whereas an HJ of arm length 12 bp is shown to be stable over all switching ion concentrations studied. The switching characteristics of this HJ have been systematically and quantitatively studied for a variety of switching ions, by measuring the required ion concentration, the sharpness of the switching transition and the fluorescent output intensity of the open and closed states. This stable HJ is shown to have favourable switch characteristics for a number of inorganic switching ions, making it a promising candidate for use in nanoscale biomolecular switch devices.

Cheminformatics methods for novel nanopore analysis of HIV DNA termini

BACKGROUND

Channel current feature extraction methods, using Hidden Markov Models (HMMs) have been designed for tracking individual-molecule conformational changes. This information is derived from observation of changes in ionic channel current blockade "signal" upon that molecule's interaction with (and occlusion of) a single nanometer-scale channel in a "nanopore detector". In effect, a nanopore detector transduces single molecule events into channel current blockades. HMM analysis tools described are used to help systematically explore DNA dinucleotide flexibility, with particular focus on HIV's highly conserved (and highly flexible/reactive) viral DNA termini. One of the most critical stages in HIV's attack is the binding between viral DNA and the retroviral integrase, which is influenced by the dynamic-coupling induced high flexibility of a CA/TG dinucleotide positioned precisely two base-pairs from the blunt terminus of the duplex viral DNA. This suggests the study of a family of such CA/TG dinucleotide molecules via nanopore measurement and cheminformatics analysis.

RESULTS

HMMs are used for level identification on the current blockades, HMM/EM with boosted variance emissions are used for level projection pre-processing, and time-domain FSAs are used to parse the level-projected waveform for kinetic information. The observed state kinetics of the DNA hairpins containing the CA/TG dinucleotide provides clear evidence for HIV's selection of a peculiarly flexible/interactive DNA terminus.

Specific Interactions of Divalent Metal Ions with a DNA Duplex Containing the d(CA)n/(GT)nTandem Repeat

Divalent metal ions are essential for maintaining functional states of the DNA molecule. Their participation in DNA structure is modulated by the base sequence and varies depending on the nature of the ion. The present investigation addresses the interaction of Ca2+ ions with a tandem repeat of two CA dinucleotides, (CA)2/(TG)2. The binding of Ca2+ to the repeat is monitored by nuclear magnetic resonance (NMR) spectroscopy using chemical shift mapping. Parallel experiments monitor binding of Mg2+ ions to the repeat as well as binding of each ion to a DNA duplex in which the (CA)2/(TG)2 repeat is eliminated. The results reveal that the direction and the magnitude of chemical shift changes induced by Ca2+ ions in the NMR spectra of the repeat are different from those induced by Mg2+ ions. The differences between the two cations are significantly diminished by the elimination of the (CA)2/(TG)2 repeat. These findings suggest a specific interaction of Ca2+ ions with the (CA)2/(TG)2 motif. The specificity of the interaction resides in the two A-T base pairs of the repeat, and it involves the major groove of the first A-T base pair and both grooves of the second A-T base pair.

Changes in drug 13C NMR chemical shifts as a tool for monitoring interactions with DNA

The antibiotic drug, netropsin, was complexed with the DNA oligonucleotide duplex [d(GGTATACC)]2 to monitor drug 13C NMR chemical shifts changes. The binding mode of netropsin to the minor groove of DNA is well-known, and served as a good model for evaluating the relative sensitivity of 13C chemical shifts to hydrogen bonding. Large downfield shifts were observed for four resonances of carbons that neighbor sites which are known to form hydrogen bond interactions with the DNA minor groove. Many of the remaining resonances of netropsin exhibit shielding or relatively smaller deshielding changes. Based on the model system presented here, large deshielding NMR shift changes of a ligand upon macromolecule binding can likely be attributed to hydrogen bond formation at nearby sites.

Anion binding to nucleic acids

Nucleic acids are generally considered as efficient cation binders. Therefore, the likelihood that negatively charged ions might intrude their first hydration shell is rarely considered. Here, we show on the basis of (i) a survey of the Nucleic Acid Database, (ii) several structures extracted from the Cambridge Structural Database, and (iii) molecular dynamics simulations, that the nucleotide electropositive edges involving mainly amino, imino, and hydroxyl groups can cast specific anion binding sites. These binding sites constitute also good locations for the binding of the negatively charged groups of the Asp and Glu residues or the nucleic acid phosphate groups. Furthermore, it is observed in several instances that anions, like water molecules and cations, do mediate protein/nucleic acid interactions. Thus, anions as well as negatively charged groups are directly involved in specific recognition and folding phenomena involving polyanionic nucleic acids.

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.

Information content and complexity in the high-order organization of DNA

Nucleic acids are characterized by a vast structural variability. Secondary structural conformations include the main polymorphs A, B, and Z, cruciforms, intrinsic curvature, and multistranded motifs. DNA secondary motifs are stabilized and regulated by the primary base sequence, contextual effects, environmental factors, as well as by high-order DNA packaging modes. The high-order modes are, in turn, affected by secondary structures and by the environment. This review is concerned with the flow of structural information among the hierarchical structural levels of DNA molecules, the intricate interplay between the various factors that affect these levels, and the regulation and physiological significance of DNA high-order structures.

The role of water in protein-DNA recognition

Is it by design or by default that water molecules are observed at the interfaces of some protein-DNA complexes? Both experimental and theoretical studies on the thermodynamics of protein-DNA binding overwhelmingly support the extended hydrophobic view that water release from interfaces favors binding. Structural and energy analyses indicate that the waters that remain at the interfaces of protein-DNA complexes ensure liquid-state packing densities, screen the electrostatic repulsions between like charges (which seems to be by design), and in a few cases act as linkers between complementary charges on the biomolecules (which may well be by default). This review presents a survey of the current literature on water in protein-DNA complexes and a critique of various interpretations of the data in the context of the role of water in protein-DNA binding and principles of protein-DNA recognition in general.

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.

Different Modes of Interaction Between Hydrated Magnesium Ion and DNA Functional Groups: Database Analysis and ab initio Studies

Role of Magnesium ion is well substantiated in DNA structure and function though the appropriate nature of DNA magnesium interaction is still not fully established. We have analyzed available DNA crystal structures in presence of magnesium ion, which show the experimental evidences for various interaction modes between DNA molecule and magnesium ion. Two preferred modes are found: direct coordinating interaction between magnesium ion and electronegative DNA atoms, and the secondary mode of interaction via formation of hydrogen bonds. This qualitative data is further supported by ab initio quantum chemical calculations using restricted Hartree-Fock and Density Functional Theory. We have analyzed the energies and partial charges of different DNA fragments and hydrated magnesium ions, following restrained and unrestrained geometry optimizations along the reaction coordinate. The restrained optimizations for the systems generally show two energy minima separated by an energy barrier, the height ranges from about 5 to 15 kcal/mol, which is in agreement with experimental observations. All these analyses suggest that both modes of interactions occur almost with equal probability, although water mediated secondary mode of interaction is preferred in most cases, which was so far neglected.

Cations as hydrogen bond donors: a view of electrostatic interactions in DNA

Cations are bound to nucleic acids in a solvated state. High-resolution X-ray diffraction studies of oligonucleotides provide a detailed view of Mg2+, and occasionally other ions bound to DNA. In a survey of several such structures, certain general observations emerge. First, cations bind preferentially to the guanine base in the major groove or to phosphate group oxygen atoms. Second, cations interact with DNA most frequently via water molecules in their primary solvation shell, direct ion-DNA contacts being only rarely observed. Thus, the solvated ions should be viewed as hydrogen bond donors in addition to point charges. Finally, ion interaction sites are readily exchangeable: The same site may be occupied by any ion, including spermine, as well as by a water molecule.

Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.

Minor groove site coordination of adenine by platinum group metal ions: effects on basicity, base pairing, and electronic structure

Dithioether- or diamine-tethered adenine derivatives react with Pt(II), Pd(II), and Rh(III) ions to give N3-coordinated complexes of the types [MCl(SSN)](+) (M = Pt or Pd), [RhCl(3)(SSN)], or [RhCl(3)(NNN)] (where SSN = 1-(N9-adenine)-3,6-dithia-heptane or 1-(N9-adenine)-4,7-dithia-octane; NNN = ethylenediamine-N,9-ethyladenine). Single-crystal X-ray analysis confirms the nature of the metal-nucleobase interaction and highlights a conserved intermolecular hydrogen-bonding motif for all the complexes, irrespective of the metal-ion geometry. Coordination significantly reduces the basicity of the adeninyl group, as indicated by a pK(a) value of -0.16 for [PtCl(N3-1-(N9-adenine)-3,6-dithia-heptane)]BF(4), compared to a pK(a) value of 4.2 for 9-ethyladenine. The site of proton binding, N1 or N7, could not be unambiguously assigned from the (1)H NMR data, because of the similar effect on the chemical shifts of the H2 and H8 protons. Density functional calculations at the BP-LACVP level suggest N1 as the site of protonation for this type of complex. This is in contrast to the N7-protonation reported for [Pt(dien)(N3-6,6',9-trimethyladenine)](2+), as reported elsewhere (Meiser et al., Chem.-Eur. J. 1997, 3, 388). However, further electronic structure calculations in the gas phase reveal that the preferred site for protonation for N3-bound complexes is conformationally dependent. N3 coordination was also found to reduce the extent of base pairing between adenine and thymine in dimethylsulfoxide for the self-complementary complex [PtCl(L3)](+) (L3 = 1-(N9-adenine)-3,6-dithia-9-(N1-thymine)nonane), compared to that for the uncomplexed ligand.

A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).

A comparative study of DNA complexation with Mg(II) and Ca(II) in aqueous solution: major and minor grooves bindings

Although structural differences for the Mg-DNA and Ca-DNA complexes are provided in the solid state, such comparative study in aqueous solution has been less investigated. The aim of this study was to examine the bindings of Mg and Ca cations with calf thymus DNA in aqueous solution at physiological pH, using constant concentration of DNA (1.25 or 12.5 mM) and various concentrations of metal ions (2 microM-650 microM). Capillary electrophoresis, UV-visible, and Fourier transform infrared spectroscopic methods were used to determine the cation-binding modes, the binding constants, and DNA structural variations in aqueous solution. Direct Ca-PO(2) binding was evident by major spectral changes (shifting and splitting) of the backbone PO(2) asymmetric stretching at 1222 cm(-1) with K = 4.80 x 10(5) M(-1), whereas an indirect Mg-phosphate interaction occurred (due to the lack of shifting and splitting of the phosphate band at 1222 cm(-1)) with K = 5.6 x 10(4) M(-1). The metal-base bindings were directly for the Mg with K = 3.20 x 10(5) M(-1) and indirectly for the Ca cation with K = 3.0 x 10(4) M(-1). Both major and minor groove bindings were observed with no alteration of the B-DNA conformation.