Charge effects for differentiation of oligodeoxynucleotide isomers containing 8-oxo-dG residues

Dissociation reactions of a series of multiply charged oligodeoxynucleotide (ODN) 12-mer anions were studied using an ion trap mass spectrometer. These mixed nucleobase 12-mers fragment first by loss of a neutral nucleobase (A, G, C, and/or 5-methyl-cytosine) followed by cleavage at 3' C-O bond of the sugar from which the base is lost to produce the complementary sequence ions, i.e., [a - B] and w type of ions. No detectable loss of 8-oxo-guanine and/or thymine from these 12-mers is observed under gentle collision conditions in the ion trap. The primary loss of a nucleobase and the subsequent backbone cleavage to generate sequence ions strongly depend on the charge state of the parent molecular ion. For low charge states (2- and 3-), product ions due to the loss of a neutral guanine base and related sequence ions are dominant in the tandem mass spectra. However, preferential loss of a neutral adenine becomes the primary reaction channel from the 5- charge state of the molecular ion. Such charge state dependent fragmentation behavior was utilized to determine the site of 8-oxo-dG residue in a series of structural isomers. The position of 8-oxo-dG residue can be simply determined from the fragmentation pattern of 3- charge state, but not of 5- charge state. It is suggested that in addition to specific modification that affects the N-glycosidic bond strength, total charge content of an ODN is an important factor for determining the differential fragmentation behavior.     

Gas-phase studies of purine 3-methyladenine DNA glycosylase II (AlkA) substrates

3-Methyladenine DNA glycosylase II (AlkA) is an enzyme that cleaves a wide range of damaged bases from DNA. The gas-phase thermochemical properties (tautomerism, acidity, and proton affinity) have been measured and calculated for a series of AlkA purine substrates (7-methyladenine, 7-methylguanine, 3-methyladenine, 3-methylguanine, purine, 6-chloropurine, xanthine) that have not been heretofore examined. The damaged nucleobases are found to be more acidic than the normal nucleobases adenine and guanine. Because of this increased acidity, the damaged bases would be expected to be more easily cleaved from DNA by AlkA (their conjugate bases should be better leaving groups). We find that the gas-phase acidity correlates to the AlkA excision rates, which lends support to an AlkA mechanism wherein the enzyme provides a nonspecific active site, and nucleobase cleavage is dependent on the intrinsic N-glycosidic bond stability.     

Chemical basis of DNA sugar-phosphate cleavage by low-energy electrons

DNA damage by low-energy electrons (LEE) was examined using a novel system in which thin solid films of oligonucleotide tetramers (CGTA and GCAT) were irradiated with monoenergetic electrons (10 eV) under ultrahigh vacuum. The products of irradiation were examined by HPLC. These analyses permitted the quantitation of 16 nonmodified nucleobase, nucleoside, and nucleotide fragments of each tetramer resulting from the cleavage of phosphodiester and N-glycosidic bonds. The distribution of nonmodified products suggests a mechanism of damage involving initial electron attachment to nucleobase moieties, followed by electron transfer to the sugar-phosphate backbone, and subsequent dissociation of the phosphodiester bond. Moreover, virtually all the nonmodified fragments contained a terminal phosphate group at the site of cleavage. These results demonstrate that the phosphodiester bond breaks by a distinct pathway in which the negative charge localizes on the phosphodiester bond giving rise to nonmodified fragments with an intact phosphate group. Conversely, the radical must localize on the sugar moiety to give as yet unidentified modifications. In summary, the reaction of LEE with simple tetramers involved dissociative electron attachment leading to phosphodiester bond cleavage and the formation of nonmodified fragments.     

A kinetic and thermodynamic study of the glycosidic bond cleavage in deoxyuridine

Density functional theory was used to study the thermodynamics and kinetics for the glycosidic bond cleavage in deoxyuridine. Two reaction pathways were characterized for the unimolecular decomposition in vacuo. However, these processes are associated with large reaction barriers and highly endothermic reaction energies, which is in agreement with experiments that suggest a (water) nucleophile is required for the nonenzymatic glycosidic bond cleavage. Two (S(N)1 and S(N)2) reaction pathways were characterized for direct hydrolysis of the glycosidic bond by a single water molecule; however, both pathways also involve very large barriers. Activation of the water nucleophile via partial proton abstraction steadily decreases the barrier and leads to a more exothermic reaction energy as the proton affinity of the molecule interacting with water increases. Indeed, our data suggests that the barrier heights and reaction energies range from that for hydrolysis by water to that for hydrolysis by the hydroxyl anion, which represents the extreme of (full) water activation (deprotonation). Hydrogen bonds between small molecules (hydrogen fluoride, water, or ammonia) and the nucleobase were found to further decrease the barrier and overall reaction energy but not to the extent that the same hydrogen-bonding interactions increase the acidity of the nucleobase. Our results suggest that the nature of the nucleophile plays a more important role in reducing the barrier to glycosidic bond cleavage than the nature of the small molecule bound, and models with more than one hydrogen fluoride molecule interacting with the nucleobase provide further support for this conclusion. Our results lead to a greater fundamental understanding of the effects of the nucleophile, activation of the nucleophile, and interactions with the nucleobase for this important biological reaction.     

Effects of ionization, metal cationization and protonation on 2'-deoxyguanosine: changes on sugar puckering and stability of the N-glycosidic bond

The influence of oxidation, protonation, and metal cationization with Cu(+) and Cu(2+) on the strength of the N-glycosidic bond in 2'-deoxyguanosine has been studied by means of quantum chemical calculations. In all cases, the N9-C1' bond distance increases (0.03-0.06 A) upon introducing positive charge in the guanine moiety, the observed variations being more important for the dicationic systems. Binding energies show that the effect of X(n)(+) in guanine hinders the homolytic dissociation, whereas it largely favors the heterolytic process. With respect to the deoxyribose ring, it has been found that metal binding, oxidation, and protonation do not significantly change the values of the phase angle of pseudorotation P. However, the glycosyl torsion angle chi varies considerably (from 242.0 degrees to 189.8 degrees) as a consequence of a stabilizing guanine-sugar (H8-O4') interaction due to the increase of acidity of guanine C8-H8 upon cationization.     

Structural characterization of isoxazolidinyl nucleosides by fast atom bombardment tandem mass spectrometry.

4'-Aza-2',3'-dideoxyerythrofuranosyl derivatives of thymine (AdT, 1) and uracil (AdU, 2) are analogues of 2',3'-dideoxyribofuranosyl thymine (ddT, 3) and uracyl (ddU, 4). Compounds 1 and 2 are representative of a new class of antiviral agents where the sugar moiety is replaced by an isoxazolidine ring. The increasing importance of isoxazolidinyl nucleosides has encouraged the exploitation of simple mass spectrometric rules for unambiguously assigning their structure. The species 1, 2, 5 and 6 were therefore synthesized in order to evaluate the role of the basic centre of the modified sugar moiety in their gas-phase chemistry. The tandem mass spectra of these compounds are similar to those of the wild-type nucleosides and display fragment ions corresponding to [B + 2H](+),[M - BH](+) and [B + 27](+) species, where B is the nucleobase. The last species derives from a retrocycloaddition process which is less evident in 2'-deoxyribosides. This behaviour is consistent with protonation of the analytes at the pyrimidine rings. Model isoxazolidines, in which the nucleobase was replaced by a phenyl or a naphthyl moiety, displayed the expected behaviour of species with a localized charge on the N-O moiety of the isoxazolidine ring.     

Gas-phase protonation thermochemistry of adenosine

The goal of this work was to obtain a detailed insight on the gas-phase protonation energetic of adenosine using both mass spectrometric experiments and quantum chemical calculations. The experimental approach used the extended kinetic method with nanoelectrospray ionization and collision-induced dissociation tandem mass spectrometry. This method provides experimental values for proton affinity, PA(adenosine) = 979 +/- 1 kJ.mol (-1), and for the "protonation entropy", Delta p S degrees (adenosine) = S degrees (adenosineH (+)) - S degrees (adenosine) = -5 +/- 5 J.mol (-1).K (-1). The corresponding gas-phase basicity is consequently equal to: GB(adenosine) = 945 +/- 2 kJ.mol (-1) at 298K. Theoretical calculations conducted at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31+G(d,p) level, including 298 K enthalpy correction, predict a proton affinity value of 974 kJ.mol (-1) after consideration of isodesmic proton transfer reactions with pyridine as the reference base. Moreover, computations clearly showed that N3 is the most favorable protonation site for adenosine, due to a strong internal hydrogen bond involving the hydroxyl group at the 2' position of the ribose sugar moiety, unlike observations for adenine and 2'-deoxyadenosine, where protonation occurs on N1. The existence of negligible protonation entropy is confirmed by calculations (theoretical Delta p S degrees (adenosine) approximately -2/-3 J.mol (-1).K (-1)) including conformational analysis and entropy of hindered rotations. Thus, the calculated protonation thermochemical properties are in good agreement with our experimental measurements. It may be noted that the new PA value is approximately 10 kJ.mol (-1) lower than the one reported in the National Institute of Standards and Technology (NIST) database, thus pointing to a correction of the tabulated protonation thermochemistry of adenosine.     

The role of nucleobase carboradical and carbanion on DNA lesions: a theoretical study

DNA base release induced by H and OH radical addition to thymine and their corresponding electron adducts is studied at the DFT B3LYP/6-31+G(d,p) level in gas phase and in solution. H atom transfer after radical formation from C2' on the sugar to the C6 site on the base is shown to be prohibited for the radical species. Their corresponding electron adducts, albeit minor events in cellular systems, show excellent capabilities to proton transfer from C2' on the sugar to the C6 site on the base. The barriers for subsequent N-glycosidic bond dissociation range from 0.1 to 1.6 kcal mol(-1) at the B3LYP level and around 5 kcal mol(-1) using the BB1K functional, implying that these reactions can serve as a source to abasic sites. Analysis of bond dissociation energies show that all the reactions are exothermic, which is consistent with the changes in N-glycosidic bond lengths during the proton-transfer reactions. Bulk solvation plays a reverse influence on proton transfer and the bond rupture reactions. Molecular orbitals, NPA charges, and electron affinities are calculated to shed further light on the properties leading up to the intramolecular reactions.     

Effects of disulfide linkages on gas-phase reactions of small multiply charged peptide ions

Multiply protonated ions of disulfide-intact and -reduced peptides were generated by electrospray ionization and studied by Fourier transform ion cyclotron resonance mass spectrometry. The effects of disulfide bonds on gas-phase deprotonation reactions and hydrogen/deuterium (H/D) exchange were investigated. Insight into conformations was gained from molecular dynamics calculations. For ions from three small peptides containing 9–14 amino acid residues, H/D exchange is more sensitive to changes in conformation than deprotonation. However, with both gas-phase reactions the more diffuse forms of the peptides (as determined by molecular modeling) react more readily. The effects of disulfide cleavage on the conformations and on the reactions were found to depend upon the sequence of the peptide. For [M + 3H]³⁺ of TGF-α (34–43), reduction of the disulfide linkage leads to a greatly extended structure and a dramatic increase in the rate and extent of H/D exchange. In contrast, [M + 2H]²⁺ of Arg⁸ -vasopressin becomes slightly more compact upon cleavage of the disulfide bond; these reduced ions are slower to react. For [M + 3H]³⁺ of somatostatin-14, reduction of the disulfide bond has little effect on conformation or gas-phase reactivity. Overall, these results indicate that there is no general rule on how cleavage of a disulfide bond will effect a peptide ion’s gas-phase reactivity.     

Transition-state analysis of the DNA repair enzyme MutY

The transition state (TS) structure of MutY-catalyzed DNA hydrolysis was solved using multiple kinetic isotope effect (KIE) measurements. MutY is a base excision repair enzyme which cleaves adenine from 8-oxo-G:A mismatches in vivo, and also from G:A mismatches in vitro. TS analysis of G:A-DNA hydrolysis revealed a stepwise S(N)1 (D(N)*A(N)(double dagger)) mechanism proceeding through a highly reactive oxacarbenium ion intermediate which would have a lifetime in solution of <10(-10) s. C-N bond cleavage is reversible; the N-glycoside bond breaks and reforms repeatedly before irreversible water attack on the oxacarbenium ion. KIEs demonstrated that MutY uses general acid catalysis by protonating N7. It enforces a 3'-exo sugar ring conformation and other sugar ring distortions to stabilize the oxacarbenium ion. Combining the experimental TS structure with the previously reported crystal structure of an abortive Michaelis complex elucidates the step-by-step catalytic sequence.     

Arrhenius activation parameters for the loss of neutral nucleobases from deprotonated oligonucleotide anions in the gas phase

Arrhenius activation parameters (E(a) and A) for the loss of neutral nucleobase from a series of doubly deprotonated oligodexoynucleotide 10-mers of the type XT(9), T(9)X, and T(5)XT(4), where X = A, C, and G, have been determined using the blackbody infrared radiative dissociation technique. At temperatures of 120 to 190 degrees C, the anions dissociate exclusively by the loss of a neutral nucleobase (XH), followed by cleavage of the sugar 3' C-O bond leading to (a-XH) and w type ions or, in the case of the T(9)X(2-) ions, the loss of H(2)O. The dissociation kinetics and energetics are sensitive to the nature and position of X. Over the temperature range investigated, the kinetics for the loss of AH and GH were similar, but approximately 100 times faster than for the loss of CH. For the loss of AH and GH, the values of E(a) are sensitive to the position of the base. The order of the E(a)s for the loss of XH from the 5' and 3' termini is: C > G > A; while for T(5)XT(4) the order is: C > A > G. The trends in the values of E(a) do not parallel the trend in deprotonation enthalpies or proton affinities of the nucleobases in the gas phase, indicating that the energetic differences do not simply reflect differences in their gas phase acidity or basicity. The pre-exponential factors (A) vary from 10(10) to 10(15) s(-1), depending on the nature and position of X. These results suggest that the reactivity of individual nucleobases is influenced by stabilizing intramolecular interactions.     

Electron attachment-induced DNA single strand breaks: C3'-O3' sigma-bond breaking of pyrimidine nucleotides predominates

A detailed understanding of DNA strand breaks induced by low energy electrons (LEE) is of crucial importance for the advancement of many areas of molecular biology and medicine. To elucidate the mechanism of DNA strand breaks by LEEs, theoretical investigations of the electron attachment-induced C3'-O3' sigma-bond breaking of the pyrimidine nucleotides have been performed. Calculations of 2'-deoxycytidine-3'-monophosphate and 2'-deoxythymidine-3'-monophosphate in their protonated form (denoted as 3'-dCMPH and 3'-dTMPH) have been carried out with the reliably calibrated B3LYP/DZP++ theoretical approach. Our results demonstrate that the transfer of the negative charge from the pi*-orbital of the radical anion of pyrimidines to the DNA backbone does not pass through the N1-glycosidic bond. Instead, the migration of the excessive negative charge through the atomic orbital overlap between the C6 of pyrimidine and the C3' of ribose most likely represents a pathway that subsequently leads to the strand breaks. The proposed mechanism of the LEE-induced single strand breaks in DNA assumes that the formation of the base-centered radical anions is the first step in this process. Subsequently, these electronically stable radical anions may undergo either C-O bond breaking or N-glycosidic bond rupture. The present investigation of 3'-dCMPH and 3'-dTMPH yields an energy barrier of 6.2-7.1 kcal/mol for the C3'-O3' sigma-bond cleavage. This is much lower than the energy barriers required for the C5'-O5' sigma-bond and the N1-glycosidic bond break. Therefore, we conclude that the C3'-O3' sigma-bond rupture dominates the LEE-induced single strand breaks of DNA.     

The chemical nature of the 2'-substituent in the pentose-sugar dictates the pseudoaromatic character of the nucleobase (pKa) in DNA/RNA

We here show that the pKa (error limit: 0.01 to 0.03 pKa unit) of a nucleobase in a nucleotide can be modulated by the chemical nature of the 2'-substituent at the sugar moiety. This has been evidenced by the measurement of nucleobase pKa in 47 different model nucleoside 3',5'-bis- and 3'-mono-ethylphosphates. The fact that the electronic character of each of the 2'-substituents (Fig. 1) alters the chemical shift of the H2' sugar proton, and also alters the pKa of the nucleobase in the nucleotides has been evidenced by a correlation plot of pKa of N3 of pyrimidine (T/C/U) or pKa of N7 of 9-guaninyl with the corresponding deltaH2' chemical shifts at the neutral pH, which shows linear correlation with high Pearson's correlation coefficients (R = 0.85-0.97). That this modulation of the pKa of the nucleobase by a 2'-substituent is a through-bond as well as through-space effect has been proven by ab initio determined pKa estimation. Interestingly, experimental pKas of nucleobases from NMR titration and the calculated pKas (by ab initio calculations utilizing closed shell HF 6-31G** basis set) are linearly correlated with R = 0.98. It has also been observed that the difference of ground and protonated/de-protonated HOMO orbital energies (DeltaHOMO, a.u.) for the nucleobases (A/G/C/T/U) are well correlated with their pK(a)s in different 2'-substituted 3',5'-bis-ethylphosphate analogs suggesting that only the orbital energy of HOMO can be successfully used to predict the modulation of the chemical reactivity of the nucleobase by the 2'-substituent. It has also been demonstrated that pKa values of nucleobases in 3',5'-bis-ethylphosphates (Table 1) are well correlated with the change in dipole moment for the respective nucleobases after protonation or de-protonation. This work thus unambiguously shows that alteration of the thermodynamic stability (Tm) of the donor-acceptor complexes [ref. 20], as found with various 2'-modified duplexes in the antisense, siRNA or in triplexes by many workers in the field, is a result of alteration of the pseudoaromatic character of the nucleobases engineered by alteration of the chemical nature of the 2'-substitution.     

Charge state effect on the zwitterion influence on stability of non-covalent interaction of single-stranded DNA with peptides

Negative ion ESI mass spectrometry was used to study the gas-phase stability and dissociation pathways of peptide-DNA complexes. We show that bradykinin and three modified peptides containing the basic residue arginine or lysine form stable interactions with single-stranded oligonucleotides. ESI-MS/MS of complexes of T(8) with PPGFSPFRR resulted in a major dissociation pathway through cleavage of the peptide covalent bond. The stability of the complex is due to electrostatic interaction between the negatively charged phosphate group and the basic side chain of the arginine and lysine residues as demonstrated by Vertes et al. and Woods et al. In fact, the present work establishes the role played by zwitterions on complex stabilisation. The presence of protons in nucleobase and/or amino acid contributes in reinforcing the strength of the salt bridge (SB) interaction. The zwitterionic form of the most basic of amino acid residues, arginine, is assumed to form a strong SB interaction to the negatively charged phosphate groups of DNA. This non-covalent complex is stable enough to withstand disruption of the non-covalent interaction and to first break the covalent bond. Moreover, the dependence of fragmentation patterns upon the complex charge state is explained by the fact that the net number of negative charges modulates the number of zwitterionic sites, which stabilise the complexes. Finally, the weak influence of the nucleobase is assumed by the existence of competition for proton addition between the nucleobase and the R/K side chain leading to a decrease in the stabilisation of the SB interaction.     

Metastable decay of negatively charged oligodeoxynucleotides analyzed with ultraviolet matrix-assisted laser desorption/ionization post-source decay and deuterium exchange

In an effort to understand the initiating step in metastable-ion decay of UV matrix-assisted laser desorption/ionization (MALDI)-produced ions, we conducted experiments in which we exchanged all active protons for deuterons of tetrameric and hexameric oligodeoxynucleotides. We wish to address the controversial proposal that in the negative-ion mode, as in the positive-ion mode, fragmentation is driven by nucleobase protonation. The results show unambiguously that proton transfer, leading to zwitterion formation, charges a nucleobase prior to its elimination. The zwitterion formation directs fragmentation of both positive and negative oligodeoxynucleotide ions. Poly-T-rich oligodeoxynucleotide tetramers show surprising differences in the negative compared to the positive-ion mode, as thymine is preferentially expelled, instead of a nucleobase with higher proton affinity. For the exceptional case of negatively charged poly-T-rich oligodeoxynucleotide tetramers generated by MALDI, we propose that zwitterion formation with positive charging of a nucleobase leads to base stabilization in the negative-ion mode through an interaction of the positive nucleobase with the excess negative charge. Moreover, backbone cleavages (accompanied by H rearrangement) of a phosophodiester bond give first-generation products that can be traced back to this tripolar complex bearing one positive and two negative charges, all of which may be interacting.     

Influence of N7 protonation on the mechanism of the N-glycosidic bond hydrolysis in 2'-deoxyguanosine. A theoretical study

The influence of N7 protonation on the mechanism of the N-glycosidic bond hydrolysis in 2'-deoxyguanosine has been studied using density functional theory (DFT) methods. For the neutral system, two different pathways (with retention and inversion of configuration at the C1' anomeric carbon) have been found, both of them consisting of two steps and involving the formation of a dihydrofurane-like intermediate. The Gibbs free energy barrier for the first step is very high in both cases (53 and 46 kcal/mol for the process with inversion and with retention, respectively). However, the N7-protonated system shows a very different mechanism which consists of two steps. The first one leads to the formation of an oxacarbenium ion intermediate, with a Gibbs free energy barrier of 27 kcal/mol, and the second one corresponds to the nucleophilic attack of the water molecule to the oxacarbenium ion and takes place with a barrier of 1.3 kcal/mol. Thus, these results agree with a stepwise SN1 mechanism (DN*AN), with a discrete intermediate formed between the leaving group and the nucleophile approach, and show that N7 protonation strongly catalyzes the hydrolysis of the N-glycosidic bond, making the guanine a better leaving group. Finally, kinetic isotope effects have been calculated for the protonated system, and the results obtained are in very good agreement with experimental data for analogous systems.     

Computational investigation and hydrogen/deuterium exchange of the fixed charge derivative tris(2,4,6-Trimethoxyphenyl) phosphonium: Implications for the aspartic acid cleavage mechanism

Aspartic acid (Asp)-containing peptides with the fixed charge derivative tris(2,4,6trimethoxyphenyl) phosphonium (tTMP-P+) were explored computationally and experimentally by hydrogen/deuterium (H/D) exchange and by fragmentation studies to probe the phenomenon of selective cleavage C-terminal to Asp in the absence of a "mobile" proton. Ab initio modeling of the tTMP-P+ electrostatic potential shows that the positive charge is distributed on the phosphonium group and therefore is not initiating or directing fragmentation as would a "mobile" proton. Geometry optimizations and vibrational analyses of different Asp conformations show that the Asp structure with a hydrogen bond between the side-chain hydroxy and backbone carbonyl lies 2.8 kcal/mol above the lowest energy conformer. In reactions with D2O, the phosphonium-derived doubly charged peptide (H+)P+LDIFSDF rapidly exchanges all 12 of its exchangeable hydrogens for deuterium and also displays a nonexchanging population. With no added proton, P+LDIFSDF exchanges a maximum of 4 of 11 exchangeable hydrogens for deuterium. No exchange is observed when all acidic groups are converted to the corresponding methyl esters. Together, these H/D exchange results indicate that the acidic hydrogens are "mobile locally" because they are able to participate in exchange even in the absence of an added proton. Fragmentation of two distinct (H+)P+LDIFSDF ion populations shows that the nonexchanging population displays selective cleavage, whereas the exchanging population fragments more evenly across the peptide backbone. This result indicates that H/D exchange can sometimes distinguish between and provide a means of separation of different protonation motifs and that these protonation motifs can have an effect on the fragmentation.     

Interaction of Tertiary Phosphines with Acetylenic Compouxds. Alkyl Migration Accompanied by C-P Bond Cleavage and Fragmentation During the Interaction of Trialkyphosphines with Phenylacetylene in the Presence of Proton Donors

Abstract

We have recently describedphosphobetaines with negative charge on β-carbon atom of vinyl group or on δ-carbon atom of l 1,3-butadienyl group prepared by interaction of trial-kylphosphines with acetylenic and vinylacetylenic compounds, respectively. In the course of developing these investigations an extra route of phosphobetaine formation has been found representing the deprotonation of corresponding phosphonium salts. It has been found that the reaction of trialkylphospliines with phenylacetylene in the presence of proton donors leads to the formation of an intermediate compound containing pentacovalent phosphorus linked with negatively charged oxygen atom. The intermediate undergoes a migration of alkyl group followed either by protonation of β-carbon atom or by unusual cleavage of P-C bond as a result of the electrons back transition:     

Metal‐stabilized rare tautomers: N4 metalated cytosine (M = Li+, Na+, K+, Rb+ and Cs+), theoretical views

Ab initio calculations indicate that metalation of the exocyclic amino group of cytosine by the elements of Group IA (Li, Na, K, Rb and Cs) induces protonation of a nucleobase ring nitrogen atom, and hence causes a proton shift from an exocyclic to an endocyclic nitrogen atom. Thus, this metal‐assisted process leads to the generation of rare nucleobase tautomers. The calculations suggest that this kind of metalation increases the protonation energies of the aromatic ring of the nucleobase. The present study reports the quantum chemistry analysis of the metal‐assisted tautomerization. The calculations clearly demonstrate that metalation of the exocyclic amino group of the nucleobase significantly increases the protonation energy of the aromatic rings of the nucleobase. Also, absolute anisotropy shift, molecular orbital and natural bond orbital calculations are compatible with these results. Copyright © 2003 John Wiley & Sons, Ltd.