pKa Calculation of Some Biologically Important Carbon Acids - An Assessment of Contemporary Theoretical Procedures

Exploring Biginelli-based scaffolds as A2B adenosine receptor antagonists: Unveiling novel structure-activity relationship trends, lead compounds, and potent colorectal anticancer agents

Antagonists of the A2B adenosine receptor have recently emerged as targeted anticancer agents and immune checkpoint inhibitors within the realm of cancer immunotherapy. This study presents a comprehensive evaluation of novel Biginelli-assembled pyrimidine chemotypes, including mono-, bi-, and tricyclic derivatives, as A2BAR antagonists. We conducted a comprehensive examination of the adenosinergic profile (both binding and functional) of a large compound library consisting of 168 compounds. This approach unveiled original lead compounds and enabled the identification of novel structure-activity relationship (SAR) trends, which were supported by extensive computational studies, including quantum mechanical calculations and free energy perturbation (FEP) analysis. In total, 25 molecules showed attractive affinity (Ki < 100 nM) and outstanding selectivity for A2BAR. From these, five molecules corresponding to the new benzothiazole scaffold were below the Ki < 10 nM threshold, in addition to a novel dual A2A/A2B antagonist. The most potent compounds, and the dual antagonist, showed enantiospecific recognition in the A2BAR. Two A2BAR selective antagonists and the dual A2AAR/A2BAR antagonist reported in this study were assessed for their impact on colorectal cancer cell lines. The results revealed a significant and dose-dependent reduction in cell proliferation. Notably, the A2BAR antagonists exhibited remarkable specificity, as they did not impede the proliferation of non-tumoral cell lines. These findings support the efficacy and potential that A2BAR antagonists as valuable candidates for cancer therapy, but also that they can effectively complement strategies involving A2AAR antagonism in the context of immune checkpoint inhibition.

Determination of Amino Acids' pKa: Importance of Cavity Scaling within Implicit Solvation Models and Choice of DFT Functionals

Protonation states of molecules significantly influence the thermodynamics and kinetics of chemical reactions. This is especially important in biochemical processes, where appropriate protonation states of amino acids control the exo/endoergicity of practically all biochemical cycles. This paper is focused on appraisal of the impact of DFT functionals and PCM solvation models on the accuracy of pKa evaluations for all proteinogenic amino acids. Eight functionals (B3LYP, PBE0, revPBE0, M06-2X, M11, M11-L, TPSSh, and ωB97X-D) and four basis sets are considered, together with four kinds of implicit solvation models when additional attention is paid to a cavity construction. An influence of nonelectrostatic contributions and Wertz's corrections on Gibbs free energy is investigated together with accuracy of provided proton solvation energy. The best model is based on the M06-2X/6-311++G**/D-PCM/UAKS computational level. The fitting procedure is utilized to improve the accuracy of the evaluated models. All of these results are also compared with values obtained from the COSMOtherm program and CCSD(T) calculations. Results for cysteine and histidine are discussed individually, as they can be found in different protonation states at neutral pH.

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.

Understanding Electron Transfer Reactions Using Constrained Density Functional Theory: Complications Due to Surface Interactions

The kinetic rates of electrochemical reactions depend on electrodes and molecules in question. In a flow battery, where the electrolyte molecules are charged and discharged on the electrodes, the efficiency of the electron transfer is of crucial importance for the performance of the device. The purpose of this work is to present a systematic atomic-level computational protocol for studying electron transfer between electrolyte and electrode. The computations are done by using constrained density functional theory (CDFT) to ensure that the electron is either on the electrode or in the electrolyte. The ab initio molecular dynamics (AIMD) is used to simulate the movement of the atoms. We use the Marcus theory to predict electron transfer rates and the combined CDFT-AIMD approach to compute the parameters for the Marcus theory where it is needed. We model the electrode with a single layer of graphene and methylviologen, 4,4'-dimethyldiquat, desalted basic red 5, 2-hydroxy-1,4-naphthaquinone, and 1,1-di(2-ethanol)-4,4-bipyridinium were selected for the electrolyte molecules. All of these molecules undergo consecutive electrochemical reactions with one electron being transferred at each stage. Because of significant electrode-molecule interactions, it is not possible to evaluate outer-sphere ET. This theoretical study contributes toward the development of a realistic-level prediction of electron transfer kinetics suitable for energy storage applications.

How to Predict the pK a of Any Compound in Any Solvent

Acid-base properties of molecules in nonaqueous solvents are of critical importance for almost all areas of chemistry. Despite this very high relevance, our knowledge is still mostly limited to the pK _a of rather few compounds in the most common solvents, and a simple yet truly general computational procedure to predict pK _a's of any compound in any solvent is still missing. In this contribution, we describe such a procedure. Our method requires only the experimental pK _a of a reference compound in water and a few standard quantum-chemical calculations. This method is tested through computing the proton solvation energy in 39 solvents and by comparing the pK _a of 142 simple compounds in 12 solvents. Our computations indicate that the method to compute the proton solvation energy is robust with respect to the detailed computational setup and the construction of the solvation model. The unscaled pK _a's computed using an implicit solvation model on the other hand differ significantly from the experimental data. These differences are partly associated with the poor quality of the experimental data and the well-known shortcomings of implicit solvation models. General linear scaling relationships to correct this error are suggested for protic and aprotic media. Using these relationships, the deviations between experiment and computations drop to a level comparable to that observed in water, which highlights the efficiency of our method.

Experimental and DFT investigation on N-functionalized biochars for enhanced removal of Cr(VI)

In this study, N-functionalized biochars with varied structural characteristics were designed by loading poplar leaf with different amounts of urea at 1:1 and 1:3 ratios through pyrolysis method. The addition of urea significantly increased the N content of biochar and facilitated the formation of amine (-NH-, -NH₂), imine (-HCNH), benzimidazole (-C₇H₅N₂), imidazole (-C₃H₃N₂), and pyrimidine (-C₄H₃N₂) groups due to substitution reaction and Maillard reaction. The effect of pH on Cr(VI) removal suggested that decrease in solution pH favored the formation of electrostatic attraction between the protonated functional groups and HCrO₄^-. And, experimental and density functional theory study were used to probe adsorption behaviors and adsorption mechanism which N-functionalized biochars interacted with Cr(VI). The protonation energy calculations indicated that N atoms in newly formed N-containing groups were better proton acceptors. Adsorption kinetics and isotherm experiments exhibited that N-functionalized biochars had greater removal rate and removal capacity for Cr(VI). The removal rate of Cr(VI) on N-functionalized biochar was 10.5-15.5 times that of untreated biochar. Meanwhile, N-functionalized biochar of NB3 with the largest number of adsorption sites for -C₇H₅N₂, -NH₂, -OH, -C₃H₃N₂, and phthalic acid (-C₈H₅O₄) exhibited the supreme adsorption capacity for Cr(VI) through H bonds and the highest adsorption energy was -5.01 kcal/mol. These mechanistic findings on the protonation and adsorption capacity are useful for better understanding the functions of N-functionalized biochars, thereby providing a guide for their use in various environmental applications.

Microscopic mechanism about the selective adsorption of Cr(VI) from salt solution on O-rich and N-rich biochars

The adsorption of Cr(VI) on biochars can be suppressed by coexisting anions, but the roles of O-containing functional groups and in particular N-containing functional groups are unclear. In this study, we combined spectroscopic and molecular simulation approaches to investigate the selective adsorption of Cr(VI) on the O-rich (PB, UB1) and N-rich (UB3, UB5) biochars under strong competition of anions. The elemental analysis and pyrolysis-gas chromatography/mass spectrometry indicated that the structures of PB and UB1 were similar, and so were the UB3 and UB5. Quantification of functional groups showed that for UB1, 75.3% of Cr(VI) removal was attributed to O-containing groups, while 53.3-72.7% of that was mediated by N-containing groups in UB3 and UB5. X-ray photoelectron spectra and density functional theory calculations confirmed that for O-rich biochars, surface complexation and strong H-bonds between carboxyl/hydroxyl and HCrO₄^- improved Cr(VI) removal in the presence of anions, while for N-rich biochars, Cr(VI) adsorption was depressed by coexisting anions in the order of Cl^->NO₃^- >SO₄^2- because of the weaker H-bond between protonated amino groups and HCrO₄^-. This study presents a novel approach for quantitative, molecular-level evaluation of the roles of biochar functional groups in the Cr(VI) removal from complex environmental systems.

Theoretical Determination of pKas of P(1)-H Phospholes and their Comparison with N(1)-H Azoles

pKas of P(1)-H deprotonation of all phosphole group compounds, phosphole (C4H4PH), 1,2 and 1,3-diphospholes (C3H3PPH), 1,2,3, and 1,2,4-triphospholes (C2H2P2PH), tetraphosphole (CHP3PH), and pentaphosphole (P4PH) are determined by DFT calculations. We have compared these to the pKas of the analogue azole group compounds and the factors that cause the differences of these values are discussed in terms of the stabilities of both the phosphole anions and azole anions.

pK-Yay: A Black-Box Method Using Density Functional Theory and Implicit Solvation Models to Compute Aqueous pKa Values of Weak and Strong Acids

The concept of pK_a is very important in chemistry and biology. Over the past two decades, electronic structure calculations have made huge strides toward becoming reliable counterparts to experiments in determining solution-phase pK_a values. However, invariably, the computation of pK_a values involves the use of methods whose error bars are intrinsically larger (definitely >2 kcal/mol with density functionals) than the accuracy desired for the estimation of pK_a values (<1 or 1.5 pK_a units). This scenario presents an ample scope for innovation in developing systematic error cancellation methods even today. In this work, we have developed the pK-Yay method. It is a user-friendly black-box method used to compute aqueous pK_a values of strong and weak acids employing routinely used and computationally inexpensive density functionals and implicit solvation models. It does not require the use of any explicit solvent molecule or modifying any other parameter in an electronic structure program. As part of evaluating the method, a comprehensive test set of 26 weak and strong organic and mineral acids covering 35 pK_a units (20 to -15) was assembled. The detailed description of the method, its performance with different functionals (ωB97X-D performed best for organic acids with a mean absolute error (MAE) of 0.8 pK_a units, and B2-PLYP performed best for strong acids with an MAE of 1.6 pK_a units), the strengths and limitations in the present version, and a future scope to improve the accuracy and reduce its empirical nature are presented herein.

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.

Accurate Estimation of pKb Values for Amino Groups from Surface Electrostatic Potential (VS,min) Calculations: The Isoelectric Points of Amino Acids as a Case Study

Theoretical calculation of equilibrium dissociation constants is a very computationally demanding and time-consuming process since it requires an extremely accurate computation of the solvation free energy changes for each of the species involved. By correlating the minimum surface electrostatic potential (V_S,min) on the nitrogen atom of several aliphatic amino groups-calculated at the density functional theory (DFT) ωB97X-D/cc-pVDZ level of theory-we obtained regression models for each kind of substitution pattern from which we interpolate their corresponding pK_b values with remarkable accuracy: primary R² = 0.9519; secondary R² = 0.9112; and tertiary R² = 0.8172 (N = 20 for each family). These models were validated with tests sets (N = 5) with mean absolute error (MAE) values of 0.1213 (primary), 0.4407 (secondary), and 0.3057 (tertiary). Combining this ansatz with another previously reported by our group to estimate pK_a values [Caballero-García, G.; et al. Molecules 2019, 24(1), 79] we are able to reproduce the isoelectric points of 13 amino acids with no titrable side chains with MAE = 0.4636 pI units.

Absolute and relative pKa predictions via a DFT approach applied to the SAMPL6 blind challenge

In this work, quantum mechanical methods were used to predict the microscopic and macroscopic pK_a values for a set of 24 molecules as a part of the SAMPL6 blind challenge. The SMD solvation model was employed with M06-2X and different basis sets to evaluate three pK_a calculation schemes (direct, vertical, and adiabatic). The adiabatic scheme is the most accurate approach (RMSE = 1.40 pK_a units) and has high correlation (R² = 0.93), with respect to experiment. This approach can be improved by applying a linear correction to yield an RMSE of 0.73 pK_a units. Additionally, we consider including explicit solvent representation and multiple lower-energy conformations to improve the predictions for outliers. Adding three water molecules explicitly can reduce the error by 2-4 pK_a units, with respect to experiment, whereas including multiple local minima conformations does not necessarily improve the pK_a prediction.

Solvation energies of the proton in methanol revisited and temperature effects

Various functionals assessing solvation free energies and enthalpies of the proton in methanol.

Development of constant-pH simulation methods in implicit solvent and applications in biomolecular systems

pH is a critical parameter for biological and technological systems directly related with electrical charges. It can give rise to peculiar electrostatic phenomena, which also makes them more challenging. Due to the quantum nature of the process, involving the forming and breaking of chemical bonds, quantum methods should ideally by employed. Nevertheless, due to the very large number of ionizable sites, different macromolecular conformations, salt conditions, and all other charged species, the CPU time cost simply becomes prohibitive for computer simulations, making this a quite complex problem. Simplified methods based on Monte Carlo sampling have been devised and will be reviewed here, highlighting the updated state-of-the-art of this field, advantages, and limitations of different theoretical protocols for biomolecular systems (proteins and nucleic acids). Following a historical perspective, the discussion will be associated with the applications to protein interactions with other proteins, polyelectrolytes, and nanoparticles.

Control of C-H Bond Activation by Mo-Oxo Complexes: pKa or Bond Dissociation Free Energy (BDFE)?

A density functional theory (DFT) study (BMK/6-31+G(d)) was initiated to investigate the activation of benzylic carbon-hydrogen bonds by a molybdenum-oxo complex with a potentially redox noninnocent supporting ligand-a simple mimic of the active species of the enzyme ethylbenzene dehydrogenase (EBDH)-through deprotonation (C-H bond heterolysis) or hydrogen atom abstraction (C-H bond homolysis) routes. Activation free-energy barriers for neutral and anionic Mo-oxo complexes were high, but lower for anionic complexes than neutral complexes. Interesting trends as a function of substituents were observed that indicated significant H^δ+ character in the transition states (TS), which was further supported by the preference for [2 + 2] addition over HAA for most complexes. Hence, it was hypothesized that C-H activation by these EBDH mimics is controlled more by the pK_a than by the bond dissociation free energy of the C-H bond being activated. Therefore, the results suggest promising pathways for designing more efficient and selective catalysts for hydrocarbon oxidation based on EBDH active-site mimics.

Microhydration of a benzoic acid molecule and its dissociation

C₆H₅COOH requires at least eight H₂O molecules for its dissociation. A linear correlation is observed between pK_a and n_H2O; an unknown pK_a can thus be predicted.

Acidity Constant (pKa ) Calculation of Large Solvated Dye Molecules: Evaluation of Two Advanced Molecular Dynamics Methods

pH‐Sensitive dyes are increasingly applied on polymer substrates for the creation of novel sensor materials. Recently, these dye molecules were modified to form a covalent bond with the polymer host. This had a large influence on the pH‐sensitive properties, in particular on the acidity constant (pK_a). Obtaining molecular control over the factors that influence the pK_a value is mandatory for the future intelligent design of sensor materials. Herein, we show that advanced molecular dynamics (MD) methods have reached the level at which the pK_a values of large solvated dye molecules can be predicted with high accuracy. Two MD methods were used in this work: steered or restrained MD and the insertion/deletion scheme. Both were first calibrated on a set of phenol derivatives and afterwards applied to the dye molecule bromothymol blue. Excellent agreement with experimental values was obtained, which opens perspectives for using these methods for designing dye molecules.

Empirically Fitted Parameters for Calculating pKa Values with Small Deviations from Experiments Using a Simple Computational Strategy

Two empirically fitted parameters have been derived for 74 levels of theory. They allow fast and reliable pKa calculations using only the Gibbs energy difference between an acid and its conjugated base in aqueous solution (ΔGs(BA)). The parameters were obtained by least-squares fits of ΔGs(BA) vs experimental pKa values for phenols, carboxylic acids, and amines using training sets of 20 molecules for each chemical family. Test sets of 10 molecules per family-completely independent from the training set-were used to verify the reliability of the fitting parameters method. It was found that, except for MP2, deviations from experiments are lower than 0.5 pKa units. Moreover, mean unsigned errors lower than 0.35 pKa units were found for the 98.6%, 98.6%, and 94.6% of the tested levels of theory for phenols, carboxylic acids and amines, respectively. The parameters estimated here are expected to facilitate computationally based estimations of pKa values of species for which this magnitude is still unknown, with uncertainties similar to the experimental ones. However, the present study deals only with molecules of modest complexity, thus the reliability of the FP method for more complex systems remains to be tested.

Prediction of pKa values using the PM6 semiempirical method

The PM6 semiempirical method and the dispersion and hydrogen bond-corrected PM6-D3H+ method are used together with the SMD and COSMO continuum solvation models to predict pKa values of pyridines, alcohols, phenols, benzoic acids, carboxylic acids, and phenols using isodesmic reactions and compared to published ab initio results. The pKa values of pyridines, alcohols, phenols, and benzoic acids considered in this study can generally be predicted with PM6 and ab initio methods to within the same overall accuracy, with average mean absolute differences (MADs) of 0.6-0.7 pH units. For carboxylic acids, the accuracy (0.7-1.0 pH units) is also comparable to ab initio results if a single outlier is removed. For primary, secondary, and tertiary amines the accuracy is, respectively, similar (0.5-0.6), slightly worse (0.5-1.0), and worse (1.0-2.5), provided that di- and tri-ethylamine are used as reference molecules for secondary and tertiary amines. When applied to a drug-like molecule where an empirical pKa predictor exhibits a large (4.9 pH unit) error, we find that the errors for PM6-based predictions are roughly the same in magnitude but opposite in sign. As a result, most of the PM6-based methods predict the correct protonation state at physiological pH, while the empirical predictor does not. The computational cost is around 2-5 min per conformer per core processor, making PM6-based pKa prediction computationally efficient enough to be used for high-throughput screening using on the order of 100 core processors.

Calculating Free Energy Changes in Continuum Solvation Models

We recently showed for a large data set of pKas and reduction potentials that free energies calculated directly within the SMD continuum model compares very well with corresponding thermodynamic cycle calculations in both aqueous and organic solvents [ Phys. Chem. Chem. Phys. 2015 , 17 , 2859 ]. In this paper, we significantly expand the scope of our study to examine the suitability of this approach for calculating general solution phase kinetics and thermodynamics, in conjunction with several commonly used solvation models (SMD-M062X, SMD-HF, CPCM-UAKS, and CPCM-UAHF) for a broad range of systems. This includes cluster-continuum schemes for pKa calculations as well as various neutral, radical, and ionic reactions such as enolization, cycloaddition, hydrogen and chlorine atom transfer, and SN2 and E2 reactions. On the basis of this benchmarking study, we conclude that the accuracies of both approaches are generally very similar-the mean errors for Gibbs free energy changes of neutral and ionic reactions are approximately 5 and 25 kJ mol(-1), respectively. In systems where there are significant structural changes due to solvation, as is the case for certain ionic transition states and amino acids, the direct approach generally afford free energy changes that are in better agreement with experiment.

Quantum mechanical based approaches for predicting pKa values of carboxylic acids: evaluating the performance of different strategies

The FP method is recommended for estimating pK_a values of carboxylic acids for which this information is still unknown, especially in combination with the PBE0 functional (MUE = 0.26) and the SMD solvent model.

Structures and spectroscopy of protonated ammonia clusters at different temperatures

Protonated ammonia clusters are all Eigen structures and the first solvation shell of the related ammonium ion core is saturated by four ammonia molecules.

Aqueous acidities of primary benzenesulfonamides: Quantum chemical predictions based on density functional theory and SMD

Aqueous pK(a) of selected primary benzenesulfonamides are predicted in a systematic manner using density functional theory methods and the SMD solvent model together with direct and proton exchange thermodynamic cycles. Some test calculations were also performed using high-level composite CBS-QB3 approach. The direct scheme generally does not yield a satisfactory agreement between calculated and measured acidities due to a severe overestimation of the Gibbs free energy changes of the gas-phase deprotonation reaction by the used exchange-correlation functionals. The relative pK(a) values calculated using proton exchange method compare to experimental data very well in both qualitative and quantitative terms, with a mean absolute error of about 0.4 pK(a) units. To achieve this accuracy, we find it mandatory to perform geometry optimization of the neutral and anionic species in the gas and solution phases separately, because different conformations are stabilized in these two cases. We have attempted to evaluate the effect of the conformer-averaged free energies in the pK(a) predictions, and the general conclusion is that this procedure is highly too costly as compared with the very small improvement we have gained.

Modulation of Phenol Oxidation in Cofacial Dyads

The presentation of two phenols on a xanthene backbone is akin to the tyrosine dyad (Y730 and Y731) of ribonucleotide reductase. X-ray crystallography reveals that the two phenol moieties are cofacially disposed at 4.35 Å. Cyclic voltammetry reveals that phenol oxidation is modulated within the dyad, which exhibits a splitting of one-electron waves with the second oxidation of the phenol dyad occurring at larger positive potential than that of a typical phenol. In contrast, a single phenol appended to a xanthene exhibits a two-electron process, consistent with reported oxidation pathways of phenols in acetonitrile. The perturbation of the phenol potential by stacking is reminiscent of a similar effect for guanines stacked within DNA base pairs.

Theoretically nanoscale study on ionization of muscimol nano drug in aqueous solution

In the present work, acid dissociation constant (pKa) values of muscimol derivatives were calculated using the Density Functional Theory (DFT) method. In this regard, free energy values of neutral, protonated and deprotonated species of muscimol were calculated in water at the B3LYP/6-31G(d) basis sets. The hydrogen bond formation of all species had been analyzed using the Tomasi's method. It was revealed that the theoretically calculated pKa values were in a good agreement with the existing experimental pKa values, which were determined from capillary electrophoresis, potentiometric titration and UV-visible spectrophotometric measurements.

Theoretical study of deuterium isotope effects on acid–base equilibria under ambient and hydrothermal conditions

The calculated difference between pK_a values in H₂O and D₂O is in excellent agreement with experiment.

Are thermodynamic cycles necessary for continuum solvent calculation of pKas and reduction potentials?

Continuum solvent calculations of pKas and reduction potentials usually entail the use of a thermodynamic cycle to express the reaction free energy in terms of gas phase energies and free energies of solvation. In this work, we present a systematic study comparing the solution phase free energy changes obtained in this manner with those directly computed within the SMD solvation model against a large test set of 117 pKas and 42 reduction potentials in water and DMSO. The inclusion of vibrational contributions in the free energy of solvation has a negligible impact on the accuracy of thermodynamic cycle predictions of pKas and reduction potentials. Additionally, when gas phase energies in the thermodynamic cycle are computed at more accurate levels of theory, very similar results (mean unsigned difference of 0.5 kcal mol(-1)) can be achieved when the high-level computations (MP2/GTMP2Large and G3(MP2)-RAD(+)) are directly carried out within the continuum model. Increasing the accuracy of the electronic structure theory may or may not improve the agreement with experiment suggesting that the error is largely in the solvation model. For amino acids where their gas and solution phase species exist as different tautomers, the direct approach provided a significant improvement in calculated pKas. These results demonstrate that direct calculation of solution phase pKas and reduction potentials within the SMD model provides a general and reliable approximation to corresponding thermodynamic cycle based protocols, and is recommended for systems where solvation induced changes in geometry are significant. Further studies are necessary to ascertain whether the results are generalisable to other continuum solvation models.