The Effect of Metal Binding to the N7 Site of Purine Nucleotides on Their Structure, Energy, and Involvement in Base Pairing

Infrared absorption detection of metal ion-deoxyguanosine monophosphate binding: experimental and theoretical study

Metal ion interactions with nucleic acids attract attention because of the environmental and biological consequences. The formation of the complex is often monitored by the vibrational spectroscopy. To identify characteristic binding patterns and marker bands on a model DNA component, infrared absorption spectra of the deoxyguanosine monophosphate complexes with Na(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), Zn(2+), and Cd(2+) cations were recorded and interpreted on the basis of density-functional computations. The aqueous environment was simulated by continuum and combined continuum-explicit solvent models. For the binding to the N7 position of the guanine base, the computation predicted a characteristic frequency upshift and splitting of the 1578 cm(-1) band, which is in accord with available experimental data. Contrary to the expectation, the modeling suggests that the binding to the carbonyl group might not be detectable, as the metal causes smaller spectral changes if compared to the hydrogen-bound water molecules. The binding to the phosphate group causes significant spectral changes in the sugar-phosphate vibrating region ( approximately 800-1200 cm(-1)), but also notable frequency shifts of the carbonyl vibrations. The Cu(2+) and Zn(2+) ions induced the largest alterations in observed vibrational absorption, which corresponds to the calculated strong interaction energies in the N7-complexes and to previous experimental experience. Additional changes in the vibrational spectra of the copper complexes were observed under high metal concentration, corresponding to the simultaneous binding to the phosphate residue. The two-step Cu(2+) binding process was also confirmed by the microcalorimetry titration curve. The computations and combination of more techniques thus help us to assign and localize spectral changes caused by the metal ion binding to nucleic acids.

General base catalysis for cleavage by the active-site cytosine of the hepatitis delta virus ribozyme: QM/MM calculations establish chemical feasibility

The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.

Interaction of the DNA bases and their mononucleotides with pyridine-2-carbaldehyde thiosemicarbazonecopper(II) complexes. Structure of the cytosine derivative

Experimental studies of the binding interactions of [CuL(NO(3))] and [{CuL'(NO(3))}(2)] (HL=pyridine-2-carbaldehyde thiosemicarbazone, and HL'=pyridine-2-carbaldehyde 4N-methylthiosemicarbazone) with adenine, guanine, cytosine, thymine and their mononucleotides (dNMP), 2-deoxyadenosine-5'-monophosphate, (dAMP), 2'-deoxyguanosine-5'-monophosphate, (dGMP), 2'-deoxycytidine-5'-monophosphate (dCMP), and thymidine-5'-monophosphate (dTMP) have been carried out in aqueous solution at pH 6.0, I=0.1M (NaClO(4)) and T=25 degrees C. The complexation constants of these compounds, calculated by Hildebrand-Benesi plots for the dye binding, D, ([CuL] or [CuL']) to the nucleobases or nucleotides (P), have shown two linear stretches in adenine, guanine, dAMP and dGMP. The data were analyzed in terms of formation of 1:1 DP and 1:2 DP(2) complexes with increasing purine base or nucleotide content. For cytosine and dCMP only 1:1 complexes have been observed, whereas for thymine and dTMP such complex structures were not observed. The [CuL(Hcyt)](ClO(4)) cytosine derivative has been isolated and characterized. The crystal structure consists of perchlorate ions and [CuL(Hcyt)](+) monomers attached by hydrogen bond, chelate pi-ring and anion-pi interactions. The Cu(2+) ions bind to the NNS chelating moiety of the thiosemicarbazone ligand and the cytosine N13 site (N3, most common notation) yielding a square-planar geometry. A pseudocoordination to the cytosine O12 site (=O2) can also be considered.

Influence of magnesium ions on spontaneous opening of DNA base pairs

A large amount of experimental evidence is available for the effects of magnesium ions on the structure and the stability of the DNA double helix. Less is known, however, on how these ions affect the dynamics of the molecule and the stability of each individual base pair. The present work addresses these questions by a study of the DNA duplex [dCGCAGATCTGCG]2, and its interactions with magnesium ions using nuclear magnetic resonance (NMR) spectroscopy and proton exchange. Two-dimensional NMR experiments indicate that binding of magnesium to this DNA duplex does not affect its structure. However, even in the absence of structural changes, magnesium ions specifically affect the exchange properties of imino protons in the four GC/CG base pairs that are located in the interior of the double helix. These specific changes do not result from alterations in the rates of spontaneous opening of these base pairs. Instead, the changes most likely reflect an enhancement in the energetic propensity for spontaneous opening of the GC/CG base pairs that is induced by the binding of magnesium ions.

Extracting stoichiometry, thermodynamics, and kinetics for the interaction of DNA with cadmium ion by capillary electrophoresis on-line coupled with electrothermal atomic absorption spectrometry

Understanding the binding of cadmium with DNA is of great importance for elucidating the mechanism of cadmium genotoxicity and carcinogenicity. In the present work, CE on-line coupled with electrothermal atomic absorption spectrometry was employed to study the binding electrophoretic behaviors, stoichiometry, thermodynamics, and kinetics for the interaction of cadmium cation (Cd(II)) with DNA. The stoichiometry (Cd(II) to DNA (as the concentration of base pairs)) for the interaction was determined to be 1:5. Two types of binding sites on DNA were observed with the binding constants of 10(6) and 10(5) L/mol, respectively, showing strong affinity of Cd(II) to DNA. The interaction of Cd(II) with both types of binding sites on DNA were driven by negative enthalpy change with a large positive entropy change. The binding of Cd(II) to DNA followed a first-order kinetics for Cd(II) with the apparent activation energy of 45.7 +/- 1.9 kJ/mol. The results obtained in present investigation would be helpful to understanding the genotoxicity and carcinogenicity of cadmium.

Cu2+/+ cation coordination to adenine--thymine base pair. Effects on intermolecular proton-transfer processes

Intermolecular proton-transfer processes in the Watson & Crick adenine-thymine Cu+ and Cu2+ cationized base pairs have been studied using the density functional theory (DFT) methods. Cationized systems subject to study are those resulting from cation coordination to the main basic sites of the base pair, N7 and N3 of adenine and O2 of thymine. For Cu+ coordinated to N7 or N3 of adenine, only the double proton-transferred product is found to be stable, similarly to the neutral system. However, when Cu+ interacts with thymine, through the O2 carbonyl atom, the single proton transfer from thymine to adenine becomes thermodynamically spontaneous, and thus rare forms of the DNA bases may spontaneously appear. For Cu2+ cation, important effects on proton-transfer processes appear due to oxidation of the base pair, which stabilizes the different single proton-transfer products. Results for hydrated systems show that the presence of the water molecules interacting with the metal cation (and their mode of coordination) can strongly influence the ability of Cu2+ to induce oxidation on the base pair.

DNA A-tracts bending: polarization effects on electrostatic interactions across their minor groove

Bending by the DNA A-tracts constitutes a contentious issue, suggesting deficiencies in the physics employed so far. Here, we inquire as to the importance in this bending of many-body polarization effects on the electrostatic interactions across their narrow minor groove. We have done this on the basis of the findings of Jarque and Buckingham who developed a procedure based on a Monte Carlo simulation for two charges of the same sign embedded in a polarizable medium. Remarkably, the present analysis reveals that for compact DNA conformations, which result from dynamic effects, an overall attractive interaction operates between the phosphate charges; this interaction is especially strong for the narrow minor groove of the A-tracts, suggesting a tendency for DNA to bend toward this groove. This tendency is in agreement with the conclusions of electrophoretic and NMR solution studies. The present analysis is also consistent with the experimental observations that the minor groove is much more easily compressible than the major groove and the bending propensity of the A-tracts is greatly reduced at "premelting" temperatures. By contrast, the dielectric screening model predicts a repulsion between the phosphate charges and is not consistent with the aforementioned bending tendency or experimental observations.

Tautomeric equilibrium, stability, and hydrogen bonding in 2'-deoxyguanosine monophosphate complexed with Mg2+

The tautomeric equilibrium and hydrogen bonding in nucleotide 2'-deoxyguanosine monophosphate that interacts with hydrated Mg2+ cation (4H2O.Mg[dGMP]) were studied at the MP2/cc-pVDZ//B3LYP/cc-pVDZ and B3LYP/aug-cc-pVTZ//B3LYP/cc-pVDZ levels of theory. The Mg2+ ion forms two inner-shell contacts with the nucleotide, similar to small phosphorylated molecules under physiological conditions. The presence of the phosphate group and the hydrated magnesium cation leads to a change in guanine tautomeric equilibrium of 4H2O.Mg[dGMP] in comparison to free guanine. The influence of the phosphate group and the magnesium cation on tautomeric equilibrium is larger in the anti conformation where the P=O-->Mg and Mg<--N7 coordinate bonds are formed. The canonical oxo form of guanine is more stable (by 6-8 kcal/mol) than the O6-hydroxo form in anti conformation. Thus, the interaction with Mg2+ ion is capable of further suppressing the likelihood of a spontaneous transient formation of the rare tautomer. In the syn conformation of 4H2O.Mg[dGMP], the interaction of the guanine nucleobase with the phosphate group and the magnesium cation is not as strong as in the anti conformation, and the relative stability of guanine tautomers is close to those in free guanine.

Molecular dynamics simulations and their application to four-stranded DNA

This review provides a critical assessment of the advantages and limitations of modeling methods available for guanine quadruplex (G-DNA) molecules. We characterize the relations of simulations to the experimental techniques and explain the actual meaning and significance of the results. The following aspects are discussed: pair-additive approximation of the empirical force fields, sampling limitations stemming from the simulation time and accuracy of description of base stacking, H-bonding, sugar-phosphate backbone and ions by force fields. Several methodological approaches complementing the classical explicit solvent molecular dynamics simulations are commented on, including enhanced sampling methods, continuum solvent methods, free energy calculations and gas phase simulations. The successes and pitfalls of recent simulation studies of G-DNA are demonstrated on selected results, including studies of cation interactions and dynamics of G-DNA stems, studies of base substitutions (inosine, thioguanine and mixed tetrads), analysis of possible kinetic intermediates in folding pathway of a G-DNA stem and analysis of loop regions of G-DNA molecules.

Structure and bonding of Ag(I)-DNA base complexes and Ag(I)-adenine-cytosine mispairs: An ab Initio study

High level ab initio calculations have been carried out to characterize the structure, bonding and energetics of Ag(I)-DNA base complexes, including adenine or cytosine, as well as Ag(I)-adenine-cytosine mispairs. The interactions of the Ag cation in all binding sites of all adenine and cytosine tautomers have been considered. The calculations show that in gas phase the canonical form of cytosine is stabilized upon metalation, whereas the lowest energy structure of Ag-adenine correspond to a rare tautomer. Interestingly, the theoretical inspection of metalated adenine-cytosine mispair reveals that the most stable structures are formed from the canonical cytosine and adenine tautomers. The lowest energy structure is planar with adenine and cytosine hydrogen-bonded. Within few kcal/mol nonplanar, conformationally very flexible structures are found, in which the Ag(I) crosslinks an endocyclic nitrogen of adenine and the oxygen of cytosine. Metalated reverse-Wobble type of structures, on the contrary, are predicted much higher in energy.

Strength of the Zn−N Coordination Bond in Zinc Porphyrins on the Basis of Experimental Thermochemistry

The compound, 5,10,15,20-tetrakis(4-methoxyphenyl)porphine zinc(II) (ZnTMPP), was prepared, and its thermochemical properties were experimentally established. The standard molar energy of combustion (Delta(c)U degrees m) was determined from oxygen rotating-bomb combustion calorimetry experiments. The standard molar enthalpies of combustion (Delta(c)H degrees m) and formation (Delta(f)H degrees m) were derived. The enthalpy of sublimation (Delta(cr)(g)H degrees m) was determined by Knudsen effusion at high temperatures. With these results, the standard molar enthalpies of formation and atomization (Delta(at)H degrees m) in the gas state were calculated. A summary of the results at T = 298.15 K (p degrees = 0.1 MPa) is shown in Table 1. Using these results and those previously obtained for the free ligand, 5,10,15,20-tetrakis(4-methoxyphenyl)porphine, the mean dissociation enthalpy for the Zn-N coordination bond is obtained as D(Zn-N) = (160 +/- 9) kJ.mol-1. This value is consistent with the results obtained using the same experimental approach in a similar system (5,10,15,20-tetraphenylporphine, TPP/ZnTPP) reported elsewhere. A discussion of the strength for the Zn-N coordination bond is made in terms of the structural and electronic features of the molecules involved.

To be dinuclear or not: Z-DNA induction by nickel complexes

The left-handed Z-DNA has been identified as a gene regulating element. Therefore the generation of Z-DNA through metal complexes might be an innovative way for the regulation of gene expression. Use of the new dinuclear complex N,N,N',N'-tetrakis-[2-(3,5-dimethylpyrazol-1-yl)ethyl]-1,3-propylenediamine-bis(nickel(II) dinitrate) (2) reversibly induced Z-DNA formation. However, when a 1:1 ratio of metal/dinucleating ligand was used as a control, the midpoint of the B- to Z-DNA transition was at the same nickel concentration as in case of the dinuclear complex. The novel mononuclear analogue, N-methyl-N,N-bis-[2-(3,5-dimethylpyrazol-1-yl)ethyl]amine-nickel(II)-dinitrate (3) was inducing the Z-DNA at a similar ratio versus nucleotides as free nickel(II) itself. For the first time, proton and nickel binding constants for the bis-[2-(pyrazol-1-yl)ethyl]amine ligand system are reported and discussed. Both nickel complexes 2 and 3 were structurally characterized by single crystal analysis. Furthermore, the synthesis of the two new ligands, N,N,N',N'-tetrakis-[2-(3,5-dimethylpyrazol-1-yl)ethyl]-1,2-propylenediamine (4) and N-methyl-N,N-bis-[2-(3,5-dimethylpyrazol-1-yl)ethyl]amine (5) is described. The two major synthetic pathways leading to polypyrazoyl amines in general are critically discussed with respect to yield, reproducibility and handling of the intermediates.

On the Bonding of First-Row Transition Metal Cations to Guanine and Adenine Nucleobases

The binding of first-row transition metal monocations (Sc+-Cu+) to N7 of guanine and N7 or N3 of adenine nucleobases has been analyzed using the hybrid B3LYP density functional theory (DFT) method. The nature of the bonding is mainly electrostatic, the electronic ground state being mainly determined by metal-ligand repulsion. M+-guanine binding energies are 18-27 kcal/mol larger than those of M+-adenine, the difference decreasing along the row. Decomposition analysis shows that differences between guanine and adenine mainly arise from Pauli repulsion and the deformation terms, which are larger for adenine. Metal cation affinity values at this level of calculation are in very good agreement with experimental data obtained by Rodgers et al. (J. Am. Chem. Soc. 2002, 124, 2678) for adenine nucleobases.

Tautomerism of Uracil and 5-Bromouracil in a Microcosmic Environment with Water and Metal Ions. What Roles Do Metal Ions Play?

The base tautomerization processes of uracil/5-bromouracil were investigated in a microcosmic environment with both H2O and Na+ (W-M environment). It was found that uracil was more stable in the W-M environment than in the microcosmic environment with only water, which suggested that the metal ions and water work cooperated to maintain the classical nucleic acid bases. However, 5-bromouracil, a chemical mutagen, was found to be less stable than uracil in the W-M environment. Why the 5-bromouracil is easier to tautomerize and therefore induce gene mutation was explained to some extent. Further research revealed that the water molecule would assist the tautomerization in the W-M environment. However, the metal ions in different regions play absolutely opposite roles: in one region, the metal ions can prevent the base from tautomerizing, whereas in another region, the metal ion can assist the tautomerization process. Furthermore, from the viewpoint of ionization of the base, it seems BrU has a stronger tendency to lose the proton at N3, which is an intrinsic consequence of the bromine atom and is not affected by the metal cation.

Polarizable model potential function for nucleic acid bases

A polarizable model potential (PMP) function for adenine (A), cytosine (C), guanine (G), thymine (T), and uracil (U) is developed on the basis of ab initio molecular orbital calculations at the MP2/6-31+G* level. The PMP function consists of Coulomb, van der Waals, and polarization terms. The permanent atomic charges of the Coulomb term are determined by using electrostatic potential (ESP) optimization. The multicenter polarizabilities of the polarization term are determined by using polarized one-electron potential (POP) optimization in which the electron density changes induced by a test charge are target. Isotropic and anisotropic polarizabilities are adopted as the multicenter polarizabilities. In the PMP calculations using the optimized parameters, the interaction energies of Watson-Crick type A-T and C-G base pairs were -15.6 and -29.4 kcal/mol, respectively. The interaction energy of Hoogsteen type A-T base pair was -17.8 kcal/mol. These results reproduce well the quantum chemistry calculations at the MP2/6-311++G(3df,2pd) level within the differences of 0.6 kcal/mol. The stacking energies of A-T and C-G were -9.7 and -10.9 kcal/mol. These reproduce well the calculation results at the MP2/6-311++G (2d,2p) level within the differences of 1.3 kcal/mol. The potential energy surfaces of the system in which a sodium ion or a chloride ion is adjacent to the nucleic acid base are calculated. The interaction energies of the PMP function reproduced well the calculation results at the MP2/6-31+G* or MP2/6-311++G(2d,2p) level. The reason why the PMP function reproduces well the high-level quantum mechanical interaction energies is addressed from the viewpoint of each energy terms.

Interguanine hydrogen-bonding patterns in adducts with water and Zn–purine complexes (purine is 9-methyladenine and 9-methylguanine). Unexpected preference of Zn(II) for adenine-N7 over guanine-N7

Guanine-guanine hydrogen bonding involving the Watson-Crick edge [N(1)H, N(2)H2] of one base and the Hoogsteen edge (N7, O6) of the other is the dominant association pattern in the solid-state structures of two hydrates of 9-ethylguanine (9-EtGH), and in adducts of 9-methylguanine (9-MeGH) with the Zn compounds [ZnCl2(H2O)(9-MeGH-N7)]*(9-MeGH) as well as [ZnCl2(H2O)(9-MeA-N7)]*2(9-MeGH) (9-MeA is 9-methyladenine). The structures of 9-EtGH*2H2O and 9-EtGH*3.5H2O are dominated by polymeric tape structures of the guanine and extended water clusters. In [ZnCl2(H2O)(9-MeGH-N7)]*(9-MeGH) the metalated guanine is involved in hydrogen bonding (GG3 motif) with a free 9-MeGH, which in turn is centrosymmetrically related to itself via hydrogen bonds involving N2H2 and N3 (GG4 motif). In [ZnCl2(H2O)(9-MeA-N7)]*2(9-MeGH) the metalated adenine base interacts via its Watson-Crick edge [N1, N(6)H2] with the sugar edge [N(2)H2, N3] of one of the guanine nucleobases of the GG pair. Crystallization of [ZnCl2(H2O)(9-MeA-N7)]*2(9-MeGH) from an aqueous solution containing 9-MeGH, 9-MeA, and ZnCl2 is fully unexpected in that the anticipated preference of Zn(II) for guanine-N7 is not realized and instead coordination to adenine-N7 is observed. The relevance of [ZnCl2(H2O)(9-MeGH-N7)]*(9-MeGH) and [ZnCl2(H2O)(9-MeA-N7)]*2(9-MeGH) for metal-containing nucleic acid triplex structures is discussed.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

Double proton transfer in the isolated and DNA-embedded guanine-cytosine base pair

The energetics and dynamics of double proton transfer (DPT) is investigated theoretically for the Watson-Crick conformation of the guanine-cytosine (GC) base pair. Using semiempirical density functional theory the isolated and DNA-embedded GC pair is considered. Differences in the energetics and dynamics of DPT thus addresses the question of how relevant studies of isolated base pairs are for the understanding of processes occurring in DNA. Two-dimensional potential energy surfaces involving the transferring hydrogen atoms and the proton donors and acceptors are presented for both systems. The DPT reaction is accompanied by a contraction of the distance between the two bases with virtually identical energetic barriers being 18.8 and 18.7 kcal/mol for the isolated and DNA-embedded system, respectively. However, the transition state for DPT in the DNA-embedded GC pair is offset by 0.1 A to larger N-H separation compared to the isolated GC pair. Using activated ab initio molecular dynamics, DPT is readily observed for the isolated base pair with a minimal amount of 21.4 kcal/mol of initial average kinetic energy along the DPT normal mode vector. On a time scale of approximately 100 fs DPT has occurred and the excess energy is redistributed. For the DNA-embedded GC pair considerably more kinetic energy is required (30.0 kcal/mol) for DPT and the process is completed within one hydrogen vibration. The relevance of studies of isolated base pairs and base pair analogs in regard of reactions or properties involving DNA is discussed.

Quantification of single-stranded nucleic acid and oligonucleotide interactions with metal ions by affinity capillary electrophoresis: part I

The interactions between oligonucleotides and inorganic cations have been measured by capillary zone electrophoresis. With increasing concentrations of divalent cations (Ca(2+), Mg(2+), Mn(2+) and Ni(2+)) in the running buffer, the migration behavior was evaluated by calculation of the binding constants. Besides these fundamental studies of binding equilibria, different buffer components, tris(hydroxymethyl)aminomethane and 3-(N-morpholino)propanesulfonic acid, have been investigated and their effects on metal ion binding quantified.

RNA unrestrained molecular dynamics ensemble improves agreement with experimental NMR data compared to single static structure: a test case

Nuclear magnetic resonance (NMR) provides structural and dynamic information reflecting an average, often non-linear, of multiple solution-state conformations. Therefore, a single optimized structure derived from NMR refinement may be misleading if the NMR data actually result from averaging of distinct conformers. It is hypothesized that a conformational ensemble generated by a valid molecular dynamics (MD) simulation should be able to improve agreement with the NMR data set compared with the single optimized starting structure. Using a model system consisting of two sequence-related self-complementary ribonucleotide octamers for which NMR data was available, 0.3 ns particle mesh Ewald MD simulations were performed in the AMBER force field in the presence of explicit water and counterions. Agreement of the averaged properties of the molecular dynamics ensembles with NMR data such as homonuclear proton nuclear Overhauser effect (NOE)-based distance constraints, homonuclear proton and heteronuclear (1)H-(31)P coupling constant (J) data, and qualitative NMR information on hydrogen bond occupancy, was systematically assessed. Despite the short length of the simulation, the ensemble generated from it agreed with the NMR experimental constraints more completely than the single optimized NMR structure. This suggests that short unrestrained MD simulations may be of utility in interpreting NMR results. As expected, a 0.5 ns simulation utilizing a distance dependent dielectric did not improve agreement with the NMR data, consistent with its inferior exploration of conformational space as assessed by 2-D RMSD plots. Thus, ability to rapidly improve agreement with NMR constraints may be a sensitive diagnostic of the MD methods themselves.

Synthesis, structure and biological activity of a new and efficient Cd(II)–uracil derivative complex system for cleavage of DNA

The new complex formed by Cd(II) and the 1:2 Schiff-base-type ligand 2,6-bis[1-(4-amino-1,2,3,6-tetrahydro-1,3-dimethyl-2,6-dioxopyrimidin-5-yl)imino]ethylpyridine (DAPDAAU) has been chemically and structurally characterized by X-ray diffraction: the ion Cd(II) is surrounded by six nitrogen atoms from two DAPDAAU ligands which coordinates each one in a tridentate fashion through the pyridine ring (N1) and both azomethine nitrogen atoms (N5). The interaction of the Cd(II) complex (compound I) with calf-thymus DNA as observed by circular dichroism spectroscopy suggests the initial unwinding of the DNA double helix strongly depends on increasing incubation times and metal-to-nucleic acid molar ratios. Electrophoretic experiments indicate that the cadmium complex induces cleavage of the plasmid pBR322 DNA to give ulterior nicking and shortening of this molecule, as a result of the complex binding to DNA, resulting in the conclusion that compound I behaves as a chemical nuclease. Cytotoxic activity of the Cd(II) complex against selected different human cancer cell lines is specific and increases with increasing concentration of the metal compound; this fact indicates the potential antitumor character of the complex. When the culture medium is supplemented with compound I, a remarkable inhibition of the growing cell is observed, important cell degeneration appears before 48 h and abundant precipitates are formed that correspond to cell residues and denatured proteins.

Stabilization of the non-canonical adenine–adeninium base pair by N(7) coordination of Zn(II)

A new zinc (II) compound with 9-ethyladenine (9-EtA) of formula [Zn(9-EtA-N7)Cl(3)](9-EtAH) has been synthesized and characterized by X-ray diffraction. Its X-structure consists of an Zn(II) anionic complex and 9-ethyladeninium as counteranion. The Zn(II) complex shows a distorted tetrahedral geometry in which three Cl and an 9-EtA coordinates through N(7) position are the ligands. An indirect chelation via intramolecular H-bond between N(6)H and an Cl ligand is present in the complex. The network of [Zn(9-EtA-N7)Cl(3)](9-EtAH) shows interesting features. Thus, self-association of coordinated adenine-adeninium takes place by H-bonding of N(6)-H...N(1) and N(6)-H...N(7), leading to a polymeric ribbon-like 1D supramolecular arrangement. Ab initio calculations have been applied in order to study the stability of the adenine-adeninium interaction due to the coordination of the Zn(II) to the N(7) position and to compare experimental and theoretical structural data.

Benzoderivatives of Nucleic Acid Bases as Modified DNA Building Blocks

The tautomeric properties of benzoderivatives of the canonical nucleic acid bases have been studied by using different computational approaches. Attention has been paid to the impact of the benzene group in altering the tautomeric preferences of the canonical bases both in the gas phase and in aqueous solution. To this end, relative solvation free energies of the tautomers determined from Self-Consistent Reaction Field continuum calculations and Monte Carlo-Free Energy Perturbation are combined with gas-phase tautomerization free energies determined from quantum mechanical calculations. The results provide a detailed picture of the tautomeric preferences of the benzoderivatives of nucleic acid bases. This information is used to examine the recognition properties of the preferred tautomers of the benzo-fused derivatives, paying particular attention to the ability to form Watson-Crick hydrogen-bonding and stacking interactions as well as to the hydrophobic nature of the modified bases. The implications of present results on the potential use of benzo-fused bases as potential building blocks in modified DNA duplexes are examined.

Theoretical study of guanine–Cu and uracil–Cu (neutral, anionic, and cationic). Is it possible to carry out a photoelectron spectroscopy experiment?

The structure and bonding of guanine-Cu and uracil-Cu (neutral, anionic, and cationic) are discussed on the basis of the calculated structures and energies. The interaction of the metal atom with guanine and uracil has been analyzed using the B3LYP density-functional approach. The removal of one electron from the neutral complexes produces the stabilization of one of the isomers, while the addition of one electron leads to a system where the metal atom is weakly bounded to guanine or uracil, according to the metal-bases bond distance that is long (2.29-2.90). For guanine-Cu and uracil-Cu, the vertical ionization energy of the anion is close to the dissociation energy of one hydrogen atom from guanine-Cu or uracil-Cu. In these cases, it could be possible to produce the detachment of one electron from the anion and also the removal of one hydrogen atom. This is important since the photoelectron spectroscopy of atomic or mixed-atomic cluster anions has proven to be a very effective tool in the study of small systems. For the analysis of copper atoms with DNA bases such as guanine and uracil, it is expected that the photoelectron spectra of the anion-bases complexes strongly resemble the spectrum of Cu(-1), just shifted to higher electron binding energies due to the product stabilization. Hopefully, this information will be useful for the experimental groups.

Structural models for the interaction of Cd(II) with DNA: trans-[Cd(9-RGH-N7)2(H2O)4]2+

Reactions of Cd(NO(3))(2) with the model nucleobases 9-alkylguanine in water at neutral pH, give the compounds trans-[Cd(9-RGH-N7)(2)(H(2)O)(4)](NO(3))(2)(R=Me, Et), with the 9-alkylguanine ligands bound to the metal cation at the N(7) position. The X-ray structures of both compounds are reported. The six-coordinate Cd(II) complexes consist of a highly regular octahedral geometry in which the two 9-alkylguanine ligands are in a trans position to each other and approximately collinear with the metal cation. In addition, the networks of both compounds show interesting features. Thus, intramolecular H-bonds between O(6) and a coordinated water molecule are present, and self-association of guanines via H-bonding of N(3)-H...N(2) take place, leading to a 1D supramolecular polymeric ribbon. Density functional theory calculations have been applied to both compounds in order to study the stability of N(7) metalated guanine-guanine associations by comparing experimental and theoretical results. The potential relevance with regard to possible Cd(II)-DNA cross-links is briefly discussed.

Charges of phosphate groups. A role in stabilization of 2'-deoxyribonucleotides. A DFT investigation

We have analyzed the relative stabilities and Gibbs tautomeric free energy for tautomeric transitions of neutral 2'-deoxyribonucleotides and its mono- and di-protonated forms. Geometry optimizations of these nucleic acid constituents have been performed at the DFT/B3LYP level using the standard 6-31G(d) basis set. The prediction of relative stabilities, Gibbs tautomeric free energy has been made at the B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) level of theory. For each nucleoside four major conformers, i.e., north/anti, north/syn, south/anti, and south/syn have been taken into consideration. We have found the substantial effect of the uncompensated charge on the relative stability of 2'-deoxyribonucleotides. In particular, when the charge of 2'-deoxyribonucleotide anions is completely compensated by protons, the syn conformations have been found to be the global minima due to stabilization provided by intramolecular hydrogen bonds. However, the negative charge that appears due to the successive removal of the protons from the phosphate group destabilizes these syn conformations and stabilizes preferably the south/anti conformations (except of 2'-deoxyguanosine phosphate). Only 2'-deoxyribonucleotides, possessing south/anti and north/anti orientations, containing guanine and cytosine can contribute significantly to the rate of spontaneous point mutations due to the formation of biologically relevant amounts of 'rare' tautomers. However, we found strong influence of uncompensated negative charge for 2'-deoxyribonucleotides which possess syn conformations. Finally we have found that the proton transfer could result in the spontaneous change of 2'-deoxyribonucleotides conformations. We conclude that this phenomenon could be considered as a new way for the stabilization of 'rare' isomers for such DNA bases as cytosine and thymine.

Marked Variations of Dissociation Energy and H-Bond Character of the Guanine-Cytosine Base Pair Induced by One-Electron Oxidation and Li+ Cation Coupling

The variation of dissociation energy and H-bond character of the G-C cation and the Li-GC cation have been investigated by employing density functional theory (B3LYP) with the 6-31+G* basis set. The one-electron oxidation and the coupling of Li(+) to the guanine-cytosine base pair can strengthen the interaction between guanine and cytosine. The interaction of the cation Li(+) with guanine is attractive and is attributed to the polarization of the H-bonds between G-C that enhances G-C interaction. The cooperativity of the three H-bonds in the GC and Li-GC cations is different from that in the neutral GC base pair. The proton-transfer process between N(1) of the guanine and N(3) of the cytosine can occur in the GC cation and the Li-GC cation. The geometries of the transition state are out of plane, especially for the transition state of the Li-GC cation. The analysis of the activation energy for the proton-transfer process shows that the GC(+) before and after proton transfer can exist simultaneously in the gas phase, but for the Li-GC(+) system, the Li-GC(+) without proton transfer is the dominating species in the gas phase.

H−Bonding Patterns in the Platinated Guanine−Cytosine Base Pair and Guanine−Cytosine−Guanine−Cytosine Base Tetrad: an Electron Density Deformation Analysis and AIM Study

The atoms in molecule theory (AIM) and electronic structure analysis are applied together to investigate H-bonding patterns in metalated nucleobase complexes. The influence of Pt on the intra GC base pair H-bonding has been found to reduce intra base pair H-bonding of N4(C)...O6(G) in the platinated GC pair and GCGC tetrad. The relaxation of geometry constrains in metalated nucleobases is found to be decisively important in the formation of novel molecular architectures from nucleobases and metal entities. The incorporation of the platinum in the GCGC tetrad benefits the formation of the unique CH...N (H5(C)...N1(G)) hydrogen bond pattern in the tetrad by offering improved geometric constraints rather than through changing the electronic properties around the H5(C) and N1(G) sites. Platination at the N7 of guanine reduces the deprotonation energy considerably.

Theoretical Calculation of the NMR Spin−Spin Coupling Constants and the NMR Shifts Allow Distinguishability between the Specific Direct and the Water-Mediated Binding of a Divalent Metal Cation to Guanine

The calculated intermolecular and intramolecular indirect NMR spin-spin coupling constants and NMR shifts were used for the discrimination between the inner-shell and the outer-shell binding motif of hydrated divalent cations Mg(2+) or Zn(2+) with a guanine base. The intermolecular coupling constants (1)J(X,O6) and (1)J(X,N7) (X = Mg(2+), Zn(2+)) can be unambiguously assigned to the specific inner-shell binding motif of the hydrated cation either with oxygen O6 or with nitrogen N7 of guanine. The calculated coupling constants (1)J(Mg,O6) and (1)J(Zn,O6) were 6.2 and -17.5 Hz, respectively, for the inner-shell complex of cation directly interacting with oxygen O6 of guanine. For the inner-shell coordination of the cation at nitrogen N7, the calculated coupling constants (1)J(Mg,N7) and (1)J(Zn,N7) were 5.6 and -36.5 Hz, respectively. When the binding of the cation is water-mediated, the coupling constant is zero. To obtain reliable shifts in NMR parameters, hydrated guanine was utilized as the reference state. The calculated change of NMR spin-spin coupling constants due to the hydration and coordination of the cation with guanine is caused mainly by the variation of Fermi-contact coupling contribution while the variation of diamagnetic spin-orbit, paramagnetic spin-orbit, and spin-dipolar coupling contributions is small. The change of s-character of guanine sigma bonding, sigma antibonding, and lone pair orbitals upon the hydration and cation coordination (calculated using the Natural Bond Orbital analysis) correlates with the variation of the Fermi-contact term. The calculated NMR shifts delta(N7) of -15.3 and -12.2 ppm upon the coordination of Mg(2+) and Zn(2+) ion are similar to the NMR shift of 19.6 ppm toward the high field measured by Tanaka for N7 of guanine upon the coordination of the Cd(2+) cation (Tanaka, Y.; Kojima, C.; Morita, E. H.; Kasai. Y.; Yamasaki, K.; Ono, A.; Kainosho, M.; Taira, K. J. Am. Chem. Soc. 2002, 124, 4595-4601). The present data indicate that measurements of NMR intermolecular coupling constants may be used to discriminate between the specific inner- and outer-shell binding of divalent cations to nucleobases in DNA and RNA.

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.

Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.