Development of an ONIOM-G3B3 method to accurately predict C-H and N-H bond dissociation enthalpies of ribonucleosides and deoxyribonucleosides

The Contribution of Theoretical Prediction Studies to the Antioxidant Activity Assessment of the Bioactive Secoiridoids Encountered in Olive Tree Products and By-Products

Assessment of the antioxidant activity of different types of natural compounds is a complex research area that encompasses various in vitro tests and in vivo studies. Sophisticated modern analytical tools permit the unambiguous characterization of the compounds present in a matrix. The contemporary researcher, knowing the chemical structure of the compounds present, can carry out quantum chemical calculations that provide important physicochemical information assisting the prediction of antioxidant potential and the mechanism behind the activity of target compounds before further experimentation. The efficiency of calculations is steadily improved due to the rapid evolution of both hardware and software. It is possible, consequently, to study compounds of medium or even larger size, incorporating also models that simulate the liquid phase (solution). This review contributes to the establishment of theoretical calculations as an inherent part of the antioxidant activity assessment process, having as a case study the complex mixtures of olive bioactive secoiridoids (oleuropein, ligstroside, and related compounds). The literature indicates great variability in theoretical approaches and models used so far for only a limited number of this group of phenolic compounds. Proposals are made for standardization of methodology (reference compounds, DFT functional, basis set size, and solvation model) to facilitate comparisons and communication of findings.

Electron-Induced Repair of 2'-Deoxyribose Sugar Radicals in DNA: A Density Functional Theory (DFT) Study

In this work, we used ωB97XD density functional and 6-31++G** basis set to study the structure, electron affinity, populations via Boltzmann distribution, and one-electron reduction potentials (E°) of 2′-deoxyribose sugar radicals in aqueous phase by considering 2′-deoxyguanosine and 2′-deoxythymidine as a model of DNA. The calculation predicted the relative stability of sugar radicals in the order C4′^• > C1′^• > C5′^• > C3′^• > C2′^•. The Boltzmann distribution populations based on the relative stability of the sugar radicals were not those found for ionizing radiation or OH-radical attack and are good evidence the kinetic mechanisms of the processes drive the products formed. The adiabatic electron affinities of these sugar radicals were in the range 2.6–3.3 eV which is higher than the canonical DNA bases. The sugar radicals reduction potentials (E°) without protonation (−1.8 to −1.2 V) were also significantly higher than the bases. Thus the sugar radicals will be far more readily reduced by solvated electrons than the DNA bases. In the aqueous phase, these one-electron reduced sugar radicals (anions) are protonated from solvent and thus are efficiently repaired via the “electron-induced proton transfer mechanism”. The calculation shows that, in comparison to efficient repair of sugar radicals by the electron-induced proton transfer mechanism, the repair of the cyclopurine lesion, 5′,8-cyclo-2′-dG, would involve a substantial barrier.

Relationship Investigation between C(sp2)-X and C(sp3)-X Bond Energies Based on Substituted Benzene and Methane

The C-X bonds of organic compounds between group X and a saturated or unsaturated carbon atom differ in bond energy. To identify the causes of variation is of great significance in terms of bond nature understanding and bond energy estimation. In this paper, the electronegativity χ[X] of group X was calculated by the "valence electron equalized electronegativity" method. Then, χ[X] and the electronic effect constant of the substituent were taken as variables to establish equations for quantitative correlation between C(sp³)-X and C(sp²)-X for the calculation of C-X bond energies. The aim is make comparison between substituted methane, Me-X, and substituted benzene, Ph-X, as well as that between Me-X and substituted ethylene, C₂H₃-X. We conducted calculation over 40 compounds that contain different X groups, and the results reveal that the C(sp³)-X and C(sp²)-X bond energies are under the influence of a number of factors. In addition to the covalent properties of C and X atoms and χ[X], the bond energies of C(sp²)-X (i.e., D[C(sp²)-X]) are under the influence of the field/inductive effect (σ_F[X]) and conjugated effect (σ_R[X]) of group X, with the former causing a decrease while the latter an increase of D[C(sp²)-X]. Using the acquired quantitative correlation equations and on the basis of a relatively rich set of measured D[Me-X] data, we estimated D[Ph-X] of Ph-X and D[C₂H₃-X] of C₂H₃-X, and the estimation accuracy is within experimental uncertainty. Employing the above method, the D[C(sp²)-X] of 33 substituted benzenes, 53 substituted ethenes, and 82 α-substituted naphthalenes was estimated with satisfactory outcomes.

Enhanced reactivity of the pyrimidine peroxyl radical towards the C–H bond in duplex DNA – a theoretical study

The peroxyl radical exhibits a much stronger reactivity towards C1′–H1′ in duplex DNA with respect to single-stranded DNA.

A new understanding towards the reactivity of DNA peroxy radicals

New radical reaction paths for the formation of DNA oxidation products are suggested, which are strongly different from the previously suggested paths with the tetraoxide intermediate.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

Predicting pKa in Implicit Solvents: Current Status and Future Directions

Computational prediction of condensed phase acidity is a topic of much interest in the field today. We introduce the methods available for predicting gas phase acidity and pKas in aqueous and non-aqueous solvents including high-level electronic structure methods, empirical linear free energy relationships (LFERs), implicit solvent methods, explicit solvent statistical free energy methods, and hybrid implicit–explicit approaches. The focus of this paper is on implicit solvent methods, and we review recent developments including new electronic structure methods, cluster-continuum schemes for calculating ionic solvation free energies, as well as address issues relating to the choice of proton solvation free energy to use with implicit solvation models, and whether thermodynamic cycles are necessary for the computation of pKas. A comparison of the scope and accuracy of implicit solvent methods with ab initio molecular dynamics free energy methods is also presented. The present status of the theory and future directions are outlined.

One-electron oxidation of neutral sugar radicals of 2'-deoxyguanosine and 2'-deoxythymidine: a density functional theory (DFT) study

One electron oxidation of neutral sugar radicals has recently been suggested to lead to important intermediates in the DNA damage process culminating in DNA strand breaks. In this work, we investigate sugar radicals in a DNA model system to understand the energetics of sugar radical formation and oxidation. The geometries of neutral sugar radicals C(1')(•), C(2')(•), C(3')(•), C(4')(•), and C(5')(•) of 2'-deoxyguanosine (dG) and 2'-deoxythymidine (dT) were optimized in the gas phase and in solution using the B3LYP and ωB97x functionals and 6-31++G(d) basis set. Their corresponding cations (C(1')(+), C(2')(+), C(3')(+), C(4')(+), and C(5')(+)) were generated by removing an electron (one-electron oxidation) from the neutral sugar radicals, and their geometries were also optimized using the same methods and basis set. The calculation predicts the relative stabilities of the neutral sugar radicals in the order C(1')(•) > C(4')(•) > C(5')(•) > C(3')(•) > C(2')(•), respectively. Of the neutral sugar radicals, C(1')(•) has the lowest vertical ionization potential (IP(vert)), ca. 6.33 eV in the gas phase and 4.71 eV in solution. C(2')(•) has the highest IP(vert), ca. 8.02 eV, in the gas phase, and the resultant C(2') cation is predicted to undergo a barrierless hydride transfer from the C(1') site to produce the C(1') cation. One electron oxidation of C(2')(•) in dG is predicted to result in a low lying triplet state consisting of G(+•) and C(2')(•). The 5',8-cyclo-2'-deoxyguanosin-7-yl radical formed by intramolecular bonding between C(5')(•) and C(8) of guanine transfers spin density from C(5') site to guanine, and this structure has IP(vert) 6.25 and 5.48 eV in the gas phase and in solution.

THEORETICAL STUDY ON OXYGEN–OXYGEN HOMOLYTIC BOND DISSOCIATION ENTHALPIES OF PEROXIDES

Seven theoretical methods (B3LYP, B3P86, X3LYP, BMK, MP2, CBS-QB3, G3B3) were employed to calculate oxygen–oxygen homolytic bond dissociation enthalpies of 13 peroxides. Comparison of the performances of these methods with the available experimental data showed that ROBMK method is suitable and economical for calculating O–O homolytic bond dissociation enthalpies of peroxides. Then ROBMK method was used to calculate O–O homolytic bond dissociation enthalpies of a series of peroxides. Meanwhile, the α-substituent effect and Hammett relationship were also discussed.

Quantum-chemical predictions of redox potentials of organic anions in dimethyl sulfoxide and reevaluation of bond dissociation enthalpies measured by the electrochemical methods

A first-principle theoretical protocol was developed that could predict the absolute pK(a) values of over 250 structurally unrelated compounds in DMSO with a precision of 1.4 pK(a) units. On this basis we developed the first theoretical protocol that could predict the standard redox potentials of over 250 structurally unrelated organic anions in DMSO with a precision of 0.11 V. Using the two new protocols we systematically reevaluated the bond dissociation enthalpies (BDEs) measured previously by the electrochemical methods. It was confirmed that for most compounds the empirical equation (BDE = 1.37 pK(HA) + 23.1E(o) + constant) was valid. The constant in this equation was determined to be 74.0 kcal/mol, compared to 73.3 kcal/mol previously reported. Nevertheless, for a few compounds the empirical equation could not be used because the solvation energy changed dramatically during the bond cleavage, which resulted from the extraordinary change of dipole moment during the reaction. In addition, we found 40 compounds (mostly oximes and amides) for which the experimental values were questionable by over 5 kcal/mol. Further analyses revealed that all these questionable BDEs could be explained by one of the three following reasons: (1) the experimental pK(a) value is questionable; (2) the experimental redox potential is questionable; (3) the solvent effect cannot be neglected. Thus, by developing practical theoretical methods and utilizing them to solve realistic problems, we hope to demonstrate that ab initio theoretical methods can now be developed to make not only reliable, but also useful, predictions for solution-phase organic chemistry.

Significant effects of phosphorylation on relative stabilities of DNA and RNA sugar radicals: remarkably high susceptibility of h-2' abstraction in RNA

The roles of nucleic acid radicals in DNA and RNA damage cannot be properly understood in the absence of knowledge of the C-H bond strengths depicting the energy cost to generate each of these radicals. However, previous theoretical studies on the relative energies of different nucleic acid radicals are not fully convincing mainly because of the use of oversimplified model compounds. In the present study we chose nucleoside 3',5'-bisphosphates as model compounds for DNA and RNA, in which the effects of both the nucleobase and phosphorylation were taken into consideration. Using the newly developed ONIOM-G3B3 methods, we calculated the gas-phase bond dissociation enthalpies and solution-phase bond dissociation free energies of all the carbohydrate C-H bonds in the model compounds. It was found that the monoanionic phosphate group (OPO3H-) was a better radical stabilization group than the OH group by 1.3 kcal/mol, whereas the neutral phosphate group (OPO3H2) was a significantly worse radical stabilization group than OH by 4.4 kcal/mol. Due to these reasons, the relative thermodynamic susceptibility of H-abstraction from deoxyribonucleotides and ribonucleotides varied considerably depending on the phosphorylation state and the charge carried by the phosphate groups. Strikingly, the bond dissociation free energy of C2'-H in ribonucleotides was dramatically lower than that of all the other C-H bonds by 5-6 kcal/mol regardless of the phosphorylation state and the charge carried by the phosphate group. This explained the previous experimental finding that radiation damage of RNA occurs mainly via H-abstraction at H-2'. A model study suggested that the strength of the hydrogen bonding interaction between the 2'-OH and 3-phosphate groups should dramatically increase from ribonucleoside 3',5'-bisphosphate to its C2' radical. The strengthened hydrogen bonding stabilized the C2' radical, rendering the C2'-H bond of RNA extraordinarily vulnerable to H-abstraction.

Theoretical Prediction of the Hydride Affinities of Various p- and o-Quinones in DMSO

The hydride affinities of 80 various p- and o-quinones in DMSO solution were predicted by using B3LYP/6-311++G (2df,p)//B3LYP/6-31+G* and MP2/6-311++G**//B3LYP/6-31+G* methods, combined with the PCM cluster continuum model for the first time. The results show that the hydride affinity scale of the 80 quinones in DMSO ranges from -47.4 kcal/mol for 9,10-anthraquinone to -124.5 kcal/mol for 3,4,5,6-tetracyano-1,2-quinone. Such a long scale of the hydride affinities (-47.4 to -124.5 kcal/mol) indicates that the 80 quinones can form a large and useful library of organic oxidants, which can provide various organic hydride acceptors that the hydride affinities are known for chemists to choose in organic syntheses. By examining the effect of substituent on the hydride affinities of quinones, it is found that the hydride affinities of quinones in DMSO are linearly dependent on the sum of the Hammett substituent parameters sigma: DeltaGH-(Q) approximately -16.0Sigmasigmai - 70.5 (kcal/mol) for p-quinones and DeltaGH-(Q) approximately -16.2Sigmasigmai - 81.5 (kcal/mol) for o-quinones only if the substituents have no large electrostatic inductive effect and large ortho-effect. Study of the effect of the aromatic properties of quinone on the hydride affinities showed that the larger the aromatic system of quinone is, the smaller the hydride affinity of the quinone is, and the decrease of the hydride affinities is linearly to take place with the increase of the number of benzene rings in the molecule of quinones, from which the hydride affinities of aromatic quinones with multiple benzene rings can be predicted. By comparing the hydride affinities of p-quinones and the corresponding o-quinones, it is found that the hydride affinities of o-quinones are generally larger than those of the corresponding p-quinones by ca. 11 kcal/mol. Analyzing the effect of solvent on the hydride affinities of quinones showed that the effects of solvent (DMSO) on the hydride affinities of quinones are mainly dependent on the electrostatic interaction of the charged hydroquinone anions (QH-) with solvent (DMSO). All the information disclosed in this work should provide some valuable clues to chemists to choose suitable quinones or hydroquinones as efficient hydride acceptors or donors in organic syntheses and to predict the thermodynamics of hydride exchange between quinones and hydroquinones in DMSO solution.

Photo-induced hole transfer from base to sugar in DNA: relationship to primary radiation damage

This work presents the hypothesis that photo-excitation of G.+ in DNA and model systems results in the same electronic states expected from direct ionization of the sugar phosphate backbone and that these states lead to specific sugar radicals on the DNA sugar phosphate backbone. As evidence we show that visible photo-excitation of guanine cation radicals (G.+) in the dinucleoside phosphate TpdG results in high yields (about 85%) of deoxyribose sugar radicals at the C1' and C3' sites. Further, we have calculated transition energies of hole transfer from G.+ in TpdG using time-dependent density functional theory (TD-DFT) at the B3LYP/6-31G(d) level in gas phase as well as in a solvated environment. These calculations clearly predict that visible excitation of G.+ in TpdG causes transitions from only inner-shell filled molecular orbitals (MOs) to the singly occupied molecular orbital (SOMO) that effectively result in hole transfer from guanine either to the sugar phosphate backbone or to the adjacent base, thymine. The hole transfer is followed by rapid deprotonation from the sugar to form C1' and C3' radicals. These experimental and theoretical results are in agreement with our previously published experimental and theoretical results that photo-excitation of G.+ results in high yields of deoxyribose sugar radicals in DNA, guanine deoxyribonucleosides and deoxyribonucleotides. Photo-excitation of G.+ therefore provides a convenient method to produce and study sugar radicals that are expected to be formed in gamma-irradiated DNA systems unencumbered by the many other pathways involved in direct ionization.

C5'- and C3'-sugar radicals produced via photo-excitation of one-electron oxidized adenine in 2'-deoxyadenosine and its derivatives

We report that photo-excitation of one-electron-oxidized adenine [A(-H)•] in dAdo and its 2′-deoxyribonucleotides leads to formation of deoxyribose sugar radicals in remarkably high yields. Illumination of A(-H)• in dAdo, 3′-dAMP and 5′-dAMP in aqueous glasses at 143 K leads to 80-100% conversion to sugar radicals at C5′ and C3′. The position of the phosphate in 5′- and 3′-dAMP is observed to deactivate radical formation at the site of substitution. In addition, the pH has a crucial influence on the site of sugar radical formation; e.g. at pH ∼5, photo-excitation of A(-H)• in dAdo at 143 K produces mainly C5′• whereas only C3′• is observed at high pH ∼12. ¹³C substitution at C5′ in dAdo yields ¹³C anisotropic couplings of (28, 28, 84) G whose isotropic component 46.7 G identifies formation of the near planar C5′•. A β-¹³C 16 G isotropic coupling from C3′• is also found. These results are found to be in accord with theoretically calculated ¹³C couplings at C5′ [DFT, B3LYP, 6-31(G) level] for C5′• and C3′•. Calculations using time-dependent density functional theory [TD-DFT B3LYP, 6-31G(d)] confirm that transitions in the near UV and visible induce hole transfer from the base radical to the sugar group leading to sugar radical formation.

Oxygen Atom Transfer Energetics: Assessment of the Effect of Method and Solvent

Several density functional methods, the semiempirical methods AM1 and PM3, Hartree-Fock, and Gaussian3 theories were applied to compute the oxygen atom transfer enthalpies for 14 X/XO couples (inorganic and organic systems, charged and neutral species, light and heavy main group element containing molecules). The calculated reaction enthalpies were compared to available experimental data. The G3 method alone was found to perform within the experimental error, while the popular B3LYP and BLYP functionals provided inadequate results. Solvent effects were estimated for 19 neutral and anionic X/XO couples by using the conductor-like polarizable continuum model and several cavity models coupled with the B3LYP/6-31++G(2d,2p) level of theory. Surprisingly, the magnitude of the aqueous solvent correction was found to vary significantly for different solute cavity models, occasionally giving larger errors than the gas-phase calculation.