Pt(II) coordination to N1 of 9-methylguanine: why it facilitates binding of additional metal ions to the purine ring.

The preparation and X-ray crystal structure analysis of {trans-[Pt(MeNH(2))(2)(9-MeG-N1)(2)]}⋅{3 K(2)[Pt(CN)(4)]}⋅6 H(2)O (3 a) (with 9-MeG being the anion of 9-methylguanine, 9-MeGH) are reported. The title compound was obtained by treating [Pt(dien)(9-MeGH-N7)](2+) (1; dien=diethylenetriamine) with trans-[Pt(MeNH(2))(2)(H(2)O)(2)](2+) at pH 9.6, 60 °C, and subsequent removal of the [(dien)Pt(II)] entities by treatment with an excess amount of KCN, which converts the latter to [Pt(CN)(4)](2-). Cocrystallization of K(2)[Pt(CN)(4)] with trans-[Pt(MeNH(2))(2)(9-MeG-N1)(2)] is a consequence of the increase in basicity of the guanine ligand following its deprotonation and Pt coordination at N1. This increase in basicity is reflected in the pK(a) values of trans-[Pt(MeNH(2))(2)(9-MeGH-N1)(2)](2+) (4.4±0.1 and 3.3±0.4). The crystal structure of 3 a reveals rare (N7,O6 chelate) and unconventional (N2,C2,N3) binding patterns of K(+) to the guaninato ligands. DFT calculations confirm that K(+) binding to the sugar edge of guanine for a N1-platinated guanine anion is a realistic option, thus ruling against a simple packing effect in the solid-state structure of 3 a. The linkage isomer of 3 a, trans-[Pt(MeNH(2))(2)(9-MeG-N7)(2)] (6 a) has likewise been isolated, and its acid-base properties determined. Compound 6 a is more basic than 3 a by more than 4 log units. Binding of metal entities to the N7 positions of 9-MeG in 3 a has been studied in detail for [(NH(3))(3)Pt(II)], trans-[(NH(3))(2)Pt(II)], and [(en)Pd(II)] (en=ethylenediamine) by using (1)H NMR spectroscopy. Without exception, binding of the second metal takes place at N7, but formation of a molecular guanine square with trans-[(Me(2)NH(2))Pt(II)] cross-linking N1 positions and trans-[(NH(3))(2)Pt(II)] cross-linking N7 positions could not be confirmed unambiguously, despite the fact that calculations are fully consistent with its existence.

Differential stabilization of adenine quartets by anions and cations.

We have investigated the structures and stabilities of four different adenine quartets with alkali and halide ions in the gas phase and in water, using dispersion-corrected density functional theory at the BLYP-D/TZ2P level. First, we examine the empty quartets and how they interact with alkali cations and halide anions with formation of adenine quartet-ion complexes. Second, we examine the interaction in a stack, in which a planar adenine quartet interacts with a cation or anion in the periphery as well as in the center of the quartet. Interestingly, for the latter situation, we find that both cations and anions can stabilize a planar adenine quartet in a stack.

A ditopic ion-pair receptor based on stacked nucleobase quartets.

Pass the salt, please! State-of-the-art computations indicate that the stacking complex of a guanine quartet and an adenine quartet (G(4)A(4)) can function as a potent ditopic receptor for NaCl in aqueous solution (see picture; Na(+), Cl(-) yellow, O red, N blue, C black, H white).

Rare tautomers of 1-methyluracil and 1-methylthymine: tuning relative stabilities through coordination to PtII complexes.

We have computationally explored how the relative stabilities of 1-methyluracil (1-MeUH) tautomers can be tuned through coordination of these tautomers to Pt(II) complexes with a particular set of ligands. This has been done using density functional theory at the BP86/TZ2P level. Thus, we have examined the water/1-MeUH exchange reactions of [Pt(II)(A)(B)(C)(OH(2))](q) + 1-MeUH to uncover: i) which tautomers are best stabilized by the Pt(II) complex, and ii) how the net charge q in the complex affects the reaction energy. The net charge q depends on the ligands A, B, and C, which can be the neutral NH(3) or anionic Cl(-). To reveal the effect of solvation, all reaction systems are studied both in the gas phase and in water. Also the stabilization of tautomers of 1-methylthymine (1-MeTH) by cisplatin is investigated. The calculations reveal that relative energies of the metal (here: Pt(II))-complexed forms of the various tautomers (here: of 1-MeUH and 1-MeTH) do not parallel those of the free tautomers. Rather, a rare nucleobase tautomer, despite its low natural abundance, may become favored over the predominant one when complexed to a metal ion.

Pi-pi stacking tackled with density functional theory.

Through comparison with ab initio reference data, we have evaluated the performance of various density functionals for describing pi-pi interactions as a function of the geometry between two stacked benzenes or benzene analogs, between two stacked DNA bases, and between two stacked Watson-Crick pairs. Our main purpose is to find a robust and computationally efficient density functional to be used specifically and only for describing pi-pi stacking interactions in DNA and other biological molecules in the framework of our recently developed QM/QM approach "QUILD". In line with previous studies, most standard density functionals recover, at best, only part of the favorable stacking interactions. An exception is the new KT1 functional, which correctly yields bound pi-stacked structures. Surprisingly, a similarly good performance is achieved with the computationally very robust and efficient local density approximation (LDA). Furthermore, we show that classical electrostatic interactions determine the shape and depth of the pi-pi stacking potential energy surface.

Nanoswitches based on DNA base pairs: why adenine-thymine is less suitable than guanine-cytosine.

Substituted Watson-Crick guanine-cytosine (GC) base pairs were recently shown to yield robust three-state nanoswitches. Here, we address the question: Can such supramolecular switches also be based on Watson-Crick adenine-thymine (AT) base pairs? We have theoretically analyzed AT 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. Thus, we explore the trend in geometrical shape and hydrogen bond strengths in AT pairs along a series of stepwise protonations of the substituents. Introducing a charge on the substituents leads to substantial and characteristic changes in the individual hydrogen bond lengths when compared to the neutral AT pair. However, the trends along the series of negative, neutral, and positive substituents are less systematic and less pronounced than for GC. In certain instances, internal proton transfer from thymine to adenine occurs. Our results suggest that AT is a less suitable candidate than GC in the quest for chemically controlled nanoswitches.

Catalytic site prediction and virtual screening of cytochrome P450 2D6 substrates by consideration of water and rescoring in automated docking.

Automated docking strategies successfully applied to binding mode predictions of ligands in Cyt P450 crystal structures in an earlier study (de Graaf et al. J. Med. Chem. 2005, 7, 2308-2318) were used for the catalytic site prediction (CSP) of 65 substrates in a CYP2D6 homology model. The consideration of water molecules at predicted positions in the active site and the rescoring of pooled docking poses from four different docking programs (AutoDock, FlexX, GOLD-Goldscore, and GOLD-Chemscore) with the SCORE scoring function enabled the successful prediction of experimentally reported sites of catalysis of more than 80% of the substrates. Three docking algorithms (FlexX, GOLD-Goldscore, and GOLD-Chemscore) were subsequently used in combination with six scoring functions (Chemscore, DOCK, FlexX, GOLD, PMF, and SCORE) to assess the ability of docking-based virtual screening methods to prioritize known CYP2D6 substrates seeded into a drug-like chemical database (in the absence and presence of active-site water molecules). Finally, the optimal docking strategy in terms of virtual screening accuracy, GOLD-Chemscore with the consideration of active-site water (60% of known substrates recovered in the top 5% of the ranked drug-like database), was verified experimentally; it was successfully used to identify high-affinity CYP2D6 ligands among a larger proprietary database.

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.