Standard state free energies, not pKas, are ideal for describing small molecule protonation and tautomeric states.

Abstract:

:The pKa is the standard measure used to describe the aqueous proton affinity of a compound, indicating the proton concentration (pH) at which two protonation states (e.g. A- and AH) have equal free energy. However, compounds can have additional protonation states (e.g. AH2+), and may assume multiple tautomeric forms, with the protons in different positions (microstates). Macroscopic pKas give the pH where the molecule changes its total number of protons, while microscopic pKas identify the tautomeric states involved. As tautomers have the same number of protons, the free energy difference between them and their relative probability is pH independent so there is no pKa connecting them. The question arises: What is the best way to describe protonation equilibria of a complex molecule in any pH range? Knowing the number of protons and the relative free energy of all microstates at a single pH, ∆G°, provides all the information needed to determine the free energy, and thus the probability of each microstate at each pH. Microstate probabilities as a function of pH generate titration curves that highlight the low energy, observable microstates, which can then be compared with experiment. A network description connecting microstates as nodes makes it straightforward to test thermodynamic consistency of microstate free energies. The utility of this analysis is illustrated by a description of one molecule from the SAMPL6 Blind pKa Prediction Challenge. Analysis of microstate ∆G°s also makes a more compact way to archive and compare the pH dependent behavior of compounds with multiple protonatable sites.

journal_name

J Comput Aided Mol Des

authors

Gunner MR,Murakami T,Rustenburg AS,Işık M,Chodera JD

doi

10.1007/s10822-020-00280-7

subject

Has Abstract

pub_date

2020-05-01 00:00:00

pages

561-573

issue

5

eissn

0920-654X

issn

1573-4951

pii

10.1007/s10822-020-00280-7

journal_volume

34

pub_type

杂志文章
  • Geometry optimization method versus predictive ability in QSPR modeling for ionic liquids.

    abstract::Computational techniques, such as Quantitative Structure-Property Relationship (QSPR) modeling, are very useful in predicting physicochemical properties of various chemicals. Building QSPR models requires calculating molecular descriptors and the proper choice of the geometry optimization method, which will be dedicat...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-016-9894-3

    authors: Rybinska A,Sosnowska A,Barycki M,Puzyn T

    更新日期:2016-02-01 00:00:00

  • An improved method to predict the entropy term with the MM/PBSA approach.

    abstract::A method is suggested to calculate improved entropies within the MM/PBSA approach (molecular mechanics combined with Poisson-Boltzmann and surface area calculations) to estimate protein-ligand binding affinities. In the conventional approach, the protein is truncated outside ~8 A from the ligand. This system is freely...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-008-9238-z

    authors: Kongsted J,Ryde U

    更新日期:2009-02-01 00:00:00

  • Classification of protein disulphide-bridge topologies.

    abstract::The preferential occurrence of certain disulphide-bridge topologies in proteins has prompted us to design a method and a program, KNOT-MATCH, for their classification. The program has been applied to a database of proteins with less than 65% homology and more than two disulphide bridges. We have investigated whether t...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1011164224144

    authors: Mas JM,Aloy P,Martí-Renom MA,Oliva B,de Llorens R,Avilés FX,Querol E

    更新日期:2001-05-01 00:00:00

  • Modelling of carbohydrate-aromatic interactions: ab initio energetics and force field performance.

    abstract::Aromatic amino acid residues are often present in carbohydrate-binding sites of proteins. These binding sites are characterized by a placement of a carbohydrate moiety in a stacking orientation to an aromatic ring. This arrangement is an example of CH/pi interactions. Ab initio interaction energies for 20 carbohydrate...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-005-9033-z

    authors: Spiwok V,Lipovová P,Skálová T,Vondrácková E,Dohnálek J,Hasek J,Králová B

    更新日期:2005-12-01 00:00:00

  • Activity cliffs in PubChem confirmatory bioassays taking inactive compounds into account.

    abstract::Activity cliffs are formed by pairs or groups of structurally similar compounds with significant differences in potency. They represent a prominent feature of activity landscapes of compound data sets and a primary source of structure-activity relationship (SAR) information. Thus far, activity cliffs have only been co...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-012-9632-4

    authors: Hu Y,Maggiora GM,Bajorath J

    更新日期:2013-02-01 00:00:00

  • Molecular surface-volume and property matching to superpose flexible dissimilar molecules.

    abstract::Steric complementarity is a prerequisite for ligand-receptor recognition; this implies that drugs with a common receptor binding site should possess sterically similar binding surfaces. This principle is used as the basis for an automatic and unbiased method that superposes molecules. One molecule is rotated and trans...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00124319

    authors: Perkins TD,Mills JE,Dean PM

    更新日期:1995-12-01 00:00:00

  • Active-site-directed 3D database searching: pharmacophore extraction and validation of hits.

    abstract::Two new computational tools, PRO_PHARMEX and PRO_SCOPE, for use in active-site-directed searching of 3D databases are described. PRO_PHARMEX is a flexible, graphics-based program facilitating the extraction of pharmacophores from the active site of a target macromolecule. These pharmacophores can then be used to searc...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00124472

    authors: Clark DE,Westhead DR,Sykes RA,Murray CW

    更新日期:1996-10-01 00:00:00

  • Prediction of the three-dimensional structure of the human Fas receptor by comparative molecular modeling.

    abstract::The Fas antigen, a cell surface receptor belonging to the tumor necrosis factor receptor (TNFR) superfamily, triggers programmed cell death (apoptosis) in the immune system. The three-dimensional structure of Fas and molecular details of the interaction between Fas and its ligand are currently unknown. A three-dimensi...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1008011024584

    authors: Bajorath J,Aruffo A

    更新日期:1997-01-01 00:00:00

  • Computational design of a Diels-Alderase from a thermophilic esterase: the importance of dynamics.

    abstract::A novel computational Diels-Alderase design, based on a relatively rare form of carboxylesterase from Geobacillus stearothermophilus, is presented and theoretically evaluated. The structure was found by mining the PDB for a suitable oxyanion hole-containing structure, followed by a combinatorial approach to find suita...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-012-9601-y

    authors: Linder M,Johansson AJ,Olsson TS,Liebeschuetz J,Brinck T

    更新日期:2012-09-01 00:00:00

  • The structure-activity relationship of inhibitors of serotonin uptake and receptor binding.

    abstract::An analysis of five different datasets of inhibitors of serotonin uptake has yielded quantitative structure/activity relationships (QSARs) which delineate the role of steric and hydrophobic properties essential for inhibition by phenylethylamine-type analogues. ...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00125664

    authors: Hansch C,Caldwell J

    更新日期:1991-10-01 00:00:00

  • Identification of family-specific residue packing motifs and their use for structure-based protein function prediction: II. Case studies and applications.

    abstract::This paper describes several case studies concerning protein function inference from its structure using our novel approach described in the accompanying paper. This approach employs family-specific motifs, i.e. three-dimensional amino acid packing patterns that are statistically prevalent within a protein family. For...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-009-9277-0

    authors: Bandyopadhyay D,Huan J,Prins J,Snoeyink J,Wang W,Tropsha A

    更新日期:2009-11-01 00:00:00

  • Molecular dynamics simulations of oligonucleotides in solution: visualization of intrinsic curvature.

    abstract::We have undertaken molecular dynamics simulations on the d(CGCAAAAAAGCG).d(CGCTTTTTTGCG) dodecamer in solution. In this study, we focus on aspects of conformation and dynamics, including the possibility of cross-strand hydrogen bonds. We compare our results with those from crystallography as well as infrared, Raman an...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00126748

    authors: de Souza ON,Goodfellow JM

    更新日期:1994-06-01 00:00:00

  • Computer Automated Structure Evaluation (CASE) of the teratogenicity of retinoids with the aid of a novel geometry index.

    abstract::The CASE (Computer Automated Structure Evaluation) program, with the aid of a geometry index for discriminating cis and trans isomers, has been used to study a set of retinoids tested for teratogenicity in hamsters. CASE identified 8 fragments, the most important representing the non-polar terminus of a retinoid with ...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00125314

    authors: Klopman G,Dimayuga ML

    更新日期:1990-06-01 00:00:00

  • Binding free energy calculations to rationalize the interactions of huprines with acetylcholinesterase.

    abstract::In the present study, the binding free energy of a family of huprines with acetylcholinesterase (AChE) is calculated by means of the free energy perturbation method, based on hybrid quantum mechanics and molecular mechanics potentials. Binding free energy calculations and the analysis of the geometrical parameters hig...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-018-0114-1

    authors: Nascimento ÉCM,Oliva M,Andrés J

    更新日期:2018-05-01 00:00:00

  • Molecular dynamics simulation of halogen bonding mimics experimental data for cathepsin L inhibition.

    abstract::A MD simulation protocol was developed to model halogen bonding in protein-ligand complexes by inclusion of a charged extra point to represent the anisotropic distribution of charge on the halogen atom. This protocol was then used to simulate the interactions of cathepsin L with a series of halogenated and non-halogen...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-014-9802-7

    authors: Celis-Barros C,Saavedra-Rivas L,Salgado JC,Cassels BK,Zapata-Torres G

    更新日期:2015-01-01 00:00:00

  • Addressing free fatty acid receptor 1 (FFAR1) activation using supervised molecular dynamics.

    abstract::The free fatty acid receptor 1 (FFAR1, formerly GPR40), is a potential G protein-coupled receptor (GPCR) target for the treatment of type 2 diabetes mellitus (T2DM), as it enhances glucose-dependent insulin secretion upon activation by endogenous long-chain free fatty acids. The presence of two allosterically communic...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-020-00338-6

    authors: Atanasio S,Deganutti G,Reynolds CA

    更新日期:2020-11-01 00:00:00

  • Improving small molecule force fields by identifying and characterizing small molecules with inconsistent parameters.

    abstract::Many molecular simulation methods use force fields to help model and simulate molecules and their behavior in various environments. Force fields are sets of functions and parameters used to calculate the potential energy of a chemical system as a function of the atomic coordinates. Despite the widespread use of force ...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-020-00367-1

    authors: Ehrman JN,Lim VT,Bannan CC,Thi N,Kyu DY,Mobley DL

    更新日期:2021-01-28 00:00:00

  • Efficient overlay of small organic molecules using 3D pharmacophores.

    abstract::Aligning and overlaying two or more bio-active molecules is one of the key tasks in computational drug discovery and bio-activity prediction. Especially chemical-functional molecule characteristics from the view point of a macromolecular target represented as a 3D pharmacophore are the most interesting similarity meas...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-006-9078-7

    authors: Wolber G,Dornhofer AA,Langer T

    更新日期:2006-12-01 00:00:00

  • Exploring sets of molecules from patents and relationships to other active compounds in chemical space networks.

    abstract::Patents from medicinal chemistry represent a rich source of novel compounds and activity data that appear only infrequently in the scientific literature. Moreover, patent information provides a primary focal point for drug discovery. Accordingly, text mining and image extraction approaches have become hot topics in pa...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-017-0061-2

    authors: Kunimoto R,Bajorath J

    更新日期:2017-09-01 00:00:00

  • The SAMPL6 challenge on predicting aqueous pKa values from EC-RISM theory.

    abstract::The "embedded cluster reference interaction site model" (EC-RISM) integral equation theory is applied to the problem of predicting aqueous pKa values for drug-like molecules based on an ensemble of tautomers. EC-RISM is based on self-consistent calculations of a solute's electronic structure and the distribution funct...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-018-0140-z

    authors: Tielker N,Eberlein L,Güssregen S,Kast SM

    更新日期:2018-10-01 00:00:00

  • Calculation of hydrophobic parameters directly from three-dimensional structures using comparative molecular field analysis.

    abstract::Capacity ratio (log k') values, which are a measure of hydrophobicity, were calculated directly from the three-dimensional structures of 17 furans and 54 triazines using the comparative molecular field analysis approach. The H2O probe and the GRID force field, including hydrogen-bond potentials, yielded excellent corr...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00125172

    authors: Kim KH

    更新日期:1995-08-01 00:00:00

  • CADD medicine: design is the potion that can cure my disease.

    abstract::The acronym "CADD" is often used interchangeably to refer to "Computer Aided Drug Discovery" and "Computer Aided Drug Design". While the former definition implies the use of a computer to impact one or more aspects of discovering a drug, in this paper we contend that computational chemists are most effective when they...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-016-0004-3

    authors: Manas ES,Green DV

    更新日期:2017-03-01 00:00:00

  • Structural organization of G-protein-coupled receptors.

    abstract::Atomic-resolution structures of the transmembrane 7-alpha-helical domains of 26 G-protein-coupled receptors (GPCRs) (including opsins, cationic amine, melatonin, purine, chemokine, opioid, and glycoprotein hormone receptors and two related proteins, retinochrome and Duffy erythrocyte antigen) were calculated by distan...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章,评审

    doi:10.1023/a:1008050821744

    authors: Lomize AL,Pogozheva ID,Mosberg HI

    更新日期:1999-07-01 00:00:00

  • Computational docking simulations of a DNA-aptamer for argininamide and related ligands.

    abstract::The binding properties of sequence-specific nucleic acids (aptamers) to low-molecular-weight ligands, macromolecules and even cells attract substantial scientific interest. These ligand-DNA complexes found different applications for sensing, nanomedicine, and DNA nanotechnology. Structural information on the aptamer-l...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-015-9844-5

    authors: Albada HB,Golub E,Willner I

    更新日期:2015-07-01 00:00:00

  • Application of a simple quantum chemical approach to ligand fragment scoring for Trypanosoma brucei pteridine reductase 1 inhibition.

    abstract::There is a need for improved and generally applicable scoring functions for fragment-based approaches to ligand design. Here, we evaluate the performance of a computationally efficient model for inhibitory activity estimation, which is composed only of multipole electrostatic energy and dispersion energy terms that ap...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-017-0035-4

    authors: Jedwabny W,Panecka-Hofman J,Dyguda-Kazimierowicz E,Wade RC,Sokalski WA

    更新日期:2017-08-01 00:00:00

  • Automated site-directed drug design: searches of the Cambridge Structural Database for bond lengths in molecular fragments to be used for automated structure assembly.

    abstract::In this paper a database of small frequently occurring molecular fragments is used for the determination of fragment bond lengths from the Cambridge Structural Database. A large number of bond types are described that have not been reported previously. ...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00125946

    authors: Chau PL,Dean PM

    更新日期:1992-08-01 00:00:00

  • Molecular modeling of the intestinal bile acid carrier: a comparative molecular field analysis study.

    abstract::A structure-binding activity relationship for the intestinal bile acid transporter has been developed using data from a series of bile acid analogs in a comparative molecular field analysis (CoMFA). The studied compounds consisted of a series of bile acid-peptide conjugates, with modifications at the 24 position of th...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1007919704457

    authors: Swaan PW,Szoka FC Jr,Oie S

    更新日期:1997-11-01 00:00:00

  • In silico probing and biological evaluation of SETDB1/ESET-targeted novel compounds that reduce tri-methylated histone H3K9 (H3K9me3) level.

    abstract::ERG-associated protein with the SET domain (ESET/SET domain bifurcated 1/SETDB1/KMT1E) is a histone lysine methyltransferase (HKMT) and it preferentially tri-methylates lysine 9 of histone H3 (H3K9me3). SETDB1/ESET leads to heterochromatin condensation and epigenetic gene silencing. These functional changes are report...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-017-0052-3

    authors: Park I,Hwang YJ,Kim T,Viswanath ANI,Londhe AM,Jung SY,Sim KM,Min SJ,Lee JE,Seong J,Kim YK,No KT,Ryu H,Pae AN

    更新日期:2017-10-01 00:00:00

  • Making priors a priority.

    abstract::When we build a predictive model of a drug property we rigorously assess its predictive accuracy, but we are rarely able to address the most important question, "How useful will the model be in making a decision in a practical context?" To answer this requires an understanding of the prior probability distribution ("t...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-010-9388-7

    authors: Segall M,Chadwick A

    更新日期:2010-12-01 00:00:00

  • Virtual screening using a conformationally flexible target protein: models for ligand binding to p38α MAPK.

    abstract::We have used virtual screening to develop models for the binding of aryl substituted heterocycles to p38α MAPK. Virtual screening was conducted on a number of p38α MAPK crystal structures using a library of 46 known p38α MAPK inhibitors containing a heterocyclic core substituted by pyridine and fluorophenyl rings (str...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-012-9569-7

    authors: Vinh NB,Simpson JS,Scammells PJ,Chalmers DK

    更新日期:2012-04-01 00:00:00