Spermine: an "invisible" component in the crystals of B-DNA. A grand canonical Monte Carlo and molecular dynamics simulation study

The association of spermine(4+) (Spm(4+)), Mg(2+) and monovalent (M(+)) ions with DNA in crystal form, have been studied using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) computer simulations. GCMC calculations were used to calculate the distribution of Spm(4+), Mg(2+), and M(+) between the equilibrating solvent and the DNA crystal under conditions mimicking the crystal-growing protocols reported in a number of recent X-ray diffraction studies of DNA oligomers. The GCMC simulations show that the composition of ions neutralizing the negative charge of DNA can vary in a broad range. The GCMC simulations were used to provide appropriate conditions for subsequent 6 ns constant pressure and temperature MD simulations of DNA in a typical crystalline environment consisting of three DNA double helix decamers in a periodic hexagonal cell, containing 1200 water molecules, eight Spm(4+), 32 Na(+) and four Cl(-) ions. Based on the simulation results, it seems possible to give an explanation why spermine molecules are usually not detected in X-ray studies in spite of their high concentration in the preparatory samples used as the crystallizing agent. It appears that this flexible polyamine molecule has several binding modes, interacting in fairly irregular manner with different sites on DNA and showing no regular ordering in the DNA crystals. Ions of Na(+) and Spm(4+) compete with each other and with water molecules in binding to bases in the minor groove and they influence the structure of the DNA hydration shell in different ways.     

A molecular dynamics simulation study of polyamine– and sodium–DNA. Interplay between polyamine binding and DNA structure

Four different molecular dynamics (MD) simulations have been performed for infinitely long ordered DNA molecules with different counterions, namely the two natural polyamines spermidine(3+) (Spd³⁺) and putrescine(2+) (Put²⁺), the synthetic polyamine diaminopropane(2+) (DAP²⁺), and the simple monovalent cation Na⁺. All systems comprised a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charge. The simulation setup mimics the DNA state in oriented DNA fibers, previously studied using NMR and other experimental methods. In this paper the interplay between polyamine binding and local DNA structure is analyzed by investigating how and if the minor groove width of DNA depends on the presence and dynamics of the counterions. The results of the MD simulations reveal principal differences in the polyamine–DNA interactions between the natural [spermine(4+), Spd³⁺, Put²⁺] and the synthetic (DAP²⁺) polyamines.Abbreviations DAP diaminopropane - DDD Drew–Dickerson dodecamer - MD molecular dynamics - Put putrescine - RDF radial distribution function - Spd spermidine - Spm spermine     

Molecular dynamics simulation study of oriented polyamine- and Na-DNA: sequence specific interactions and effects on DNA structure

Molecular dynamics (MD) computer simulations have been carried out on four systems that correspond to an infinite array of parallel ordered B-DNA, mimicking the state in oriented DNA fibers and also being relevant for crystals of B-DNA oligonucleotides. The systems were all comprised of a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charges. The sequence of the double helical DNA decamer was d(5'-ATGCAGTCAG)xd(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine(3+) (Spd(3+)) and putrescine(2+) (Put(2+)), the synthetic polyamine diaminopropane(2+) (DAP(2+)), and the simple monovalent cation Na(+). This work compares the specific structures of the polyamine- and Na-DNA systems and how they are affected by counterion interactions. It also describes sequence-specific hydration and interaction of the cations with DNA. The local DNA structure is dependent on the nature of the counterion. Even the very similar polyamines, Put(2+) and DAP(2+), show clear differences in binding to DNA and in effect on hydration and local structure. Generally, the polyamines disorder the hydration of the DNA around their binding sites whereas Na(+) being bound to DNA attracts and organizes water in its vicinity. Cation binding at the selected sites in the minor and in the major groove is compared for the different polyamines and Na(+). We conclude that the synthetic polyamine (DAP(2+)) binds specifically to several structural and sequence-specific motifs on B-DNA, unlike the natural polyamines, Spd(3+) and Put(2+). This specificity of DAP(2+) compared to the more dynamic behavior of Spd(3+) and Put(2+) may explain why the latter polyamines are naturally occurring in cells.     

Binding of cations of group IA and IIA to bovine serum amine oxidase: effect on the activity

In this paper, we report on the presence of cation binding areas on bovine serum amine oxidase, where metal ions of the groups IA and IIA, such as Na(+), K(+), Cs(+), Mg(2+), and Ca(2+), bind with various affinities. We found a cation-binding area that influences the enzyme activity if occupied, so that the catalytic reaction may be altered by some physiologically relevant cations, such as Ca(2+) and K(+). This binding area appears to be localized inside the enzyme active site, because some of these cations act as competitive inhibitors when highly charged amines, such as spermine and spermidine, are used as substrates. In particular, dissociation constant values (K(d)) of 23 and 27 mM were measured for Cs(+) and Ca(2+), respectively, using, as substrate, spermine, a polyamine of plasma. An additional cation-binding area, where metal ions such as Cs(+) (K(d) congruent with 0.1 mM) and Na(+) (K(d) congruent with 54 mM) bind without affecting the enzyme activity, was found by NMR.     

Binding sites of the polyamines putrescine, cadaverine, spermidine and spermine on A- and B-DNA located by simulated annealing

Molecular dynamics simulations with simulated annealing are performed on polyamine-DNA systems in order to determine the binding sites of putrescine, cadaverine, spermidine and spermine on A- and B-DNA. The simulations either contain no additional counterions or sufficient Na+ ions, together with the charge on the polyamine, to provide 73% neutralisation of the charges on the DNA phosphates. The stabilisation energies of the complexes indicate that all four polyamines should stabilise A-DNA in preference to B-DNA, which is in agreement with experiment in the case of spermine and spermidine, but not in the case of putrescine or cadaverine. The major groove is the preferred binding site on A-DNA of all the polyamines. Putrescine and cadaverine tend to bind to the sugar-phosphate backbone of B-DNA, whereas spermidine and spermine occupy more varied sites, including binding along the backbone and bridging both the major and minor grooves.     

Molecular dynamics simulations of DNA in solutions with different counter-ions

Molecular dynamics simulations of the [d(ATGCAGTCAG]2 fragment of DNA, in water and in the presence of three different counter-ions (Li+, Na+ and Cs+) are reported. Three-dimensional hydration structure and ion distribution have been calculated using spatial distribution functions for a detailed picture of local concentrations of ions and water molecules around DNA. According to the simulations, Cs+ ions bind directly to the bases in the minor groove, Na+ ions bind prevailing to the bases in the minor groove through one water molecule, whereas Li+ ions bind directly to the phosphate oxygens. The different behavior of the counter-ions is explained by specific hydration structures around the DNA and the ions. It is proposed how the observed differences in the ion binding to DNA may explain different conformational behavior of DNA. Calculated self-diffusion coefficients for the ions agree well with the available NMR data.     

Conformational and hydration effects of site-selective sodium,calcium and strontium ion binding to the DNA Holliday junction structure d(TCGGTACCGA)(4)

The role of metal ions in determining the solution conformation of the Holliday junction is well established, but to date the picture of metal ion binding from structural studies of the four-way DNA junction is very incomplete. Here we present two refined structures of the Holliday junction formed by the sequence d(TCGGTACCGA) in the presence of Na(+) and Ca(2+), and separately with Sr(2+) to resolutions of 1.85A and 1.65A, respectively. This sequence includes the ACC core found to promote spontaneous junction formation, but its structure has not previously been reported. Almost complete hydration spheres can be defined for each metal cation. The Na(+) sites, the most convincing observation of such sites in junctions to date, are one on either face of the junction crossover region, and stabilise the ordered hydration inside the junction arms. The four Ca(2+) sites in the same structure are at the CG/CG steps in the minor groove. The Sr(2+) ions occupy the TC/AG, GG/CC, and TA/TA sites in the minor groove, giving ten positions forming two spines of ions, spiralling through the minor grooves within each arm of the stacked-X structure. The two structures were solved in the two different C2 lattices previously observed, with the Sr(2+) derivative crystallising in the more highly symmetrical form with two-fold symmetry at its centre. Both structures show an opening of the minor groove face of the junction of 8.4 degrees in the Ca(2+) and Na(+) containing structure, and 13.4 degrees in the Sr(2+) containing structure. The crossover angles at the junction are 39.3 degrees and 43.3 degrees, respectively. In addition to this, a relative shift in the base pair stack alignment of the arms of 2.3A is observed for the Sr(2+) containing structure only. Overall these results provide an insight into the so-far elusive stabilising ion structure for the DNA Holliday junction.     

Molecular dynamics simulations of A. T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A. T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na(+) intrusions in the first layer of hydration of the minor groove, with long residence times, up to approximately 3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove.     

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.     

Localization of ammonium ions in the minor groove of DNA duplexes in solution and the origin of DNA A-tract bending

Monovalent cation binding sites on nucleic acids in solution can be localized using the isotopically labeled ammonium ion (15NH4+) as a probe in high resolution NMR spectroscopy experiments. The application of this technique to a series of DNA duplexes reveals a preference for the binding of ammonium cations in the minor groove of A-tract sequences. These results are consistent with a recent report which indicates that some solvent electron densities previously identified as water molecules in DNA X-ray crystal structures are partially occupied by sodium ions. The sequence-specific nature of monovalent cation binding sites demonstrated here for A-tract DNA provides an explanation for the origin of sequence-directed bending.     

Binding Sites of the Polyamines Putrescine,Cadaverine, Spermidine and Spermine on A- and B-DNA Located by Simulated Annealing

Abstract

Molecular dynamics simulations with simulated annealing are performed on polyamine-DNA systems in order to determine the binding sites of putrescine, cadaverine, spermidine and spermine on A- and B-DNA. The simulations either contain no additional counterions or sufficient Na⁺ ions, together with the charge on the polyamine, to provide 73% neutralisation of the charges on the DNA phosphates. The stabilisation energies of the complexes indicate that all four polyamines should stabilise A-DNA in preference to B-DNA, which is in agreement with experiment in the case of spermine and spermidine, but not in the case of putrescine or cadaverine. The major groove is the preferred binding site on A-DNA of all the polyamines. Putrescine and cadaverine tend to bind to the sugar-phosphate backbone of B-DNA, whereas spermidine and spermine occupy more varied sites, including binding along the backbone and bridging both the major and minor grooves.     

Liquid crystalline phase behavior of high molecular weight DNA: a comparative study of the influence of metal ions of different size, charge and binding mode

The ability of Li(+), Na(+), K(+), Rb(+), Cs(+), Mg(2+), Ca(2+), Sr(2+), Ba(2+), Cu(2+), Cd(2+), Al(3+), V(4+), Hg(2+), Pd(2+), Au(3+), and Pt(4+) to provoke liquid crystalline (LC) phases in high molecular weight DNA was investigated. The alkali and alkaline earth metal ions provoked typical cholesteric/columnar structures, whereas transition metal ions precipitated DNA into solid/translucent gel-like aggregates. Heavy metal ions reduced viscosity of DNA solution, disrupting rigid, rod-like DNA structure necessary for LC textures. Three-layer quantum mechanical-molecular mechanical (QM/MM) studies of Li(+), Na(+), K(+), Mg(2+), and Ca(2+) binding DNA fragment suggested several possible binding modes of these ions to the phosphate groups. The dianion mode of metal binding, involving the phosphate groups of both strands of DNA, allowed for higher DNA binding affinity of the alkaline earth metal ions. These results have implications in understanding the biological role of metal ions and developing DNA-based sensors and nanoelectronic devices.     

Critical amount of oligovalent ion binding required for the B-Z transition of poly (dG-m5dC)

By working at very low Na+ concentrations (1mM and less), the number of bound Mg(2+), cobalt hexamine (3+), and spermine (4+) necessary to induce the B-Z transition of poly (dG-m5dC) has been directly measured. The results show that if as little as 1 cobalt hexamine(3+) or spermine(4+) is bound per 40-50 nucleotides the transition will occur. A greater fraction of bound Mg(2+) is required, 1 bound Mg(2+)/10 nucleotides. The dependence of the transition midpoint concentrations of oligovalent ions on Na+ concentrations suggests that specific ion binding energies, not included in counterion condensation theory, are responsible for the transition.     

Raman microspectroscopic study of effects of Na(I) and Mg(II) ions on low pH induced DNA structural changes

In this work a confocal Raman microspectrometer is used to investigate the influence of Na(+) and Mg(2+) ions on the DNA structural changes induced by low pH. Measurements are carried out on calf thymus DNA at neutral pH (7) and pH 3 in the presence of low and high concentrations of Na(+) and Mg(2+) ions, respectively. It is found that low concentrations of Na(+) ions do not protect DNA against binding of H(+). High concentrations of monovalent ions can prevent protonation of the DNA double helix. Our Raman spectra show that low concentrations of Mg(2+) ions partly protect DNA against protonation of cytosine (line at 1262 cm(-1)) but do not protect adenine and guanine N(7) against binding of H(+) (characteristic lines at 1304 and 1488 cm(-1), respectively). High concentrations of Mg(2+) can prevent protonation of cytosine and protonation of adenine (disruption of AT pairs). By analyzing the line at 1488 cm(-1), which obtains most of its intensity from a guanine vibration, high magnesium salt protect the N(7) of guanine against protonation. A high salt concentration can prevent protonation of guanine, cytosine, and adenine in DNA. Higher salt concentrations cause less DNA protonation than lower salt concentrations. Magnesium ions are found to be more effective in protecting DNA against binding of H(+) as compared with calcium ions presented in a previous study. Divalent metal cations (Mg(2+), Ca(2+)) are more effective in protecting DNA against protonation than monovalent ions (Na(+)).     

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

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

Use of the parmbsc0 force field and trajectory analysis to study the binding of netropsin to the DNA fragment (5'CCAATTGG)(2) in the presence of excess NaCl salt in aqueous solution

The parmbsc0 force field was applied to study in detail the binding of netropsin, at a salt concentration of 0.28M Na(+), to the minor groove of an 8-mer (5'CCAATTGG)(2) DNA duplex forming a netropsin·DNA complex which previously has been characterized by X-ray crystallography, albeit with the use of closely related DNA duplexes. The X-ray structure revealed that the terminal guanidinium and amidinium groups of netropsin interact with the extreme ends of the palindromic AATT sequence of the receptor DNA. The parmbsc0 parameters of B-DNA and AMBER v9 parameters of netropsin generated a stable 6ns molecular dynamics (MD) trajectory for a 1:1 class I binding motif of this complex. Trajectory analysis for the salt and hydration effects on the binding of netropsin to the 8-mer DNA duplex revealed that 18 water molecules and 2 Na(+) are displaced from the DNA upon netropsin binding. A hydration density map of the complex parallels the X-ray data showing that two structured water molecules are localized near the netropsin guanidinium and amidinium groups forming H-bond bridges between the receptor and the ligand.     

DNA structure control by polycationic species: polyamine, cobalt ammines, and di-metallo transition metal chelates

Many polycationic species bind to DNA and induce structural changes. The work reported here is the first phase of a program whose long-term aim is to create a class of simple and inexpensive sequence-selective compounds that will enable enhanced DNA structure control for a wide range of applications. Three classes of molecule have been included in this work: the polyamine spermine (charge: 4(+)) and spermidine (charge: 3(+)) (which are known to induce a wide range of DNA conformational changes but whose binding modes are still not well understood); cobalt (III) amine transition metal complexes as potential polyamine mimics and [Fe(H(2)O)(6)](3+); and the first member of a new class of di-metallo tris-chelated cylinders of helical structure (charge 4(+)). Temperature-dependent absorption, circular dichroism, linear dichroism, gel electrophoresis, and molecular modeling data are presented. The cobalt amines prove to be effective polyamine mimics, although their binding appears to be restricted to backbone and major groove. All the ligands stabilize the DNA, but the 4(+) di-iron tris-chelate does so comparatively weakly and seems to have a preference for single-stranded DNA. All the molecules studied bend the DNA, with the di-iron tris-chelate having a particularly dramatic effect even at very low drug load.     

Structural analysis of DNA interactions with biogenic polyamines and cobalt(III)hexamine studied by Fourier transform infrared and capillary electrophoresis

Biogenic polyamines, such as putrescine, spermidine, and spermine are small organic polycations involved in numerous diverse biological processes. These compounds play an important role in nucleic acid function due to their binding to DNA and RNA. It has been shown that biogenic polyamines cause DNA condensation and aggregation similar to that of inorganic cobalt(III)hexamine cation, which has the ability to induce DNA conformational changes. However, the nature of the polyamine.DNA binding at the molecular level is not clearly established and is the subject of much controversy. In the present study the effects of spermine, spermidine, putrescine, and cobalt(III)hexamine on the solution structure of calf-thymus DNA were investigated using affinity capillary electrophoresis, Fourier transform infrared, and circular dichroism spectroscopic methods. At low polycation concentrations, putrescine binds preferentially through the minor and major grooves of double strand DNA, whereas spermine, spermidine, and cobalt(III)hexamine bind to the major groove. At high polycation concentrations, putrescine interaction with the bases is weak, whereas strong base binding occurred for spermidine in the major and minor grooves of DNA duplex. However, major groove binding is preferred by spermine and cobalt(III)hexamine cations. Electrostatic attractions between polycation and the backbone phosphate group were also observed. No major alterations of B-DNA were observed for biogenic polyamines, whereas cobalt(III)hexamine induced a partial B --> A transition. DNA condensation was also observed for cobalt(III)hexamine cation, whereas organic polyamines induced duplex stabilization. The binding constants calculated for biogenic polyamines are K(Spm) = 2.3 x 10(5) M(-1), K(Spd) = 1.4 x 10(5) M(-1), and K(Put) = 1.02 x 10(5) M(-1). Two binding constants have been found for cobalt(III)hexamine with K(1) = 1.8 x 10(5) M(-1) and K(2) = 9.2 x 10(4) M(-1). The Hill coefficients indicate a positive cooperativity binding for biogenic polyamines and a negative cooperativity for cobalt(III)hexamine.     

EPR spectroscopic analysis of TAR RNA-metal ion interactions

Metal ion-induced changes in HIV-1 TAR RNA internal dynamics were determined by the changes in EPR spectral width for TAR RNAs containing spin-labeled nucleotides (U23, U25, U38, and U40). This gave a dynamic signature for each of 10 metal ions studied, which fell into one of three distinct groups. While Li(+) and K(+) had little effect on TAR RNA internal dynamics, Na(+) unexpectedly had a dynamic signature that was similar to Ca(2+) and Sr(2+), with a decrease in mobility at U23 and U38, little or no change at U25, and an increase in mobility at U40. Mg(2+), Co(2+), Ni(2+), Zn(2+), and Ba(2+) had similar effects on U23, U38, and U40, but the mobility of U25 was markedly increased. Our results show that RNA dynamics change upon metal binding to the TAR RNA bulge, indicating that RNA structure adapts to accommodate metal ions of different size and coordination properties.