Calculations of pKa Values for a Series of Naturally Occurring Modified Nucleobases.

Modified nucleobases are found in functionally important regions of RNA and are often responsible for essential structural roles. Many of these nucleobase modifications are dynamically regulated in nature, with each modification having a different biological role in RNA. Despite the high abundance of modifications, many of their characteristics are still poorly understood. One important property of a nucleobase is its pKa value, which has been widely studied for unmodified nucleobases, but not for the modified versions. In this study, the pKa values of modified nucleobases were determined by performing ab initio quantum mechanical calculations using a B3LYP density functional with the 6-31+G(d,p) basis set and a combination of implicit-explicit solvation systems. This method, which was previously employed to determine the pKa values of unmodified nucleobases, is applicable to a variety of modified nucleobases. Comparisons of the pKa values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.

Computational Investigation into the Oxidation of Guanine to Form Imidazolone (Iz) and Related Degradation Products.

Imidazolone (Iz) is one of the many products resulting from oxidative damage to DNA. Three pathways for the formation of Iz and related degradation products have been studied by density functional theory using the ωB97XD functional with the 6-31+G(d,p) basis set and SMD implicit water solvation plus a small number of explicit water molecules positioned to help stabilize charged species and facilitate reaction steps. The first pathway starts with guanine radical and the addition of superoxide at C5. Endoperoxide formation was calculated to have slightly lower barriers than diol formation. The next steps are pyrimidine ring opening and decarboxylation. Ring migration then proceeds via an acyclic intermediate rather than a bicyclic intermediate and is followed by formamide loss to yield Iz. The second pathway starts with 8oxoG and proceeds via C5 superoxide addition and diol formation to a relatively stable intermediate, oxidized guanidinohydantoin (Ghox). The barriers for hydroxide ion addition to Ghox are much lower than for water addition and should yield more Iz and parabanic acid at higher pH. The third pathway starts with 8-hydroxy guanine radical formed by hydroxyl radical addition to C8 of guanine or water addition to C8 of guanine radical. Superoxide addition at C5 is followed by diol formation, ring opening and decarboxylation similar to pathways 1 and 2, subsequently leading to Iz formation. The calculated pathways are in good agreement with experimental observations.

Computational Study of the Oxidation of Guanine To Form 5-Carboxyamido-5-formamido-2-iminohydantoin (2Ih).

Oxidative damage to DNA leads to a number of two-electron oxidation products of guanine such as 8-oxo-7,8-dihydroguanine (8oxoG). 5-Carboxyamido-5-formamido-2-iminohydantoin (2Ih) is another two-electron oxidation product that forms in competition with 8oxoG. The pathways for the formation of 2Ih have been studied by density functional theory using the ωB97XD functional with the 6-31+G(d,p) basis set and SMD implicit water solvation plus a small number of explicit water molecules positioned to help stabilize charged species and facilitate reaction steps. For oxidative conditions that produce hydroxyl radical, such as Fenton chemistry, hydroxy radical can add at C4, C5, or C8. Addition at C4 or C5 followed by loss of H2O produces guanine radical. Guanine radical can also be produced directly by oxidation of guanine by reactive oxygen species (ROS). A C5-OH intermediate can be formed by addition of superoxide to C5 of guanine radical followed by reduction. Alternatively, the C5-OH intermediate can be formed by hydroxy radical addition at C5 and oxidation by 3O2. The competition between oxidative and reductive pathways depends on the reaction conditions. Acyl migration of the C5-OH intermediate yields reduced spiroiminodihydantoin (Spred). Subsequent water addition at C8 of Spred and N7-C8 ring opening produces 2Ih. Hydroxy radical addition at C8 can lead to a number of products. Oxidation and tautomerization produces 8oxoG. Alternatively, addition of superoxide at C5 and reduction results in a C5, C8 dihydroxy intermediate. For this species, the low energy pathway to 2Ih is N7-C8 ring opening followed by acyl migration. Ring opening occurs more easily at C8-N9 but leads to a higher energy analogue of 2Ih. Thus, the dominant pathway for the production of 2Ih depends on the nature of the reactive oxygen species and on the presence or absence of reducing agents.

Computational Study of the Formation of C8, C5, and C4 Guanine:Lysine Adducts via Oxidation of Guanine by Sulfate Radical Anion.

Oxidative damage to DNA can lead to DNA-protein cross-links which can interfere with DNA transcription, replication, and repair. In experimental studies modeling oxidative damage to DNA, oxidation of guanosine by sulfate radical anion in the presence of lysine produced a mixture of lysine (Lys)-substituted spiroiminodihydantoins (Sp): ∼65% 5-Lys-Sp, ∼30% 8-Lys-Sp, and ∼5% 5,8-diLys-Sp. Pathways for formation of the lysine adducts during the oxidation of guanine by sulfate radical anions have been mapped out using B3LYP density functional theory and the SMD solvation model. Methylamine was used as a model for lysine, and imidazole served as a proton acceptor. The lowest barrier for methylamine reaction with guanine radical is addition at C8, yielding mainly 8-NHR-Sp and some 5,8-diNR-Sp. This is in good agreement with the cross-link ratios for mild oxidations mediated by type I photosensitizers such as benzophenone, but this is not in agreement with the product ratios for strong oxidants such as sulfate radical anion. The calculations explored pathways for oxidation of guanine by sulfate radical anion that produced guanine radical and radical cation and doubly oxidized guanine (Gox) and its cation. Sulfate radical anion can also oxidize methylamine to produce neutral methylamine radical (CH3NH•) after deprotonation. The calculations qualitatively reproduced the observed product ratio at pH 7 via a pathway involving the barrierless addition of methylamine radical at C5 and C8 of guanine radical. After C5 addition of methylamine radical, the lowest barrier is for H2O addition at C8 leading exclusively to 5-NHR-Sp. After C8 addition of methylamine radical, H2O and methylamine addition to C5 lead to 8-NHR-Sp and some 5,8-diNR-Sp.

Computational Study of the pH-Dependent Competition between Carbonate and Thymine Addition to the Guanine Radical.

When oligonucleotides are oxidized by carbonate radical, thymine and carbonate can add to guanine radical, yielding either a guanine-thymine cross-link product (G∧T) or 8-oxo-7,8-dehydroguanine (8oxoG) and its further oxidation products such as spiroiminodihydantoin (Sp) and guanidinohydantoin (Gh). The ratio of thymine addition to carbonate addition depends strongly on the pH. Details of the mechanism have been explored by density functional calculations using the ωB97XD/6-31+G(d,p) level of theory with the SMD implicit solvation method, augmented with a few explicit waters. Free energies of intermediates and transition states in aqueous solution have been calculated along the pathways for addition of thymine, CO32-/HCO3- and carbonate radical to guanine radical. The pH dependence was examined by using appropriate explicit proton donors/acceptors as computational models for buffers at pH 2.5, 7, and 10. Deprotonation of thymine is required for nucleophilic addition at C8 of guanine radical, and thus is favored at higher pH. The barrier for carbonate radical addition is lower than for bicarbonate or carbonate dianion addition; however, for low concentrations of carbonate radical, the reaction may proceed by addition of bicarbonate/carbonate dianion to guanine radical. Thymine and bicarbonate/carbonate dianion addition are followed by oxidation by O2, loss of a proton from C8 and decarboxylation of the carbonate adduct. At pH 2.5, guanine radical cation can be formed by oxidization with sulfate radical. Water addition to guanine radical cation is the preferred path for forming 8oxoG at pH 2.5.

Synthetic and Computational Study of Tin-Free Reductive Tandem Cyclizations of Neutral Aminyl Radicals.

5- exo, 5- exo Cyclizations of conformationally unbiased propargylic aminyl radicals proceed with excellent yield, chemoselectivity, and diastereoselectivity under tin-free reductive cyclization conditions, regardless of the electronic environments and intermediate radical stabilization resulting from various olefin substituents. These conditions avoid the need for slow addition of initiator and reductant. By contrast, analogous 6- exo, 5- exo cyclizations require substituents capable of intermediate radical stabilization to avoid premature reduction products. These experimental results are corroborated by computations that further establish the reactivity of these aminyl radicals upon exposure to tin-free cyclization conditions.

Investigation into 9(S)-HPODE-derived allene oxide to cyclopentenone cyclization mechanism via diradical oxyallyl intermediates.

The cyclopentane core is ubiquitous among a large number of biologically relevant natural products. Cyclopentenones have been shown to be versatile intermediates for the stereoselective preparation of highly substituted cyclopentane derivatives. Allene oxides are oxygenated fatty acids which are involved in the pathways of cyclopentenone biosynthesis in plants and marine invertebrates; however, their cyclization behavior is not well understood. Recent work by Brash and co-workers (J. Biol. Chem., 2013, 288, 20797) revealed an unusual cyclization property of the 9(S)-HPODE-derived allene oxides: the previously unreported 10Z-isomer cyclizes to a cis-dialkylcyclopentenone in hexane/isopropyl alcohol (100 : 3, v/v), but the known 10E-isomer does not yield cis-cyclopentenone under the same conditions. The mechanism for cyclization has been investigated for unsubstituted and methyl substituted vinyl allene oxide using a variety of methods including CASSCF, ωB97xD, and CCSD(T) and basis sets up to cc-pVTZ. The lowest energy pathway proceeds via homolytic cleavage of the epoxide ring, formation of an oxyallyl diradical, which closes readily to a cyclopropanone intermediate. The cyclopropanone opens to the requisite oxyallyl which closes to the experimentally observed product, cis-cyclopentenone. The calculations show that the open shell, diradical pathway is lower in energy than the closed shell reactions of allene oxide to cyclopropanone, and cyclopropanone to cyclopentenone.

Intramolecular charge-assisted hydrogen bond strength in pseudochair carboxyphosphate.

Carboxyphosphate, a suspected intermediate in ATP-dependent carboxylases, has not been isolated nor observed directly by experiment. Consequently, little is known concerning its structure, stability, and ionization state. Recently, carboxyphosphate as either a monoanion or dianion has been shown computationally to adopt a novel pseudochair conformation featuring an intramolecular charge-assisted hydrogen bond (CAHB). In this work, additive and subtractive correction schemes to the commonly employed open-closed method are used to estimate the strength of the CAHB. Truhlar's Minnesota M06-2X functional with Dunning's aug-cc-pVTZ basis set has been used for geometry optimization, energy evaluation, and frequency analysis. The CHARMM force field has been used to approximate the Pauli repulsive terms in the closed and open forms of carboxyphosphate. From our additive correction scheme, differential Pauli repulsion contributions between the pseudochair (closed) and open conformations of carboxyphosphate are found to be significant in determining the CAHB strength. The additive correction modifies the CAHB prediction (ΔEclosed-open) of -14 kcal/mol for the monoanion and -12 kcal/mol for the dianion to -22.9 and -18.4 kcal/mol, respectively. Results from the subtractive technique reinforce those from our additive procedure, where the predicted CAHB strength ranges from -17.8 to -25.4 kcal/mol for the monoanion and from -15.7 to -20.9 kcal/mol for the dianion. Ultimately, we find that the CAHB in carboxyphosphate meets the criteria for short-strong hydrogen bonds. However, carboxyphosphate has a unique energy profile that does not result in the symmetric double-well behavior of low-barrier hydrogen bonds. These findings provide deeper insight into the pseudochair conformation of carboxyphosphate, and lead to an improved mechanistic understanding of this intermediate in ATP-dependent carboxylases.