Determination of acidity constants of pyridines,imidazoles, and oximes by capillary electrophoresis

The acidity constant in the form of pK_a is one of the most important physicochemical quantities. There are prediction tools available for calculating the pK_a, but they only deliver precise calculated values for a relatively small set of chemicals. For complex structures with multiple functional groups in particular, the error in the predicted pK_a is high due to the application domain of the corresponding models. Thus, we aim to enlarge the dataset of experimentally determined pK_a values using capillary electrophoresis. We, therefore, selected various pyridines, imidazoles, and oximes to determine the pK_a values using the internal standard approach and the classical method. Especially oximes were not investigated in the past, and predictions for them include larger errors. Thus, our experimentally determined values could contribute to an improved understanding of various functional groups impacting the pK_a values and serve as additional datasets to develop improved pK_a prediction tools.     

Avoiding gas-phase calculations in theoretical p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> predictions

CBS-QB3, two simplified and less computationally demanding versions of CBS-QB3, DFT-B3LYP, and HF quantum chemistry methods have been used in conjunction with the CPCM continuum solvent model to calculate the free energies of proton exchange reactions in water solution following an isodesmic reaction approach. According to our results, the precision of the predicted pK _a values when compared to experiment is equivalent to that of the thermodynamic cycles that combine gas-phase and solution-phase calculations. However, in the aqueous isodesmic reaction schema, the accuracy of the results is less sensitive to the presence of explicit water molecules and to the global charges of the involved species since the free energies of solvation are not required. In addition, this procedure makes easier the prediction of pK _a values for molecules that undergo large conformational changes in solvation process and makes possible the pK _a prediction of unstable species in gas-phase such as some zwitterionic tautomers. The successive pK _a values of few amino acids corresponding to the ionization of the α-carboxylic acid and α-amine groups, which is one of the problematic cases for thermodynamic cycles, were successfully calculated by employing the aqueous isodesmic reaction yielding mean absolute deviations of 0.22 and 0.19 pK _a units for the first and second ionization processes, respectively.     

Theoretical calculation of the pKa values of some drugs in aqueous solution

In this work, calculations of pK_a values have been performed on benzoic acid and its para‐substituted derivatives and some drugs by using Gaussian 98 software package. Gas‐phase energies were calculated with HF/6‐31 G** and B3LYP/6‐31 G** levels of theory. Free energies of solvation have been computed using the polarizable continuum model (PCM), conductor‐like PCM (CPCM), and the integral equation formalism‐PCM at the same levels which have been used for geometry determination in the gas‐phase. The results that show the calculated pK_a values using the B3LYP are better than those using the corresponding HF. In comparison to the other models, the results obtained indicate that the PCM model is a suitable solvation model for calculating pK_a values. For the investigated compounds, a good agreement between the experimental and the calculated pK_a values was also observed. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Microscopic Protonation/Deprotonation Equilibria of the Anti-Inflammatory Agent Piroxicam

The microscopic ionization behavior of piroxicam was investigated using two different approaches, i.e., direct UV spectroscopy and an indirect analogue approach (deductive method). The best microscopic pK_a values (pK_a12 = 4.60, pK_a21 = 5.40, pK_a22 = 2.72, and pK_a11 = 1.92) were obtained by the deductive method using as pK_a22 the pK_a of the enolic O-methylated piroxicam 2 . The results show remarkable electrostatic effects in the protonation/deprotonation equilibria, a marked increase in the acidity of the enolic function (2.68 pK_a units) being caused by the pyridinium group. The electronic structure of piroxicam was studied based on ¹H-NMR chemical shifts at various ionization states, indicating an extended electron conjugation through the molecule. The partition measurements in octan-1-ol/H₂O of zwitterionic compound 3 (the pyridyl N-methyl derivative of piroxicam ( 1 )) suggest that the two opposite charges in zwitterionic piroxicam are indeed in a close intramolecular proximity.     

Effect of charge distribution on electrostatic chromophore—protein interactions in Bacteriorhodopsin

Charge distributions of a protonated and unprotonated Schiff base model compound are determined using different quantum chemical methods. After fitting the model molecule onto the protonated retinal Schiff base in Bacteriorhodopsin, electrostatic interaction energies between the model molecule and protein are calculated. Interaction energies as well as the calculated pK_1/2 values of the model molecule are shown to depend considerably on the chosen charge distribution. Electrostatic potential derived partial charges determined at different ab initio levels reveal interaction energies between the model molecule and nearby residues such as ARG-82, ASP-85, and ASP-212, which are relatively method independent. Consequently, such charge distributions also result in pK_1/2 values for the model molecule that are very similar. Larger deviations in the electrostatic interaction energies, however, are found in the case of charge distributions derived according to the Mulliken population analysis. Nevertheless, some sets of Mulliken derived partial charges predicted pK_1/2 values for the model molecule that are close to those determined with electrostatic potential derived partial charges. This agreement, however, is only achieved because the individual errors of the contributing terms are approximately compensated. The use of the extended atom model is shown to be problematic. Although potential derived charges can correctly describe electrostatic interaction energies, they fail to predict pK_1/2 values. On the basis of the present investigation a new set of partial charges for the protonated and unprotonated retinal Schiff base is proposed to be used in molecular dynamics simulations and electrostatics calculations. © 1997 by John Wiley & Sons, Inc.     

Importance of electrostatic polarizability in calculating cysteine acidity constants and copper(I) binding energy of Bacillus subtilis CopZ

CopZ is a copper chaperone from Bacillus subtilis. It is an important part of Cu(I) trafficking. We have calculated pK_a values for the CXXC motif of this protein, which is responsible for the Cu(I) binding, and the Cu(I) binding constants. Polarizable and fixed‐charges formalisms were used, and solvation parameters for the both models have been refitted. We had to partially redevelop parameters for the protonated and deprotonated cysteine residues. We have discovered that the polarizable force field (PFF) is qualitatively superior and allows a uniformly better level of energetic results. The PFF pK_a values for cysteine are within about 0.8–2.8 pH units of the experimental data, while the fixed‐charges OPLS formalism yields errors of up to tens of units. The PFF magnitude of the copper binding energy is about 10 kcal/mol or 50% higher than the experimental value, while the using the refitted OPLS parameters leads to an overall positive binding energy, thus predicting no thermodynamically stable complex. At the same time, the agreement of the polarizable S···Cu(I) distances with the experimental results is within 0.08 Å range, and the nonpolarizable calculations lead to an error of about 0.4 Å. Moreover, the accuracy of the PFF has been achieved without any explicit fitting to either pK_a or CopZ···Cu(I) binding energies. We believe that this makes our polarizable technique a choice method in reproducing protein—copper binding and further supports the notion that explicit treatment of electrostatic polarization is crucial in many biologically relevant studies, especially ion binding and transport. © 2012 Wiley Periodicals, Inc.     

pKa prediction for acidic phosphorus‐containing compounds using multiple linear regression with computational descriptors

Ninety‐six acidic phosphorus‐containing molecules with pK_a 1.88 to 6.26 were collected and divided into training and test sets by random sampling. Structural parameters were obtained by density functional theory calculation of the molecules. The relationship between the experimental pK_a values and structural parameters was obtained by multiple linear regression fitting for the training set, and tested with the test set; the R² values were 0.974 and 0.966 for the training and test sets, respectively. This regression equation, which quantitatively describes the influence of structural parameters on pK_a, and can be used to predict pK_a values of similar structures, is significant for the design of new acidic phosphorus‐containing extractants. © 2016 Wiley Periodicals, Inc.     

Accurate pKa Calculation of the Conjugate Acids of Alkanolamines,Alkaloids and Nucleotide Bases by Quantum Chemical Methods

The pK_a of the conjugate acids of alkanolamines, neurotransmitters, alkaloid drugs and nucleotide bases are calculated with density functional methods (B3LYP, M08‐HX and M11‐L) and ab initio methods (SCS‐MP2, G3). Implicit solvent effects are included with a conductor‐like polarizable continuum model (CPCM) and universal solvation models (SMD, SM8). G3, SCS‐MP2 and M11‐L methods coupled with SMD and SM8 solvation models perform well for alkanolamines with mean unsigned errors below 0.20 pK_a units, in all cases. Extending this method to the pK_a calculation of 35 nitrogen‐containing compounds spanning 12 pK_a units showed an excellent correlation between experimental and computational pK_a values of these 35 amines with the computationally low‐cost SM8/M11‐L density functional approach.     

Computation of the influence of chemical substitution on the pK a of pyridine using semiempirical and ab initio methods

The physical properties of chemicals are strongly influenced by their protonation state, affecting, for example, solubility or hydrogen-bonding characteristics. The ability to accurately calculate protonation states in the form of pK _as is, therefore, desirable. Calculations of pK _a changes in a series of substituted pyridines are presented. Computations were performed using both ab initio and semiempirical approaches, including free energies of solvation via reaction-field models. The selected methods are readily accessible with respect to both software and computational feasibility. Comparison of calculated and experimental pK _as shows the experimental trends to be reasonably reproduced by the computations with root-mean-square differences ranging from 1.22 to 4.14 pK _a units. Of the theoretical methods applied the best agreement occurred using the second-order M?ller–Plesset/6-31G(d)/isodensity surface polarized continuum solvation model, while the more computationally accessible Austin model 1/Solvent model 2 (SM2) approach yielded results similar to the ab initio methods. Analysis of component contributions to the calculated pK _as indicates the largest source of error to be associated with the free energies of solvation of the protonated species followed by the gas-phase protonation energies; while the latter may be improved via the use of higher levels of theory, enhancements in the former require improvements in the solvation models. The inclusion of alternate minimum in the computation of pK _as is also indicated to contribute to differences between experimental and calculated pK _a values. Received: 27 April 1999 / Accepted: 27 July 1999 / Published online: 2 November 1999     

Thermodynamic Estimate of pKa Values of the Carboxylic Acids in Aqueous Solution with the Density Functional Theory

The 96 pK_a values of 85 carboxylic acids in aqueous solution were calculated with the density functional theory method at the level of B3LYP/6‐31+G(d,p) and the polarizable continuum model (PCM) was used to describe the solvent. In the calculations of pK_a values, the dissociation Gibbs free energies were directly calculated using carboxylic acid dissociation reactions in aqueous solution, i. e., no thermodynamic cycle was employed, which is different from the previous literatures. A highly significant correlation of R²=0.95 with a standard deviation (SD) of 0.36 between the experimental pK_a values and the calculated dissociation Gibbs free energies [ΔG(calc.)] was found. The slope of pK_a vs. (G(calc.)/(20303RT) is only 47.6% of the theoretically expected value, which implies that the ΔG(calc.) value from the theoretical calculation is larger than the actual one for all 85 carboxylic acids studied. Thus, by adding the 0.476 scaling‐factor into the slope, we can derive a reliably procedure that can reproduce the experimental pK_a values of carboxylic acids. The pK_a values furnished by this procedure are in good agreement with the experimental results for carboxylic acids in aqueous solution.     

SAMPL6: calculation of macroscopic p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> values from ab initio quantum mechanical free energies

Macroscopic pK_a values were calculated for all compounds in the SAMPL6 blind prediction challenge, based on quantum chemical calculations with a continuum solvation model and a linear correction derived from a small training set. Microscopic pK_a values were derived from the gas-phase free energy difference between protonated and deprotonated forms together with the Conductor-like Polarizable Continuum Solvation Model and the experimental solvation free energy of the proton. pH-dependent microstate free energies were obtained from the microscopic pK_as with a maximum likelihood estimator and appropriately summed to yield macroscopic pK_a values or microstate populations as function of pH. We assessed the accuracy of three approaches to calculate the microscopic pK_as: direct use of the quantum mechanical free energy differences and correction of the direct values for short-comings in the QM solvation model with two different linear models that we independently derived from a small training set of 38 compounds with known pK_a. The predictions that were corrected with the linear models had much better accuracy [root-mean-square error (RMSE) 2.04 and 1.95 pK_a units] than the direct calculation (RMSE 3.74). Statistical measures indicate that some systematic errors remain, likely due to differences in the SAMPL6 data set and the small training set with respect to their interactions with water. Overall, the current approach provides a viable physics-based route to estimate macroscopic pK_a values for novel compounds with reasonable accuracy.     

Negative influence of pKa on activation energy barrier: A case study for double proton transfer reaction in inorganic acid dimers

Strength of acid can be determined by means of pK_a value. Attempts have been made to find a relationship between pK_a and activation energy barrier for a double proton transfer (DPT) reaction in inorganic acid dimers. Negative influence of pK_a is observed on activation energy (E_a) which is contrary to the general convention of pK_a. Four different levels of theories with two different basis sets have been used to calculate the activation energy barrier of the DPT reaction in inorganic acid dimers. A model based on first and second order polynomial has been created to find the relationship between activation energy for DPT reaction. © 2018 Wiley Periodicals, Inc.     

Deprotometalation-Iodolysis and Direct Iodination of 1-Arylated 7-Azaindoles: Reactivity Studies and Molecule Properties

Five protocols were first compared for the copper-catalyzed C-N bond formation between 7-azaindole and aryl/heteroaryl iodides/bromides. The 1-arylated 7-azaindoles thus obtained were subjected to deprotometalation-iodolysis sequences using lithium 2,2,6,6-tetramethylpiperidide as the base and the corresponding zinc diamide as an in situ trap. The reactivity of the substrate was discussed in light of the calculated atomic charges and the pK_a values. The behavior of the 1-arylated 7-azaindoles in direct iodination was then studied, and the results explained by considering the HOMO orbital coefficients and the atomic charges. Finally, some of the iodides generated, generally original, were involved in the N-arylation of indole. While crystallographic data were collected for fifteen of the synthesized compounds, biological properties (antimicrobial, antifungal and antioxidant activity) were evaluated for others.     

MCCE2: Improving protein pKa calculations with extensive side chain rotamer sampling

Multiconformation continuum electrostatics (MCCE) explores different conformational degrees of freedom in Monte Carlo calculations of protein residue and ligand pK_as. Explicit changes in side chain conformations throughout a titration create a position dependent, heterogeneous dielectric response giving a more accurate picture of coupled ionization and position changes. The MCCE2 methods for choosing a group of input heavy atom and proton positions are described. The pK_as calculated with different isosteric conformers, heavy atom rotamers and proton positions, with different degrees of optimization are tested against a curated group of 305 experimental pK_as in 33 proteins. QUICK calculations, with rotation around Asn and Gln termini, sampling His tautomers and torsion minimum hydroxyls yield an RMSD of 1.34 with 84% of the errors being <1.5 pH units. FULL calculations adding heavy atom rotamers and side chain optimization yield an RMSD of 0.90 with 90% of the errors <1.5 pH unit. Good results are also found for pK_as in the membrane protein bacteriorhodopsin. The inclusion of extra side chain positions distorts the dielectric boundary and also biases the calculated pK_as by creating more neutral than ionized conformers. Methods for correcting these errors are introduced. Calculations are compared with multiple X‐ray and NMR derived structures in 36 soluble proteins. Calculations with X‐ray structures give significantly better pK_as. Results with the default protein dielectric constant of 4 are as good as those using a value of 8. The MCCE2 program can be downloaded from http://www.sci.ccny.cuny.edu/～mcce . © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Computing ionization states of proteins with a detailed charge model

A convenient computational approach for the calculation of the p K_as of ionizable groups in a protein is described. The method uses detailed models of the charges in both the neutral and ionized form of each ionizable group. A full derivation of the theoretical framework is presented, as are details of its implementation in the UHBD program. Application to four proteins whose crystal structures are known shows that the detailed charge model improves agreement with experimentally determined pK_as when a low protein dielectric constant is assumed, relative to the results with a simpler single-site ionization model. It is also found that use of the detailed charge model increases the sensitivity of the computed pK_as to the details of proton placement. © 1996 by John Wiley & Sons, Inc.     

Determination of chromatographic dissociation constants of some carbapenem group antibiotics and quantification of these compounds in human urine

The dissociation constant values (_s^spK_a) of some carbapenem group drugs (ertapenem, meropenem, doripenem) in different percentages of methanol–water binary mixtures (18, 20 and 22%, v/v) were determined from the mobile phase pH dependence of their retention factor. Evaluation of these data was performed using the NLREG program. From calculated pK_a values, the aqueous pK_a values of these subtances were calculated by different approaches. Moreover, the correlation established between retention factor and the pH of the water–methanol mobile phase was used to determine the optimum separation conditions. In order to validate the optimized conditions, these drugs were studied in human urine. The chromatographic separation was realized using a Gemini NX C₁₈ column (250 × 4.6 mm i.d., 5 µm particles) and UV detector set at 220 and 295 nm. Copyright © 2013 John Wiley & Sons, Ltd.     

Relative pKa of some anilinium derivatives in methanol,acetonitrile, and tetrahydrofurane solvents

In this study, the relative pK_a values of nine anilinium derivatives in methanol (MeOH), acetonitrile (AN), and tetrahydrofurane (THF) solutions were successfully calculated with mean absolute deviations of 0.63, 0.68, and 0.75 pK_a units, respectively. To this aim, their gas‐phase basicities were computed using the CBS‐QB3 composite method. Also, conductor‐like polarizable continuum model (CPCM) with UAHF, UAKS and UA0 cavities and SM8 solvation models at HF/6‐31+G(d) level of theory were applied for the calculation of the solvation Gibbs free energies. The obtained results indicate that there is reliable correlation between the experimental and computed pK_a values in the studied solutions. Therefore, to extend the pK_a database for anilines, correlation equations were used to predict the pK_a values in the investigated solvents.     

Combined Experimental and Computational Studies on the Physical and Chemical Properties of the Renewable Amide, 3‐Acetamido‐5‐acetylfuran

The pK_a of 3‐acetamido‐5‐acetylfuran (3A5AF) was predicted to be in the range 18.5–21.5 by using the B3LYP/6‐311+G(2d,p) method and several amides as references. The experimental pK_a value, 20.7, was determined through UV/Vis titrations. Its solubility was measured in methanol‐modified supercritical CO₂ (mole fraction, 3.23×10^?4, cloud points 40–80 °C) and it was shown to be less soluble than 5‐hydroxymethylfurfural (5‐HMF). Dimerization energies were calculated for 3A5AF and 5‐HMF to compare hydrogen bonding, as such interactions will affect their solubility. Infrared and ¹H nuclear magnetic resonance spectra of 3A5AF samples support the existence of intermolecular hydrogen bonding. The highest occupied molecular orbital, lowest unoccupied molecular orbital, and electrostatic potential of 3A5AF were determined through molecular orbital calculations using B3LYP/6‐311+G(2d,p). The π–π* transition energy (time‐dependent density functional theory study) was compared with UV/Vis data. Calculated atomic charges were used in an attempt to predict the reactivity of 3A5AF. A reaction between 3A5AF and CH₃MgBr was conducted. As 3A5AF is a recently developed renewable compound that has previously not been studied extensively, these studies will be helpful in designing future reactions and processes involving this molecule.