Direct determination of tautomerism in purine derivatives by low-temperature NMR spectroscopy

Improved synthesis of 6-bromo-7-[11C]methylpurine for clinical use

BACKGROUND

Multidrug resistance-associated protein 1 (MRP1), an energy-dependent efflux pump, is expressed widely in various tissues and contributes to many physiological and pathophysiological processes. 6-Bromo-7-[11C]methylpurine ([11C]7m6BP) is expected to be useful for the assessment of MRP1 activity in the human brain and lungs. However, the radiochemical yield (RCY) in the synthesis of [11C]7m6BP was low, limiting its clinical application, because the methylation of the precursor with [11C]CH3I provided primarily the undesired isomer, 6-bromo-9-[11C]methylpurine ([11C]9m6BP). To increase the RCY of [11C]7m6BP, we investigated conditions for improving the [11C]7m6BP/[11C]9m6BP selectivity of the methylation reaction.

RESULTS

[11C]7m6BP was manually synthesized via the methylation of 6-bromopurine with [11C]CH3I in various solvents and at different temperatures in the presence of potassium carbonate for 5 min. Several less polar solvents, including tetrahydrofuran (THF), 2-methyltetrahydrofuran (2-MeTHF), and ethyl acetate (AcOEt) improved the [11C]7m6BP/[11C]9m6BP selectivity from 1:1 to 2:1, compared with the conventionally used solvents for the alkylation of 6-halopurines, acetone, acetonitrile, and N,N-dimethylformamide. However, a higher temperature (140 °C or 180 °C) was needed to progress the 11C-methylation in the less polar solvents, and the manual conditions could not be directly translated to an automated synthesis. [11C]Methyl triflate ([11C]CH3OTf) was thus used as a methylating agent to increase the conversion at a lower temperature. The 11C-methylation using [11C]CH3OTf at 100 °C proceeded efficiently in THF, 2-MeTHF, and AcOEt with maintenance of the improved selectivity. Starting from 28 to 34 GBq [11C]CO2, [11C]7m6BP was produced with 2.3-2.6 GBq for THF, 2.7-3.3 GBq for AcOEt, and 2.8-3.9 GBq for 2-MeTHF at approximately 30 min after the end of bombardment (n = 3 per solvent). The isolated RCYs (decay corrected) for THF, 2-MeTHF, and AcOEt were 24-28%, 29-35%, and 22-31% (n = 3), respectively.

CONCLUSIONS

The use of THF, 2-MeTHF, and AcOEt improved the [11C]7m6BP/[11C]9m6BP selectivity in the methylation reaction, and the improved method provided [11C]7m6BP with sufficient radioactivity for clinical use.

Influence of the Solvent on the Stability of Aminopurine Tautomers and Properties of the Amino Group

Amino derivatives of purine (2-, 6-, 8-, and N-NH2) have found many applications in biochemistry. This paper presents the results of a systematic computational study of the substituent and solvent effects in these systems. The issues considered are the electron-donating properties of NH2, its geometry, π-electron delocalization in purine rings and tautomeric stability. Calculations were performed in ten environments, with 1 < ε < 109, using the polarizable continuum model of solvation. Electron-donating properties were quantitatively described by cSAR (charge of the substituent active region) parameter and π-electron delocalization by using the HOMA (harmonic oscillator model of aromaticity) index. In aminopurines, NH2 proximity interactions depend on its position and the tautomer. The results show that they are the main factor determining how solvation affects the electron-donating strength and geometry of NH2. Proximity with the NH∙∙∙HN repulsive interaction between the NH2 and endocyclic NH group results in stronger solvent effects than the proximity with two attractive NH∙∙∙N interactions. The effect of amino and nitro (previously studied) substitution on aromaticity was compared; these two groups have, in most cases, the opposite effect, with the largest being in N1H and N3H purine tautomers. The amino group has a smaller effect on the tautomeric preferences of purine than the nitro group. Only in 8-aminopurine do tautomeric preferences change: N7H is more stable than N9H in H2O.

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.

The thermochemical, structural, and spectroscopic analyses of the tautomers of sulfur and selenium modified emissive nucleobases

In this study, density functional theory (DFT) and time dependent density functional theory (TD-DFT) are used to investigate the stabilities and spectral properties [IR, UV–vis, and two-photon absorption (2PA)] of two sets of modified RNA nucleobase tautomers. The modifications introduce either a sulfur or selenium atom to form an isothiazolo[4,3-d]pyrimidine or isoselenazolo[4,3-d]pyrimidine heterocyclic core, respectively. The relative stabilities of both sets of modified tautomers determined with B3LYP/6-31++G(d,p) reveal that in water (with a polarizable continuum model), the 6-keto-2-amino tautomer of guanine and the rare 4-imino-2-keto tautomer of cytosine may be present at significant populations, whereas the 6-enol-2-amino tautomer of guanine is more common in the gas phase. The identification of these modified tautomers is possible due to the natural differences in their vibrational modes and hence IR spectra. Furthermore, the photophysical properties of both these sets of modified tautomers indicate that excitation and emission energies are shifted relative to their more abundant form in both one photon absorption and emission and 2PA spectra, as determined at the B3LYP/6-31++G(d,p) and CAM-B3LYP/aug-cc-pVDZ levels of theory, respectively. Even though the 2PA cross sections in water for all of the species are small (0.3–2.3 GM), the modified cytosine tautomer shows promise, as its cross section is larger than the more dominant form. The spectral separation between the dominant form and the tautomers of isoselenazole and isothiazole modified cytosine and guanine are relatively similar, suggesting both modifications could be useful in elucidating the tautomers from their more abundant counterparts.

Deconvolution of fast exchange equilibrium states in NMR spectroscopy using virtual reference standards and probability theory

A methodology for deconvolution of fast exchange equilibrium states in NMR spectroscopy (DFEQNMR) was developed based on DFT-GIAO NMR chemical shift prediction and a probability theory algorithm. Proof-of-concept studies were performed to estimate the protonation state of N-containing organic molecules involving fast proton exchange equilibrium and evaluate the solution tautomerism of a purine derivative. DFT-GIAO calculations were optimized to achieve good accuracy in ¹³C, ¹H and ¹⁵N chemical shift prediction for protonated species. The probability theory algorithm enabled the determination of solution species ratios and yielded 95% confidence regions by comparing experimental and simulated chemical shift data sets. The calculation showed good accuracy for model partial salts with various functionalities and application in structure elucidation of complex natural product partial salts was also demonstrated. This method showed promising potential in acquisition of important insight into fast exchange equilibrium systems with only one experimental NMR chemical shift data set.

NMR-Based Investigation of Hydrogen Bonding in a Dihydroanthracen-1(4 H)one from Rubia philippinensis and Its Soluble Epoxide Hydrolase Inhibitory Potential

Hydrogen bonding is a vital feature of a large ensemble of chemical structures. Soluble epoxide hydrolase (sEH) has been targeted for development of the treatment for inflammation-associated diseases. Compounds 1 and 2 were purified from Rubia philippinensis, and their structures were established via physical data analysis. Compound 1 possesses intramolecular hydrogen bonding, sufficiently robust to transfer heteronuclear magnetization via a nonbonded interaction. The bonding strength was assessed using the ¹H NMR chemical shift temperature coefficients (-1.8 ppb/K), and the heteronuclear coupling constants were measured. The stereochemical details were investigated using interproton distance analysis and ECD. Purified compounds displayed moderate sEH-inhibitory activity.

Tautomerism provides a molecular explanation for the mutagenic properties of the anti-HIV nucleoside 5-aza-5,6-dihydro-2'-deoxycytidine

Viral lethal mutagenesis is a strategy whereby the innate immune system or mutagenic pool nucleotides increase the error rate of viral replication above the error catastrophe limit. Lethal mutagenesis has been proposed as a mechanism for several antiviral compounds, including the drug candidate 5-aza-5,6-dihydro-2'-deoxycytidine (KP1212), which causes A-to-G and G-to-A mutations in the HIV genome, both in tissue culture and in HIV positive patients undergoing KP1212 monotherapy. This work explored the molecular mechanism(s) underlying the mutagenicity of KP1212, and specifically whether tautomerism, a previously proposed hypothesis, could explain the biological consequences of this nucleoside analog. Establishing tautomerism of nucleic acid bases under physiological conditions has been challenging because of the lack of sensitive methods. This study investigated tautomerism using an array of spectroscopic, theoretical, and chemical biology approaches. Variable temperature NMR and 2D infrared spectroscopic methods demonstrated that KP1212 existed as a broad ensemble of interconverting tautomers, among which enolic forms dominated. The mutagenic properties of KP1212 were determined empirically by in vitro and in vivo replication of a single-stranded vector containing a single KP1212. It was found that KP1212 paired with both A (10%) and G (90%), which is in accord with clinical observations. Moreover, this mutation frequency is sufficient for pushing a viral population over its error catastrophe limit, as observed before in cell culture studies. Finally, a model is proposed that correlates the mutagenicity of KP1212 with its tautomeric distribution in solution.

Stability of N-glycosidic bond of (5'S)-8,5'-cyclo-2'-deoxyguanosine

8,5'-Cyclopurine deoxynucleosides are unique tandem lesions containing an additional covalent bond between the base and the sugar. These mutagenic and genotoxic lesions are repaired only by nucleotide excision repair. The N-glycosidic (or C1'-N9) bond of 2'-deoxyguanosine (dG) derivatives is usually susceptible to acid hydrolysis, but even after cleavage of this bond of the cyclopurine lesions, the base would remain attached to the sugar. Here, the stability of the N-glycosidic bond and the products formed by formic acid hydrolysis of (5'S)-8,5'-cyclo-2'-deoxyguanosine (S-cdG) were investigated. For comparison, the stability of the N-glycosidic bond of 8,5'-cyclo-2',5'-dideoxyguanosine (ddcdG), 8-methyl-2'-deoxyguanosine (8-Me-dG), 7,8-dihydro-8-oxo-2'-deoxyguanosine (8-Oxo-dG), and dG was also studied. In various acid conditions, S-cdG and ddcdG exhibited similar stability to hydrolysis. Likewise, 8-Me-dG and dG showed comparable stability, but the half-lives of the cyclic dG lesions were at least 5-fold higher than those of dG or 8-Me-dG. NMR studies were carried out to investigate the products formed after the cleavage of the C1'-N9 bond. 2-Deoxyribose generated α and β anomers of deoxyribopyranose and deoxyribopyranose oligomers following acid treatment. S-cdG gave α- and β-deoxyribopyranose linked guanine as the major products, but α and β anomers of deoxyribofuranose linked guanine and other products were also detected. The N-glycosidic bond of 8-Oxo-dG was found exceptionally stable in acid. Computational studies determined that both the protonation of the N7 atom and the rate constant in the bond breaking step control the overall kinetics of hydrolysis, but both varied for the molecules studied indicating a delicate balance between the two steps. Nevertheless, the computational approach successfully predicted the trend observed experimentally. For 8-Oxo-dG, the low pK(a) of O(8) and N3 prevented appreciable protonation, making the free energy for N-glycosidic bond cleavage in the subsequent step very high.

1H-, 13C-, and 15N-NMR chemical shifts for selected glucosides and ribosides of aromatic cytokinins

The (1)H and (13)C NMR resonances of 16 purine glucosides were assigned by a combination of one- and two-dimensional NMR experiments, including gs-COSY, gs-HSQC, and gs-HMBC, in order to characterize the effect of substituent and the position of glucose unit on the NMR chemical shifts. In addition, (15)N NMR chemical shifts for selected derivatives were investigated by using (1)H-(15)N chemical shift correlation techniques. To map the influence of sugar moiety on the directly bonded nitrogen atom, selected N(9)-glucosides and their ribose analogs were compared. Characteristic long-range (1)H-(15)N coupling constants, measured by using (1)H-(15)N gradient-selected single-quantum multiple bond correlation (GSQMBC), are also reported and discussed. All compounds investigated here belong to cytokinins, an important group of plant hormones.

NMR study of 5-substituted pyrazolo[3,4-c]pyridine derivatives

Substituted pyrazolopyridines are potent inhibitors of phosphodiesterases and cyclin-dependent kinases. In this study, NMR was used to investigate the potential N1-H and N2-H tautomerism of 5-substituted pyrazolo[3,4-c]pyridine derivatives. Six compounds were fully characterized by using (1)H, (13)C, and (15)N chemical shifts and indirect (1)H--(13)C and (1)H--(15)N coupling constants. The (1)H NMR spectra were measured over a broad range of temperatures. All of the compounds were shown to exist predominantly in the N1-H tautomeric form. Complementary quantum-chemical calculations of the chemical shieldings and indirect spin-spin couplings support the structural conclusions drawn.

Adenine tautomers: relative stabilities, ionization energies, and mismatch with cytosine

In this study, we have investigated 12 tautomers of the DNA base adenine at the BP86/TZ2P and BP86/QZ4P levels of density functional theory. The vertical and adiabatic ionization energies of all tautomers were determined as the difference in energy between the radical cation and the corresponding neutral system. Furthermore, an evaluation is made for the eigenvalue spectra calculated with the SAOP functional, which is shown to lead to substantial improvements for orbital energies compared to BP86. We have also explored the correlations between the Kohn-Sham orbitals of the different tautomers at the BP86/QZ4P and SAOP/QZ4P levels. Finally, we discuss implications of the existence of the tautomeric forms of adenine for the DNA replication.