Next Article in Journal
Preparation and Electrochemical Performance of Three-Dimensional Vertically Aligned Graphene by Unidirectional Freezing Method
Next Article in Special Issue
Scaled in Cartesian Coordinates Ab Initio Molecular Force Fields of DNA Bases: Application to Canonical Pairs
Previous Article in Journal
Structural Identification of Metalloproteomes in Marine Diatoms, an Efficient Algae Model in Toxic Metals Bioremediation
Previous Article in Special Issue
Thermodynamically Stable Cationic Dimers in Carboxyl-Functionalized Ionic Liquids: The Paradoxical Case of “Anti-Electrostatic” Hydrogen Bonding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Estimation of the Acidities of Pyrimidines and Related Compounds †

Department of Chemistry, Wright State University, Dayton, OH 45435, USA
*
Author to whom correspondence should be addressed.
Dedicated to Prof. Frank Weinhold on the occasion of his 80th birthday.
Molecules 2022, 27(2), 385; https://doi.org/10.3390/molecules27020385
Submission received: 25 November 2021 / Revised: 20 December 2021 / Accepted: 1 January 2022 / Published: 7 January 2022

Abstract

:
Pyrimidines are key components in the genetic code of living organisms and the pyrimidine scaffold is also found in many bioactive and medicinal compounds. The acidities of these compounds, as represented by their pKas, are of special interest since they determine the species that will prevail under different pH conditions. Here, a quantum chemical quantitative structure–activity relationship (QSAR) approach was employed to estimate these acidities. Density-functional theory calculations at the B3LYP/6-31+G(d,p) level and the SM8 aqueous solvent model were employed, and the energy difference ∆EH2O between the parent compound and its dissociation product was used as a variation parameter. Excellent estimates for both the cation → neutral (pKa1, R2 = 0.965) and neutral → anion (pKa2, R2 = 0.962) dissociations were obtained. A commercial package from Advanced Chemical Design also yielded excellent results for these acidities.

1. Introduction

Pyrimidines play a central role in the terrestrial genetic code and the pyrimidine framework is also found in many bioactive compounds and medicines. As measurements of molecular properties are frequently difficult or expensive, it is of interest to develop theoretical means for estimating these properties. Among the most important and interesting properties of the pyrimidines are their acidities, as represented by their pKa values, which determine the forms of the compounds that will prevail in solution under different pH conditions. As a result, there has been a long-standing interest in estimating the pKas of chemical compounds using theoretical approaches [1,2]. This interest continues, as demonstrated by the broad range of methods employed in recent pKa studies [3,4,5,6,7,8,9,10]. In an earlier study by our group [11], we presented computational estimates of the pKas of the biologically related purines and indoles. In the present work we develop estimates for the acidities of pyrimidines and related compounds.
Three main approaches have been used in studies estimating compound acidities [1,12]. In the first approach, first-principles or absolute pKas are determined in a straight-forward manner using direct calculations, often relying on a thermodynamic cycle to separate different hypothetical stages [13,14]. The advantage of this approach is that it follows standard procedures and does not depend on prior knowledge of experimental results. Its disadvantage is that it normally requires a high level of computational effort to achieve reasonable accuracy. A second approach employs the development of an appropriate Quantitative Structure–Activity Relationship (QSAR) for the acidities of a set of compounds [2,15]. This approach requires first the collection of measured experimental pKa values for the set and then the identification of some suitable molecular parameter that is closely related to the acidity. An advantage of this approach is that it can generally achieve high accuracy while employing more modest computational efforts, and it accordingly allows the estimation of the pKas of related, unreported compounds using the same modest computational effort. Furthermore, digressions from the derived regression algorithms (i.e., outliers) can alert practitioners to compounds that may be exhibiting behaviors different from the others in the sample [16,17,18]. In a third approach, a number of commercial programs use algorithms derived from large acidity databases and employ empirical parameters such as Hammett constants to estimate pKas [19,20]. The latter two approaches were employed in the present study.
As noted, the QSAR approach relies on the discovery of some parameter or property of the compounds examined that correlates with the activity of interest, in the present case the pKa. In many cases, partial atomic charges have been employed as variation parameters for acidities and other properties. However, because the notion of an “atomic charge” in a molecule is not a proper quantum chemical observable, a variety of different approximate schemes have been developed to represent these atomic charges. Of these proposed schemes, our own group [21,22,23,24] and others [5,25,26,27] have found the natural population analysis (NPA) orbitals and charges developed by Weinhold et al. [28,29,30] to be especially useful in this role. More recently, we have also used another parameter, the energy difference ∆EH2O between the parent compound and its main dissociation product for this purpose. It is this parameter that was employed in the present study.
A particular difficulty arises in a study of the pyrimidines because these compounds do not typically appear as a single species in gas phase or solution, but rather are present as a collection of related tautomers, a situation that also prevailed in our earlier study of purines [11]. Accordingly, some accommodation must be made for this condition, as will be described in the following section.

2. Methods

Our aim in a recent series of reports has been to design and evaluate a QSAR protocol suitable for producing accurate pKa estimates for selected classes of compounds while employing relatively modest and commonly available computational tools.
The initial step in a pKa QSAR study involves the collection of reported experimental pKa values from the literature for the class of compounds of interest. In the present case, values for both the cation-to-neutral dissociation, AH2+ → AH + H+, which we will designate pKa1, and the neutral-to-anion dissociation, AH → A + H+, which we will call pKa2, were available for a number of pyrimidines and related compounds. These values are tabulated in Table 1. Given in Table 1 also are computed pKa values obtained from the Advanced Chemical Development (ACD) [31] commercial software package.
Calculations were carried out using the Spartan’10 software package (Wavefunction, Inc., Irvine, CA). In earlier pKa studies, we found that density-functional computations at the B3LYP/6-31+G(d,p) level provided accurate accounts of molecular properties while still requiring only modest computational demands. After testing this assumption (vide infra) against available gas-phase experimental results for pyrimidines in the NIST database [40], this level of theory was adopted in the present work. For the studies in aqueous solution, the SM8 aqueous solvent model of Marenich et al. [41,42] was used. This same solvent model was also employed in our earlier study of purines and indoles [11] and other studies [16,17,43], where it was shown to perform well in helping to reproduce the experimental data.
As noted above, many of the compounds examined here exist in several tautomeric forms in solution. (For example, uracil has three cationic tautomers, six neutral tautomers, and two anionic tautomers.) Accordingly, the relative stabilities of the tautomeric forms—cationic, neutral, and anionic—of each compound were evaluated within the SM8 aqueous solvent model, and the most stable tautomer for each condition was taken as representative of that compound for computational purposes [44]. Fortunately, there frequently exists a substantial energy gap (>25 kJ/mol) between the most stable tautomer and the remaining tautomers of that species; however, ultimately the validity of this approximation must await justification from the results of the subsequent analysis. In earlier studies, we found that for neutral → anion dissociations (pKa2) the value ∆EH2O = EH2O(A) − EH2O(AH) for the energy difference between the parent compound AH and its dissociation product A- in aqueous solution provides an excellent regression parameter for QSAR pKa studies, and after testing several other parameters this descriptor was used in the present studies. For the cation → neutral dissociation (pKa1), the analogous expression ∆EH2O = EH2O(AH) − EH2O(AH2+) was used.
We also provide a note of caution regarding directly comparing numerical results found here using Spartan’10 with those obtained using the popular Gaussian computational package (Gaussian, Inc., Wallingford, CT 06492, USA), since these programs use different basis sets for some atoms [18,45].

3. Results and Discussion

We first wished to assure that the B3LYP/6-31+G(d,p) level of computation was sufficient to provide accurate results for the dissociations in question. For this, we turned to gas-phase reaction results reported in the NIST chemical database [24]. This database contains gas-phase thermodynamic data for the Gibbs energy change ∆rG° for the anion + H+ → neutral reaction, for six of the compounds studied here. These experimental data are compared with our computed ∆rG° values in Table 2 (note that G°(H+) is −26.3 kJ/mol at 298.15 K [1]). As can be seen, there is very close agreement between the experimental and calculated ΔrG° values. With the exception of succinimide, all of the calculated values fall within the estimated experimental errors given in the NIST database. This, and the coefficient of determination of R2 = 0.998 between the computed and experimental values, suggest that the level of computation described should provide a good account of the reactions to be considered.
As noted above, in previous studies we have found that the energy difference ΔEH2O between the parent compound and its dissociation product provides an excellent regression variable for pKa QSAR estimations. We first examined employment of the gas-phase ΔE values for this purpose. As expected, the use of ΔEgas values yielded good, but not exceptional, correlations for both the pKa1 (cation → neutral, R2 = 0.707) and pKa2 (neutral → anion, R2 = 0.874) dissociations.
We next optimized the compounds within the SM8 aqueous solvent model. Using the solvent-optimized structures and the calculated ∆EH2O values for the appropriate reactions, we obtained the following QSAR models:
pKa1 (calc.) = −0.131 (± 0.008) · ∆EH2O − 151.54 (± 10)
n = 17, R2 = 0.965, s = 1.25, F = 145
pKa2 (calc.) = −0.141 (± 0.008) · ∆EH2O − 159.42 (± 10)
n = 12, R2 = 0.962, s = 0.612, F = 304
In these equations, n = the number of compounds in the sample, R2 is the coefficient of determination (the fraction of the variance in the data accounted for by the model), s is the standard error of the estimate, and F is the Fisher statistic. It is evident that optimization of the structures within the solvent model significantly increases the accuracy of the model, as was also shown in earlier work [1,7,8,9].
The results for pKa1 and pKa2 are plotted in Figure 1 and Figure 2, and the calculated values are compared with the experimental values in Table 3 and Table 4. We note that several of the pKas for the cation → neutral dissociation (pKa1) fall into the difficult-to-measure negative value range and carry large uncertainties. Accordingly, these values are also not well characterized for use in this range in forming the regression Equation (1), and we prefer to recognize this uncertainty by simply indicating modestly negative (<0) or significantly negative (<<0) for the pKas of these compounds.
We also tested the ability of a commercial software program, Advanced Chemical Development, Inc.’s ACD/Labs PhysChem Percepta Suite, to estimate these pKas. The results for pKa1 showed an excellent correlation:
pka1 (exp.) = 1.07 (±0.03) ∙ pka1 (ACD) + 0.34 (±0.2)
where n = 19, R2 = 0.991, s = 0.334, and F = 1583.
The results for pKa2 were also very good:
pKa2 (exp.) = 0.95 (±0.05) ∙ pKa1 (ACD) + 0.92 (±0.52)
where n = 13, R2 = 0.966, s = 0.612, and F = 316.
These results encourage use of this software for studies of the pKas of these compounds.

4. Conclusions

A primary endeavor of scientific studies is to develop models of physical, chemical, and biological systems for the purpose of understanding these systems better. All models are by their very nature approximate. However, as Gauch has noted [46], in some cases—counterintuitively—a model can be more accurate than the data from which it is constructed “because it amplifies hidden patterns and discards unwanted noise” inherent in the system being examined. The QSAR equations used here take advantage of this property by “averaging through” the noise, or random errors, in the experimental pKa data. It is evident that both the QSAR Equations (1) and (2) above provide relatively simple means, via mathematical models, for estimating the pKas of the pyrimidines and related compounds. For example, in order to estimate the pKas of unmeasured compounds in this class or to check the reported pKas of measured compounds, using the QSAR equations one merely needs to determine ∆EH2O for the compound and then evaluate the pKa from the appropriate regression equation. Therefore, the equations provided should allow reasonable estimations for the pKas of other pyrimidines and compounds similar to the pyrimidines. Use of the commercial ACD/Labs program can provide a further independent and very useful check on the pKa estimates.

Author Contributions

Conceptualization, P.G.S.; formal analysis, R.A.H.; investigation, R.A.H.; supervision, P.G.S.; validation, R.A.H.; writing—original draft, P.G.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

See Ref. [44] for available data.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Shields, G.C.; Seybold, P.G. Computational Approaches for the Prediction of pKa Values; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  2. Ho, J. Predicting pKa in Implicit Solvents: Current Status and Future Directions. Aust. J. Chem. 2014, 67, 1441–1460. [Google Scholar] [CrossRef]
  3. Pracht, P.; Grimme, S. Efficient Quantum-Chemical Calculations of Acid Dissociation Constants from Free-Energy Relationships. J. Phys. Chem. A 2021, 125, 5681–5692. [Google Scholar] [CrossRef]
  4. Morency, M.; Néron, S.; Iftimie, R.; Wuest, J.D. Predicting pKa Values of Quinols and Related Aromatic Compounds with Multiple OH Groups. J. Org. Chem. 2021, 86, 14444–14460. [Google Scholar] [CrossRef] [PubMed]
  5. Haslak, Z.P.; Zareb, S.; Dogan, I.; Aviyente, V.; Monard, G. Using Atomic Charges to Describe the pKa of Carboxylic Acids. J. Chem. Inf. Model. 2021, 61, 2733–2743. [Google Scholar] [CrossRef]
  6. Thapa, B.; Raghavachari, K. Accurate pKa Evaluations for Complex Bio-Organic Molecules in Aqueous Media. J. Chem. Theory Comput. 2019, 15, 6025–6035. [Google Scholar] [CrossRef]
  7. Hunt, P.; Hosseini-Gerami, L.; Chrien, T.; Plante, J.; Ponting, D.J.; Segall, M. Predicting pKa Using a Combination of Semi-Empirical Quantum Mechanics and Radial Basis Function Methods. J. Chem. Inf. Model. 2019, 35, 7896–7904. [Google Scholar] [CrossRef] [PubMed]
  8. Dutra, F.R.; de Souza Silva, C.; Custodio, R. On the Accuracy of the Direct Method to Calculate pKa from Electronic Structure Calculations. J. Phys. Chem. A 2021, 125, 65–73. [Google Scholar] [CrossRef]
  9. Pereira, R.W.; Ramabhadran, R.O. pK-Yay: A Black-Box Method Using Density Functional Theory and Implicit Solvation Models to Compute Aqueous pKa Values of Weak and Strong Acids. J. Phys. Chem. A 2020, 124, 9061–9074. [Google Scholar] [CrossRef] [PubMed]
  10. Yang, Q.; Liu, Y.; Zhang, L.; Luo, S.; Cheng, J.-P. Holistic Prediction of the pKa in Diverse Solvents Based on a Machine Learning Approach. Angew. Chem. Int. Ed. 2020, 59, 19282–19291. [Google Scholar] [CrossRef]
  11. Geremia, K.L.; Seybold, P.G. Computational Estimation of the Acidities of Purines and Indoles. J. Mol. Model. 2019, 25, 12. [Google Scholar] [CrossRef]
  12. Seybold, P.G.; Shields, G.C. Computational estimation of pKa values. WIREs Comput. Mol. Sci. 2015, 5, 290–297. [Google Scholar] [CrossRef]
  13. Alongi, K.S.; Shields, G.C. Theoretical Calculations of Acid Dissociation Constants: A Review Article. Ann. Rep. Comput. Chem. 2010, 6, 113–138. [Google Scholar]
  14. Liptak, M.D.; Gross, K.C.; Seybold, P.G.; Feldgus, S.; Shields, G.C. Absolute pKa determinations for substituted phenols. J. Am. Chem. Soc. 2002, 124, 6421–6427. [Google Scholar] [CrossRef]
  15. Seybold, P.G. Quantum Chemical QSPR Estimation of the Acidities and Basicities of Organic Compounds. Adv. Quantum Chem. 2012, 64, 83–104. [Google Scholar]
  16. Seybold, P.G. Quantum Chemical Estimation of the Acidities of Some Inorganic Oxoacids. Mol. Phys. 2015, 113, 232–236. [Google Scholar] [CrossRef]
  17. Seybold, P.G. Computational Estimation of the Acidities of Some Inorganic Nitrogen Acids. Mol. Phys. 2016, 114, 389–393. [Google Scholar] [CrossRef]
  18. Baldasare, C.A.; Seybold, P.G. Computational Estimation of the Gas-Phase and Aqueous Acidities of Carbon Acids. J. Phys. Chem. A 2020, 125, 2152–2159. [Google Scholar] [CrossRef] [PubMed]
  19. Liao, C.; Nicklaus, M.C. Comparison of nine programs predicting pKa values of pharmaceutical substances. J. Chem. Inf. Model. 2009, 49, 2801–2812. [Google Scholar] [CrossRef]
  20. Dearden, J.C.; Rotureau, P.; Fayet, G. QSPR prediction of physico-chemical properties for REACH. SAR QSAR Environ. Res. 2013, 34, 279–318. [Google Scholar] [CrossRef]
  21. Gross, K.C.; Seybold, P.G. Substituent Effects on the Physical Properties and pKa of Aniline. Int. J. Quantum Chem. 2000, 80, 1107–1115. [Google Scholar] [CrossRef]
  22. Gross, K.C.; Seybold, P.G.; Hadad, C.M. Comparison of Different Atomic Charge Schemes for Predicting pKa Variations in Substituted Anilines and Phenols. Int. J. Quantum Chem. 2002, 90, 445–458. [Google Scholar] [CrossRef]
  23. Gross, K.C.; Hadad, C.M.; Seybold, P.G. Charge Competition in Halogenated Hydrocarbons. Int. J. Quantum Chem. 2011, 112, 219–229. [Google Scholar] [CrossRef]
  24. Scheiner, S.; Seybold, P.G. Quantum chemical analysis of the energetics of the anti and gauche conformers of ethanol. Struct. Chem. 2009, 20, 43–48. [Google Scholar] [CrossRef]
  25. Varekova, R.S.; Geidl, S.; Ionescu, C.-M.; Skrehota, O.; Kudera, M.; Sehnal, D.; Bouchal, T.; Abagyan, R.; Huber, H.J.; Koca, J. Predicting pKa Values of Substituted Phenols from Atomic Charges: Comparison of Different Quantum Mechanical Methods and Charge Distribution Schemes. J. Chem. Inf. Model. 2011, 51, 1795–1806. [Google Scholar] [CrossRef]
  26. Ugur, I.; Marion, A.; Parant, S.; Jensen, J.H.; Monard, G. Rationalization of the pKa Values of Alcohols and Thiols Using Atomic Charge Descriptors and Its Application to the Prediction of Amino Acid pKa’s. J. Chem. Inf. Model. 2014, 54, 2200–2213. [Google Scholar] [CrossRef]
  27. Mull, H.F.; Turney, J.M.; Douberly, G.E.; Schaefer, H.F., III. Kinetic Stability of Pentazole. J. Phys. Chem. A 2021, 125, 9092–9098. [Google Scholar] [CrossRef]
  28. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Natural population analysis. J. Chem. Phys. 1985, 83, 735–746. [Google Scholar] [CrossRef]
  29. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Intermolecular Interactions from a Natural Bond Orbital, Donor-Acceptor Viewpoint. Chem. Rev. 1988, 88, 899–926. [Google Scholar] [CrossRef]
  30. Glendening, E.D.; Landis, C.R.; Weinhold, F. Natural Bond Orbital Methods. WIREs Comput. Mol. Sci. 2012, 2, 1–42. [Google Scholar] [CrossRef]
  31. Advanced Chemistry Development, Inc.: Toronto, ON, Canada. Available online: https://www.acdlabs.com/ (accessed on 24 November 2021).
  32. Katritzky, A.R.; Ramsden, C.A.; Joule, J.A.; Zhdankin, V.V. Handbook of Heterocyclic Chemistry, 3nd ed.; Elsevier: New York, NY, USA, 2010. [Google Scholar]
  33. Serjeant, E.P.; Dempsey, B. Ionization Constants of Organic Acids in Aqueous Solution; UPAC Chemical Data Series-No. 23; International Union of Pure and Applied Chemistry: Oxford, UK, 1979. [Google Scholar]
  34. CRC Handbook of Chemistry and Physics, 89th ed.; CRC Press: Boca Raton, FL, USA, 2008.
  35. Dean, J.A. Lange’s Handbook of Chemistry, 13th ed.; McGraw-Hill: New York, NY, USA, 1985. [Google Scholar]
  36. Levene, P.A.; Bass, L.W.; Simms, H.S. The ionization of pyrimidines in relation to the structure of pyrimidine nucleosides. J. Biol. Chem. 1926, 70, 229–241. [Google Scholar] [CrossRef]
  37. Darnall, K.R.; Townsend, L.B.; Robins, R.K. The structure of showdomycin, a novel carbon-linked nucleoside antibiotic related to uridine. Proc. Nat. Acad. Sci. USA 1967, 57, 548–553. [Google Scholar] [CrossRef] [Green Version]
  38. Privat, E.J.; Sowers, L.C. A proposed mechanism for the mutagenicity of 5-formyluracil. Mutat. Res. 1996, 354, 151–156. [Google Scholar] [CrossRef]
  39. Handschmacher, R.E.; Welch, A.D. Microbial studies of 6-azauracil, an antagonist of uracil. Cancer Res. 1956, 16, 965–969. [Google Scholar]
  40. NIST Chemistry WebBook. Available online: http://webbook.nist.gov/chemistry/ (accessed on 13 July 2016).
  41. Marenich, A.V.; Olson, R.M.; Kelly, C.P.; Cramer, C.J.; Truhlar, D.G. Self-consistent reaction field model for aqueous and nonaqueous solutions based on accurate polarized partial charges. J. Chem. Theory Comput. 2007, 3, 2011–2033. [Google Scholar] [CrossRef] [PubMed]
  42. Cramer, C.J.; Truhlar, D.G. A Universal Approach to Solvation Modeling. Accts. Chem. Res. 2008, 41, 760–768. [Google Scholar] [CrossRef] [PubMed]
  43. Hunter, N.E.; Seybold, P.G. Theoretical estimation of the aqueous pKas of thiols. Mol. Phys. 2014, 112, 340–348. [Google Scholar] [CrossRef]
  44. Wessner, R.A. Theoretical Estimation of pKas of Pyrimidines and Related Heterocycles. Master’s Thesis, Wright State University, Dayton, OH, USA, 2016. [Google Scholar]
  45. Scheiner, S.; Seybold, P.G. (Utah State University, Logan, Utah). Unpublished work. 2009.
  46. Gauch, H.G., Jr. Prediction, Parsimony and Noise. Am. Sci. 1993, 81, 468–478. [Google Scholar]
Figure 1. Plot of experimental pKa1s vs. calculated ΔEH2Os.
Figure 1. Plot of experimental pKa1s vs. calculated ΔEH2Os.
Molecules 27 00385 g001
Figure 2. Plot of experimental pKa2s vs. calculated ΔEH2Os.
Figure 2. Plot of experimental pKa2s vs. calculated ΔEH2Os.
Molecules 27 00385 g002
Table 1. Reported experimental pKa1 and pKa2 values and ACD computed values for the pyrimidines and related heterocycles studied.
Table 1. Reported experimental pKa1 and pKa2 values and ACD computed values for the pyrimidines and related heterocycles studied.
No.CompoundFormulapKa1ACD pKa1pKa2ACD pKa2
1azauracilC3H3N3O2-−4.4 ± 0.2-7.8 ± 0.2
2aziridineC2H5N7.98 [32], 8.05 [33]8.1 ± 0.2--
3creatinineC4H7N3O4.8 [34]4.8 ± 0.1--
4cytosineC4H5N3O4.32 [30], 4.58 [35], 4.6 [31]4.4 ± 0.113 [33], 12.15 [36]12.3 ± 0.1
5flucytosineC4H3FN2O23.26 [34]2.6 ± 0.1-10.5 ± 0.1
6imidazoleC3H4N27.15 [33], 6.99 [34,35]7.2 ± 0.614.44 [33]13.9 ± 0.1
71-methylimidazoleC4H6N26.95 [34]7.0 ± 0.1N/A-
84-methylimidazoleC4H6N27.55 [35]7.7 ± 0.6-14.3 ± 0.1
9isocytosineC4H5N3O4.01 [34]3.4 ± 0.59.42 [36]9.6 ± 0.4
10isoxazoleC3H6NO−2.0 [34]−2.0 ± 0.5N/A
11maleimideC4H3NO2-−5.7 ± 0.29.46 [37] 8.5 ± 0.2
12morpholineC4H9NO8.492 [35]9.0 ± 0.2N/A-
13piperidineC5H11N11.12 [35]10.4 ± 0.1--
14piperazineC4H10N29.78 [35], 9.73 [34]9.6 ± 0.1--
151-methylpiperazineC5H12N210.19 [35]9.3 ± 0.1--
16oxazoleC3H6NO0.8 [34]1.0 ± 0.1N/A-
17pyrazineC4H4N20.65 [34]1.2 ± 0.1N/A-
18pyrazoleC3H4N22.61 [35]2.8 ± 0.114.21 [33]14.0 ± 0.5
19pyridazineC4H43N22.3 [32]3.1 ± 0.1N/A-
20pyridineC5H5N5.23 [34]5.2 ± 0.1N/A-
21pyrimidineC4H4N21.3 [32]1.8 ± 0.1N/A-
22pyrroleC4H5N−3.8 [34]−0.3 ± 0.517.0 [33]17.0 ± 0.5
23pyrrolidineC4H9N12.10 [33], 11.31 [34,35]10.5 ± 0.1--
24succinimideC4H5NO2-−4.4 ± 0.29.62 [34,35], 9.68 [36]9.6 ± 0.1
25thymineC5H6N2O2-−4.1 ± 0.49.9 [33], 9.79 [35], 9.44 [34]9.2 ± 0.1
26uracilC4H4N2O2-−4.2 ± 0.19.43 [38], 9.45 [34,36]8.9 ± 0.1
275-bromouracilC4H3BrN2O2--7.91 [38]6.8 ± 0.1
285-chlorouracilC4H3CIN2O2--7.92 [38]6.8 ± 0.1
29fluorouracilC4H4FN3O--8.04 [30], 8.00 [39], 7.93 [38]6.7 ± 0.1
305-formyluracilC5H4N2O3--6.84 [38]7.3 ± 0.1
315-nitrouracilC4H3N3O4--5.3 [38]5.2 ± 0.1
Table 2. Experimental and calculated gas-phase ΔrG° values (kJ/mol) for the reaction A + H+ → AH.
Table 2. Experimental and calculated gas-phase ΔrG° values (kJ/mol) for the reaction A + H+ → AH.
CompoundExp. ΔraCalc. ΔrbCalc. ΔE b
pyridine 1601 16051648
pyrazine160516051643
pyrimidine157715791614
pyridazine156515621601
imidazole143314321466
succinimide141414011436
a From the NIST database, Ref. [40]; b B3LYP/6-31+G(d,p).
Table 3. Literature pKa1s and estimated pKa1s.
Table 3. Literature pKa1s and estimated pKa1s.
CompoundΔE kJ/molExp. pKa1Calc. pKa1 aResidual
azauracil−1084-<<0-
aziridine−12168.017.760.25
creatinine−11884.84.090.71
cytosine−11834.53.431.07
flucytosine−11693.261.601.66
imidazole−12137.077.36−0.29
1-methylimidazole−12167.957.760.19
4-methylimidazole −12157.557.63−0.08
isocytosine−11834.013.430.58
isoxazole−1115−2<0-
maleimide−1035-<<0-
oxazole−8720.8<<0-
piperidine −123911.1210.770.35
piperazine−12269.769.070.69
1-methylpiperazine−123110.199.720.47
pyrazine−11660.651.21−0.56
pyrazole−11742.612.250.36
pyridazine −11862.33.83−1.53
pyridine−12005.235.66−0.43
pyrimidine−11721.31.99−0.69
pyrrolidine−123911.7110.770.94
succinimide−1007-<<0-
a pKa1s were estimated using Equation (1).
Table 4. Literature pKa2s and estimated pKa2s.
Table 4. Literature pKa2s and estimated pKa2s.
CompoundΔE kJ/molExp. pKa2Calc. pKa2 aResiduals
azauracil−1179-6.82-
aziridine−1325-27.41-
creatinine−1217-12.18-
cytosine−121612.5712.040.53
isocytosine−11969.429.220.20
flucytosine−1254-17.39-
imidazole−121614.412.042.36
4-methylimidazole−1224-13.16-
maleimide−11969.59.220.28
piperazine−1313-25.71-
1-methylpiperazine−1395-37.28-
piperidine−1310-25.29-
pyrrole−12431715.841.16
pyrrolidine−1390-36.57-
thymine−11999.719.640.07
uracil−12009.449.78−0.34
5-bromouracil−11847.917.520.39
5-chlorouracil−11827.927.240.68
fluorouracil−11937.998.79−0.80
5-formyluracil−11806.846.96−0.12
5-nitrouracil−11595.34.001.30
a pKa2s were estimated using Equation (2).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Holt, R.A.; Seybold, P.G. Computational Estimation of the Acidities of Pyrimidines and Related Compounds. Molecules 2022, 27, 385. https://doi.org/10.3390/molecules27020385

AMA Style

Holt RA, Seybold PG. Computational Estimation of the Acidities of Pyrimidines and Related Compounds. Molecules. 2022; 27(2):385. https://doi.org/10.3390/molecules27020385

Chicago/Turabian Style

Holt, Rachael A., and Paul G. Seybold. 2022. "Computational Estimation of the Acidities of Pyrimidines and Related Compounds" Molecules 27, no. 2: 385. https://doi.org/10.3390/molecules27020385

APA Style

Holt, R. A., & Seybold, P. G. (2022). Computational Estimation of the Acidities of Pyrimidines and Related Compounds. Molecules, 27(2), 385. https://doi.org/10.3390/molecules27020385

Article Metrics

Back to TopTop