The effects of conformation on the acidity of ascorbic acid: a density functional study

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.

Calculating the Aqueous pKa of Phenols: Predictions for Antioxidants and Cannabinoids

We aim to develop a theoretical methodology for the accurate aqueous pK_a prediction of structurally complex phenolic antioxidants and cannabinoids. In this study, five functionals (M06-2X, B3LYP, BHandHLYP, PBE0, and TPSS) and two solvent models (SMD and PCM) were combined with the 6-311++G(d,p) basis set to predict pK_a values for twenty structurally simple phenols. None of the direct calculations produced good results. However, the correlations between the calculated Gibbs energy difference of each acid and its conjugate base, ΔGaq(BA)°=ΔGaqA-°-ΔGaq(HA)°, and the experimental aqueous pK_a values had superior predictive accuracy, which was also tested relative to an independent set of ten molecules of which six were structurally complex phenols. New correlations were built with twenty-seven phenols (including the phenols with experimental pK_a values from the test set), which were used to make predictions. The best correlation equations used the PCM method and produced mean absolute errors of 0.26-0.27 pK_a units and R² values of 0.957-0.960. The average range of predictions for the potential antioxidants (cannabinoids) was 0.15 (0.25) pK_a units, which indicates good agreement between our methodologies. The new correlation equations could be used to make pK_a predictions for other phenols in water and potentially in other solvents where they might be more soluble.

A theoretical study of ascorbic acid oxidation and HOO˙/O2˙- radical scavenging

Ascorbic acid is a well-known antioxidant and radical scavenger. It can be oxidized by losing two protons and two electrons, but normally loses only one electron at a time. The reactivity of the ascorbate radical is unusual, in that it can either disproportionate or react with other radicals, but it reacts poorly with non-radical species. To explore the oxidation mechanism of ascorbic acid, the pK_a's and reduction potentials have been calculated using the B3LYP/6-31+G(d,p) and CBS-QB3 levels of theory with the SMD implicit solvent model and explicit waters. Calculations show that the most stable form of dehydroascorbic acid in water is the bicyclic hydrated structure, in agreement with NMR studies. The possible oxidation reactions at different pH conditions can be understood by constructing a potential-pH (Pourbaix) diagram from the calculated pK_a's and standard reduction potentials. At physiological pH disproportionation of the intermediate radical is thermodynamically favored. The calculations show that disproportionation proceeds via dimerization of ascorbate radical and internal electron transfer, as suggested by Bielski. In the dimer, one of the ascorbate units cyclizes. Then protonation and dissociation yields the fully reduced and bicyclic fully oxidized structures. Calculations show that this mechanism also explains the reaction of the ascorbic acid radical with other radical species such as superoxide. Ascorbate radical combines with the radical, and intramolecular electron transfer followed by cyclization and hydrolysis yields dehydroascorbic acid and converts the radical to its reduced form.

Frontier orbitals and transition states in the oxidation and degradation of L-ascorbic acid: a DFT study

DFT calculations were carried out to investigate reaction paths of L-ascorbic acid (AAH2), hydroxyl radicals and water clusters. Frontier-orbital analyses were also performed to examine the regioselectivity of the OH˙ addition. Transition states of the electrolytic dissociation of AAH2 and intermediate carboxylic acids were found to have very small activation energies through proton transfers along hydrogen bonds. The ionized species (anions) are subject to the electrophilic attack of OH˙. The elementary processes of AAH2 → A˙(-) → dehydroascorbic acid → diketogulonic acid → threonic, oxalic, xylonic and lyxonic acids were investigated and discussed. The processes involved in the conversion of dehydroascorbic acid into a bicyclic hemiketal were also examined as a side-chain participating reaction. The oxidation and degradation of vitamin C up to threonic acid were described mainly as a donor (AAH2)-acceptor (OH˙) reaction.

Ab initio and DFT study of hydrogen bond interactions between ascorbic acid and dimethylsulfoxide based on FT-IR and FT-Raman spectra

The hydrogen bonding of 1:1 complexes formed between l-ascorbic acid (LAA) and dimethylsulfoxide (DMSO) has been studied by means of ab initio and density functional theory (DFT) calculations. Solutions of l-ascorbic acid (AA) in dimethylsulfoxide (DMSO) have been studied by means of both FT-IR (4000-220 cm(-1)) and FT-Raman spectroscopy. Ab initio Hartree-Fock (HF) and DFT methods have been used to determine the structure and energies of stable conformers of various types of L-AA/DMSO complexes in gas phase and solution. The basis sets 6-31++G(**) and 6-311+G(*) were used to describe the structure, energy, charges and vibrational frequencies of interacting complexes in the gas phase. The optimized geometric parameters and interaction energies for various complexes at different theories have been estimated. Binding energies have been corrected for basis set superposition error (BSSE) and harmonic vibrational frequencies of the structures have been calculated to obtain the stable forms of the complexes. The self-consistent reaction field (SCRF) has been used to calculate the effect of DMSO as the solvent on the geometry, energy and charges of complexes. The solvent effect has been studied using the Onsager models. It is shown that the polarity of the solvent plays an important role on the structures and relative stabilities of different complexes. The results obtained show that there is a satisfactory correlation between experimental and theoretical predictions.

Understanding the antioxidant behavior of some vitamin molecules: a first-principles density functional approach

The structures, energetics, vertical and adiabatic ionization potentials, electron affinities, and global reactivity descriptors of antioxidant vitamins (both water- and fat-soluble) in neutral, positively charged, and negatively charged states were investigated theoretically. We worked within the framework of first-principles density functional theory (DFT), using the B3LYP functional and both localized (6-311G+(d,p) and plane-wave basis sets. Solvent effects were modeled via the polarizable continuum model (PCM), using the integral equation formalism variant (IEFPCM). From the computed structural parameters, ionization potentials, electron affinities, and spin densities, we deduced that these vitamins prefer to lose electrons to neutral reactive oxygen species (·OH and ·OOH), making them good antioxidants. Conceptual DFT was used to determine global chemical reactivity parameters. The computed chemical hardnesses showed that these antioxidant vitamins are more reactive than neutral reactive oxygen species (ROS), thus supporting their antioxidant character towards these species. However, in the neutral state, these vitamins do not act as antioxidants for [Formula: see text]. The reactivity of vitamins towards ROS depends on the nature of the solvent. Amongst the ROS, ·OH has the greatest propensity to attract electrons from a generic donor. The reactivities of fat-soluble vitamins towards neutral reactive oxygen species were found to be larger than those of water-soluble vitamins towards these species, showing that the former are better antioxidants.