Computational estimation of the acidities of purines and indoles

Accurate Computation of Aqueous pKas of Biologically Relevant Organic Acids: Overcoming the Challenges Posed by Multiple Conformers, Tautomeric Equilibria, and Disparate Functional Groups with the Fully Black-Box pK-Yay Method

The use of static electronic structure calculations to compute solution-phase pKas offers a great advantage in that a macroscopic bulk property could be computed via microscopic computations involving very few molecules. There are various sources of errors in the quantum chemical calculations though. Overcoming these errors to accurately compute pKas of a plethora of acids is an active area of research in physical chemistry pursued by both computational as well as experimental chemists. We recently developed the pK-Yay method in our attempt to accurately compute aqueous pKas of strong and weak acids. The method is fully black-box, computationally inexpensive, and is very easy for even a nonexpert to use. However, the method was thus far tested on very few molecules (only 16 in all). Herein, in order to assess the future applicability of pK-Yay, we study the effect of multiple conformers, the presence of tautomers under equilibrium, and the impact of a wide variety of functional groups (derivatives of acetic acid with substituents at various positions, dicarboxylic acids, aromatic carboxylic acids, amines and amides, phenols and thiols, and fluorine bearing organic acids). Starting with more than 1000 conformers and tautomers, this study establishes that overall errors of ∼ 1.0 pKa units are routinely obtained for a majority of the molecules. Larger errors are noted in cases where multiple charges, intramolecular hydrogen bonding, and several ionizable functional groups are simultaneously present. An important conclusion to emerge from this work is that, the computed pKas are insensitive (difference <0.5) to whether we consider multiple conformers/tautomers or only choose the most stable conformer/tautomer. Further, pK-Yay captures the stereoelectronic effects arising due to differing axial vs equatorial pattern, and is useful to predict the dominant acid-base equilibrium in a system featuring several equilibria. Overall, pK-Yay may be employed in several chemical applications featuring organic molecules and biomonomers.

Mechanistic Studies of a Skatole-Forming Glycyl Radical Enzyme Suggest Reaction Initiation via Hydrogen Atom Transfer

Gut microbial decarboxylation of amino acid-derived arylacetates is a chemically challenging enzymatic transformation which generates small molecules that impact host physiology. The glycyl radical enzyme (GRE) indoleacetate decarboxylase from Olsenella uli (Ou IAD) performs the non-oxidative radical decarboxylation of indole-3-acetate (I3A) to yield skatole, a disease-associated metabolite produced in the guts of swine and ruminants. Despite the importance of IAD, our understanding of its mechanism is limited. Here, we characterize the mechanism of Ou IAD, evaluating previously proposed hypotheses of: (1) a Kolbe-type decarboxylation reaction involving an initial 1-e^– oxidation of the carboxylate of I3A or (2) a hydrogen atom abstraction from the α-carbon of I3A to generate an initial carbon-centered radical. Site-directed mutagenesis, kinetic isotope effect experiments, analysis of reactions performed in D₂O, and computational modeling are consistent with a mechanism involving initial hydrogen atom transfer. This finding expands the types of radical mechanisms employed by GRE decarboxylases and non-oxidative decarboxylases, more broadly. Elucidating the mechanism of IAD decarboxylation enhances our understanding of radical enzymes and may inform downstream efforts to modulate this disease-associated metabolism.

Computational Estimation of the Acidities of Pyrimidines and Related Compounds

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.

Computational Estimation of the Aqueous Acidities of Alcohols, Hydrates, and Enols

Alcohols are ubiquitous chemical species, and their acid dissociation constants are critical properties in many applications. It is useful to have a reliable and convenient computational procedure for estimating these acidities. Here we describe a quantitative structure-activity relationship (QSAR) for the pK_as of alcohols in aqueous solution employing density functional theory computations at the B3LYP/6-31+G(d,p) level with the CPCM implicit solvent model for the aqueous solutions. For the pK_as, the energy difference in solution between the parent compound and the anion, ΔE(H₂O), was used as a single descriptor. High quality QSAR equations are obtained both for the gas-phase Gibbs energy changes (ΔG_r⁰s, R² = 0.9764) and the aqueous pK_as (R² = 0.9594). It is shown that the aqueous equation can also be used to estimate the pK_as of fluorinated alcohols, hydrates formed by carbonyl compounds in aqueous solution, and enols that appear as minor components in keto-enol equilibria.

Simultaneous Quantification of Antioxidants Paraxanthine and Caffeine in Human Saliva by Electrochemical Sensing for CYP1A2 Phenotyping

The enzyme CYP1A2 is responsible for the metabolism of numerous antioxidants in the body, including caffeine, which is transformed into paraxanthine, its main primary metabolite. Both molecules are known for their antioxidant and pro-oxidant characteristics, and the paraxanthine-to-caffeine molar ratio is a widely accepted metric for CYP1A2 phenotyping, to optimize dose–response effects in individual patients. We developed a simple, cheap and fast electrochemical based method for the simultaneous quantification of paraxanthine and caffeine in human saliva, by differential pulse voltammetry, using an anodically pretreated glassy carbon electrode. Cyclic voltammetry experiments revealed for the first time that the oxidation of paraxanthine is diffusion controlled with an irreversible peak at ca. +1.24 V (vs. Ag/AgCl) in a 0.1 M H₂SO₄ solution, and that the mechanism occurs via the transfer of two electrons and two protons. The simultaneous quantification of paraxanthine and caffeine was demonstrated in 0.1 M H₂SO₄ and spiked human saliva samples. In the latter case, limits of detection of 2.89 μM for paraxanthine and 5.80 μM for caffeine were obtained, respectively. The sensor is reliable, providing a relative standard deviation within 7% (n = 6). Potential applicability of the sensing platform was demonstrated by running a small scale trial on five healthy volunteers, with simultaneous quantification by differential pulse voltammetry (DPV) of paraxanthine and caffeine in saliva samples collected at 1, 3 and 6 h postdose administration. The results were validated by ultra-high pressure liquid chromatography and shown to have a high correlation factor (r = 0.994).

Purine tautomeric preferences and bond-length alternation in relation with protonation-deprotonation and alkali metal cationization

Quantum chemical calculations were carried out for deprotonated (P-) and protonated purine (PH+) and for adducts with one alkali metal cation (P-M+ and PM+, where M+ is Li+ or Na+) in the gas phase {B3LYP/6-311+G(d,p)}, a model of perfectly apolar environment, and for selected structures in aqueous solution {PCM(water)//B3LYP/6-311+G(d,p)}, a reference polar medium for biological studies. All potential isomers of purine derivatives were considered, the favored structures indicated, and the preferred sites for protonation/deprotonation and cationization reactions determined. Proton and metal cation basicities of purine in the gas phase were discussed and compared with those of imidazole and pyrimidine. Bond-length alternations in the P, PH+, P-M+, and PM+ forms were quantitatively measured using the harmonic oscillator model of electron delocalization (HOMED) indices and compared with those for P. Variations of the HOMED values when proceeding from the purine structural building blocks, pyrimidine and imidazole, to the bicyclic purine system were also examined. Generally, the isolated NH isomers exhibit a strongly delocalized π-system (HOMED > 0.8). Deprotonation slightly increases the HOMED values, whereas protonation and cationization change the HOMED indices in different way. For bidentate M+-adducts, the HOMED values are larger than 0.9 like for the largely delocalized P-. The HOMED values correlate well in a comprehensive relationship with the relative Gibbs energies (ΔG) calculated for individual isomers whatever the purine form is, neutral, protonated, or cationized. When PCM-DFT model was utilized for P-, PH+, PM+, and P-M+ (M+ = Li+) both electron delocalization and relative stability are different from those for the molecules in vacuo. The solvation effects cause a slight increase in HOMEDs, whereas the ΔEs decrease, but in different ways. Hence, contribution of particular isomers in the isomeric mixtures of PH+, PM+, and P-M+ also varies. HOMED variations for the favored neutral, deprotonated, protonated, and lithiated forms of purine in the gas phase and aqueous solution.

Computational Estimation of the Gas-Phase and Aqueous Acidities of Carbon Acids

Carbon acids are compounds which ionize by dissociating along carbon-hydrogen bonds. Although commonly noted for the extremely low acidities of some members of this class, these compounds in fact display a wide range of pK_as and upon proper substitution can even form strong acids. This study employs density functional theory to estimate the gas-phase acidities (ΔG°s) of these compounds and applies a quantitative structure-activity relationship (QSAR) method at the B3LYP/6-31+G** level with the CPCM aqueous solvent model to estimate their aqueous pK_as. For the latter study, the energy difference ΔE(H₂O) in water between the parent compounds and their dissociation products was used as a single parameter. Good quality QSAR regression equations were obtained for both the gas-phase (R² = 0.9905) and aqueous (R² = 0.9647) dissociations. These equations should be useful for the estimation of missing pK_as for compounds in this class. A general discussion of the features affecting the pK_as of these compounds is also given.