Atomic-resolution crystal structures of B-DNA reveal specific influences of divalent metal ions on conformation and packing

Probing hydration of monovalent cations condensed around polymeric nucleic acids

We report the relative molar sound velocity increments, [U], partial molar volumes, V(o), and partial molar adiabatic compressibilities, K(S)(o), of the Li(+), Na(+), K(+), Rb(+), Cs(+), NH(4)(+), and N(CH(3))(4)(+) salts of poly(dAdT)poly(dAdT), poly(dGdC)poly(dGdC), poly(dIdC)poly(dIdC), poly(rA)poly(rU), poly(rG)poly(rC), poly(rI)poly(rC), and poly(rU) at 25 degrees C. When analyzing these data, we take into account the Donnan membrane equilibrium effect. Comparison between the values of [U], V(o), and K(S)(o) exhibited by the nucleic acid salts and respective chlorides (LiCl, NaCl, KCl, RbCl, CsCl, NH(4)Cl, and N(CH(3))(4)Cl) yields information about the state of counterion hydration in the vicinity of each nucleic acid structure studied here. Our analysis reveals that the poly(dGdC)poly(dGdC), poly(dIdC)poly(dIdC), and poly(rI)poly(rC) duplexes and single-stranded poly(rU) do not significantly influence the hydration properties of their condensed counterions. In the vicinity of these polymers, counterions retain their full hydration shells (within +/-15%). By contrast, counterions condensed around the poly(dAdT)poly(dAdT), poly(rA)poly(rU), and poly(rG)poly(rC) duplexes are significantly dehydrated and retain, respectively, only 65(+/-18)%, 34(+/-21)%, and 33(+/-9)% of their original hydration shells. Taken together, the volumetric data reported here provide important new information that ultimately may help us understand the central role that hydration and counterions play in modulating the conformational preferences of nucleic acids and the energetics of DNA recognition events.

RNA Adducts with Chlorophyll and Chlorophyllin: Stability and Structural Features

Porphyrins and their metal derivatives are strong nucleic acids binders. Some of these compounds have been used for radiation sensitization therapy of cancer and are targeted to interact with cellular DNA. Chlorophyll (Chl) binds DNA via guanine N-7 atom (major groove) and the backbone phosphate group (Neault and Tajmir-Riahi. Biophys. J. 76, 2177, 1999), whereas chlorophyllin (Chln) intercalates into A-T and G-C regions (Neault and Tajmir-Riahi. J. Phys. Chem. B. 102, 1610, 1998). This study was designed to examine the interaction of RNA with chlorophyll a and chlorophyllin in aqueous solution at physiological pH with pigment/RNA(phosphate) ratios (r) of 1/80 to 1/2. Fourier transform infrared (FTIR) and UV-visible difference spectroscopic methods were used to characterize the nature of pigment-RNA interaction and to establish correlation between spectral changes and the pigment binding mode, binding constant, RNA secondary structure and structural variations of pigment-RNA complexes in aqueous solution. Spectroscopic results showed that Chl and Chln bind RNA through G-C and A-U bases and the backbone phosphate group with overall binding constants of KChl = 1.95 x 10(5) M(-1) and KChln = 1.61 x 10(5) M(-1). The larger K value obtained for Chl-RNA complexes is attributed to the formation of more stable five or six-coordinate Mg cation in the RNA adducts, while the four-coordination Cu(II) in Chln can be more stable than that of the five or six-coordinated copper ion in the Chln-RNA complexes. Aggregation of pigment-RNA complexes occurs at high metalloporphyrin concentrations. No biopolymer secondary structural changes were observed upon pigment interaction and RNA remains in the A-family structure in these pigment complexes.

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.

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.

Face-to-face and edge-to-face pi-pi interactions in a synthetic DNA hairpin with a stilbenediether linker

Synthetic conjugates possessing bis(2-hydroxyethyl)stilbene-4,4'-diether linkers (Sd2) form the most stable DNA hairpins reported to date. Factors that affect stability are length and flexibility of the linkers and pi-stacking of the stilbene moiety on the adjacent base pair. The crystal structure of the hairpin d(GT(4)G)-Sd2-d(CA(4)C) was determined at 1.5 A resolution. The conformations of the two molecules in the asymmetric unit differ both in the linker and the stem portions. One of them shows a planar stilbene that is stacked on the adjacent G:C base pair. The other displays considerable rotation between the phenyl rings and an unprecedented edge-to-face orientation of stilbene and base pair. The observation of considerable variations in the conformation of the Sd moiety in the crystal structure allows us to exclude restriction of motion as the reason for the absence of Sd photoisomerization in the hairpins. Conformational differences in the stem portion of the two hairpin molecules go along with different Mg(2+) binding modes. Most remarkable among them is the sequence-specific coordination of a metal ion in the narrow A-tract minor groove. The crystal structure provides unequivocal evidence that a fully hydrated Mg(2+) ion can penetrate the narrow A-tract minor groove, causing the groove to further contract. Overall, the structural data provide a better understanding of the origins of hairpin stability and their photochemical behavior in solution.

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.

Influence of copper(II) and magnesium(II) ions on the ciprofloxacin binding to DNA

The influence of magnesium(II) and copper(II) ions on the binding of ciprofloxacin to double stranded calf thymus DNA was studied by fluorescence emission spectroscopy, ultraviolet- and circular dichroism (CD) spectroscopy. The interaction of ciprofloxacin and copper(II) ions was followed by strong fluorescence quenching which was almost unaffected by the presence of DNA. On the other hand, only a slight decrease in fluorescence emission intensity, which was enhanced in the presence of DNA, was observed for ciprofloxacin interaction with magnesium(II) ions. Furthermore, magnesium(II) ions increase the thermal stability of the DNA, while, in the presence of ciprofloxacin, the degree of stabilisation is smaller. In contrast, copper(II) ions destabilise double helical DNA to heat, while ciprofloxacin slightly affects only the second transition of the biphasic melting curve of calf thymus DNA. Magnesium(II) ions at 25 degrees C induce conformational transitions of DNA at concentrations of 1.5 mM and 2.5 M, as monitored by CD. On the other hand copper(II) ions induce only one conformational transition, at a concentration of 12.7 microM. At higher concentrations of copper(II) ions (c>700 microM) DNA starts to precipitate. Significant changes in the CD spectra of DNA were observed after addition of ciprofloxacin to a solution containing DNA and copper(II) ions, but not to DNA and magnesium(II) ions. Based on our spectroscopic results, we propose that copper(II) ions are not directly involved into ciprofloxacin binding to DNA via phosphate groups as it has been suggested for magnesium(II) ions.

Surprising roles of electrostatic interactions in DNA-ligand complexes

The positions of cations in x-ray structures are modulated by sequence, conformation, and ligand interactions. The goal here is to use x-ray diffraction to help resolve structural and thermodynamic roles of specifically localized cations in DNA-anthracycline complexes. We describe a 1.34 A resolution structure of a CGATCG(2)-adriamycin(2) complex obtained from crystals grown in the presence of thallium (I) ions. Tl(+) can substitute for biological monovalent cations, but is readily detected by distinctive x-ray scattering, obviating analysis of subtle differences in coordination geometry and x-ray scattering of water, sodium, potassium, and ammonium. Six localized Tl(+) sites are observable adjacent to each CGATCG(2)-adriamycin(2) complex. Each of these localized monovalent cations are found within the G-tract major groove of the intercalated DNA-drug complex. Adriamycin appears to be designed by nature to interact favorably with the electrostatic landscape of DNA, and to conserve the distribution of localized cationic charge. Localized inorganic cations in the major groove are conserved upon binding of adriamycin. In the minor groove, inorganic cations are substituted by a cationic functional group of adriamycin. This partitioning of cationic charge by adriamycin into the major groove of CG base pairs and the minor groove of AT base pairs may be a general feature of sequence-specific DNA-small molecule interactions and a potentially useful important factor in ligand design.

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.

DNA polymorphism: a comparison of force fields for nucleic acids

The improvements of the force fields and the more accurate treatment of long-range interactions are providing more reliable molecular dynamics simulations of nucleic acids. The abilities of certain nucleic acid force fields to represent the structural and conformational properties of nucleic acids in solution are compared. The force fields are AMBER 4.1, BMS, CHARMM22, and CHARMM27; the comparison of the latter two is the primary focus of this paper. The performance of each force field is evaluated first on its ability to reproduce the B-DNA decamer d(CGATTAATCG)(2) in solution with simulations in which the long-range electrostatics were treated by the particle mesh Ewald method; the crystal structure determined by Quintana et al. (1992) is used as the starting point for all simulations. A detailed analysis of the structural and solvation properties shows how well the different force fields can reproduce sequence-specific features. The results are compared with data from experimental and previous theoretical studies.

QSAR and molecular modelling studies on B-DNA recognition of minor groove binders

Aromatic bisamidines have been proved to be efficient compounds against Leishmania spp. and Pneumocystis carinii. Although the mode of action is still not known, these molecules are supposed to be DNA minor groove binders (MGBs). This paper describes a molecular modelling study for a set of MGBs in order to rank them through their complementarity to the Dickerson Drew Dodecamer (DDD) according to their interaction energies with B-DNA. A comparative molecular field analysis (CoMFA) has shown the importance of relatively bulky positively charged groups attached to the MGB aromatic rings, and small and negatively charged substituents into the middle chain. Models were obtained for DNA denaturation related to H-bonding processes of binding modes. Validation of the model demonstrated the robustness of CoMFA in terms of independent test set of similar MGBs. GRID results allotted bioisosteric substitution of z.sbnd;Oz.sbnd; by z.sbnd;NHz.sbnd; in furan ring of furamidine and related compounds as being capable to enhance the binding to DDD.

Influence of netropsin's charges on the minor groove width of d(CGCGAATTCGCG)2

The exact understanding of the interaction of minor groove binding drugs with DNA is of interest due to their importance as transcription controlling drugs. In this study we performed four molecular dynamics simulations, one of the uncomplexed d(CGCGAATTCGCG)(2) dodecamer and three simulations of the DNA complexed with the minor groove binder netropsin. The charged guanidinium and amidinium ends of the small ligand were in one simulation formally uncharged, in the second one normally charged, and in the third simulation we doubled the charges of the two ends. So we are able to filter out the influence the charges exert on the DNA structure. The positive charges reduce the width of the minor groove showing that charges are able to modify the groove width by charge neutralization of the negative phosphate groups. The quality of the used force field was successfully tested by comparing the results of the uncomplexed dodecamer with already reported NMR and x-ray studies. Thus our simulations should be able to describe the minor groove width of DNA in a correct manner underlying the validity of the results.

DNA-dependent divalent cation binding in the nucleosome core particle

Sequence-specific binding of divalent cations to nucleosomal DNA can potentially influence nucleosome position and mobility, as well as modulate interactions with nuclear factors. We define the bonding and specificity of divalent cation interaction with nucleosomal DNA by characterizing Mn2+ binding in the x-ray structure of the nucleosome core particle at 1.9-A resolution. Manganese ions are found ligated with high occupancy in the major groove to 12 of 22 GG and GC base pair steps. The specific location and mode of metal binding is the consequence of unambiguous conformational differences between dinucleotide sites, owing to their sequence context and orientation in the nucleosome core.

Characterization of divalent cation localization in the minor groove of the A(n)T(n) and T(n)A(n) DNA sequence elements by (1)H NMR spectroscopy and manganese(II)

The localization of Mn(2+) in A-tract DNA has been studied by (1)H NMR spectroscopy using a series of self-complementary dodecamer oligonucleotides that contain the sequence motifs A(n)(n) and T(n)A(n), where n = 2, 3, or 4. Mn(2+) localization in the minor groove is observed for all the sequences that have been studied, with the position and degree of localization being highly sequence-dependent. The site most favored for Mn(2+) localization in the minor groove is near the 5'-most ApA step for both the T(n)A(n) and the A(n)T(n) series. For the T(n)A(n) series, this results in two closely spaced symmetry-related Mn(2+) localization sites near the center of each duplex, while for the A(n)T(n) series, the two symmetry-related sites are separated by as much as one half-helical turn. The degree of Mn(2+) localization in the minor groove of the T(n)A(n) series decreases substantially as the AT sequence element is shortened from T(4)A(4) to T(2)A(2). The A(n)T(n) series also exhibits length-dependent Mn(2+) localization; however, the degree of minor groove occupancy by Mn(2+) is significantly less than that observed for the T(n)A(n) series. For both A(n)T(n) and T(n)A(n) sequences, the 3'-most AH2 resonance is the least broadened of the AH2 resonances. This is consistent with the observation that the minor groove of A-tract DNA narrows in the 5' to 3' direction, apparently becoming too narrow after two base pairs for the entry of a fully hydrated divalent cation. The results that are reported illustrate the delicate interplay that exists between DNA nucleotide sequence, minor groove width, and divalent cation localization. The proposed role of cation localization in helical axis bending by A-tracts is also discussed.

Thymine-methyl/pi interaction implicated in the sequence-dependent deformability of DNA

The crystal structures of deoxy-oligonucleotides were retrieved from the Nucleic Acid Database and analyzed with the use of our program CHPI. The structure of 5'-ApTpApT-3' has been shown to be stabilized by the 5-methyl group in the thymidine moiety that favorably interacts with the adenine pi-ring preceding it. H2' of the deoxyribose in adenine also interacts with the thymine ring next to it. Since a 5'-ApT-3' sequence is accompanied by another 5'-ApT-3' in the complementary strand, the interaction is duplicated, thus forming a 'twin A/T-Me interaction'. Coordinates of oligonucleotides with A-T rich sequences were retrieved and analyzed. In almost every case, the thymidine 5-methyl group favorably interacts with an adenine ring in the same strand. The structure of duplexes incorporating A-tracts was also analyzed. The 5-methyl group in the thymidine moiety has been found to interact favorably with the base pi-ring before it. Since an A-tract is lined with an oligo-T sequence in the complementary strand, a successive N/T-Me stacking may contribute in making the A-tracts robust and straight. The possible involvement of the N/T-Me and the twin A/T-Me motif in the deformability of DNA has been suggested. The role of methyl groups in modified DNA has been discussed on a similar basis.

Metal ions and flexibility in a viral RNA pseudoknot at atomic resolution

Many pathogenic viruses use programmed -1 ribosomal frameshifting to regulate translation of their structural and enzymatic proteins from polycistronic mRNAs. Frameshifting is commonly stimulated by a pseudoknot located downstream from a slippery sequence, the latter positioned at the ribosomal A and P sites. We report here the structures of two crystal forms of the frameshifting RNA pseudoknot from beet western yellow virus at resolutions of 1.25 and 2.85 A. Because of the very high resolution of 1.25 A, ten mono- and divalent metal ions per asymmetric unit could be identified, giving insight into potential roles of metal ions in stabilizing the pseudoknot. A magnesium ion located at the junction of the two pseudoknot stems appears to play a crucial role in stabilizing the structure. Because the two crystal forms exhibit mostly unrelated packing interactions and local crystallographic disorder in the high-resolution form was resolvable, the two structures offer the most detailed view yet of the conformational preference and flexibility of an RNA pseudoknot.

Water mediated Dickerson-Drew-type crystal of DNA dodecamer containing 2'-deoxy-5-formyluridine

To investigate the role of divalent cations in crystal packing, a Dickerson-Drew-type dodecamer with the sequence d(CGCGAATXCGCG), containing 2'-deoxy-5-formyluridine at X, was crystallized under several conditions with Ba(2+) ion instead of Mg(2+) ion. The crystal structure is isomorphous with the original Dickerson-type crystal containing Mg(2+) ion. In the Mg(2+)-free crystals, however, a five-membered ring of water molecules occupies the same position as the magnesium site found in the Mg(2+)-containing crystals, and connects the two duplexes similarly to the hydrated Mg(2+) ion. It has been concluded that the five-membered water molecules can take the place of the hydrated magnesium cation in crystallization. The 5-formyluracil residues form the canonical Watson-Crick pair with the opposite adenine residues.

DNA-cation interactions: quo vadis?

The interactions between double helical DNA and cations, specifically mono- and divalent metal ions, have recently received increased attention. Molecular dynamics simulations, solution NMR, and X-ray crystallography have all shed light on the coordination of ions in the major and minor grooves of DNA. Metal ion interactions may play key roles in the control of DNA conformation and topology, but despite progress in locating the ions and determining their precise binding modes, it remains difficult to figure out just how important ions really are. What have we learned and what remains to be done?

Cation-chromatin binding as shown by ion microscopy is essential for the structural integrity of chromosomes

Mammalian interphase and mitotic cells were analyzed for their cation composition using a three-dimensional high resolution scanning ion microprobe. This instrument maps the distribution of bound and unbound cations by secondary ion mass spectrometry (SIMS). SIMS analysis of cryofractured interphase and mitotic cells revealed a cell cycle dynamics of Ca2+, Mg2+, Na+, and K+. Direct analytical images showed that all four, but no other cations, were detected on mitotic chromosomes. SIMS measurements of the total cation content for diploid chromosomes imply that one Ca2+ binds to every 12.5-20 nucleotides and one Mg2+ to every 20-30 nucleotides. Only Ca2+ was enriched at the chromosomal DNA axis and colocalized with topoisomerase IIalpha (Topo II) and scaffold protein II (ScII). Cells depleted of Ca2+ and Mg2+ showed partially decondensed chromosomes and a loss of Topo II and ScII, but not hCAP-C and histones. The Ca2+-induced inhibition of Topo II catalytic activity and direct binding of Ca2+ to Topo II by a fluorescent filter-binding assay supports a regulatory role of Ca2+ during mitosis in promoting solely the structural function of Topo II. Our study directly implicates Ca2+, Mg2+, Na+, and K+ in higher order chromosome structure through electrostatic neutralization and a functional interaction with nonhistone proteins.

Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts

Introduction of a T-A or pyrimidine-purine step into a straight and rigid A-tract can cause a positive roll deformation that kinks the DNA helix at that step. In CCTTTAAAGG, the central T-A step has an 8.6 degrees bend toward the major groove. We report the structural analysis of CCTTTAAAGG and a comparison with 25 other representative crystal structures from the NDB containing at least four consecutive A or T bases. On average, more local bending occurs at the disruptive T-A step (8.21 degrees ) than at an A-T step (5.71 degrees ). In addition, A-tracts containing an A-T step are more bent than are pure A-tracts, and hence A-A and A-T steps are not equivalent. All T-A steps examined exhibit positive roll, bending towards the major groove, while A-T steps display negative roll and bend slightly towards the minor groove. This illustrates how inherent negative and positive roll are, respectively, at A-T and T-A steps within A-tracts. T-A steps are more deformable, showing larger and more variable deformations of minor groove width, rise, cup, twist, and buckle. Standard deviations of twist, rise, and cup for T-A steps are 6.66 degrees, 0.55 A, and 15.90 degrees, versus 2.28 degrees, 0.21 A, and 2.99 degrees for A-T steps. Packing constraints determine which local values of these helical parameters an individual T-A step will adopt. For instance, with CCTTTAAAGG and three isomorphous structures, CGATTAATCG, CGATATATCG, and CGATCGATCG, crystal packing forces lead to a series of correlated changes: widened minor groove, large slide, low twist, and large rise. The difference in helical parameters between A-T steps lying within A-tracts, versus A-T steps within alternating AT sequences, demonstrates the importance of neighboring steps on the conformation of a given dinucleotide step.

Locating monovalent cations in the grooves of B-DNA

Here we demonstrate that monovalent cations can localize around B-DNA in geometrically regular, sequence-specific sites in oligonucleotide crystals. Positions of monovalent ions were determined from high-resolution X-ray diffraction of DNA crystals grown in the presence of thallium(I) cations (Tl(+)). Tl(+) has previously been shown to be a useful K(+) mimic. Tl(+) positions determined by refinement of model to data are consistent with positions determined using isomorphous F(Tl) - F(K) difference Fouriers and anomalous difference Fouriers. None of the observed Tl(+) sites surrounding CGCGAATTCGCG are fully occupied by Tl(+) ions. The most highly occupied sites, located within the G-tract major groove, have estimated occupancies ranging from 20% to 35%. The occupancies of the minor groove sites are estimated to be around 10%. The Tl(+) positions in general are not in direct proximity to phosphate groups. The A-tract major groove appears devoid of localized cations. The majority of the observed Tl(+) ions interact with a single duplex and so are not engaged in lattice interactions or crystal packing. The locations of the cation sites are dictated by coordination geometry, electronegative potential, avoidance of electropositive amino groups, and cation-pi interactions. It appears that partially dehydrated monovalent cations, hydrated divalent cations, and polyamines compete for a common binding region on the floor of the G-tract major groove.

Magnesium increases the curvature of duplex DNA that contains dA tracts

Distinct structural features of DNA, such as the curvature of dA tracts, are important in the recognition, packaging, and regulation of DNA. Physiologically relevant concentrations of magnesium have been found to enhance the curvature of dA tract DNAs, as monitored by solution-state NMR, indicating that the structure of DNA depends on the cations present in solution. A model is presented which accounts for the sequence-dependent effects of magnesium on DNA curvature as well as for the previously known sequence-independent effect on DNA flexibility.

Re-examination of the intrinsic, dynamic and hydration properties of phosphoramidate DNA

Intrinsic energetic and solvation factors contributing to the unusual structural and biochemical properties of N3'-phosphoramidate DNA analogs have been re-examined using a combination of quantum mechanical and molecular dynamics methods. Evaluation of the impact of the N3'-H substitution was performed via comparison of N3'-phosphoramidate DNA starting from both A- and B-form structures, B-form DNA and A-form RNA. The N3'-H group is shown to be flexible, undergoing reversible inversion transitions associated with motion of the hydrogen atom attached to the N3' atom. The inversion process is correlated with both sugar pucker characteristics as well as other local backbone torsional dynamics, yielding increased dihedral flexibility over DNA. Solvation of N3'-phosphoramidate DNA is shown to be similar to RNA, consistent with thermodynamic data on the two species. A previously unobserved intrinsic conformational perturbation caused by the N5'-phosphoramidate substitution is identified and suggested to be linked to the differences in the properties of N3'- and N5'-phosphoramidate oligonucleotide analogs.

Molecular modeling of nucleic acid structure: electrostatics and solvation

This unit presents an overview of computer simulation techniques as applied to nucleic acid systems, ranging from simple in vacuo molecular modeling techniques to more complete all-atom molecular dynamics treatments that include an explicit representation of the environment. The third in a series of four units, this unit focuses on critical issues in solvation and the treatment of electrostatics.

DNA bending by an adenine--thymine tract and its role in gene regulation

To gain insight into the structural basis of DNA bending by adenine-thymine tracts (A-tracts) and their role in DNA recognition by gene-regulatory proteins, we have determined the crystal structure of the high-affinity DNA target of the cancer-associated human papillomavirus E2 protein. The three independent B-DNA molecules of the crystal structure determined at 2.2-A resolution are examples of A-tract-containing helices where the global direction and magnitude of curvature are in accord with solution data, thereby providing insights, at the base pair level, into the mechanism of DNA bending by such sequence motifs. A comparative analysis of E2-DNA conformations with respect to other structural and biochemical studies demonstrates that (i) the A-tract structure of the core region, which is not contacted by the protein, is critical for the formation of the high-affinity sequence-specific protein-DNA complex, and (ii) differential binding affinity is regulated by the intrinsic structure and deformability encoded in the base sequence of the DNA target.

DNA-cation interactions: The major and minor grooves are flexible ionophores

Several crystallographic, solution-state and theoretical studies carried out this past year provide new support for the sequence-specific nature of monovalent and divalent cation coordination within the DNA major and minor grooves. Correlations observed between groove width and cation coordination indicate that the grooves are flexible and respond to cation binding.

Detection of alkali metal ions in DNA crystals using state-of-the-art X-ray diffraction experiments

The observation of light metal ions in nucleic acids crystals is generally a fortuitous event. Sodium ions in particular are notoriously difficult to detect because their X-ray scattering contributions are virtually identical to those of water and Na(+.)O distances are only slightly shorter than strong hydrogen bonds between well-ordered water molecules. We demonstrate here that replacement of Na(+) by K(+), Rb(+) or Cs(+) and precise measurements of anomalous differences in intensities provide a particularly sensitive method for detecting alkali metal ion-binding sites in nucleic acid crystals. Not only can alkali metal ions be readily located in such structures, but the presence of Rb(+) or Cs(+) also allows structure determination by the single wavelength anomalous diffraction technique. Besides allowing identification of high occupancy binding sites, the combination of high resolution and anomalous diffraction data established here can also pinpoint binding sites that feature only partial occupancy. Conversely, high resolution of the data alone does not necessarily allow differentiation between water and partially ordered metal ions, as demonstrated with the crystal structure of a DNA duplex determined to a resolution of 0.6 A.

DNA structure: what's in charge?

DNA structure is well known to be sensitive to hydration and ionic strength. Recent theoretical predictions and experimental observations have raised the idea of the intrusion of monovalent cations into the minor groove spine of hydration in B-form DNA. To investigate this further, extensions and further analysis of molecular dynamics (MD) simulations on d(CGCCGAATTCGCG), d(ATAGGCAAAAAATAGGCAAAAATGG) and d(G(5)-(GA(4)T(4)C)(2)-C(5)), including counterions and water, have been performed. To examine the effective of minor groove ions on structure, we analyzed the MD snapshots from a 15 ns trajectory on d(CGCGAATTCGCG) as two subsets: those exhibiting a minor groove water spine and those with groove-bound ions. The results indicate that Na(+) at the ApT step of the minor groove of d(CGCCGAATTCGCG) makes only small local changes in the DNA structure, and these changes are well within the thermal fluctuations calculated from the MD. To examine the effect of ions on the differential stability of a B-form helix, further analysis was performed on two longer oligonucleotides, which exhibit A-tract-induced axis bending localized around the CpG step in the major groove. Plots of axis bending and proximity of ions to the bending locus were generated as a function of time and revealed a strong linear correlation, supporting the idea that mobile cations play a key role in local helix deformations of DNA and indicating ion proximity just precedes the bending event. To address the issue of "what's in charge?" of DNA structure more generally, the relative free energy of A and B-form d(CGCGAATTCGCG) structures from MD simulations under various environmental circumstances were estimated using the free energy component method. The results indicate that the dominant effects on conformational stability come from the electrostatic free energy, but not exclusively from groove bound ions per se, but from a balance of competing factors in the electrostatic free energy, including phosphate repulsions internal to the DNA, the electrostatic component of hydration (i.e. solvent polarization), and electrostatic effects of the counterion atmosphere. In summary, free energy calculations indicate that the electrostatic component is dominant, MD shows temporal proximity of mobile counterions to be correlated with A-track-induced bending, and thus the mobile ion component of electrostatics is a significant contributor. However, the MD structure of the dodecamer d(CGCGAATTCGCG) is not highly sensitive to whether there is a sodium ion in the minor groove.

Local conformational variations observed in B-DNA crystals do not improve base stacking: computational analysis of base stacking in a d(CATGGGCCCATG)(2) B<-->A intermediate crystal structure

The crystal structure of d(CATGGGCCCATG)(2) shows unique stacking patterns of a stable B<-->A-DNA intermediate. We evaluated intrinsic base stacking energies in this crystal structure using an ab initio quantum mechanical method. We found that all crystal base pair steps have stacking energies close to their values in the standard and crystal B-DNA geometries. Thus, naturally occurring stacking geometries were essentially isoenergetic while individual base pair steps differed substantially in the balance of intra-strand and inter-strand stacking terms. Also, relative dispersion, electrostatic and polarization contributions to the stability of different base pair steps were very sensitive to base composition and sequence context. A large stacking flexibility is most apparent for the CpA step, while the GpG step is characterized by weak intra-strand stacking. Hydration effects were estimated using the Langevin dipoles solvation model. These calculations showed that an aqueous environment efficiently compensates for electrostatic stacking contributions. Finally, we have carried out explicit solvent molecular dynamics simulation of the d(CATGGGCCCATG)(2) duplex in water. Here the DNA conformation did not retain the initial crystal geometry, but moved from the B<-->A intermediate towards the B-DNA structure. The base stacking energy improved in the course of this simulation. Our findings indicate that intrinsic base stacking interactions are not sufficient to stabilize the local conformational variations in crystals.

A novel form of intercalation involving four DNA duplexes in an acridine-4-carboxamide complex of d(CGTACG)(2)

The structures of the complexes formed between 9-amino-[N:-(2-dimethyl-amino)butyl]acridine-4-carboxamide and d(CG(5Br)UACG)(2) and d(CGTACG)(2) have been solved by X-ray crystallography using MAD phasing methodology and refined to a resolution of 1.6 A. The complexes crystallised in space group C222. An asymmetric unit in the brominated complex comprises two strands of DNA, one disordered drug molecule, two cobalt (II) ions and 19 water molecules (31 in the native complex). Asymmetric units in the native complex also contain a sodium ion. The structures exhibit novel features not previously observed in crystals of DNA/drug complexes. The DNA helices stack in continuous columns with their central 4 bp adopting a B-like motif. However, despite being a palindromic sequence, the terminal GC base pairs engage in quite different interactions. At one end of the duplex there is a CpG dinucleotide overlap modified by ligand intercalation and terminal cytosine exchange between symmetry-related duplexes. A novel intercalation complex is formed involving four DNA duplexes, four ligand molecules and two pairs of base tetrads. The other end of the DNA is frayed with the terminal guanine lying in the minor groove of the next duplex in the column. The structure is stabilised by guanine N7/cobalt (II) coordination. We discuss our findings with respect to the effects of packing forces on DNA crystal structure, and the potential effects of intercalating agents on biochemical processes involving DNA quadruplexes and strand exchanges. NDB accession numbers: DD0032 (brominated) and DD0033 (native).

1 A crystal structures of B-DNA reveal sequence-specific binding and groove-specific bending of DNA by magnesium and calcium

The 1 A resolution X-ray crystal structures of Mg(2+) and Ca(2+) salts of the B-DNA decamers CCAACGTTGG and CCAGCGCTGG reveal sequence-specific binding of Mg(2+) and Ca(2+) to the major and minor grooves of DNA, as well as non-specific binding to backbone phosphate oxygen atoms. Minor groove binding involves H-bond interactions between cross-strand DNA base atoms of adjacent base-pairs and the cations' water ligands. In the major groove the cations' water ligands can interact through H-bonds with O and N atoms from either one base or adjacent bases, and in addition the softer Ca(2+) can form polar covalent bonds bridging adjacent N7 and O6 atoms at GG bases. For reasons outlined earlier, localized monovalent cations are neither expected nor found.Ultra-high atomic resolution gives an unprecedented view of hydration in both grooves of DNA, permits an analysis of individual anisotropic displacement parameters, and reveals up to 22 divalent cations per DNA duplex. Each DNA helix is quite anisotropic, and alternate conformations, with motion in the direction of opening and closing the minor groove, are observed for the sugar-phosphate backbone. Taking into consideration the variability of experimental parameters and crystal packing environments among these four helices, and 24 other Mg(2+) and Ca(2+) bound B-DNA structures, we conclude that sequence-specific and strand-specific binding of Mg(2+) and Ca(2+) to the major groove causes DNA bending by base-roll compression towards the major groove, while sequence-specific binding of Mg(2+) and Ca(2+) in the minor groove has a negligible effect on helix curvature. The minor groove opens and closes to accommodate Mg(2+) and Ca(2+) without the necessity for significant bending of the overall helix. The program Shelxdna was written to facilitate refinement and analysis of X-ray crystal structures by Shelxl-97 and to plot and analyze one or more Curves and Freehelix output files.

Solvent organization in an oligonucleotide crystal. The structure of d(GCGAATTCG)2 at atomic resolution

We describe the crystal structure of d(GCGAATTCG) determined by x-ray diffraction at atomic resolution level (0.89 A). The duplex structure is practically identical to that described at 2.05 A resolution (Van Meervelt, L., Vlieghe, D., Dautant, A., Gallois, B., Précigoux, G., and Kennard, O. (1995) Nature 374, 742-744), however about half of the phosphate groups show multiple conformations. The crystal has three regions with different solvent structure. One of them contains several ordered Mg(+2) ions and can be considered as an ionic crystal. A second region is formed by a network of ordered water molecules with a polygonal organization that binds three duplexes. The third region is formed by channels of solvent in which very few ordered solvent molecules are visible. The less ordered phosphates are found facing this channel. The latter region provides a view of DNA with highly movable charges, both negative phosphates and counterions, without a precise location.

Influence of divalent cations on the conformation of phosphorothioate oligodeoxynucleotides: a circular dichroism study

Phosphorothioate oligodeoxynucleotides (ODNs) have been extensively investigated in vivo and in vitro for antisense control of gene expression. It has been shown that cellular uptake of phosphorothioate ODNs in some in vitro cell systems increases in the presence of divalent cations. In this work, we analyze the conformation of phosphorothioate ODNs and specific changes induced in it by various divalent cations using circular dichroism (CD) spectroscopy. CD data were obtained with several phosphorothioate ODNs in the absence and presence of the divalent cations Mg(2+), Ca(2+), Sr(2+), Ba(2+) and Mn(2+). All CD spectra indicated stable conformations of the ODNs in solution. The spectra were strongly dependent on ODN sequence and composition. Some ODNs such as T(23) and another with 'random' distribution of bases showed CD spectra characteristic of B-form DNA. Other ODNs which had at least three consecutive guanines in their sequences exhibited spectra characteristic of parallel G-tetraplexes. CD spectra of antisense ODNs exhibited specific responses to divalent cations. Changes in the conformation were not simply due to ionic strength effects. Mn(2+) diminished secondary structure in some ODNs. Group II divalent ions stabilized the parallel G-tetraplexes, and Mg(2+) generally had the weakest stabilizing efficiency. Each sequence/ion combination had a specific response so these effects cannot be generalized. These sequence-dependent, divalent ion-sensitive, and structurally unique solution conformations may be related to ion-mediated ODN uptake.

Interaction of alkaline-earth metal ions with calf thymus DNA. Volume and compressibility effects in diluted aqueous solutions

The binding of Mg2+, Ca2+, Sr2+ and Ba2+ ions to calf thymus DNA in solutions has been investigated by ultrasonic and densimetric techniques. The obtained parameters, the apparent molar volume, phiV, and the apparent molar adiabatic compressibility, phiK(S), are very sensitive to hydration of investigated molecules. The interaction between the cations and DNA is accompanied by overlapping their hydration shells and consequently releasing the water molecules from hydration shells to bulk state. The change in the hydration is reflected in the measured parameters, phiV and phiK(S). The magnitude of these hydration changes is determined by the position of the cation relative to DNA atomic groups involved in the binding, and thus can characterize the structure of cation-DNA complexes. The values of the dehydration effects of the binding, deltaphiV and deltaphiK(S), correspond to two direct or higher number of indirect contacts between calf thymus DNA and the cations.

Intermolecular interactions and water structure in a condensed phase B-DNA crystal

By controlled dehydration, the unit cells of dodecamer DNA-drug crystals have been shrunk from 68,000 (normal state) to 60,000 (partially dehydrated intermediate state) to 51,000 A(3) (fully dehydrated state), beyond which no further solvent loss occurs. The total solvent content in the normal crystals is approximately 40% by volume, reducing to approximately 20% in the fully dehydrated phase. The 25% reduction in cell volume induced a dramatic enhancement in the resolution of the X-ray diffraction data (from 2. 6 to beyond 1.5 A). We have determined the structures of the normal, partially dehydrated and fully dehydrated crystals. Details of the ligand binding have been presented in the preceding article. The present paper describes the unique features of the structure of the fully dehydrated phase. This structure was refined with 9,015 unique observed reflections to R = 14.9%, making it one of the most reliable models of B -form DNA available. The crystals exist as infinite polymeric networks, in which neighbouring dodecamer duplexes are crosslinked through phosphate oxygens via direct bonding to magnesium cations. The DNA is packed so tightly that there is essentially only a single layer of solvent between adjacent molecules. The details of the crystal packing, magnesium bridging, DNA hydration and DNA conformation are described and compared with other experimental evidence related to DNA condensation.

Conformational flexibility of B-DNA at 0.74 A resolution: d(CCAGTACTGG)(2)

The affinity and specificity of a ligand for its DNA site is a function of the conformational changes between the isolated and complexed states. Although the structures of a hydroxypyrrole-imidazole-pyrrole polyamide dimer with 5'-CCAGTACTGG-3' and the trp repressor recognizing the sequence 5'-GTACT-3' are known, the baseline conformation of the DNA site would contribute to our understanding of DNA recognition by these ligands. The 0.74 A resolution structure of a B-DNA double helix, 5'-CCAGTACTGG-3', has been determined by X-ray crystallography. Six of the nine phosphates, two of four bound calcium ions and networks of water molecules hydrating the oligonucleotide have alternate conformations. By contrast, nine of the ten bases have a single, unique conformation with hydrogen atoms visible in most cases. The polyamide molecules alter the geometry of the phosphodiester backbone, and the water molecules mediating contacts in the trp repressor/operator complex are conserved in the unliganded DNA. Furthermore, the multiple conformational states, ions and hydration revealed by this ultrahigh resolution structure of a B-form oligonucleotide are potentially general considerations for understanding DNA-binding affinity and specificity by ligands.