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 function of surrounding water. Following-up on the workflow developed after the SAMPL5 challenge on cyclohexane-water distribution coefficients we extended and improved the methodology by taking into account exact electrostatic solute-solvent interactions taken from the wave function in solution. As before, the model is calibrated against Gibbs energies of hydration from the "Minnesota Solvation Database" and a public dataset of acidity constants of organic acids and bases by adjusting in total 4 parameters, among which only 3 are relevant for predicting pKa values. While the best-performing training model yields a root-mean-square error (RMSE) of 1 pK unit, the corresponding test set prediction on the full SAMPL6 dataset of macroscopic pKa values using the same level of theory exhibits slightly larger error (1.7 pK units) than the best test set model submitted (1.7 pK units for corresponding training set vs. test set performance of 1.6). Post-submission analysis revealed a number of physical optimization options regarding the numerical treatment of electrostatic interactions and conformational sampling. While the experimental test set data revealed after submission was not used for reparametrizing the methodology, the best physically optimized models consequentially result in RMSEs of 1.5 if only improved electrostatic interactions are considered and of 1.1 if, in addition, conformational sampling accounts for quantum-chemically derived rankings. We conclude that these numbers are probably near the ultimate accuracy achievable with the simple 3-parameter model using a single or the two best-ranking conformations per tautomer or microstate. Finally, relations of the present macrostate approach to microstate pKa results are discussed and some illustrative results for microstate populations are presented.

journal_name

J Comput Aided Mol Des

authors

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

doi

10.1007/s10822-018-0140-z

subject

Has Abstract

pub_date

2018-10-01 00:00:00

pages

1151-1163

issue

10

eissn

0920-654X

issn

1573-4951

pii

10.1007/s10822-018-0140-z

journal_volume

32

pub_type

杂志文章
  • Molecular mechanics calculations on deaminooxytocin and on deamino-arginine-vasopressin and its analogues.

    abstract::The backbone conformations of the cyclic moieties of 1-[beta-mercaptopropionic acid]-oxytocin [( Mpa1]-OT), [1-beta-mercaptopropionic acid]-arginine-vasopressin [( Mpa1]-AVP), [1-(beta'-mercapto-beta,beta-cyclopentamethylene)propionic acid]-arginine-vasopressin [( Cpp1]-AVP), and [1-thiosalicylic acid]-arginine-vasopr...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF01532991

    authors: Liwo A,Tempczyk A,Grzonka Z

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

  • Can we separate active from inactive conformations?

    abstract::Molecular modeling methodologies such as molecular docking, pharmacophore modeling, and 3D-QSAR, rely on conformational searches of small molecules as a starting point. All of these methodologies seek conformations of the small molecules as they bind to target proteins, i.e., their active conformations. Thus the quest...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1016320106741

    authors: Diller DJ,Merz KM Jr

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

  • Comparative molecular field analysis and energy interaction studies of thrombin-inhibitor complexes.

    abstract::A Comparative Molecular Field Analysis (CoMFA) and an interaction energy-based method were applied on a database holding the 3D structures of 29 thrombin-inhibitor complexes. Several parameters were optimized in both methods in order to obtain the best correlation between theoretical and experimentally determined bind...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1008010016362

    authors: Bursi R,Grootenhuis PD

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

  • Quantitative surface field analysis: learning causal models to predict ligand binding affinity and pose.

    abstract::We introduce the QuanSA method for inducing physically meaningful field-based models of ligand binding pockets based on structure-activity data alone. The method is closely related to the QMOD approach, substituting a learned scoring field for a pocket constructed of molecular fragments. The problem of mutual ligand a...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-018-0126-x

    authors: Cleves AE,Jain AN

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

  • Protein-ligand docking using FFT based sampling: D3R case study.

    abstract::Fast Fourier transform (FFT) based approaches have been successful in application to modeling of relatively rigid protein-protein complexes. Recently, we have been able to adapt the FFT methodology to treatment of flexible protein-peptide interactions. Here, we report our latest attempt to expand the capabilities of t...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-017-0069-7

    authors: Padhorny D,Hall DR,Mirzaei H,Mamonov AB,Moghadasi M,Alekseenko A,Beglov D,Kozakov D

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

  • Solvent effect on the synthesis of clarithromycin: a molecular dynamics study.

    abstract::Clarithromycin (6-O-methylerythromycin A) is a 14-membered macrolide antibiotic which is active in vitro against clinically important gram-positive and gram-negative bacteria. The selectivity of the methylation of the C-6 OH group is studied on erythromycin A derivatives. To understand the effect of the solvent on the...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/b:jcam.0000030037.67742.cb

    authors: Duran D,Aviyente V,Baysa C

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

  • Molecular modeling of the human vasopressin V2 receptor/agonist complex.

    abstract::The V2 vasopressin renal receptor (V2R), which controls antidiuresis in mammals, is a member of the large family of heptahelical transmembrane (7TM) G protein-coupled receptors (GPCRs). Using the automated GPCR modeling facility available via Internet (http:/(/)expasy.hcuge.ch/swissmod/SWISS-MODEL.+ ++html) for constr...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1007969526447

    authors: Czaplewski C,Kaźmierkiewicz R,Ciarkowski J

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

  • AstexViewer: a visualisation aid for structure-based drug design.

    abstract::AstexViewer is a Java molecular graphics program that can be used for visualisation in many aspects of structure-based drug design. This paper describes its functionality, implementation and examples of its use. The program can run as an Applet in a web browser allowing structures to be displayed without installing ad...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1023813504011

    authors: Hartshorn MJ

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

  • The IUPAC aqueous and non-aqueous experimental pKa data repositories of organic acids and bases.

    abstract::Accurate and well-curated experimental pKa data of organic acids and bases in both aqueous and non-aqueous media are invaluable in many areas of chemical research, including pharmaceutical, agrochemical, specialty chemical and property prediction research. In pharmaceutical research, pKa data are relevant in ligand de...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-014-9764-9

    authors: Slater AM

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

  • In silico molecular docking analysis of the human Argonaute 2 PAZ domain reveals insights into RNA interference.

    abstract::RNA interference (RNAi) is a critical cellular pathway activated by double stranded RNA and regulates the gene expression of target mRNA. During RNAi, the 3' end of siRNA binds with the PAZ domain, followed by release and rebinding in a cyclic manner, which deemed essential for proper gene silencing. Recently, we prov...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-013-9665-3

    authors: Kandeel M,Kitade Y

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

  • Identification of ligand efficient, fragment-like hits from an HTS library: structure-based virtual screening and docking investigations of 2H- and 3H-pyrazolo tautomers for Aurora kinase A selectivity.

    abstract::Furanopyrimidine 1 (IC50 = 273 nM, LE = 0.36, LELP = 10.28) was recently identified by high-throughput screening (HTS) of an in-house library (125,000 compounds) as an Aurora kinase inhibitor. Structure-based hit optimization resulted in lead molecules with in vivo efficacy in a mouse tumour xenograft model, but no or...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-014-9807-2

    authors: Sarvagalla S,Singh VK,Ke YY,Shiao HY,Lin WH,Hsieh HP,Hsu JT,Coumar MS

    更新日期:2015-01-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

  • Computer simulation of the binding of amonafide and azonafide to DNA.

    abstract::Intercalative binding of the antitumor drugs amonafide and azonafide to the oligonucleotide duplex d(GGCCGGCCGG).d(CCGGCCGGCC) was compared using molecular dynamics in vacuum with the AMBER force field. A number of reasonable possible binding conformations were obtained, with the azonafide complexes favored over the a...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00402824

    authors: Bear S,Remers WA

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

  • An improved scoring function for suboptimal polar ligand complexes.

    abstract::Learning strategies can be used to improve the efficiency of virtual screening of very large databases. In these strategies new compounds to be screened are selected on the basis of the results obtained in previous stages, even if truly good ligands have not yet been identified. This approach requires that the scoring...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-008-9246-z

    authors: Cincilla G,Vidal D,Pons M

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

  • New designs for MRI contrast agents.

    abstract::New designs for Magnetic Resonance Imaging contrast agents are presented. Essentially, they all are host-guest inclusion complexes between y-cyclodextrins and polyazamacrocycles of gadolinium (III) ion. Substitutions have been made to the host to optimise the host-guest association. Molecular mechanics calculations ha...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1027347527385

    authors: Fernandes PA,Carvalho AT,Marques AT,Pereira AL,Madeira AP,Ribeiro AS,Carvalho AF,Ricardo ET,Pinto FJ,Santos HA,Mangericão HD,Martins HM,Pinto HD,Santos HR,Moreira IS,Azeredo MJ,Abreu RP,Oliveira RM,Sousa SF,Silva RJ

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

  • Blind prediction of cyclohexane-water distribution coefficients from the SAMPL5 challenge.

    abstract::In the recent SAMPL5 challenge, participants submitted predictions for cyclohexane/water distribution coefficients for a set of 53 small molecules. Distribution coefficients (log D) replace the hydration free energies that were a central part of the past five SAMPL challenges. A wide variety of computational methods w...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-016-9954-8

    authors: Bannan CC,Burley KH,Chiu M,Shirts MR,Gilson MK,Mobley DL

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

  • Performance of multiple docking and refinement methods in the pose prediction D3R prospective Grand Challenge 2016.

    abstract::We describe the performance of multiple pose prediction methods for the D3R 2016 Grand Challenge. The pose prediction challenge includes 36 ligands, which represent 4 chemotypes and some miscellaneous structures against the FXR ligand binding domain. In this study we use a mix of fully automated methods as well as hum...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-017-0053-2

    authors: Fradera X,Verras A,Hu Y,Wang D,Wang H,Fells JI,Armacost KA,Crespo A,Sherborne B,Wang H,Peng Z,Gao YD

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

  • A novel view of modelling interactions between synthetic and biological polymers via docking.

    abstract::Multipoint interactions between synthetic and natural polymers provide a promising platform for many topical applications, including therapeutic blockage of virus-specific targets. Docking may become a useful tool for modelling of such interactions. However, the rigid docking cannot be correctly applied to synthetic p...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-012-9621-7

    authors: Tsvetkov VB,Serbin AV

    更新日期:2012-12-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

  • An automated method for predicting the positions of hydrogen-bonding atoms in binding sites.

    abstract::Hydrogen bonds are the most specific, and therefore predictable of the intermolecular interactions involved in ligand-protein binding. Given the structure of a molecule, it is possible to estimate the positions at which complementary hydrogen-bonding atoms could be found. Crystal-survey data are used in the design of ...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1007900527102

    authors: Mills JE,Perkins TD,Dean PM

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

  • LASSO-ligand activity by surface similarity order: a new tool for ligand based virtual screening.

    abstract::Virtual Ligand Screening (VLS) has become an integral part of the drug discovery process for many pharmaceutical companies. Ligand similarity searches provide a very powerful method of screening large databases of ligands to identify possible hits. If these hits belong to new chemotypes the method is deemed even more ...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-007-9164-5

    authors: Reid D,Sadjad BS,Zsoldos Z,Simon A

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

  • Cavity search: an algorithm for the isolation and display of cavity-like binding regions.

    abstract::A set of algorithms designed to enhance the display of protein binding cavities is presented. These algorithms, collectively entitled CAVITY SEARCH, allow the user to isolate and fully define the extent of a particular cavity. Solid modeling techniques are employed to produce a detailed cast of the active site region,...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00117400

    authors: Ho CM,Marshall GR

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

  • SAMPL6: calculation of macroscopic pKa values from ab initio quantum mechanical free energies.

    abstract::Macroscopic pKa 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 pKa values were derived from the gas-phase free energy difference between pro...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-018-0138-6

    authors: Selwa E,Kenney IM,Beckstein O,Iorga BI

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

  • Hydrophobic molecular similarity from MST fractional contributions to the octanol/water partition coefficient.

    abstract::The use of a recently proposed hydrophobic similarity index for the alignment of molecules and the prediction of their differences in biological activity is described. The hydrophobic similarity index exploits atomic contributions to the octanol/water transfer free energy, which are evaluated by means of the fractiona...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-005-7928-3

    authors: Muñoz-Muriedas J,Perspicace S,Bech N,Guccione S,Orozco M,Luque FJ

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

  • Structure-toxicity relationships of polycyclic aromatic hydrocarbons using molecular quantum similarity.

    abstract::The establishment of quantitative structure-activity relationship (QSAR) models for the toxicity of polycylic aromatic hydrocarbons (PAHs) is described. Two properties, in vitro percutaneous absorption in rat skin and discrete levels of carcinogenic activity, are examined using molecular quantum similarity measures (M...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/a:1011150003086

    authors: Gallegos A,Robert D,Gironés X,Carbó-Dorca R

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

  • Conformational behaviour of the antineoplastic peptide dolastatin-10 and of two mutated derivatives.

    abstract::The three-dimensional structure of dolastatin-10, an extremely potent cytostatic and antineoplastic peptide extracted from the mollusc Dolabella auricularia, has not yet been fully characterized in an experimental way. By means of a systematic conformational search of the natural peptide and of two mutated analogs, ca...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/BF00124000

    authors: Fantucci P,Marino T,Russo N,Villa AM

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

  • SAMPL7 blind predictions using nonequilibrium alchemical approaches.

    abstract::In the context of the SAMPL7 challenge, we computed, employing a non-equilibrium (NE) alchemical technique, the standard binding free energy of two series of host-guest systems, involving as a host the Isaac's TrimerTrip, a Cucurbituril-like open cavitand, and the Gilson's Cyclodextrin derivatives. The adopted NE alch...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-020-00365-3

    authors: Procacci P,Guarnieri G

    更新日期:2021-01-04 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

  • eFindSite: improved prediction of ligand binding sites in protein models using meta-threading, machine learning and auxiliary ligands.

    abstract::Molecular structures and functions of the majority of proteins across different species are yet to be identified. Much needed functional annotation of these gene products often benefits from the knowledge of protein-ligand interactions. Towards this goal, we developed eFindSite, an improved version of FINDSITE, design...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1007/s10822-013-9663-5

    authors: Brylinski M,Feinstein WP

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

  • A promising enantioselective Diels-Alder dienophile by computer-assisted rational design: 2,5-diphenyl-1-vinyl-borolane.

    abstract::Several chiral vinylboranes have been theoretically evaluated as enantioselective Diels-Alder dienophiles. Theoretical results suggest that optically pure 2,5-diphenyl-1-vinyl-borolane should be the reagent of choice for performing such transformations efficiently. ...

    journal_title:Journal of computer-aided molecular design

    pub_type: 杂志文章

    doi:10.1023/b:jcam.0000035200.30340.41

    authors: Pellegrinet SC,Silva MA,Goodman JM

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