Theoretical Study of the Interaction between Thymine and Water. Protonation and Deprotonation Enthalpies and Comparison with Uracil

Uracil-water interaction revisited - in search of single H-bonded secondary minima

Monohydrated uracil (UW) complexes are stabilized by both O⋯HO and NH⋯O hydrogen bonds (H-bonds), simultaneously participating in forming three stable cyclic structures. The role and contribution of these individual H-bonds (O⋯HO and NH⋯O) to the stability of the three UW complexes are still not understood, because of the technical problems in obtaining their optimized structures by standard geometry optimization. The present study explores a non-standard approach to identify three single H-bonded local minima structures without imaginary frequency using DFT (M06-2X, B3LYP and B3LYP-D3), MP2 and CCSD(T) theories and Dunning's correlation-consistent aug-cc-pVTZ basis set, in both vacuum and aqueous media (CPCM method). The results reveal that these new structures are very shallow secondary minima between two deep wells or next to a deep well of primary minima (double H-bonded structures) in the potential energy surface. The H-bond energy of these single H-bonded complexes is found to be less sensitive to a wide range (about 15-20 degrees) of O⋯HO and NH⋯O angles, and the linearity is preferred in the stable three single H-bonded structures. The technical method used to locate such a shallow minimum is described in detail and may be useful for identifying local minima in other cases where consecutive multiple H-bonded structures are global minima. Energy decomposition (using symmetry adapted perturbation theory (SAPT)) of interaction energy, electron redistribution, and relevant vibrational modes are discussed.

Density Functional Method Study on the Cooperativity of Intermolecular H-bonding and π-π+ Stacking Interactions in Thymine-[Cnmim]Br (n = 2, 4, 6, 8, 10) Microhydrates

The exploration of the ionic liquids’ mechanism of action on nucleobase’s structure and properties is still limited. In this work, the binding model of the 1-alkyl-3-methylimidazolium bromide ([Cnmim]Br, n = 2, 4, 6, 8, 10) ionic liquids to the thymine (T) was studied in a water environment (PCM) and a microhydrated surroundings (PCM + wH2O). Geometries of the mono-, di-, tri-, and tetra-ionic thymine (T-wH2O-y[Cnmim]+-xBr−, w = 5~1 and x + y = 0~4) complexes were optimized at the M06-2X/6-311++G(2d, p) level. The IR and UV-Vis spectra, QTAIM, and NBO analysis for the most stable T-4H2O-Br−-1, T-3H2O-[Cnmim]+-Br−-1, T-2H2O-[Cnmim]+-2Br−-1, and T-1H2O-2[Cnmim]+-2Br−-1 hydrates were presented in great detail. The results show that the order of the arrangement stability of thymine with the cations (T-[Cnmim]+) by PCM is stacking > perpendicular > coplanar, and with the anion (T-Br−) is front > top. The stability order for the different microhydrates is following T-5H2O-1 < T-4H2O-Br−-1 < T-3H2O-[Cnmim]+-Br−-1 < T-2H2O-[Cnmim]+-2Br−-1 < T-1H2O-2[Cnmim]+-2Br−-1. A good linear relationship between binding EB values and the increasing number (x + y) of ions has been found, which indicates that the cooperativity of interactions for the H-bonding and π-π+ stacking is varying incrementally in the growing ionic clusters. The stacking model between thymine and [Cnmim]+ cations is accompanied by weaker hydrogen bonds which are always much less favorable than those in T-xBr− complexes; the same trend holds when the clusters in size grow and the length of alkyl chains in the imidazolium cations increase. QTAIM and NBO analytical methods support the existence of mutually reinforcing hydrogen bonds and π-π cooperativity in the systems.

Acid-base and lipophilic properties of peptide nucleic acid derivatives

The first combined experimental and theoretical study on the ionization and lipophilic properties of peptide nucleic acid (PNA) derivatives, including eleven PNA monomers and two PNA decamers, is described. The acidity constants (pK_a) of individual acidic and basic centers of PNA monomers were measured by automated potentiometric pH titrations in water/methanol solution, and these values were found to be in agreement with those obtained by MoKa software. These results indicate that single nucleobases do not change their pK_a values when included in PNA monomers and oligomers. In addition, immobilized artificial membrane chromatography was employed to evaluate the lipophilic properties of PNA monomers and oligomers, which showed the PNA derivatives had poor affinity towards membrane phospholipids, and confirmed their scarce cell penetrating ability. Overall, our study not only is of potential relevance to evaluate the pharmacokinetic properties of PNA, but also constitutes a reliable basis to properly modify PNA to obtain mimics with enhanced cell penetration properties.

•

The first study on acid-base and lipophilic properties of peptide nucleic acids (PNA).

•

pK_a of acid-base centers of PNA evaluated by potentiometric method and MoKa prediction.

•

NMR experiments provide additional information on the protonation of PNA monomers.

•

Lipophilicity of PNA monomers and oligomers is investigated by IAM chromatography.

•

This study can lay the basis of evaluating the pharmacokinetic properties of PNA.

Reactions of Cytosine, Thymine and Uracil toward 2,4,6-trinitro-chlorobenzene

The title reactions are an example of direct introduction of an aryl group into DNA bases.

Computational chemistry applied to vibrational spectroscopy: A tool for characterization of nucleic acid bases and some of their 5-substituted derivatives

Computational chemistry can be applied to vibrational spectroscopy in different ways, such as for a better characterization and assignment of all the bands of the experimental spectra, as a tool in the identification of the tautomers present in the gas phase and in the solid state through their spectra and for the simulation of the solid and liquid phase of a compound and the consequent simulation and interpretation of their spectra. In the present study, as an example of the applicability of computational chemistry, the structure and spectra of cytosine and uracil nucleic acid bases and two cytosine derivatives are shown. The FTIR and Raman spectra were analysed with the support of ab initio (Hartree-Fock (HF), MP2) and density functional theory (DFT) (B3LYP, PBE, B-P, etc.) calculations using several basis sets and several scaling equations. The calculations predict an easier tautomerization of cytosine than uracil molecule, but the tautomerization is hindered in the 5-bromocytosine molecule. Thus, in the solid state, this molecule only exists in the amino-oxo tautomeric form.

Quantitative Single-Molecule Surface-Enhanced Raman Scattering by Optothermal Tuning of DNA Origami-Assembled Plasmonic Nanoantennas

DNA origami is a powerful approach for assembling plasmonic nanoparticle dimers and Raman dyes with high yields and excellent positioning control. Here we show how optothermal-induced shrinking of a DNA origami template can be employed to control the gap sizes between two 40 nm gold nanoparticles in a range from 1 to 2 nm. The high field confinement achieved with this optothermal approach was demonstrated by detection of surface-enhanced Raman spectroscopy (SERS) signals from single molecules that are precisely placed within the DNA origami template that spans the nanoparticle gap. By comparing the SERS intensity with respect to the field enhancement in the plasmonic hot-spot region, we found good agreement between measurement and theory. Our straightforward approach for the fabrication of addressable plasmonic nanosensors by DNA origami demonstrates a path toward future sensing applications with single-molecule resolution.

Microhydration of Deprotonated Nucleobases

Hydration reactions of deprotonated nucleobases (uracil, thymine, 5-fluorouracil,2-thiouracil, cytosine, adenine, and hypoxanthine) produced by electrospray have been experimentally studied in the gas phase at 10 mbar using a pulsed ion-beam high-pressure mass spectrometer. The thermochemical data, ΔH (o) , ΔS (o) , and ΔG (o) , for the monohydrated systems were determined. The hydration enthalpies were found to be similar for all studied systems and varied between 39.4 and 44.8 kJ/mol. A linear correlation was found between water binding energies in the hydrated complexes and the corresponding acidities of the most acidic site of nucleobases. The structural and energetic aspects of the precursors for the hydrated complexes are discussed in conjunction with available literature data. Graphical Abstract ᅟ.

Simulated Raman correlation spectroscopy for quantifying nucleic acid-silver composites

Plasmonic devices are of great interest due to their ability to confine light to the nanoscale level and dramatically increase the intensity of the electromagnetic field, functioning as high performance platforms for Raman signal enhancement. While Raman spectroscopy has been proposed as a tool to identify the preferential binding sites and adsorption configurations of molecules to nanoparticles, the results have been limited by the assumption that a single binding site is responsible for molecular adsorption. Here, we develop the simulated Raman correlation spectroscopy (SRCS) process to determine which binding sites of a molecule preferentially bind to a plasmonic material and in what capacity. We apply the method to the case of nucleic acids binding to silver, discovering that multiple atoms are responsible for adsorption kinetics. This method can be applied to future systems, such as to study the molecular orientation of adsorbates to films or protein conformation upon adsorption.

Microhydration effects on the structures and electrophilic properties of cytidine

Adiabatic electron affinities (AEAs) for cytidine hydrates with up to four water molecules.

Gas-phase acidity and dynamics of the protonation processes of carbidopa and levodopa. A QM/QD study

The present work details, our efforts to investigate and compare the acid-base properties of levodopa (LD) and carbidopa (CD), a drug combination being used in the treatment of Parkinson's disease. Protonation and deprotonation were examined for all possible sites in the two drugs. Proton affinity and proton detachment enthalpies were computed at the B3LYP/6-311++G** level of theory. Results of the present work reveal that CD is more basic and can abstract protons in solution much more efficiently than LD. Furthermore, at all deportation sites considered, CD is more acidic than LD. DFT-based ADMP, dynamic simulations have been performed to explore the dynamics of the protonation processes in LD and CD. Thus, while the dynamics of the protonation process of LD is very straightforward leading to the formation of the corresponding cation, the protonation process in CD is very complex involving a major geometry change and rearrangement. Results of the present work reveal that the active species in acid medium is not CD in its normal geometry but a carbonyl hydrazine form instead. The presence of the carbonyl group β to the hydrazine group may very well underlie its enhanced activity which allows it to bind to the active site of the DDC enzyme. The relative stabilities of various water-water-CD complexes have been computed and compared.

Structure of protonated thymidine characterized by infrared multiple photon dissociation and quantum calculations

RATIONALE

Many fundamental studies are motivated by the probable relationship between the presence of rare enol tautomers of nucleobases and point mutation developing during nucleic acid replication. The evaluation of the tautomeric behaviour of nucleobases is therefore of fundamental importance. This can be probed in the gas phase by combining action spectroscopy and mass spectrometry.

METHODS

Experimental Infrared Multiple Photon Dissociation spectra in the fingerprint region of electrospray-generated and subsequently selected ions were recorded at the CLIO free electron laser (FEL) facility, by coupling the FEL to a quadrupole ion trap, and compared to calculated harmonic vibrational infrared spectra of the different low-lying isomers computed at the B3LYP/6-31++G(d,p) level. Relative energies were refined using the extended basis set 6-311++G(3df,2p).

RESULTS

The Density Functional Theory (DFT) study shows that, as for protonated thymine, the global energy minimum of protonated thymidine corresponds to an enol tautomer, whose infrared absorption spectrum is found to be in very good agreement with the experimental IRMPD spectrum. A very weak IRMPD signal observed at ~1780 cm(-1) is very likely the signature of an oxo tautomer. Consequently, as for thymine, protonated thymidine generated by electrospray corresponds to a mixture of at least two tautomeric forms.

CONCLUSIONS

Tautomerization can be characterized by IRMPD spectroscopy. Interestingly, the dominant enolic tautomeric form(s) presently observed cannot be directly generated from the most stable neutral tautomer of the thymine residue.

Structural Characterization of Single-Stranded DNA Monolayers Using Two-Dimensional Sum Frequency Generation Spectroscopy

DNA-covered materials are important in technological applications such as biosensors and microarrays, but obtaining structural information on surface-bound biomolecules is experimentally challenging. In this paper, we structurally characterize single-stranded DNA monolayers of poly(thymine) from 10 to 25 bases in length with an emerging surface technique called two-dimensional sum frequency generation (2D SFG) spectroscopy. These experiments are carried out by adding a mid-IR pulse shaper to a femtosecond broad-band SFG spectrometer. Cross peaks and 2D line shapes in the 2D SFG spectra provide information about structure and dynamics. Because the 2D SFG spectra are heterodyne detected, the monolayer spectra can be directly compared to 2D infrared (2D IR) spectra of poly(thymine) in solution, which aids interpretation. We simulate the 2D SFG spectra using DFT calculations and an excitonic Hamiltonian that relates the molecular geometry to the vibrational coupling. Intrabase cross peaks help define the orientation of the bases and interbase cross peaks, created by coupling between bases, and resolves features not observed in 1D SFG spectra that constrain the relative geometries of stacked bases. We present a structure for the poly(T) oligomer that is consistent with the 2D SFG data. These experiments provide insight into the DNA monolayer structure and set precedent for studying complex biomolecules on surfaces with 2D SFG spectroscopy.

An experimental and theoretical vibrational study of interaction of adenine and thymine with artificial seawaters: A prebiotic chemistry experiment

Nucleic acid bases play important roles in living beings. Thus, their interaction with salts the prebiotic Earth could be an important issue for the understanding of origin of life. In this study, the effect of pH and artificial seawaters on the structure of adenine and thymine was studied via parallel determinations using FT-IR, Raman spectroscopy and theoretical calculations. Thymine and adenine lyophilized in solutions at basic and acidic conditions showed characteristic bands of the enol-imino tautomer due to the deprotonation and the hydrochloride form due to protonation, respectively. The interaction of thymine and adenine with different seawaters representative of different geological periods on Earth was also studied. In the case of thymine a strong interaction with Sr(2+) promoted changes in the Raman and infrared spectra. For adenine changes in infrared and Raman spectra were observed in the presence of salts from all seawaters tested. The experimental results were compared to theoretical calculations, which showed structural changes due to the presence of ions Na(+), Mg(2+), Ca(2+) and Sr(2+) of artificial seawaters. For thymine the bands arising from C4=C5 and C6=O stretching were shifted to lower values, and for adenine, a new band at 1310cm(-1) was observed. The reactivity of adenine and thymine was studied by comparing changes in nucleophilicity and energy of the HOMO orbital.

Structures, electrophilic properties, and hydrogen bonds of cytidine, uridine, and their radical anions: Microhydration effects

Structures, electrophilic properties, and hydrogen bonds of the neutral and anionic monohydrated nucleoside, (cytidine)H2O, and (uridine)H2O have been systematically investigated using density functional theory. Various water-binding sites were predicted by explicitly considering the optimized monohydrated structures. Meanwhile, predictions of electron affinities and vertical detachment energies were also carried out to investigate their electrophilic properties. By examining the singly occupied molecular orbital and natural population analysis, we found the excess negative charge is localized on the cytidine and uridine moiety in anionic monohydrates. This may be the reason why the strength of hydrogen bonding undergoes an obvious change upon the extra electron attachment. Based on the electron density (ρ) and reduced density gradient (RDG), we present an approach to map and analyze the weak interaction (especially hydrogen bond) in monohydrated cytidine and uridine. The scatter plots of RDG versus ρ allow us to identify the different type interactions. Meanwhile, the maps of the gradient isosurfaces show a rich visualization of hydrogen bond, van der Waals interaction, and steric effect.

Polarization in the structures of uracil and thiouracils: Implication for binding with orotidine 5'-monophosphate decarboxylase

The structures of the uracil and thiouracils were examined using NMR spectroscopy and crystal structure data when available. The relationships between the extent of polarization and the C5-C6 bond length as well as the H5-H6 coupling constants were probed. It was found that the bond length and coupling constants correlate well with the proton affinities at the carbonyl or thiocarbonyl groups at C4 but not C2. The possible implication in the tighter binding of thiouracil based nucleotides to orotidine-5'-monophosphate decarboxylase was discussed.

Specific photodynamics in thymine clusters: the role of hydrogen bonding

A photoionization detected IR study of thymine and 1-methylthymine monohydrates and of their homodimers was carried out to shed some light on the structure of the thymine clusters whose complex photodynamics has recently been the subject of great interest. Under supersonic jet conditions, thymine forms doubly H-bonded cyclic clusters with water or another base preferentially via its N1-H group and the adjacent carbonyl group. This hydrate is of no biological relevance since the N1-H group is the sugar binding site in thymidine. On the other hand, 1-methylthymine forms the donor H-bonds only via the N3-H group. Hence, properties of the N1-H and the N3-H bound clusters of thymine can be studied using thymine and 1-methylthymine molecules, respectively. No biologically relevant conformations of the dimers and hydrates of thymine, contrary to those of 1-methylthymine, are observed under supersonic jet conditions. Thymine homodimer, which extensively fragments upon UV ionization by formation of a protonated monomer, exhibits two N1-H···O═C2 hydrogen bonds. The photodynamics of hydrated thymines is found to be extremely sensitive to the hydration site: ranging from an ultrafast relaxation in less than 100 fs up to formation of a dark state with the lifetime on the microsecond time scale.

Adsorption of thymine and uracil on 1:1 clay mineral surfaces: comprehensive ab initio study on influence of sodium cation and water

This computational study performed using the density functional theory shows that hydrated and non-hydrated tetrahedral and octahedral kaolinite mineral surfaces in the presence of a cation adsorb the nucleic acid bases thymine and uracil well. Differences in the structure and chemistry of specific clay mineral surfaces led to a variety of DNA bases adsorption mechanisms. The energetically most predisposed positions for an adsorbate molecule on the mineral surface were revealed. The target molecule binding with the surface can be characterized as physisorption, which occurs mainly due to a cation-molecular oxygen interaction, with hydrogen bonds providing an additional stabilization. The adsorption strength is proportional to the number of intermolecular interactions formed between the target molecule and the surface. From the Atoms in Molecules analysis and comparison of binding energy values of studied systems it is concluded that the sorption activity of kaolinite minerals for thymine and uracil depends on various factors, among which are the structure and accessibility of the organic compounds. The adsorption is governed mostly by the surface type, its properties and presence of cation, which cause a selective binding of the nucleobase. Adsorbate stabilization on the mineral surface increases only slightly with explicit addition of water. Comparison of activity of different studied kaolinite mineral models reveals the following order for stabilization: octahedral-Na-water > octahedral-Na > tetrahedral-Na > tetrahedral-Na-water. Further investigation of the electrostatic potentials helps understanding of the adsorption process and confirmation of the active sites on the kaolinite mineral surfaces. Based on the conclusions that clay mineral affinity for DNA and RNA bases can vary due to different structural and chemical properties of the surface, a hypothesis on possible role of clays in the origin of life was made.

Dielectric properties of fully hydrated nucleotides in the terahertz frequency range

We use terahertz time domain spectroscopy (THz-TDS) to determine the complex frequency-dependent dielectric response of all four nucleotides at different dilute concentrations. In addition, the suspension model's ability to extract the dielectric response of just the nucleotide with the hydration shell epsilon(b) excluding the dielectric information relating to the bulk will be verified. The suspension model enables us to make the determination that the nucleotides have influences on the water molecules out to the fourth hydration shell. We use a two Debye relaxation fit model for water, all concentrations and all epsilon(b) values. We observed how the nucleotides affect the relaxation parameters in relation to that of pure bulk water. With this information, we notice a transition between purines and pyrimidines, where one is a hydrogen-bond network structure building type material with a low concentration increment and the other is a structure breaking type material with a low concentration decrement. Due to conductivity measurements, we determine that kinetic depolarization is a negligible affect compared to that of dielectric saturation, which we find to dominate where a decrement is found.

An ab initio and AIM investigation into the hydration of 2-thioxanthine

BACKGROUND

Hydration is a universal phenomenon in nature. The interactions between biomolecules and water of hydration play a pivotal role in molecular biology. 2-Thioxanthine (2TX), a thio-modified nucleic acid base, is of significant interest as a DNA inhibitor yet its interactions with hydration water have not been investigated either computationally or experimentally. Here in, we reported an ab initio study of the hydration of 2TX, revealing water can form seven hydrated complexes.

RESULTS

Hydrogen-bond (H-bond) interactions in 1:1 complexes of 2TX with water are studied at the MP2/6-311G(d, p) and B3LYP/6-311G(d, p) levels. Seven 2TX...H2O hydrogen bonded complexes have been theoretically identified and reported for the first time. The proton affinities (PAs) of the O, S, and N atoms and deprotonantion enthalpies (DPEs) of different N-H bonds in 2TX are calculated, factors surrounding why the seven complexes have different hydrogen bond energies are discussed. The theoretical infrared and NMR spectra of hydrated 2TX complexes are reported to probe the characteristics of the proposed H-bonds. An improper blue-shifting H-bond with a shortened C-H bond was found in one case. NBO and AIM analysis were carried out to explain the formation of improper blue-shifting H-bonds, and the H-bonding characteristics are discussed.

CONCLUSION

2TX can interact with water by five different H-bonding regimes, N-H...O, O-H...N, O-H...O, O-H...S and C-H...O, all of which are medium strength hydrogen bonds. The most stable H-bond complex has a closed structure with two hydrogen bonds (N(7)-H...O and O-H...O), whereas the least stable one has an open structure with one H-bond. The interaction energies of the studied complexes are correlated to the PA and DPE involved in H-bond formation. After formation of H-bonds, the calculated IR and NMR spectra of the 2TX-water complexes change greatly, which serves to identify the hydration of 2TX.

Microhydration of protonated nucleic acid bases and protonated nucleosides in the gas phase

Thermochemical data, DeltaH(o)(n), DeltaS(o)(n), and DeltaG(o)(n), for the hydration of protonated nucleic acid bases and protonated nucleosides have been experimentally studied by equilibrium measurements using an electrospray high-pressure mass spectrometer equipped with a pulsed ion-beam reaction chamber. For protonated nucleobases the hydration enthalpies were found to be similar for all studied systems and varied between 12.4-13.1 kcal/mol for the first and 11.2-11.5 kcal/mol for the second water molecule. While for protonated nucleosides the water binding enthalpies (11.7-13.3 kcal/mol) are very close to those for protonated nucleobases, the entropy values are "more negative." The structural and energetic aspects of hydrated ions are discussed in conjunction with the available theoretical data.

A DFT study of thymine and its tautomers

The structures and the relative energies of six possible tautomers of the thymine base have been studied by density functional theory (DFT) using the B3LYP and BP86 functionals. The keto-thymine (T1) is predicted to be the most stable thymine tautomer, which is consistent with the other theoretical results and experimental data. The corresponding thymine cations and anions are studied using the same level of theory with double-ζ plus polarization and diffuse functions (DZP++) basis sets. The ionization potentials (IPs), the electron affinities (EAs), and proton affinities (PAs) for different protonation sites in thymine base are obtained. T1 has the largest ionization potential and the lowest proton affinity among all the considered thymine tautomers.

Infrared spectra of protonated uracil, thymine and cytosine

The gas-phase structures of protonated uracil, thymine, and cytosine are probed by using mid-infrared multiple-photon dissociation (IRMPD) spectroscopy performed at the Free Electron Laser facility of the Centre Laser Infrarouge d'Orsay (CLIO), France. Experimental infrared (IR) spectra are recorded for ions that were generated by electrospray ionization, isolated, and then irradiated in a quadrupole ion trap; the results are compared to the calculated infrared absorption spectra of the different low-lying isomers (computed at the B3LYP/6-31++G(d,p) level). For each protonated base, the global energy minimum corresponds to an enolic tautomer, whose infrared absorption spectrum matched very well with the experimental IRMPD spectrum, with the exception of a very weak IRMPD signal observed at about 1800 cm(-1) in the case of the three protonated bases. This signal is likely to be the signature of the second-energy-lying oxo tautomer. We thus conclude that within our experimental conditions, two tautomeric ions are formed which coexist in the quadrupole ion trap.

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.

Effects of OH Radical Addition on Proton Transfer in the Guanine−Cytosine Base Pair

Double proton transfer (PT) reactions in guanine-cytosine OH radical adducts are studied by the hybrid density functional B3LYP approach. Concerted and stepwise proton-transfer processes are explored between N1(H) on guanine (G) and N3 on cytosine (C), and between N4(H) on C and O6 on G. All systems except GC6OH display a concerted mechanism. 8OHGC has the highest dissociation energy and is 1.2 kcal/mol more stable than the nonradical GC base pair. The origin of the interactions are investigated through the estimation of intrinsic acid-basic properties of the *OH-X monomer (X = G or C). Solvent effects play a significant role in reducing the dissociation energy. The reactions including *OH-C adducts have significantly lower PT barriers than both the nonradical GC pair and the *OH-G adducts. All reactions are endothermic, with the GC6OH --> GC6OHPT reaction has the lowest reaction energy (4.6 kcal/mol). In accordance with earlier results, the estimated NBO charges show that the G moiety carries a slight negative charge (and C a corresponding positive one) in each adduct. The formation of a partial ion pair may be a potential factor leading to the PT reactions being thermodynamically unfavored.