Predicting protein pKa by environment similarity

A statistical method to predict protein pK_a has been developed by using the 3D structure of a protein and a database of 434 experimental protein pK_a values. Each pK_a in the database is associated with a fingerprint that describes the chemical environment around an ionizable residue. A computational tool, MoKaBio, has been developed to identify automatically ionizable residues in a protein, generate fingerprints that describe the chemical environment around such residues, and predict pK_a from the experimental pK_a values in the database by using a similarity metric. The method, which retrieved the pK_a of 429 of the 434 ionizable sites in the database correctly, was crossvalidated by leave‐one‐out and yielded root mean square error (RMSE) = 0.95, a result that is superior to that obtained by using the Null Model (RMSE 1.07) and other well‐established protein pK_a prediction tools. This novel approach is suitable to rationalize protein pK_a by comparing the region around the ionizable site with similar regions whose ionizable site pK_a is known. The pK_a of residues that have a unique environment not represented in the training set cannot be predicted accurately, however, the method offers the advantage of being trainable to increase its predictive power. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

The antimalarial activity of Ru–chloroquine complexes against resistant Plasmodium falciparum is related to lipophilicity, basicity, and heme aggregation inhibition ability near water/n-octanol interfaces

We have measured water/n-octanol partition coefficients, pK _a values, heme binding constants, and heme aggregation inhibition activity of a series of ruthenium–π-arene–chloroquine (CQ) complexes recently reported to be active against CQ-resistant strains of Plasmodium falciparum. Measurements of heme aggregation inhibition activity of the metal complexes near water/n-octanol interfaces qualitatively predict their superior antiplasmodial action against resistant parasites, in relation to CQ; we conclude that this modified method may be a better predictor of antimalarial potency than standard tests in aqueous acidic buffer. Some interesting tendencies emerge from our data, indicating that the antiplasmodial activity is related to a balance of effects associated with the lipophilicity, basicity, and structural details of the compounds studied. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Chemically accurate protein structures: Validation of protein NMR structures by comparison of measured and predicted pK a values

A new method is presented for evaluating the quality of protein structures obtained by NMR. This method exploits the dependence between measurable chemical properties of a protein, namely pK _a values of acidic residues, and protein structure. The accurate and fast empirical computational method employed by the PROPKA program () allows the user to test the ability of a given structure to reproduce known pK _a values, which in turn can be used as a criterion for the selection of more accurate structures. We demonstrate the feasibility of this novel idea for a series of proteins for which both␣NMR and X-ray structures, as well as pK _a values of all ionizable residues, have been determined. For the 17 NMR ensembles used in this study, this criterion is shown effective in the elimination of a large number of NMR structure ensemble members.     

Improving Compound Quality through in vitro and in silico Physicochemical Profiling

Many compounds entering clinical studies do not survive the numerous hurdles for a good pharmacological lead to a drug on the market. The reasons for attrition have been widely studied which resulted in more early attention to compound quality related to physical chemistry, drug metabolism and pharmacokinetics (DMPK), and toxicology/safety. This paper will briefly review current physicochemical in vitro assays and in silico predictions to support compound and library design through to lead optimization. The most important physicochemical properties include lipophilicity (log P/D), pK_a, solubility, and permeability. These drive key ADMET properties such as absorption, cell penetration, access to the brain, volume of distribution, plasma protein binding, metabolism, and toxicity, as well as biopharmaceutical behavior. Much data are now available from medium‐ to high‐throughput physchem and ADMET in vitro assays, either in the public domain (see, e.g., PubChem, PubMed) or in drug companies' in‐house databases. Such data are increasingly being computer‐modelled and used in predictive chemistry. New pipelining technology makes it easier to build and update QSAR models so that such models can use the latest available data to produce robust local and global predictive in silico ADMET models.     

Assessment of antidotal efficacy of cholinesterase reactivators against paraoxon: In vitro reactivation kinetics and physicochemical properties

The search of proficient oximes as reactivators of irreversibly inhibited-AChE by organophosphate poisoning necessitates an appropriate assessment of their physicochemical properties and reactivation kinetics. Therefore, herein acid dissociation constant; pK_a, lipophilicity; log P, polar surface area, hydrogen bond donor and acceptor counts of structurally different oximes (two tertiary oximes and thirteen pyridinium aldoxime derivatives) have been evaluated. The experimentally obtained data for pK_a has been comparatively analyzed by using non-linear regression. Further the tested oximes were screened through in vitro reactivation kinetics against paraoxon-inhibited AChE. The pK_a values of all the examined oximes were within the range of 7.50–9.53. pK_a values of uncharged and mono-pyridinium oximes were in good correlation with their reactivation potency. The high negative log P values of pyridinium oxime reactivators indicate their high hydrophilic character; hence oximes with improved lipophilicity should be designed for the development of novel and more potent antidotes. Propane and butane linked oximes were superior reactivators than xylene linked bis-oxime reactivators. It is concluded from the present study that pK_a value is not only ruled by the position of oximino functionality in the pyridinium ring, but also by the position of linker. Although, pyridinium oximes are proved to be better reactivators but their lipophilicity has to be improved.     

Influence of p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> Shifts on the Calculated Dipole Moments of Proteins

The protein dipole moment is a low-resolution parameter that characterizes the second-order charge organization of a biomolecule. Theoretical approaches to calculate protein dipole moments rely on pK _a values, which are either computed individually for each ionizable residue or obtained from model compounds. The influence of pK _a shifts are evaluated first by comparing calculated and measured dipole moments of β-lactoglobulin. Second, calculations are made on a dataset of 66 proteins from the Protein Data Bank, and average differences are determined between dipole moments calculated with model pK _as, pK _as derived using a Poisson–Boltzmann approach, and empirically-calculated pK _as. Dipole moment predictions that neglect pK _a shifts are consistently larger than predictions in which they are included. The importance of pK _a shifts are observed to vary with protein size, internal permittivity, and solution pH.     

Functional tuning of the catalytic residue pKa in a de novo designed esterase

AlleyCatE is a de novo designed esterase that can be allosterically regulated by calcium ions. This artificial enzyme has been shown to hydrolyze p‐nitrophenyl acetate (pNPA) and 4‐nitrophenyl‐(2‐phenyl)‐propanoate (pNPP) with high catalytic efficiency. AlleyCatE was created by introducing a single‐histidine residue (His¹⁴⁴) into a hydrophobic pocket of calmodulin. In this work, we explore the determinants of catalytic properties of AlleyCatE. We obtained the pK_a value of the catalytic histidine using experimental measurements by NMR and pH rate profile and compared these values to those predicted from electrostatics pK_a calculations (from both empirical and continuum electrostatics calculations). Surprisingly, the pK_a value of the catalytic histidine inside the hydrophobic pocket of calmodulin is elevated as compared to the model compound pK_a value of this residue in water. We determined that a short‐range favorable interaction with Glu¹²⁷ contributes to the elevated pK_a of His¹⁴⁴. We have rationally modulated local electrostatic potential in AlleyCatE to decrease the pK_a of its active nucleophile, His¹⁴⁴, by 0.7 units. As a direct result of the decrease in the His¹⁴⁴ pK_a value, catalytic efficiency of the enzyme increased by 45% at pH 6. This work shows that a series of simple NMR experiments that can be performed using low field spectrometers, combined with straightforward computational analysis, provide rapid and accurate guidance to rationally improve catalytic efficiency of histidine‐promoted catalysis. Proteins 2017; 85:1656–1665. © 2017 Wiley Periodicals, Inc.     

Bayesian model aggregation for ensemble‐based estimates of protein pKa values

This article investigates an ensemble‐based technique called Bayesian Model Averaging (BMA) to improve the performance of protein amino acid pK_a predictions. Structure‐based pK_a calculations play an important role in the mechanistic interpretation of protein structure and are also used to determine a wide range of protein properties. A diverse set of methods currently exist for pK_a prediction, ranging from empirical statistical models to ab initio quantum mechanical approaches. However, each of these methods are based on a set of conceptual assumptions that can effect a model's accuracy and generalizability for pK_a prediction in complicated biomolecular systems. We use BMA to combine eleven diverse prediction methods that each estimate pK_a values of amino acids in staphylococcal nuclease. These methods are based on work conducted for the pK_a Cooperative and the pK_a measurements are based on experimental work conducted by the García‐Moreno lab. Our cross‐validation study demonstrates that the aggregated estimate obtained from BMA outperforms all individual prediction methods with improvements ranging from 45 to 73% over other method classes. This study also compares BMA's predictive performance to other ensemble‐based techniques and demonstrates that BMA can outperform these approaches with improvements ranging from 27 to 60%. This work illustrates a new possible mechanism for improving the accuracy of pK_a prediction and lays the foundation for future work on aggregate models that balance computational cost with prediction accuracy. Proteins 2014; 82:354–363. © 2013 Wiley Periodicals, Inc.     

In silico Prediction of Human Oral Absorption Based on QSAR Analyses of PAMPA Permeability

The parallel artificial membrane permeation assay (PAMPA) was developed as a model for the prediction of transcellular permeation in the process of drug absorption. Our research group has measured the PAMPA permeability of peptide‐related compounds, diverse drugs, and agrochemicals. This work led to a classical quantitative structure–activity relationship (QSAR) equation for PAMPA permeability coefficients of structurally diverse compounds based on simple physicochemical parameters such as lipophilicity at a particular pH (log P_oct and |pK_a?pH|), H‐bond acceptor ability (SA_HA), and H‐bond donor ability (SA_HD). Since the PAMPA permeability of lipophilic compounds decreased with their apparent lipophilicity due to the unstirred water layer (UWL) barrier on membrane surfaces and to membrane retention, a bilinear QSAR model was introduced to explain the permeability of a broader set of compounds using the same physicochemical parameters as those used for the linear model. We also compared PAMPA and Caco‐2 cell permeability coefficients of compounds transported by various absorption mechanisms. The compounds were classified according to their absorption pathway (passively transported compounds, actively transported compounds, and compounds excreted by efflux systems) in the plot of Caco‐2 vs. PAMPA permeability. Finally, based on the QSAR analyses of PAMPA permeability, an in silico prediction model of human oral absorption for possibly transported compounds was proposed, and the usefulness of the model was examined.     

A predictive model for the selective accumulation of chemicals in tumor cells

Cationic lipophilic dyes can accumulate in mitochondria, and especially in mitochondria of tumor cells. We investigated the chemical properties and the processes allowing selective uptake into tumor cells using the Fick–Nernst–Planck equation. The model simulates uptake into cytoplasm and mitochondria and is valid for neutral molecules and ions, and thus also for weak electrolytes. The differential equation system was analytically solved for the steady-state and the dynamic case. The parameterization was for a generic human cell, with a 60 mV more negative potential at the inner mitochondrial membrane of generic tumor cells. The chemical input data were the lipophilicity (logK_OW), the acid/base dissociation constant (pK_a) and the electric charge (z). Accumulation in mitochondria occurred for polar acids with pK_a between 5 and 9 owing to the ion trap, and for lipophilic bases with pK_a>11 or permanent cations owing to electrical attraction. Selective accumulation in tumor cells was found for monovalent cations or strong bases with logK_OW of the cation between –2 and 2, with the optimum near 0. The results are in agreement with experimental results for rhodamine 123, a series of cationic triarylmethane dyes, F16 and MKT-077, an anticancer drug targeting tumor mitochondria.     

Mutation of non-conserved amino acids surrounding catalytic site to shift pH optimum of Bacillus circulans xylanase

Improvement of enzyme function by engineering pH dependence of enzymatic activity is of importance for industrial application of Bacillus circulans xylanases. Target mutation sites were selected by structural alignment between B. circulans xylanase and other xylanases having different pH optima. We selected non-conserved mutant sites within 8 Å from the catalytic residues, to see whether these residues have some role in modulating pK_as of the catalytic residues. We hypothesized that the non-conserved residues which may not have any role in enzyme catalysis might perturb pK_as of the catalytic residues. Change in pK_a of a titratable group due to change in electrostatic potential of a mutation was calculated and the change in pH optimum was predicted from the change in pK_a of the catalytic residues. Our strategy is proved to be useful in selection of promising mutants to shift the pH optimum of the xylanases towards desired side.     

Electrostatics of cysteine residues in proteins: Parameterization and validation of a simple model

One of the most popular and simple models for the calculation of pK_as from a protein structure is the semi‐macroscopic electrostatic model MEAD. This model requires empirical parameters for each residue to calculate pK_as. Analysis of current, widely used empirical parameters for cysteine residues showed that they did not reproduce expected cysteine pK_as; thus, we set out to identify parameters consistent with the CHARMM27 force field that capture both the behavior of typical cysteines in proteins and the behavior of cysteines which have perturbed pK_as. The new parameters were validated in three ways: (1) calculation across a large set of typical cysteines in proteins (where the calculations are expected to reproduce expected ensemble behavior); (2) calculation across a set of perturbed cysteines in proteins (where the calculations are expected to reproduce the shifted ensemble behavior); and (3) comparison to experimentally determined pK_a values (where the calculation should reproduce the pK_a within experimental error). Both the general behavior of cysteines in proteins and the perturbed pK_a in some proteins can be predicted reasonably well using the newly determined empirical parameters within the MEAD model for protein electrostatics. This study provides the first general analysis of the electrostatics of cysteines in proteins, with specific attention paid to capturing both the behavior of typical cysteines in a protein and the behavior of cysteines whose pK_a should be shifted, and validation of force field parameters for cysteine residues. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

The pKa Values of Acidic and Basic Residues Buried at the Same Internal Location in a Protein Are Governed by Different Factors

The pK_a values of internal ionizable groups are usually very different from the normal pK_a values of ionizable groups in water. To examine the molecular determinants of pK_a values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK_a, whereas Asp38 and Glu38 titrate with elevated pK_a values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK_a of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK_a values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK_a values of internal groups in proteins. This study suggests that improvements to computational methods for pK_a calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.     

Effect of methylation on the side‐chain pKa value of arginine

Arginine methylation is important in biological systems. Recent studies link the deregulation of protein arginine methyltransferases with certain cancers. To assess the impact of methylation on interaction with other biomolecules, the pK_a values of methylated arginine variants were determined using NMR data. The pK_a values of monomethylated, symmetrically dimethylated, and asymmetrically dimethylated arginine are similar to the unmodified arginine (14.2 ± 0.4). Although the pK_a value has not been significantly affected by methylation, consequences of methylation include changes in charge distribution and steric effects, suggesting alternative mechanisms for recognition.     

A genetic component of extinction risk in mammals

Genetic factors may play an important role in species extinction but their actual effect remains poorly understood, particularly because of a strong and potentially masking effect expected from ecological traits. We investigated the role of genetics in mammal extinction taking both ecological and genetic factors into account. As a proxy for the role of genetics we used the ratio of the rates of nonsynonymous (amino acid changing) to synonymous (leaving the amino acid unchanged) nucleotide substitutions, K_a / K_s. Because most nonsynonymous substitutions are likely to be slightly deleterious and thus selected against, this ratio is a measure of the inefficiency of selection: if large (but less than 1), it implies a low efficiency of selection against nonsynonymous mutations. As a result, nonsynonymous mutations may accumulate and thus contribute to extinction. As a proxy for the role of ecology we used body mass W, with which most extinction‐related ecological traits strongly correlate. As a measure of extinction risk we used species’ affiliation with the five levels of extinction threat according to the IUCN Red List of Threatened Species. We calculated K_a / K_s for mitochondrial protein‐coding genes of 211 mammalian species, each of which was characterized by body mass and the level of threat. Using logistic regression analysis, we then constructed a set of logistic regression models of extinction risk on ln(K_a / K_s) and lnW. We found that K_a / K_s and body mass are responsible for a 38% and a 62% increase in extinction risk, respectively. Given that the standard error of these values is 13%, the contribution of genetic factors to extinction risk in mammals is estimated to be one‐quarter to one‐half of the total of ecological and genetic effects. We conclude that the effect of genetics on extinction is significant, though it is almost certainly smaller than the effect of ecological traits. Synthesis Mutation provides the material for evolution. However, most mutations that play a role in evolution are slightly deleterious and thus may contribute to extinction. We assess the role of mitochondrial DNA mutations in mammalian extinction risk and find it to be one‐quarter to one‐half of the total of mutation and body mass effects, where body mass represents an integral measure of extinction‐related ecological traits. Genetic factors may be all the more important, because ecological traits associated with large body mass would both promote and protect from extinction, while mutation accumulation caused by low effective population size seems to have no counterbalance.     

Conformational sampling and polarization of Asp26 in pKa calculations of thioredoxin

Thioredoxin is a protein that has been used as model system by various computational methods to predict the pK_a of aspartate residue Asp26 which is 3.5 units higher than a solvent exposed one (eg, Asp20). Here, we use extensive atomistic molecular dynamics simulations of two different protonation states of Asp26 in combination with conformational analysis based on RMSD clustering and principle component analysis to identify representative conformations of the protein in solution. For each conformation, the Gibbs free energy of proton transfer between Asp26 and Asp20, which is fully solvated in a loop region of the protein, is calculated with the Amber99sb force field in alchemical transformations. The varying polarization of the two residues in different molecular environments and protonation states is described by Hirshfeld-I (HI) atomic charges obtained from the averaged polarized electron density. Our results show that the Gibbs free energy of proton transfer is dependent on the protein conformation, the proper sampling of the neighboring Lys57 residue orientations and on water molecules entering the hydrophobic cavity upon deprotonating Asp26. The inclusion of the polarization of both aspartate residues in the free energy cycle by HI atomic charges corrects the results from the non-polarizable force field and reproduces the experimental ΔpK_a value of Asp26.     

The Unstirrable Water Layer Is Unstirrable because It Does Not Exist

A number of membrane‐permeation models require the incorporation of an unstirred or unstirrable water layer (UWL). An example occurs in PAMPA models when the effective permeation rate of lipophilic acids and bases, P_e, falls behind the expected permeation rate, P_m, at pH values providing a high concentration of unionized species in the donor phase. In such cases, the compound has an apparent pK_a of a weaker acid or base. The explanation is that an UWL adjacent to the membrane provides a rate‐limiting diffusion barrier for such compounds. The thickness of the UWL is correlated with the difference between the aqueous pK_a and the apparent pK_a (pK ). Here, we provide an explanation for the pK term that requires no UWL. It comes from the fact that, in the process of passing into a membrane, an ionizable compound undergoes a change in pK_a. At some point along its path into the membrane, the compound attains a maximum free energy, at which point it is as likely to continue into the membrane, as it is to return to the donor phase. This is the transition state for absorption. The pK is the pK_a of the compound at the transition state. This is a testable hypothesis (see text). The relevance of absorption to permeation depends on the rate‐limiting step of permeation.     

The concerted action of a positive charge and hydrogen bonds dynamically regulates the pKa of the nucleophilic cysteine in the NrdH‐redoxin family

NrdH‐redoxins shuffle electrons from the NADPH pool in the cell to Class Ib ribonucleotide reductases, which in turn provide the precursors for DNA replication and repair. NrdH‐redoxins have a CVQC active site motif and belong to the thioredoxin‐fold protein family. As for other thioredoxin‐fold proteins, the pK_a of the nucleophilic cysteine of NrdH‐redoxins is of particular interest since it affects the catalytic reaction rate of the enzymes. Recently, the pK_a value of this cysteine in Corynebacterium glutamicum and Mycobacterium tuberculosis NrdH‐redoxins were determined, but structural insights explaining the relatively low pK_a remained elusive. We subjected C. glutamicum NrdH‐redoxin to an extensive molecular dynamics simulation to expose the factors regulating the pK_a of the nucleophilic cysteine. We found that the nucleophilic cysteine receives three hydrogen bonds from residues within the CVQC active site motif. Additionally, a fourth hydrogen bond with a lysine located N‐terminal of the active site further lowers the cysteine pK_a. However, site‐directed mutagenesis data show that the major contribution to the lowering of the cysteine pK_a comes from the positive charge of the lysine and not from the additional Lys‐Cys hydrogen bond. In 12% of the NrdH‐redoxin family, this lysine is replaced by an arginine that also lowers the cysteine pK_a. All together, the four hydrogen bonds and the electrostatic effect of a lysine or an arginine located N‐terminally of the active site dynamically regulate the pK_a of the nucleophilic cysteine in NrdH‐redoxins.     

Multidimensional optimization of promising antitumor xanthone derivatives

A promising antitumor xanthone derivative was optimized following a multidimensional approach that involved the synthesis of 17 analogues, the study of their lipophilicity and solubility, and the evaluation of their growth inhibitory activity on four human tumor cell lines. A new synthetic route for the hit xanthone derivative was also developed and applied for the synthesis of its analogues. Among the used cell lines, the HL-60 showed to be in general more sensitive to the compounds tested, with the most potent compound having a GI₅₀ of 5.1 μM, lower than the hit compound. Lipophilicity was evaluated by the partition coefficient (K_p) of a solute between buffer and two membrane models, namely liposomes and micelles. The compounds showed a log K_p between 3 and 5 and the two membrane models showed a good correlation (r² = 0.916) between each other. Studies concerning relationship between solubility and structure were developed for the hit compound and 5 of its analogues.     

Differentiation of various types of biological oxidation of nitrogen in organic compounds

The various types of nitrogen which occur in organic compounds and which are susceptible to biological oxidation are clearly divided into groups depending upon the pK_a, of the constituent nitrogen. The enzymatic processes which give rise to the N-oxidation products are reviewed by a consideration of species differences, age of animal, pH optima, influence of inducing agents, inhibitors and microsomal pretreatments, as well as the stereochemistry of the nitrogen atom.From the data collected, a concept is developed which suggests that all basic amines (group I) are oxidised by a flavine adenine nucleotide (FAD)-dependent enzyme system, whereas non-basic nitrogen-containing compounds (group III) are oxidised by a cytochrome P450-dependent system.It is further suggested that compounds of intermediary pK_a,i.e. between 1 and 7 (group II), may be substrates for both enzyme systems, which would yield the same products, but by different processes. The extent to which N-oxidation occurs in a species would therefore depend on the pK_a of the substrate and the amounts and ratio of the two enzymes present, a lower pK_a favouring oxidation by the cytochrome P450 system and a higher pK_a favouring oxidation by the FAD system.In a similar manner, it is suggested that the oxidation of aromatic heterocyclic amines depends upon the pK_a of the nitrogen, compounds having a low pK_a being preferentially metabolised by nitrogen oxidation.