The Influence of N7Guanine Modifications on the Strength of Watson−Crick Base Pairing and Guanine N1Acidity: Comparison of Gas-Phase and Condensed-Phase Trends

Intermolecular hydrogen bonding in DNA base pairs interacting with different numbers of bare and hydrated Li+: NBO, QTAIM, and computational spectroscopic studies

In this study, the effect of different numbers of Li+ interacting with various sites of DNA base pairs (adenine-thymine (AT) and cytosine-guanine (GC)) on the base pair structures, the strength of hydrogen bonding between the bases, and spectroscopic properties (IR and absorption spectra) of the base pairs was investigated. Two quantum computational analyses, the natural bonding orbitals (NBO) and the quantum theory of atoms in molecules (QTAIM), were used to follow the change in the strength of hydrogen bonds between the bases in each pair. The type of base pair's site interacting with Li+ showed different effects on the change in the strength of the hydrogen bonds between the bases. The IR and absorption spectra of the lithiated base pairs were calculated and compared with those of bare base pairs. This comparison provided the changes in the spectra as a fingerprint for the structural identification of different lithiated base pairs. Also, the determination of the change in the strength of hydrogen bonds in the lithiated base pairs compared to their bare base pairs. In the other part of this study, the effect of the hydration of the attached Li+ in the structure of lithiated base pairs on the strength of their hydrogen bonds and spectra was investigated.

Cytotoxic sub-nanometer aqueous platinum clusters as potential antitumoral agents

Ligand-free sub-nanometer metal clusters (MCs) of Pt, Ir, Rh, Au and Cu, are prepared here in neat water and used as extremely active (nM) antitumoral agents for HeLa and A2870 cells. The preparation just consists of adding the biocompatible polymer ethylene-vinyl alcohol (EVOH) to an aqueous solution of the corresponding metal salt, to give liters of a MC solution after filtration of the polymer. Since the MC solution is composed of just neat metal atoms and water, the intrinsic antitumoral activity of the different sub-nanometer metal clusters can now fairly be evaluated. Pt clusters show an IC₅₀ of 0.48 μM for HeLa and A2870 cancer cells, 23 times higher than that of cisplatin and 1000 times higher than that of Pt NPs, and this extremely high cytotoxicity also occurs for cisplatin-resistant (A2870 cis) cells, with a resistance factor of 1.4 (IC₅₀ = 0.68 μM). Rh and Ir clusters showed an IC₅₀ ∼ 1 μM. Combined experimental and computational studies support an enhanced internalization and cytotoxic activation.

On the reactivity of mRNA Cap0: C-H oxidative addition of 7-methylguanosine to Pt0 and base pairing studies

7-Methylguanosine, whose lability to form an ylide/NHC has been known for decades, reacts with [Pt(PPh3)4] via C-H oxidative addition to yield a hydrido-PtII carbene complex. 1H NMR studies on Watson-Crick base-pairs showed no significant effect of metallation.

Assessing the importance of proton transfer reactions in DNA

Although engineered by millions of years of evolution, the cellular machinery is not flawless, and errors regularly appear during DNA replication. The subsequent alteration of the stored genetic message results in a mutation and might be the starting point of important health disorders. The question therefore is what causes DNA mutations? All living organisms are constantly exposed to a number of external agents such as free radicals and to radiation, which may lead to induced mutations. There are also mutations happening without invoking the action of any exogenous element, the so-called spontaneous mutations. The former can be partially controlled by avoiding exposure to high-risk environments, while the latter are more intriguing because their origin is unclear and difficult to determine. As noted by Watson and Crick when they first discovered the DNA structure, the correct replication of DNA rests on the assumption that the base pairs remain in their most stable, canonical form. However, protons along the interbase hydrogen-bond network are not static entities. They can in fact interchange their positions in DNA bases through proton transfer (PT) reactions before strands unwind, giving rise to noncanonical structures defined as rare tautomers. The importance of these rare tautomers was also cleverly anticipated by Watson and Crick and some years later claimed by Löwdin to be a source of spontaneous mutations. In Watson and Crick's words: "It would be of interest to know the precise difference in free energy between the various tautomeric forms under physiological conditions." Unfortunately, rare tautomeric forms are very difficult to detect, so no direct and accurate free energy measure has been discerned. In contrast, theoretical chemistry is making good progress toward the quantification of PT reactions in DNA and their biological consequences. This Account touches upon the theoretical studies devoted to appraising the importance of rare tautomers as promoters of spontaneous mutations. We focus in particular on the crucial role played by the biological environment on DNA stability. It has now been demonstrated that valuable macroscopic predictions require not only highly accurate theories but also refined chemical models. Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations performed on short but complete DNA sequence fragments emerge in this context as the most adequate tools. In addition, these methods can be used to quantify the effect of different external agents on the PT tautomeric equilibria and, eventually, to conveniently handle them. This is the case for the possible alteration of the naturally observed mutation rate by exposure to intense electric fields. Theoretical predictions envision in this respect promising applications of ultrashort electric pulses in medicine to selectively modify the mutated/canonical ratio in DNA.

Mechanism of the cis-[Pt(1R,2R-DACH)(H2O)2]2+ intrastrand binding to the double-stranded (pGpG)·(CpC) dinucleotide in aqueous solution: a computational DFT study

A mechanism of the intrastrand 1,2-cross-link formation between the double-stranded pGpG·CpC dinucleotide (ds(pGpG)) and fully aquated oxaliplatin cis-[Pt(DACH)(H2O)2](2+) (DACH = cyclohexane-1R,2R-diamine) is presented. All structures of the reaction pathways including the transition states (TSs) were fully optimized in water solvent using DFT methodology with dispersion corrections. Both 5' → 3' and 3' → 5' binding directions were considered. In the first step there is a slight kinetic preference for 5'-guanine (5'G) monoadduct formation with an activation Gibbs free energy of 18.7 kcal/mol since the N7 center of the 5'G base is fully exposed to the solvent. On the other hand, the N7 atom of 3'-guanine (3'G) is sterically shielded by 5'G. The lowest energy path for formation of the 3'G monoadduct with an activation barrier of 19.3 kcal/mol is connected with a disruption of the 'DNA-like' structure of ds(pGpG). Monoadduct formation is the rate-determining process. The second step, chelate formation, is kinetically preferred in the 3' → 5' direction. The whole process of the platination is exergonic by up to -18.8 kcal/mol. Structural changes of ds(pGpG), charge transfer effects, and the influence of platination on the G·C base pair interaction strengths are also discussed in detail.

Cisplatin interaction with amino acids cysteine and methionine from gas phase to solutions with constant pH

This work is focused on the computational studies of reactions of hydrated forms of cisplatin with sulphur-containing amino acids cysteine and methionine. First, the appropriate model for solvation of the examined complexes was searched for. The suggested procedure employs the B3LYP density functional, 6-311++G(2df,2pd) basis set with Stuttgart-Dresden pseudopotentials on heavy atoms, the D-PCM solvation model and the UAKS cavity which uses more realistic NPA partial charges instead of formal partial charges for platinum ligands. In the second part this model is applied to the evaluation of the Legendre transformed reaction Gibbs free energy of cisplatin with cysteine and methionine in solution at constant pH.

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.

Comparison of the electronic properties, and thermodynamic and kinetic parameters of the aquation of selected platinum(II) derivatives with their anticancer IC50 indexes

Three potential anticancer agents {trans-[PtCl(2)(NH(3))(thiazole)], cis-[PtCl(2)(NH(3))(piperidine)], and PtCl(2)(NH(3))(cyclohexylamine) (JM118)} were explored and compared with cisplatin and the inactive [PtCl(dien)](+) complex. Basic electronic properties, bonding and stabilization energies were determined, and thermodynamic and kinetic parameters for the aquation reaction were estimated at the B3LYP/6-311++G(2df,2pd) level of theory. Since the aquation process represents activation of these agents, the obtained rate constants were compared with the experimental IC(50) values for several tumor cells. Despite the fact that the processes in which these drugs are involved and the way in which they affect cells are very complex, some correlations can be deduced.

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.

Multiple Proton-Transfer Reactions in DNA Base Pairs by Coordination of Pt Complex

The possibility of multiple proton-transfer reactions in DNA base pairs because of coordination of cisplatin is theoretically elucidated by density functional theory (DFT) and by quantum mechanics/molecular mechanics (QM/MM) methods with an ONIOM method. From the energetics of two base pairs with the cisplatin, it is theoretically confirmed that the Pt complex is likely to bind in the order cis-(CG)-Pt-(GC), cis-(CG)-Pt-(AT), cis-(TA)-Pt-(AT), where G, C, A, and T are guanine, cytosine, adenine, and thymine, respectively, and the Pt atom bonds to the N7 site of G and A. This result supports the experimental evidence, where the structure cis-A-Pt-A is seldom observed at room temperature. The single proton-transfer reaction occurs in one of the two GC pairs. No simultaneous single proton-transfer reaction can occur in both base pairs. Two different single proton-transferred structures (cis-(CG*)(d)-Pt-(GC)(p) and cis-(CG)(d)-Pt-(G*C)(p), where the asterisk means a proton donor of G) are as stable as the original structures (CG)(d)-Pt-(GC)(p). The same tendency was observed with cis-(CG*)-Pt-(AT). In contrast to cisplatin, multiple single proton-transfer reactions occur in the system consisting of two base pairs with transplatin. The optimized structure agrees with the experimental data for Pt-G coordination except for the hydrogen-bonding length.

Pt-bridges in various single-strand and double-helix DNA sequences. DFT and MP2 study of the cisplatin coordination with guanine, adenine, and cytosine

In this study, various platinum cross-links in DNA bases were explored. Some of these structures occur in many cis/trans-platinated double-helixes or single-stranded adducts. However, in the models studied, no steric hindrance from sugar-phosphate backbone or other surroundings is considered. Such restrictions can change the bonding picture partially but hopefully the basic energy characteristics will not be changed substantially. The optimization of the structures explored was performed at the DFT level with the B3LYP functional and the 6-31G(d) basis set. Perturbation theory at the MP2/6-31++G(2df,2pd) level was used for the single-point energy and 6-31+G(d) basis set for the electron-property analyses. It was found that the most stable structures are the diguanine complexes followed by guanine-cytosine Pt-cross-links, ca 5 kcal mol(-1) less stable. The adenine-containing complexes are about 15 kcal mol(-1) below the stability of diguanine structures. This stability order was also confirmed by the BE of Pt-N bonds. For a detailed view on dative and electrostatic contributions to Pt-N bonds, Natural Population Analysis, determination of electrostatic potentials, and canonical Molecular Orbitals description of the examined systems were used.

Loss of amine from platinum(II) complexes: implications for cisplatin inactivation, storage, and resistance

Potential consequences of the binding of the anticancer drug cisplatin to various biomolecules in the cell have been investigated by using a combined density functional theory and continuum dielectric model approach. Since the amine ligands remain coordinated at the metal upon formation of the most frequent DNA adducts, whereas they were found to be displaced from the metal upon formation of drug metabolites, we have analyzed the factors governing amine loss from platinum(II) complexes as a possible pathway of cisplatin inactivation. The calculations systematically show the effect of 1) the trans ligand, 2) the charge of complex, 3) the nucleophile, and 4) the environment on the thermodynamic instability and kinetic lability of the platinum-amine bonds. After initial binding of cisplatin hydrolysis products to thioethers or thiols, loss of the amine trans to this sulfur ligand rather than replacement of the sulfur ligand itself by other nucleophiles like guanine-N7 is predicted to be the predominant reaction. The results of this study contribute to an understanding of the modes of cisplatin inactivation prior to DNA binding, for example, by elevated glutathione levels in cisplatin-resistant cancer cells.

Supramolecular Switches Based on the Guanine–Cytosine (GC) Watson–Crick Pair: Effect of Neutral and Ionic Substituents

We have theoretically analyzed Watson-Crick guanine-cytosine (GC) base pairs in which purine-C8 and/or pyrimidine-C6 positions carry a substituent X = NH(-), NH(2), NH(3) (+) (N series), O(-), OH, or OH(2) (+) (O series), using the generalized gradient approximation (GGA) of density functional theory at the BP86/TZ2P level. The purpose is to study the effects on structure and hydrogen-bond strength if X= H is substituted by an anionic, neutral, or cationic substituent. We found that replacing X = H by a neutral substituent has relatively small effects. Introducing a charged substituent, on the other hand, led to substantial and characteristic changes in hydrogen-bond lengths, strengths, and hydrogen-bonding mechanism. In general, introducing an anionic substituent reduces the hydrogen-bond-donating and increases the hydrogen-bond-accepting capabilities of a DNA base, and vice versa for a cationic substituent. Thus, along both the N and O series of substituents, the geometric shape and bond strength of our DNA base pair can be chemically switched between three states, thus yielding a chemically controlled supramolecular switch. Interestingly, the orbital-interaction component in some of these hydrogen bonds was found to contribute to more than 49 % of the attractive interactions and is thus virtually equal in magnitude to the electrostatic component, which provides the other (somewhat less than) 51 % of the attraction.

Cisplatin interaction with cysteine and methionine, a theoretical DFT study

Interactions of hydrated cisplatin complexes with sulphur-containing amino acids cysteine and methionine were explored. The square-planar cis-[Pt(NH3)2(H2O)X]+ complexes (where X=Cl- and OH-) were chosen as mono- and dihydrated reactants. Calculations using density functional theory (DFT) techniques with B3LYP functional were performed. The isolated molecules and the supermolecular approaches were employed for the determination of the reaction energies. Bond dissociation energies (BDE) were estimated in the model of isolated molecules and supermolecules were used for the determination of the association energies between the two interacting parts. Formation of monodentate complexes by replacing the aqua-ligand with the S, N, and O-sites of both amino acids represents an exothermic process. The highest BDE was found in cysteine structures for the Pt-S coordination. The bonding energy is about 114 kcal/mol, which is comparable with cisplatin-guanine adducts. Analogous BDE for methionine complexes is smaller by about 40 kcal/mol. This correlates well with the known fact that cysteine forms irreversible cisplatin adducts while similar adducts in the methionine case are reversible. The formation of chelate structures is an exothermic reaction only for the hydroxo-form of reactants in the supermolecular approach where additional association interactions between the released water and chelate molecules sufficiently stabilize the final product.

Vibrational Spectroscopy of the G· · ·C Base Pair: Experiment, Harmonic and Anharmonic Calculations, and the Nature of the Anharmonic Couplings

The results of harmonic and anharmonic frequency calculations on a guanine-cytosine complex with an enolic structure (a tautomeric form with cytosine in the enol form and with a hydrogen at the 7-position on guanine) are presented and compared to gas-phase IR-UV double resonance spectral data. Harmonic frequencies were obtained at the RI-MP2/cc-pVDZ, RI-MP2/TZVPP, and semiempirical PM3 levels of electronic structure theory. Anharmonic frequencies were obtained by the CC-VSCF method with improved PM3 potential surfaces; the improved PM3 potential surfaces are obtained from standard PM3 theory by coordinate scaling such that the improved PM3 harmonic frequencies are the same as those computed at the RI-MP2/cc-pVDZ level. Comparison of the data with experimental results indicates that the average absolute percentage deviation for the methods is 2.6% for harmonic RI-MP2/cc-pVDZ (3.0% with the inclusion of a 0.956 scaling factor that compensates for anharmonicity), 2.5% for harmonic RI-MP2/TZVPP (2.9% with a 0.956 anharmonicity factor included), and 2.3% for adapted PM3 CC-VSCF; the empirical scaling factor for the ab initio harmonic calculations improves the stretching frequencies but decreases the accuracy of the other mode frequencies. The agreement with experiment supports the adequacy of the improved PM3 potentials for describing the anharmonic force field of the G...C base pair in the spectroscopically probed region. These results may be useful for the prediction of the pathways of vibrational energy flow upon excitation of this system. The anharmonic calculations indicate that anharmonicity along single mode coordinates can be significant for simple stretching modes. For several other cases, coupling between different vibrational modes provides the main contribution to anharmonicity. Examples of strongly anharmonically coupled modes are the symmetric stretch and group torsion of the hydrogen-bonded NH2 group on guanine, the OH stretch and torsion of the enol group on cytosine, and the NH stretch and NH out-of-plane bend of the non-hydrogen-bonded NH group on guanine.

Hydration process as an activation of trans- and cisplatin complexes in anticancer treatment. DFT andab initio computational study of thermodynamic and kinetic parameters

The thermodynamic and kinetic aspects of hydration reactions of cis-/transplatin were explored. The polarizable continuum model was used for estimation of solvent effects. Using the B3LYP/6-31+G(d) method, the structures were optimized and vibrational frequencies estimated. Interaction energies and activation barriers were determined at the CCSD(T)/6-31++G(d,p) level within the COSMO approach. An associative mechanism was assumed with a trigonal-bipyramidal structure of the transition state. Within the applied model, all the hydration reactions are slightly endothermic. The Gibbs energies of cisplatin hydration amount to 7.0 and 14.2 kcal/mol for the chloride and ammonium replacement, respectively. Analogous values for the transplatin reactions are 6.8 and 11.9 kcal/mol. The determined rate constants are by several (three to four) orders of magnitude larger for the dechlorination process than for deammination. The cisplatin dechlorination rate constant was established as 1.3 x 10(-4) s(-1) in excellent accord with the experiment.

Hydrogen Bonding and Covalent Effects in Binding of Cisplatin to Purine Bases: Ab Initio and Atoms in Molecules Studies

Ab initio and density functional calculations are employed to investigate the role of hydrogen bonding in the binding of cisplatin to the purine bases guanine and adenine. Through the use of the theory of atoms in molecules (AIM), it is shown that hydrogen bonds are ubiquitous in such systems, with N-H...N and N-H...Cl interactions present in addition to the expected N-H...O. This in turn means that the known stability of cisplatin-guanine complexes cannot be ascribed solely to hydrogen bonding and allows decomposition of total binding energy into contributions from covalent and hydrogen bonds. To do so, a new method for predicting hydrogen bond energies from bond critical point properties is proposed, employing partial least-squares analysis to remove the family dependence of simple models. Still more hydrogen bond motifs are found in bifunctional complexes of the general type purine-[Pt(NH(3))(2)](2+)-purine, including purine...purine contacts, though again the energetics of these are insufficient to explain the observed trends in stability. Finally, the effect of platination on the pairing of guanine with cytosine is studied in a similar manner, revealing large redistributions of hydrogen bonding but surprisingly small overall changes in pairing energy.

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.

Hydrogen bonding, solvation, and hydrolysis of cisplatin: A theoretical study

Density functional calculations on a range of hydrogen bonded clusters of cisplatin are reported. A systematic search of 1:1 cisplatin:water complexes reveals only three stable minima, which contain a number of common, recurring interactions, such as an N-H...O-H...Cl bridging mode. Expanding these clusters by adding water molecules leads to a model of the first solvation shell of cisplatin, which contains the above motifs along with several strong water-water interactions. The strengths of such interactions are rationalized on the basis of electrostatic potentials, and quantified by use of Atoms in Molecules properties. This analysis also allows us to estimate cisplatin's position on Abraham's hydrogen bond acidity and basicity scales, indicating that cisplatin is a strong donor and acceptor of hydrogen bonds due to the dominance of hard, electrostatic interactions. The effects of this explicit solvation on the barrier to hydrolysis, and hence activation, of cisplatin are explored, indicating a slightly higher barrier than in the gas phase, leading to better agreement with experiment than either gas phase or continuum solvation calculations.

Activation barriers and rate constants for hydration of platinum and palladium square-planar complexes: An ab initio study

In the present work, an ab initio study on hydration (a metal-ligand replacement by water molecule or OH- group) of cis- and transplatin and their palladium analogs was performed within a neutral pseudomolecule approach (e.g., metal-complex+water as reactant complex). Subsequent replacement of the second ligand was considered. Optimizations were performed at the MP2/6-31+G(d) level with single-point energy evaluation using the CCSD(T)/6-31++G(d,p) approach. For the obtained structures of reactants, transition states (TS's), and products, both thermodynamic (reaction energies and Gibbs energies) and kinetic (rate constants) characteristics were estimated. It was found that all the hydration processes are mildly endothermic reactions-in the first step they require 8.7 and 10.2 kcal/mol for ammonium and chloride replacement in cisplatin and 13.8 and 17.8 kcal/mol in the transplatin case, respectively. Corresponding energies for cispalladium amount to 5.2 and 9.8 kcal/mol, and 11.0 and 17.7 kcal/mol for transpalladium. Based on vibrational analyses at MP2/6-31+G(d) level, transition state theory rate constants were computed for all the hydration reactions. A qualitative agreement between the predicted and known experimental data was achieved. It was also found that the close similarities in reaction thermodynamics of both Pd(II) and Pt(II) complexes (average difference for all the hydration reactions are approximately 1.8 kcal/mol) do not correspond to the TS characteristics. The TS energies for examined Pd(II) complexes are about 9.7 kcal/mol lower in comparison with the Pt analogs. This leads to 10(6) times faster reaction course in the Pd cases. This is by 1 or 2 orders of magnitude more than the results based on experimental measurements.

cis-{Pt(NH3)2(L)}2+/+ (L = Cl, H2O, NH3) Binding to Purines and CO: Does π-Back-Donation Play a Role?

The ability of cis-[Pt(NH(3))(2)(L)](2+/+), a molecular fragment of the anticancer drug cisplatin, to bind to purines and CO by pi-back-donation from Pt to the ligand was examined computationally. Optimized geometries and computed vibrational frequencies suggest that cis-[Pt(NH(3))(2)(L)](2+/+) (L = Cl, H(2)O, NH(3)) is a poor pi-donor and that pi-back-donation does not play an important role for Pt(II)-ligand interactions in general.

Theoretical Study of Cisplatin Binding to Purine Bases: Why Does Cisplatin Prefer Guanine over Adenine?

The thermodynamics and kinetics for the monofunctional binding of the antitumor drug cisplatin, cis-diamminedichloroplatinum(II), to a purine base site of DNA were studied computationally using guanine and adenine as model reactants. A dominating preference for initial attack at the N7-position of guanine is established experimentally, which is a crucial first step for the formation of a 1,2-intrastrand cross-link of adjacent guanine bases that leads to bending and unwinding of DNA. These structural distortions are proposed ultimately to be responsible for the anticancer activity of cisplatin. Utilizing density functional theory in combination with a continuum solvation model, we developed a concept for the initial Pt-N7 bond formation to atomic detail. In good agreement with experiments that suggested DeltaG++ = approximately 23 kcal/mol for the monofunctional platination of guanine, our model gives DeltaG++ = 24.6 kcal/mol for guanine, whereas 30.2 kcal/mol is computed when adenine is used. This result predicts that guanine is 3-4 orders of magnitude more reactive toward cisplatin than adenine. A detailed energy decomposition and molecular orbital analysis was conducted to explain the different barrier heights. Two effects are equally important to give the preference for guanine over adenine: First, the transition state is characterized by a strong hydrogen bond between the ammine-hydrogen of cisplatin and the O=C6 moiety of guanine in addition to a stronger electrostatic interaction between the two reacting fragments. When adenine binds, only a weak hydrogen bond forms between the chloride ligand of cisplatin and the H(2)N-C6 group of adenine. Second, a significantly stronger molecular orbital interaction is identified for guanine compared to adenine. A detailed MO analysis is presented to provide an intuitive view into the different electronic features governing the character of the Pt-N7 bond in platinated purine bases.

Molecular Interactions of Nucleic Acid Bases. A Review of Quantum-Chemical Studies

Ab initio quantum-chemical calculations with inclusion of electron correlation significantly contributed to our understanding of molecular interactions of DNA and RNA bases. Some of the most important findings are introduced in the present overview: structures and energies of hydrogen bonded base pairs, nature of base stacking, interactions between metal cations and nucleobases, nonplanarity of isolated nucleobases and other monomer properties, tautomeric equilibria of nucleobases, out-of-plane hydrogen bonds and amino acceptor interactions. The role of selected molecular interactions in nucleic acids is discussed and representative examples where these interactions occur are given. Also, accuracy of density functional theory, semiempirical methods, distributed multipole analysis and empirical potentials is commented on. Special attention is given to our very recent reference calculations on base stacking and H-bonding. Finally, we briefly comment on the relationship between advanced ab initio quantum-chemical methods and large-scale explicit solvent molecular dynamics simulations of nucleic acids.