Ab initio MP2 and DFT calculations of geometry and solution tautomerism of purine and some purine derivatives

Intramolecular Charge Transfer in the Azathioprine Prodrug Quenches Intersystem Crossing to the Reactive Triplet State in 6-Mercaptopurine†

The thiopurine prodrugs 6-mercaptopurine and azathioprine are among the world's essential medications for acute lymphoblastic leukemia, immunosuppression and several autoimmune conditions. Thiopurine prodrugs are efficient UVA absorbers and singlet oxygen generators and the long-term treatment with these prodrugs correlates with a high incidence of sunlight-induced skin cancer in patients. In this contribution, we show that the electronic relaxation mechanisms and photochemical properties of azathioprine are remarkably different from those of 6-mercaptopurine upon absorption of UVA radiation. UVA excitation of 6-mercaptopurine results in nearly 100% triplet yield and up to 30% singlet oxygen generation, whereas excitation of azathioprine with UVA leads to triplet yields of 15-3% depending on pH of the aqueous solution and <1% singlet oxygen generation. While photoexcitation of 6-mercaptopurine and other thiopurine prodrugs can facilitate oxidatively generated cell damage, azathioprine's poor photosensitization ability reveals the use of interchromophoric charge-transfer interactions for the molecular design of photostable prodrugs exhibiting a remarkable reduction in photocytotoxic side effects before drug metabolization.

Solvent Effect on the Stability and Reverse Substituent Effect in Nitropurine Tautomers

The solvent effect on the stability and electron-accepting properties (EA) of the nitro group attached to the C2, C6, or C8 position of nitropurine NH tautomers is investigated. For this purpose, the density functional theory (DFT) and the polarizable continuum model (PCM) of solvation in a wide range of solvents (1 < ε < 109) are used. We show that the EA properties of the NO2 group, described by the charge of the substituent active region (cSAR) model, are linearly dependent on the reciprocal of the solvent dielectric constant; in all cases, solvation enhances the EA properties of this group. Furthermore, the sensitivity of EA properties of the nitro group to the solvent effect depends on the proximity effects. It has been shown that the proximity of two endocyclic N atoms (two repulsive interactions) results in higher sensitivity than the asymmetric proximity of the endocyclic N atom and NH group (one repulsive and one attractive interaction). To explain this phenomenon, the geometry of the nitro group in coplanar form and after forcing its rotation around the CN bond is discussed. Relative stabilities of nitropurine tautomers in different solvents are also presented. Differences in the stabilities and solvation energies are explained by aromaticity, electronic structure, and intramolecular interactions of the nitropurine tautomers.

Structural and Electronic Characterization of a Photoresponsive Lanthanum(III) Complex Incorporated into Electrospun Fibers for Phosphate Ester Catalysis

Herein, we present the light-induced synthesis and characterization of a La³⁺/spiropyran derivative complex (LaMC) and its application as a catalyst when incorporated into electrospun polycaprolactone (PCL) fibers. In addition to experimental methods, computational calculations were also essential to better understand the structure and electronic characteristics of LaMC. The LaMC complex was identified as a 10-coordinated structure with the La³⁺ ion coordinated by four oxygens from the phenolate and the carbonyl of the carboxyl acid group from both MC ligands and by six oxygens from three nitrate ligands. In addition, LaMC was capable of getting reversibly isomerized by UV or visible light cycling. All PCL fibers were successively obtained, and their morphologies, surface properties, and catalytic behavior were studied. Results showed that PCL/LaMC fibers were capable of catalyzing bis(2,4-dinitrophenyl)phosphate degradation efficiently. Complete hydrolysis was accomplished in only 1.5 days relative to the half-life time of 35 days for the uncatalyzed hydrolysis at pH 8.1 and 25 °C.

Purine tautomeric preferences and bond-length alternation in relation with protonation-deprotonation and alkali metal cationization

Quantum chemical calculations were carried out for deprotonated (P-) and protonated purine (PH+) and for adducts with one alkali metal cation (P-M+ and PM+, where M+ is Li+ or Na+) in the gas phase {B3LYP/6-311+G(d,p)}, a model of perfectly apolar environment, and for selected structures in aqueous solution {PCM(water)//B3LYP/6-311+G(d,p)}, a reference polar medium for biological studies. All potential isomers of purine derivatives were considered, the favored structures indicated, and the preferred sites for protonation/deprotonation and cationization reactions determined. Proton and metal cation basicities of purine in the gas phase were discussed and compared with those of imidazole and pyrimidine. Bond-length alternations in the P, PH+, P-M+, and PM+ forms were quantitatively measured using the harmonic oscillator model of electron delocalization (HOMED) indices and compared with those for P. Variations of the HOMED values when proceeding from the purine structural building blocks, pyrimidine and imidazole, to the bicyclic purine system were also examined. Generally, the isolated NH isomers exhibit a strongly delocalized π-system (HOMED > 0.8). Deprotonation slightly increases the HOMED values, whereas protonation and cationization change the HOMED indices in different way. For bidentate M+-adducts, the HOMED values are larger than 0.9 like for the largely delocalized P-. The HOMED values correlate well in a comprehensive relationship with the relative Gibbs energies (ΔG) calculated for individual isomers whatever the purine form is, neutral, protonated, or cationized. When PCM-DFT model was utilized for P-, PH+, PM+, and P-M+ (M+ = Li+) both electron delocalization and relative stability are different from those for the molecules in vacuo. The solvation effects cause a slight increase in HOMEDs, whereas the ΔEs decrease, but in different ways. Hence, contribution of particular isomers in the isomeric mixtures of PH+, PM+, and P-M+ also varies. HOMED variations for the favored neutral, deprotonated, protonated, and lithiated forms of purine in the gas phase and aqueous solution.

Direct observation of an intramolecular charge transfer state in epigenetic nucleobase N6-methyladenine

N6-Methyladenine (6MeAde), the most abundant internal modification in mRNA, has proved to be an important epigenetic biomarker for gene regulation just like 5-methylcytosine in DNA. Recently, a unique UV-induced response of 6MeAde was reported, which makes it instructive and intriguing to reveal the excited state relaxation mechanism in this methylated adenine and its derivatives. In this work, we investigated 6MeAde and its ribose species N6-methyladenosine (6MeAdo) by using femtosecond time-resolved fluorescence up-conversion (FUC) and broadband transient absorption (TA) spectroscopy. Both 6MeAde and 6MeAdo exhibit a hundreds of femtoseconds lifetime, which originates from the efficient depletion of the ππ* (La) state. A several picoseconds lifetime is also observed and it should be attributed to the ππ* (Lb) state. Surprisingly, dual peak fluorescence emission is observed in 6MeAde and the long wavelength emission is ascribed to an intramolecular charge transfer (ICT) state. The lifetime of this ICT state is determined to be 107 ps. The kinetic isotope effect shows that the ICT state is closely associated with the solute-solvent H-bonding in aqueous solution. In 6MeAdo, the ICT state is apparently quenched and adenine-like excited state dynamics suggests that DNA/RNA containing such modification could still possess excellent photostability under UV irradiation. Our results contain an important insight for understanding excited state properties in epigenetic modified DNA/RNA.

Physico-chemical profiles of the wobble ↔ Watson-Crick G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) tautomerisations: a QM/QTAIM comprehensive survey

This study is intended to clarify in detail the tautomeric transformations of the wobble (w) G*·2AP(w) and A·2AP(w) nucleobase mispairs involving 2-aminopurine (2AP) into the Watson-Crick (WC) G·2AP(WC) and A*·2AP(WC) base mispairs (asterisks denote mutagenic tautomers of the DNA bases), respectively, by quantum-mechanical methods and Bader's Quantum Theory of Atoms in Molecules. Our previously reported methodology has been used, which allows the evolution of the physico-chemical parameters to be tracked along the entire internal reaction coordinate (IRC), not exclusively in the stationary states of these reactions. These biologically important G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) w ↔ WC tautomerisations, which are involved in mutagenic tautomerically-conformational pathways, determine the origin of the transitions and transversions induced by 2AP. In addition, it is established that they proceed through planar, highly stable, zwitterionic transition states and they exhibit similar physico-chemical profiles and stages of sequential intrapair proton transfer, followed by spatial rearrangement of the nucleobases relative to each other within the base pairs. These w ↔ WC tautomerisations occur non-dissociatively and are accompanied by a significant alteration in geometry (from wobble to Watson-Crick and vice versa) and redistribution of the specific intermolecular interactions, which can be divided into 10 patterns including AHB H-bonds and loosened A-H-B covalent bridges along the IRC of tautomerisation. Based on the redistribution of the geometrical and electron-topological parameters of the intrapair hydrogen bonds, exactly 9 key points have been allocated to characterize the evolution of these reactions.

Efficient intersystem crossing in 2-aminopurine riboside probed by femtosecond time-resolved transient vibrational absorption spectroscopy

The photophysical dynamics of 2-aminopurine, a fluorescent analogue of the canonical nucleobase adenine, has been studied by femtosecond transient vibrational absorption spectroscopy.

Whether the amino–imino tautomerism of 2-aminopurine is involved into its mutagenicity? Results of a thorough QM investigation

2AP* mutagenic tautomer is able to induce only one incorporation error – transversion – by pairing through the H-bonds into the G·2AP* mispair.

Structural grounds for the 2-aminopurine mutagenicity: a novel insight into the old problem of the replication errors

Mutagenic pressure of the 2AP molecule on DNA during its replication is realized via the more intensive generation of the T* mutagenic tautomers through the reaction 2AP·T(WC) → 2AP·T*(w).

Electronic and structural elements that regulate the excited-state dynamics in purine nucleobase derivatives

The excited-state dynamics of the purine free base and 9-methylpurine are investigated using experimental and theoretical methods. Femtosecond broadband transient absorption experiments reveal that excitation of these purine derivatives in aqueous solution at 266 nm results primarily in ultrafast conversion of the S2(ππ*) state to the vibrationally excited (1)nπ* state. Following vibrational and conformational relaxation, the (1)nπ* state acts as a doorway state in the efficient population of the triplet manifold with an intersystem crossing lifetime of hundreds of picoseconds. Experiments show an almost 2-fold increase in the intersystem crossing rate on going from polar aprotic to nonpolar solvents, suggesting that a solvent-dependent energy barrier must be surmounted to access the singlet-to-triplet crossing region. Ab initio static and surface-hopping dynamics simulations lend strong support to the proposed relaxation mechanism. Collectively, the experimental and computational results demonstrate that the accessibility of the nπ* states and the topology of the potential energy surfaces in the vicinity of conical intersections are key elements in controlling the excited-state dynamics of the purine derivatives. From a structural perspective, it is shown that the purine chromophore is not responsible for the ultrafast internal conversion in the adenine and guanine monomers. Instead, C6 functionalization plays an important role in regulating the rates of radiative and nonradiative relaxation. C6 functionalization inhibits access to the (1)nπ* state while simultaneously facilitating access to the (1)ππ*(La)/S0 conical intersection, such that population of the (1)nπ* state cannot compete with the relaxation pathways to the ground state involving ring puckering at the C2 position.

Geometric and energetic consequences of prototropy for adenine and its structural models – a review

Prototropy for adenine and its convenient models causes parallel changes of geometric (HOMED) and energetic (ΔE) parameters for neutral tautomers.

Quantum chemical and spectroscopic investigations of 3-methyladenine

FTIR, FT-Raman and UV-Vis spectra of 3-methyladenine have been recorded and investigated using quantum chemical calculations. The molecular geometry and vibrational spectra of 3-methyladenine in the ground state are computed by using HF and DFT methods with 6-311G(d,p) basis set. VSCF, CC-VSCF methods based on 2MR-QFF and PT2 (Barone method) have been utilized for computing anharmonic vibrational frequencies. These methods yield results that are in remarkable agreement with the experimental data. The magnitudes of coupling between pair of modes have been also computed. Vibrational modes are assigned with the help of visual inspection of atomic displacements. The electronic spectra, simulated at TD-B3LYP/6-311++G(d,p) level of theory, are compared to the experiment. The global quantities: electronic chemical potential, electrophilicity index, chemical hardness and softness based on HOMO and LUMO energy eigenvalues are also computed at B3LYP/6-311++G(d,p) level of theory.

Variations of the tautomeric preferences and π-electron delocalization for the neutral and redox forms of purine when proceeding from the gas phase (DFT) to water (PCM)

Quantum-chemical calculations were performed for all possible nine neutral tautomers of purine and their oxidized and reduced forms in water {PCM//DFT(B3LYP)/6-311+G(d,p)} and compared to those in the gas phase {DFT(B3LYP)/6-311+G(d,p)}. PCM hydration influences geometries, π-electron delocalization, and relative energies of purine tautomers in different ways. Generally, the harmonic oscillator model of electron delocalization (HOMED) indices increase when proceeding from the gas phase to aequeous solution for the neutral and redox forms of purine. Their changes for the neutral and oxidized tautomers are almost parallel to the relative energies showing that aromaticity plays an important role in the tautomeric preferences. Tautomeric stabilities and tautomeric preferences vary when proceeding from the gas phase to water indicating additionally that intra- and intermolecular interactions affect tautomeric equilibria. The tautomeric mixture of neutral purine in the gas phase consists mainly of the N9H tautomer, whereas two tautomers (N9H and N7H) dominate in water. For oxidized purine, N9H is favored in the gas phase, whereas N1H in water. A gain of one electron dramatically changes the relative stabilities of the CH and NH tautomers that C6H and C8H dominate in the tautomeric mixture in the gas phase, whereas N3H in water. These variations show exceptional sensitivity of the tautomeric purine system on environment in the electron-transfer reactions.

Effect of the H-bonding on aromaticity of purine tautomers

Four tautomers of purine (1-H, 3-H, 7-H, and 9-H) and their equilibrium H-bonded complexes with F(-) and HF for acidic and basic centers, respectively, were optimized by means of the B3LYP/6-311++G(d,p) level of theory. Purine tautomer stability increases in the following series: 1-H < 3-H < 7-H < 9-H, consistent with increasing aromaticity. Furthermore, the presence of a hydrogen bond with HF does not change this order. For neutral H-bonded complexes, the strongest and the weakest intermolecular interactions occur (-14.12 and -10.49 kcal/mol) for less stable purine tautomers when the proton acceptor is located in the five- and six-membered rings, respectively. For 9-H and 7-H tautomers the order is reversed. The H-bond energy for the imidazole complex with HF amounts to -14.03 kcal/mol; hence, in the latter case, the fusion of imidazole to pyrimidine decreases its basicity. The ionic H-bonds of N(-)···HF type are stronger by ~10 kcal/mol than the neutral N···HF intermolecular interactions. The hydrogen bond N(-)···HF energies in pyrrole and imidazole are -32.28 and -30.03 kcal/mol, respectively, and are substantially stronger than those observed in purine complexes. The aromaticity of each individual ring and of the whole molecule for all tautomers in ionic complexes is very similar to that observed for the anion of purine. This is not the case for neutral complexes and purine as a reference. The N···HF bonds perturb much more the π-electron structure of five-membered rings than that of the six-membered ones. The H-bonding complexes for 7-H and 9-H tautomers are characterized by higher aromaticity and a much lower range of HOMA variability.

One-electron reduction of 2-aminopurine in the aqueous phase. A DFT and pulse radiolysis study

The electron affinity and the subsequent proton affinity of the electron adducts of 2-aminopurine (abbreviated 2AP) and adenine (A) are calculated with density functional theory (DFT). According to these calculations 2AP*- and A*- have similar thermochemical parameters leading to the conclusion that their reaction pathways should be close to analogous. Using the pulse radiolysis technique 2AP*- is formed by reaction with the hydrated electron (e(-)aq) and the resulting transient absorption spectrum is assigned to 2AP(NH)*. Additionally, it was found, employing the laser flash photolysis technique, that the excited singlet state of 2AP is incapable of oxidizing guanine in the aqueous phase. However, the one-electron oxidized 2AP (2AP*+) has sufficient energy to ionize guanine according to the DFT results in agreement with experimental data from the literature.

Time-resolved fluorescent properties of 8-vinyl-deoxyadenosine and 2-amino-deoxyribosylpurine exhibit different sensitivity to their opposite base in duplexes

8-Vinyl-deoxyadenosine (8VA) has been recently introduced as a fluorescent analogue of adenosine that is less perturbing and less quenched than the well-established 2-amino-deoxyribosylpurine (2AP) probe when inserted in oligonucleotides. To further validate 8VA as a fluorescent substitute of A, we compared the ability of 8VA and 2AP in sequences of the type d(CGT TTT XNX TTT TGC) (with N=8VA or 2AP and X=T and C) to discriminate the nature of the opposite base (Y) in duplexes. For both probes, systematic variations in the amplitudes of the short- and long-lived lifetimes of the fluorescence intensity decays as well as in the amplitude of the fast rotational correlation time of the fluorescence anisotropy decays were observed as a function of the nature of Y. From these parameters, we inferred a stability order 8VA-T > 8VA-G > 8VA-A > 8VA-C, similar to the stability order with the native A base, but different from the stability order with 2AP. Using a combination of molecular mechanics and ab initio calculations, we found that the time-resolved parameters of 8VA, but not the 2AP ones, correlate well with the geometry and the strength of the A-Y base-pairing interaction. This may be rationalized by the smaller structural and electronic perturbations induced by the vinyl group in position 8 as compared to the amino group at position 2. As a consequence, substitution of A by 8VA in a base pair was found to only minimally modify the structure and interaction energy of the base pair. Thus, 8VA can be used as a fluorescent substitute of the natural A, to straightforwardly discriminate the nature of the opposite base. This may find interesting applications notably in the elucidation of the mechanisms and dynamics of the DNA mismatch repair system.

The mid-IR spectra of 9-ethyl guanine, guanosine, and 2-deoxyguanosine

We present the mid-IR (400-1800 cm(-1)) spectra of 9-ethyl guanine, guanosine, and 2-deoxyguanosine measured by IR-UV double-resonance spectroscopy. We compare the recorded mid-IR spectra with the spectra of the most stable structures obtained from RI-MP2 and RI-DFT-D calculations. The results confirm the enol form for all structures and demonstrate the efficacy of a new approach to DFT calculations that includes dispersion interactions.

Three-state conical intersections in nucleic acid bases

The involvement of three-state conical intersections in the photophysics and radiationless decay processes of the nucleobases has been investigated using multireference configuration interaction methods. Three-state conical intersections have been located for the pyrimidine base, uracil, and the purine base, adenine. In uracil, a three-state degeneracy between the S(0), S(1), and S(2) states has been located at 6.2 eV above the ground-state minimum energy. This energy is 0.4 eV higher than vertical excitation to S(2) and at least 1.3 eV higher than the two-state conical intersections found previously. In adenine, two different three-state degeneracies between the S(1), S(2), and S(3) states have been located at energies close to the vertical excitation energies. The energetics of these three-state conical intersections suggest they can play a role in a radiationless decay pathway present in adenine. The existence of two different seams of three-state conical intersections indicates that these features are common and complicate the potential energy surfaces of adenine and possibly many other aromatic molecules.

Photophysics and X-ray Structure of Crystalline 2-Aminopurine

To explore the effect of intermolecular interactions on the photophysics of 2-aminopurine (2AP) in a well-defined environment, we have investigated the fluorescence properties of single 2AP crystals, having determined their X-ray structure. In the crystal, 2AP is subject to base-stacking and hydrogen-bonding interactions similar to those found in DNA. The crystal shows dual fluorescence: pi-stacked molecules in the bulk of the lattice have redshifted excitation and emission spectra, while molecules at defect sites have spectra similar to those of 2AP in solution or in DNA. Heterogeneous intermolecular interactions in the crystal give rise to multiexponential fluorescence decay characteristics similar to those observed for 2AP-labelled DNA. The presence of about 13 % of the 7H tautomer in the crystal confirms that 9H-7H tautomerisation of 2AP occurs in the ground state. Long-wavelength excitation of a 2AP-labelled oligonucleotide duplex produced redshifted emission similar to that observed in the crystal, indicating that pi-stacking interaction of 2AP with nucleobases gives rise to a low energy excited state.

Ab Initio Studies on the Radiationless Decay Mechanisms of the Lowest Excited Singlet States of 9H-Adenine

The mechanisms that are responsible for the rapid deactivation of the (1)npi and( 1)pipi excited singlet states of the 9H isomer of adenine have been investigated with multireference ab initio methods (complete-active-space self-consistent-field (CASSCF) method and second-order perturbation theory based on the CASSCF reference (CASPT2)). Two novel photochemical pathways, which lead to conical intersections of the S(1) excited potential-energy surface with the electronic ground-state surface, have been identified. They involve out-of-plane deformations of the six-membered aromatic ring via the twisting of the N(3)C(2) and N(1)C(6) bonds. These low-lying conical intersections are separated from the minimum energy of the lowest ((1)npi) excited state in the Franck-Condon region by very low energy barriers (of the order of 0.1 eV). These properties of the S(1) and S(0) potential-energy surfaces explain the unusual laser-induced fluorescence spectrum of jet-cooled 9H-adenine, showing sharp structures only in a narrow energy interval near the origin, as well as the extreme excess-energy dependence of the lifetime of the singlet excited states. It is suggested that internal-conversion processes via conical intersections, which are accessed by out-of-plane deformation of the six-membered ring, dominate the photophysics of the lowest vibronic levels of adenine in the gas phase, while hydrogen-abstraction photochemistry driven by repulsive (1)pisigma states may become competitive at higher excitation energies. These ultrafast excited-state deactivation processes provide adenine with a high degree of intrinsic photostability.

N-Glucosylation of Cytokinins by Glycosyltransferases of Arabidopsis thaliana

Cytokinins are plant hormones that can be glucosylated to form O-glucosides and N-glucosides. The glycoconjugates are inactive and are thought to play a role in homeostasis of the hormones. Although O-glucosyltransferases have been identified that recognize cytokinins, the enzymes involved in N-glucosylation have not been identified even though the process has been recognized for many years. This study utilizes a screening strategy in which 105 recombinant glycosyltransferases (UGTs) of Arabidopsis have been analyzed for catalytic activity toward the classical cytokinins: trans-zeatin, dihydrozeatin, N(6)-benzyladenine, N(6)-isopentenyladenine, and kinetin. Five UGTs were identified in the screen. UGT76C1 and UGT76C2 recognized all cytokinins and glucosylated the hormones at the N(7) and N(9) positions. UGT85A1, UGT73C5, and UGT73C1 recognized trans-zeatin and dihydrozeatin, which have an available hydroxyl group for glucosylation and formed the O-glucosides. The biochemical characteristics of the N-glucosyltransferases were analyzed, and highly effective inhibitors of their activities were identified. Constitutive overexpression of UGT76C1 in transgenic Arabidopsis confirmed that the recombinant enzyme functioned in vivo to glucosylate cytokinin applied to the plant. The role of the N-glucosyltransferases in cytokinin metabolism is discussed.