Pka Measurements for the SAMPL6 Prediction Challenge for A

Total Page:16

File Type:pdf, Size:1020Kb

Pka Measurements for the SAMPL6 Prediction Challenge for A bioRxiv preprint doi: https://doi.org/10.1101/368787; this version posted July 13, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. PrEPRINT AHEAD OF SUBMISSION — July 13, 2018 P MEASUREMENTS FOR THE SAMPL6 1 K A PREDICTION CHALLENGE FOR A SET OF KINASE 2 inhibitor-likE FRAGMENTS 3 1,2 3 1,4 5 4 Mehtap Işık , DorOTHY LeVORSE , Ariën S. RustenburG , IKENNA E. Ndukwe , Heather 6 5 5 5 6 5 WANG , Xiao WANG , Mikhail Reibarkh , Gary E. Martin , AleXEY A. MakarOV , David L. 7 3* 1* 6 MobleY , Timothy Rhodes , John D. ChoderA Computational AND Systems Biology PrOGRam, Sloan Kettering Institute, Memorial Sloan 7 1 Kettering Cancer Center, NeW York, NY 10065, United States; Tri-INSTITUTIONAL PhD PrOGRAM IN 8 2 Chemical Biology, WEILL Cornell GrADUATE School OF Medical Sciences, Cornell University, NeW York, 9 NY 10065, United States; Pharmaceutical Sciences, MRL, MerCK & Co., Inc., 126 East Lincoln 10 3 Avenue, Rahway, NeW JerseY 07065, United States; GrADUATE PrOGRAM IN Physiology, Biophysics, 11 4 AND Systems Biology, WEILL Cornell Medical College, NeW York, NY 10065, United States; PrOCESS 12 5 AND Analytical ResearCH AND Development, MerCK & Co., Inc., Rahway, NJ 07065, United States; 13 Analytical ResearCH & Development, MRL, MerCK & Co., Inc., MRL, 126 East Lincoln Avenue, 14 6 Rahway, NeW JerseY 07065, United States; Department OF Pharmaceutical Sciences AND 15 7 Department OF Chemistry, University OF California, Irvine, Irvine, California 92697, United States 16 17 *For CORRespondence: 18 [email protected] (TR); [email protected] (JDC) 19 20 AbstrACT Determining THE NET CHARGE AND PROTONATION STATES POPULATED BY A SMALL MOLECULE IN AN envi- 21 RONMENT OF INTEREST—SUCH AS solvent, A PROTEIN BINDING site, OR A LIPID BILAYER—OR THE COST OF ALTERING THOSE 22 PROTONATION STATES UPON TRANSFER TO ANOTHER ENVIRONMENT IS A PREREQUISITE FOR PREDICTING ITS PHYSICOCHEMICAL 23 AND PHARMACEUTICAL PRoperties, AS WELL AS INTERACTIONS WITH BIOLOGICAL MACROMOLECULES USING COMPUTATIONAL 24 models. INCORRECTLY MODELING THE DOMINANT PROTONATION state, SHIFTS IN DOMINANT PROTONATION state, OR THE 25 POPULATION OF SIGNIfiCANT MIXTURES OF PROTONATION STATES CAN LEAD TO LARGE MODELING ERRORS THAT DEGRADE THE 26 ACCURACY OF PHYSICAL MODELING AND HINDER THE ABILITY TO USE PHYSICAL MODELING APPROACHES FOR MOLECULAR 27 design. For SMALL molecules, THE ACID DISSOCIATION CONSTANT (pK A) IS THE PRIMARY QUANTITY NEEDED TO deter- 28 MINE THE IONIC STATES POPULATED BY A MOLECULE IN AN AQUEOUS SOLUTION AT A GIVEN pH. As A PART OF SAMPL6 29 COMMUNITY challenge, WE ORGANIZED A BLIND PK A PREDICTION COMPONENT TO ASSESS THE ACCURACY WITH WHICH 30 CONTEMPORARY PK A PREDICTION METHODS CAN PREDICT THIS QUANTITY, WITH THE ULTIMATE AIM OF ASSESSING THE 31 EXPECTED IMPACT ON MODELING ERRORS THIS WOULD induce. While A MULTITUDE OF APPROACHES FOR PREDICTING PK A 32 VALUES CURRENTLY Exist, PREDICTING THE PK AS OF drug-likE MOLECULES CAN BE DIffiCULT DUE TO CHALLENGING PROPERTIES 33 SUCH AS MULTIPLE TITRATABLE sites, HETERocycles, AND tautomerization. For THIS challenge, WE FOCUSED ON SET 34 OF 24 SMALL MOLECULES SELECTED TO RESEMBLE SELECTIVE KINASE INHIBITORS—AN IMPORTANT CLASS OF THERAPEUTICS 35 REPLETE WITH TITRATABLE moieties. Using A Sirius T3 INSTRUMENT THAT PERFORMS AUTOMATED acid-base TITRations, 36 WE USED UV absorbance-based PK A MEASUREMENTS TO CONSTRUCT A high-quality EXPERIMENTAL REFERENCE DATASET 37 OF MACROSCOPIC PK AS FOR THE EVALUATION OF COMPUTATIONAL PK A PREDICTION METHODOLOGIES THAT WAS UTILIZED IN 38 THE SAMPL6 PK A challenge. For SEVERAL COMPOUNDS IN WHICH THE MICROSCOPIC PROTONATION STATES ASSOCIATED 39 WITH MACROSCOPIC PK AS WERE ambiguous, WE PERFORMED FOLLOw-up NMR EXPERIMENTS TO DISAMBIGUATE THE 40 MICROSTATES INVOLVED IN THE TRansition. This DATASET PROVIDES A USEFUL STANDARD BENCHMARK DATASET FOR THE 41 EVALUATION OF PK A PREDICTION METHODOLOGIES ON KINASE inhibitor-likE compounds. 1 OF 27 bioRxiv preprint doi: https://doi.org/10.1101/368787; this version posted July 13, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. PrEPRINT AHEAD OF SUBMISSION — July 13, 2018 42 43 KeYWORDS ⋅ ⋅ ⋅ ⋅ 44 ACID DISSOCIATION CONSTANTS SPECTROPHOTOMETRIC PK A MEASUREMENT BLIND PREDICTION CHALLENGE SAMPL ⋅ ⋅ ⋅ 45 MACROSCOPIC PK A MICROSCOPIC PK A MACROSCOPIC PROTONATION STATE MICROSCOPIC PROTONATION STATE 46 AbbrEVIATIONS 47 SAMPL Statistical Assessment OF THE Modeling OF PrOTEINS AND Ligands log 48 PK A – 10 ACID DISSOCIATION EQUILIBRIUM CONSTANT log 49 PSK A – 10 APPARENT ACID DISSOCIATION EQUILIBRIUM CONSTANT IN COSOLVENT 50 DMSO Dimethyl SULFOXIDE 51 ISA Ionic-strENGTH ADJUSTED 52 SEM StandarD ERROR OF THE MEAN 53 TFA TARGET FACTOR ANALYSIS 54 LC-MS Liquid CHROMATOGRAPHY - MASS SPECTROMETRY 55 NMR Nuclear MAGNETIC RESONANCE SPECTROSCOPY 56 HMBC HeterONUCLEAR Multiple-Bond CorrELATION 57 TFA-D DEUTERo-triflUOROACETIC ACID INTRODUCTION 58 59 SAMPL (Statistical Assessment OF THE Modeling OF PrOTEINS AND Ligands) IS A RECURRING SERIES OF BLIND PREDICTION 60 CHALLENGES FOR THE COMPUTATIONAL CHEMISTRY COMMUNITY [1, 2]. ThrOUGH THESE challenges, SAMPL AIMS TO 61 EVALUATE AND ADVANCE COMPUTATIONAL TOOLS FOR RATIONAL DRUG design. By FOCUSING THE COMMUNITY ON SPECIfiC 62 PHENOMENA RELEVANT TO DRUG DISCOVERY POORLY PREDICTED BY CURRENT models, ISOLATING THAT PHENOMENON 63 FROM OTHER CONFOUNDING FACTORS IN well-designed TEST systems, EVALUATING TOOLS PROSPECTIVELY, ENABLING DATA 64 SHARING TO LEARN FROM FAILURes, AND RELEASING THE RESULTING high-quality DATASETS INTO THE COMMUNITY AS 65 BENCHMARK sets, SAMPL HAS DRIVEN PROGRESS IN A NUMBER OF AREAS OVER SEVEN PREVIOUS ROUNDS OF CHALLENGE 66 CYCLES [3–7, 7–15]. 67 As A STEPPING STONE TO ENABLING THE ACCURATE PREDICTION OF PRotein-ligand BINDING Affinities, SAMPL 68 HAS FOCUSED ON EVALUATING HOW WELL PHYSICAL AND EMPIRICAL MODELING METHODOLOGIES CAN PREDICT VARIOUS 69 PHYSICOCHEMICAL PROPERTIES RELEVANT TO BINDING AND DRUG DISCOVERY, SUCH AS HYDRATION FREE ENERGIES (which 70 MODEL ASPECTS OF DESOLVATION IN isolation), DISTRIBUTION COEffiCIENTS (which MODEL TRANSFER FROM RELATIVELY 71 HOMOGENEOUS AQUEOUS TO NONPOLAR ENVIRonments), AND host-guest BINDING AffiNITIES (which MODEL high- 72 AffiNITY ASSOCIATION WITHOUT THE COMPLICATION OF SLOW PROTEIN dynamics). These PHYSICOCHEMICAL PROPERTY 73 PREDICTION CHALLENGES—IN ADDITION TO ASSESSING THE PREDICTIVE ACCURACY OF QUANTITIES THAT ARE USEFUL IN VARIOUS 74 STAGES OF DRUG DISCOVERY IN THEIR OWN RIGHT—HAVE BEEN HELPFUL IN PINPOINTING DEfiCIENCIES IN COMPUTATIONAL 75 MODELS THAT CAN LEAD TO SUBSTANTIAL ERRORS IN AffiNITY PRedictions. Neglect OF PROTONATION STATE EffECTS CAN LEAD TO LARGE MODELING ERRORS 76 77 As PART OF THE SAMPL5 CHALLENGE series, A NEW CYCLOHExane-water DISTRIBUTION CONSTANT (log D) PREDICTION 78 CHALLENGE WAS INTRoduced, WHERE PARTICIPANTS PREDICTED THE TRANSFER FREE ENERGY OF SMALL drug-likE MOLECULES 79 BETWEEN AN AQUEOUS BUffER PHASE AT pH 7.4 AND A NONAQUEOUS CYCLOHEXANE PHASE [16, 17]. While octanol- 80 WATER DISTRIBUTION COEffiCIENT MEASUREMENTS ARE MORE common, CYCLOHEXANE WAS SELECTED FOR THE SIMPLICITY 81 OF ITS LIQUID PHASE AND RELATIVE DRYNESS COMPARED TO WET OCTANOL phases. While THE EXPECTATION WAS THAT 82 THIS CHALLENGE WOULD BE RELATIVELY STRAIGHTFORWARD GIVEN THE LACK OF COMPLEXITY OF CYCLOHEXANE phases, 83 ANALYSIS OF PARTICIPANT PERFORMANCE REVEALED THAT MULTIPLE FACTORS CONTRIBUTED TO SIGNIfiCANT PREDICTION 84 FAILURes: POOR CONFORMATIONAL SAMPLING OF flEXIBLE SOLUTE molecules, MISPREDICTION OF RELEVANT PROTONATION 85 AND TAUTOMERIC STATES (or FAILURE TO ACCOMMODATE SHIFTS IN THEIR populations), AND FORCE fiELD INACCURACIES 86 RESULTING IN BIAS TOWARDS THE CYCLOHEXANE phase. While THESE fiNDINGS JUSTIfiED THE BENEfiT OF FUTURE ITERATIONS 87 OF BLIND DISTRIBUTION OR PARTITION COEffiCIENT challenges, THE MOST SURPRISING OBSERVATION FROM THIS INITIAL LOG D 88 CHALLENGE WAS THAT PARTICIPANTS ALMOST UNIFORMLY NEGLECTED TO ACCURATELY MODEL PROTONATION STATE Effects, 2 OF 27 bioRxiv preprint doi: https://doi.org/10.1101/368787; this version posted July 13, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. PrEPRINT AHEAD OF SUBMISSION — July 13, 2018 89 AND THAT NEGLECT OF THESE EffECTS LED TO SURPRISINGLY LARGE ERRORS IN TRANSFER FREE ENERGIES [16–18]. CarEFUL 90 QUANTUM CHEMICAL ASSESSMENTS OF THE MAGNITUDE OF THESE PROTONATION STATE EffECTS FOUND THAT THEIR NEGLECT 91 COULD INTRODUCE ERRORS UP TO 6–8 Kcal/mol FOR SOME COMPOUNDS [18]. This EffECT STEMS FROM THE NEED TO 92 ACCOUNT FOR THE FREE ENERGY DIffERENCE BETWEEN THE MAJOR IONIZATION STATE IN CYCLOHEXANE (most LIKELY NEUTRAL 93 state) AND IN WATER PHASE (which COULD BE NEUTRAL OR CHARged). 94 TO ISOLATE THESE SURPRISINGLY LARGE PROTONATION STATE MODELING ERRORS FROM DIffiCULTIES RELATED TO lipophilic- 95 ITY (log P AND LOG D) PREDICTION methods, WE DECIDED TO ORGANIZE A SET OF STAGED PHYSICOCHEMICAL PROPERTY 96 CHALLENGES USING A CONSISTENT SET OF MOLECULES THAT RESEMBLE SMALL MOLECULE KINASE INHIBITORS—AN IMPORTANT 97 DRUG CLASS REPLETE
Recommended publications
  • Evaluating Small Molecule Microscopic and Macroscopic Pka
    bioRxiv preprint doi: https://doi.org/10.1101/2020.10.15.341792; this version posted October 15, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. 1 Overview of the SAMPL6 pK a Challenge: 2 Evaluating small molecule microscopic and 3 macroscopic pK a predictions 4 Mehtap Işık (ORCID: 0000-0002-6789-952X)1,2*, Ariën S. Rustenburg (ORCID: 0000-0002-3422-0613)1,3, Andrea 5 Rizzi (ORCID: 0000-0001-7693-2013)1,4, M. R. Gunner (ORCID: 0000-0003-1120-5776)6, David L. Mobley (ORCID: 6 0000-0002-1083-5533)5, John D. Chodera (ORCID: 0000-0003-0542-119X)1 7 1Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, 8 New York, NY 10065, United States; 2Tri-Institutional PhD Program in Chemical Biology, Weill Cornell Graduate 9 School of Medical Sciences, Cornell University, New York, NY 10065, United States; 3Graduate Program in 10 Physiology, Biophysics, and Systems Biology, Weill Cornell Medical College, New York, NY 10065, United States; 11 4Tri-Institutional PhD Program in Computational Biology and Medicine, Weill Cornell Graduate School of Medical 12 Sciences, Cornell University, New York, NY 10065, United States; 5Department of Pharmaceutical Sciences and 13 Department of Chemistry, University of California, Irvine, Irvine, California 92697, United States; 6Department of 14 Physics, City College of New York, New York NY 10031 15 *For correspondence: 16 [email protected] (MI) 17 18 Abstract 19 K The prediction of acid dissociation constants (p a) is a prerequisite for predicting many other properties of a small molecule, 20 such as its protein-ligand binding affinity, distribution coefficient (log D), membrane permeability, and solubility.
    [Show full text]
  • Accurate Receptor-Ligand Binding Free Energies from Fast QM Conformational Chemical Space Sampling
    International Journal of Molecular Sciences Article Accurate Receptor-Ligand Binding Free Energies from Fast QM Conformational Chemical Space Sampling Esra Boz and Matthias Stein * Max Planck Institute for Dynamics of Complex Technical Systems, Molecular Simulations and Design Group, 39106 Magdeburg, Germany; [email protected] * Correspondence: [email protected] Abstract: Small molecule receptor-binding is dominated by weak, non-covalent interactions such as van-der-Waals hydrogen bonding or electrostatics. Calculating these non-covalent ligand-receptor interactions is a challenge to computational means in terms of accuracy and efficacy since the ligand may bind in a number of thermally accessible conformations. The conformational rotamer ensemble sampling tool (CREST) uses an iterative scheme to efficiently sample the conformational space and calculates energies using the semi-empirical ‘Geometry, Frequency, Noncovalent, eXtended Tight Binding’ (GFN2-xTB) method. This combined approach is applied to blind predictions of the modes and free energies of binding for a set of 10 drug molecule ligands to the cucurbit[n]urils CB[8] receptor from the recent ‘Statistical Assessment of the Modeling of Proteins and Ligands’ (SAMPL) challenge including morphine, hydromorphine, cocaine, fentanyl, and ketamine. For each system, the conformational space was sufficiently sampled for the free ligand and the ligand-receptor complexes using the quantum chemical Hamiltonian. A multitude of structures makes up the final conformer-rotamer ensemble, for which then free energies of binding are calculated. For those large and complex molecules, the results are in good agreement with experimental values with a mean Citation: Boz, E.; Stein, M. Accurate error of 3 kcal/mol.
    [Show full text]
  • 20180921 Michel SAMPL6-Ms-Revised
    Edinburgh Research Explorer Blinded predictions of standard binding free energies: lessons learned from the SAMPL6 challenge Citation for published version: Papadourakis, M, Bosisio, S & Michel, J 2018, 'Blinded predictions of standard binding free energies: lessons learned from the SAMPL6 challenge', Journal of computer-Aided molecular design. https://doi.org/10.1007/s10822-018-0154-6 Digital Object Identifier (DOI): 10.1007/s10822-018-0154-6 Link: Link to publication record in Edinburgh Research Explorer Document Version: Peer reviewed version Published In: Journal of computer-Aided molecular design General rights Copyright for the publications made accessible via the Edinburgh Research Explorer is retained by the author(s) and / or other copyright owners and it is a condition of accessing these publications that users recognise and abide by the legal requirements associated with these rights. Take down policy The University of Edinburgh has made every reasonable effort to ensure that Edinburgh Research Explorer content complies with UK legislation. If you believe that the public display of this file breaches copyright please contact [email protected] providing details, and we will remove access to the work immediately and investigate your claim. Download date: 08. Oct. 2021 Noname manuscript No. (will be inserted by the editor) Blinded predictions of standard binding free energies: lessons learned from the SAMPL6 challenge Michail Papadourakis, Stefano Bosisio, Julien Michel* EaStCHEM School of Chemistry, Joseph Black Building, King’s Buildings, Edinburgh EH9 3FJ, UK. Abstract In the context of the SAMPL6 challenges, series of blinded predictions of standard binding free energies were made with the SOMD software for a dataset of 27 host-guest systems featuring two octa-acids hosts (OA and TEMOA) and a cucurbituril ring (CB8) host.
    [Show full text]
  • Anew Drug Design Strategy in the Liht of Molecular Hybridization Concept
    www.ijcrt.org © 2020 IJCRT | Volume 8, Issue 12 December 2020 | ISSN: 2320-2882 “Drug Design strategy and chemical process maximization in the light of Molecular Hybridization Concept.” Subhasis Basu, Ph D Registration No: VB 1198 of 2018-2019. Department Of Chemistry, Visva-Bharati University A Draft Thesis is submitted for the partial fulfilment of PhD in Chemistry Thesis/Degree proceeding. DECLARATION I Certify that a. The Work contained in this thesis is original and has been done by me under the guidance of my supervisor. b. The work has not been submitted to any other Institute for any degree or diploma. c. I have followed the guidelines provided by the Institute in preparing the thesis. d. I have conformed to the norms and guidelines given in the Ethical Code of Conduct of the Institute. e. Whenever I have used materials (data, theoretical analysis, figures and text) from other sources, I have given due credit to them by citing them in the text of the thesis and giving their details in the references. Further, I have taken permission from the copyright owners of the sources, whenever necessary. IJCRT2012039 International Journal of Creative Research Thoughts (IJCRT) www.ijcrt.org 284 www.ijcrt.org © 2020 IJCRT | Volume 8, Issue 12 December 2020 | ISSN: 2320-2882 f. Whenever I have quoted written materials from other sources I have put them under quotation marks and given due credit to the sources by citing them and giving required details in the references. (Subhasis Basu) ACKNOWLEDGEMENT This preface is to extend an appreciation to all those individuals who with their generous co- operation guided us in every aspect to make this design and drawing successful.
    [Show full text]
  • Understanding TRPV1 Activation by Ligands: Insights from the Binding
    Understanding TRPV1 activation by ligands: Insights PNAS PLUS from the binding modes of capsaicin and resiniferatoxin Khaled Elokelya,b,c, Phanindra Velisettyd, Lucie Delemottea, Eugene Palovcaka, Michael L. Kleina,b,1, Tibor Rohacsd, and Vincenzo Carnevalea,1 aInstitute for Computational Molecular Science, Temple University, Philadelphia, PA 19122; bDepartment of Chemistry, Temple University, Philadelphia, PA 19122; cDepartment of Pharmaceutical Chemistry, Tanta University, 31527 Tanta, Egypt; and dDepartment of Pharmacology, Physiology and Neuroscience, Rutgers–New Jersey Medical School, Newark, NJ 07103 Contributed by Michael L. Klein, December 7, 2015 (sent for review November 6, 2015; reviewed by Kenton J. Swartz and Vladimir Yarov-Yarovoy) The transient receptor potential cation channel subfamily V member TRPV1 is known to be the target of capsaicin (CAPS), the active 1 (TRPV1) or vanilloid receptor 1 is a nonselective cation channel that component of chili peppers, and it can also be referred to as the is involved in the detection and transduction of nociceptive stimuli. capsaicin receptor (18). Resiniferatoxin (RTX), a phorbol ester Inflammation and nerve damage result in the up-regulation of TRPV1 isolated from the irritant lattices of the Moroccan cactus, shows a transcription, and, therefore, modulators of TRPV1 channels are much higher affinity for TRPV1 than CAPS (19). Both compounds potentially useful in the treatment of inflammatory and neuropathic activate TRPV1, causing the channel to be more permeable to pain. Understanding the binding modes of known ligands would cations, ultimately resulting in an analgesic effect due to channel significantly contribute to the success of TRPV1 modulator drug desensitization. CAPS can be subdivided into three structural re- design programs.
    [Show full text]
  • Survey Results
    Results of the 2017 Roadmap survey of the Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) challenge community David L. Mobley1, John D. Chodera2 and Michael K. Gilson3 1Departments of Pharmaceutical Sciences and Chemistry, University of California, Irvine 2Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center 3Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego Abstract: The Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) series of blind prediction challenges provide unbiased, prospective tests of computational methods, serving as a crowdsourcing mechanism to drive innovation. These challenges solicit submissions from the computational community to predict properties of interest, such as solvation, partition, or binding of drug-like molecules. These challenges provided substantial benefit to the community, and have led to roughly 100 publications, many of which are broadly cited (see attached bibliography). We are currently seeking funding from the NIH and surveyed the community concerning experiences with SAMPL and how our future plans for SAMPL can best align with the community’s interests and needs. This document summarizes the results of this survey and describes our findings. On the whole, the community enthusiastically supports our plans for the future of SAMPL, and provided modest suggestions to further strengthen our plans. For up-to-date info please see the SAMPL website. Survey methods and results Here, we reproduce the results of the survey in full except for anonymizing respondents. Questions 1 and 2 dealt with identifying information (name and e-mail) and are thus bypassed here. There were 44 respondents, though not all respondents answered all questions.
    [Show full text]
  • Multipose Binding in Molecular Docking
    Int. J. Mol. Sci. 2014, 15, 2622-2645; doi:10.3390/ijms15022622 OPEN ACCESS International Journal of Molecular Sciences ISSN 1422-0067 www.mdpi.com/journal/ijms Article Multipose Binding in Molecular Docking Kalina Atkovska 1,†, Sergey A. Samsonov 1,†,*, Maciej Paszkowski-Rogacz 2 and M. Teresa Pisabarro 1,* 1 Structural Bioinformatics, BIOTEC TU Dresden, Tatzberg 47-51, Dresden 01307, Germany; E-Mail: [email protected] 2 Department of Medical Systems Biology, TU Dresden, Medical Faculty, University Hospital Carl Gustav Carus, Dresden 01307, Germany; E-Mail: [email protected] † These authors contributed equally to this work. * Authors to whom correspondence should be addressed; E-Mails: [email protected] (M.T.P.); [email protected] (S.A.S.); Tel.: +49-351-463-40071 (M.T.P. & S.A.S.); Fax: +49-351-463-40087 (M.T.P. & S.A.S.). Received: 18 December 2013; in revised form: 17 January 2014 / Accepted: 21 January 2014 / Published: 14 February 2014 Abstract: Molecular docking has been extensively applied in virtual screening of small molecule libraries for lead identification and optimization. A necessary prerequisite for successful differentiation between active and non-active ligands is the accurate prediction of their binding affinities in the complex by use of docking scoring functions. However, many studies have shown rather poor correlations between docking scores and experimental binding affinities. Our work aimed to improve this correlation by implementing a multipose binding concept in the docking scoring scheme. Multipose binding, i.e., the property of certain protein-ligand complexes to exhibit different ligand binding modes, has been shown to occur in nature for a variety of molecules.
    [Show full text]
  • Overview of the SAMPL5 Host-Guest Challenge: Are We Doing Better? Jian Yin†, Niel M
    Overview of the SAMPL5 Host-Guest Challenge: Are We Doing Better? Jian Yin†, Niel M. Henriksen†, David R. Slochower†, Michael R. Shirts‡, Michael W. ¶ † * Chiu , David L. Mobleyѱ and Michael K. Gilson †Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, La Jolla, CA 92093, USA ‡Department of Chemical Engineering, University of Virginia, Charlottesville, Virginia 22904, USA ¶Qualcomm Institute, University of California, San Diego, La Jolla, CA 92093, USA ѱDepartments of Pharmaceutical Sciences and Chemistry, University of California Irvine, Irvine, CA 92697, USA Keywords Host-guest, molecular recognition, computer-aided drug design, Blind challenge, Binding affinity *Corresponding author: Prof. Michael K. Gilson Email: [email protected] Telephone: 858-822-0622 Fax: 858-822-7726 Abstract The ability to computationally predict protein-small molecule Binding affinities with high accuracy would accelerate drug discovery and reduce its cost By eliminating rounds of trial-and- error synthesis and experimental evaluation of candidate ligands. As academic and industrial groups work toward this capability, there is an ongoing need for datasets that can Be used to rigorously test new computational methods. Although protein-ligand data are clearly important for this purpose, their size and complexity make it difficult to oBtain well-converged results and to trouBleshoot computational methods. Host-guest systems offer a valuable alternative class of test cases, as they exemplify noncovalent molecular recognition But are far smaller and simpler. As a consequence, host-guest systems have Been part of the prior two rounds of SAMPL prediction exercises, and they also figure in the present SAMPL5 round. In addition to Being Blinded, and thus avoiding Biases that may arise in retrospective studies, the SAMPL challenges have the merit of focusing multiple researchers on a common set of molecular systems, so that methods may Be compared and ideas exchanged.
    [Show full text]
  • Biomolecular Modeling and Simulation: a Prospering Multidisciplinary Field
    Annual Review of Biophysics Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field Tamar Schlick,1,2,3 Stephanie Portillo-Ledesma,1 Christopher G. Myers,1 Lauren Beljak,4 Justin Chen,4 Sami Dakhel,4 Daniel Darling,4 Sayak Ghosh,4 Joseph Hall,4 Mikaeel Jan,4 Emily Liang,4 Sera Saju,4 Mackenzie Vohr,4 Chris Wu,4 Yifan Xu,4 and Eva Xue4 1Department of Chemistry, New York University, New York, New York 10003, USA; email: [email protected] 2Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA 3New York University–East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200122, China 4College of Arts and Science, New York University, New York, New York 10003, USA Annu. Rev. Biophys. 2021. 50:267–301 Keywords First published as a Review in Advance on biomolecular modeling, biomolecular simulation, biomolecular dynamics, February 19, 2021 structure prediction, protein folding, DNA folding, RNA folding, The Annual Review of Biophysics is online at multiscale modeling, citizen science projects, machine learning, artificial biophys.annualreviews.org intelligence https://doi.org/10.1146/annurev-biophys-091720- 102019 Abstract Copyright © 2021 by Annual Reviews. We reassess progress in the field of biomolecular modeling and simulation, Annu. Rev. Biophys. 2021.50:267-301. Downloaded from www.annualreviews.org All rights reserved following up on our perspective published in 2011. By reviewing metrics Access provided by New York University - Bobst Library on 05/12/21. For personal use only. for the field’s productivity and providing examples of success, we under- score the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated.
    [Show full text]
  • Octanol-Water Partition Coefficient Measurements for the SAMPL6
    bioRxiv preprint doi: https://doi.org/10.1101/757393; this version posted September 8, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. PrEPRINT AHEAD OF SUBMISSION — September 8, 2019 Octanol-water PARTITION COEffiCIENT 1 MEASUREMENTS FOR THE SAMPL6 Blind PrEDICTION 2 Challenge 3 1,2 3 4 Mehtap Işık (ORCID: 0000-0002-6789-952X) , DorOTHY LeVORSE , David L. MobleY (ORCID: 4 3* 5 0000-0002-1083-5533) , Timothy Rhodes (ORCID: 0000-0001-7534-9221) , John D. ChoderA (ORCID: 1* 6 0000-0003-0542-119X) Computational AND Systems Biology PrOGRam, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, 7 1 NeW York, NY 10065, United States; Tri-INSTITUTIONAL PhD PrOGRAM IN Chemical Biology, WEILL Cornell GrADUATE School 8 2 OF Medical Sciences, Cornell University, NeW York, NY 10065, United States; Pharmaceutical Sciences, MRL, MerCK & 9 3 Co., Inc., 126 East Lincoln Avenue, Rahway, NeW JerseY 07065, United States; Department OF Pharmaceutical 10 4 Sciences AND Department OF Chemistry, University OF California, Irvine, Irvine, California 92697, United States 11 12 *For CORRespondence: 13 [email protected] (TR); [email protected] (JDC) 14 15 AbstrACT Partition COEffiCIENTS DESCRIBE THE EQUILIBRIUM PARTITIONING OF A single, DEfiNED CHARGE STATE OF A SOLUTE BETWEEN 16 TWO IMMISCIBLE phases, TYPICALLY A NEUTRAL solute. Octanol-water PARTITION COEffiCIENTS (Kow), OR THEIR LOGARITHMS (log P), ARE 17 FREQUENTLY USED AS A MEASURE OF LIPOPHILICITY IN DRUG DISCOVERY.
    [Show full text]
  • UC Irvine Mobley Lab Datasets
    UC Irvine Mobley Lab Datasets Title Advancing predictive modeling through focused development of model systems to drive new modeling innovations Permalink https://escholarship.org/uc/item/7cf8c6cr Authors Mobley, David L Chodera, John D Isaacs, Lyle et al. Publication Date 2016-10-05 License https://creativecommons.org/licenses/by-nc-nd/4.0/ 4.0 eScholarship.org Powered by the California Digital Library University of California Advancing predictive modeling through focused development of model systems to drive new modeling innovations A funding proposal recently submitted to the NIH David L. Mobley Associate Professor, Departments of Pharmaceutical Sciences and Chemistry, University of California, Irvine John D. Chodera Assistant Faculty Member, Computational Biology Program, Memorial Sloan-Kettering Cancer Center Lyle Isaacs Professor, Department of Chemistry and Biochemistry, University of Maryland Bruce C. Gibb Professor, Department of Chemistry, Tulane University ABSTRACT This work seeks to advance quantitative methods for biomolecular design, especially for predicting biomolecular interactions, via a focused series of community blind prediction challenges. Physical methods for predicting binding free energies, or “free energy methods”, are poised to dramatically reshape early stage drug discovery, and are already finding applications in pharmaceutical lead optimization. However, performance is unreliable, the domain of applicability is limited, and failures in pharmaceutical applications are often hard to understand and fix. On the other hand, these methods can now typically predict a variety of simple physical properties such as solvation free energies or relative solubilities, though there is still clear room for improvement in accuracy. In recent years, competitions and crowdsourcing have proven an effective model for driving innovations in diverse fields.
    [Show full text]
  • Overview of the SAMPL6 Host-Guest Binding Affinity Prediction Challenge
    bioRxiv preprint doi: https://doi.org/10.1101/371724; this version posted July 19, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. PrEPRINT AHEAD OF SUBMISSION — July 18, 2018 OvervieW OF THE SAMPL6 host-guest 1 BINDING AffiNITY PREDICTION CHALLENGE 2 1,2 3 3 4 4 5 3 AndrEA Rizzi , SteVEN Murkli , John N. McNeill , WEI YAO , MattheW Sullivan , Michael K. Gilson , Michael W. 6 3 4 7* 1* 4 Chiu , Lyle ISAACS , Bruce C. Gibb , David L. MobleY , John D. ChoderA 1 5 Computational AND Systems Biology PrOGRam, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, NeW 2 6 York, NY 10065, USA; Tri-INSTITUTIONAL TRAINING PrOGRAM IN Computational Biology AND Medicine, NeW York, NY 10065, USA; 3 4 7 Department OF Chemistry AND Biochemistry, University OF Maryland, College Park, MD 20742, USA; Department OF 5 8 Chemistry, TULANE University, Louisiana, LA 70118, USA; Skaggs School OF Pharmacy AND Pharmaceutical Sciences, 6 9 University OF California, San Diego, La Jolla, CA 92093, USA; Qualcomm Institute, University OF California, San Diego, La 7 10 Jolla, CA 92093, USA; Department OF Pharmaceutical Sciences AND Department OF Chemistry, University OF California, 11 Irvine, California 92697, USA 12 *For CORRespondence: 13 [email protected] (DLM); [email protected] (JDC) 14 15 AbstrACT The ABILITY TO ACCURATELY PREDICT THE BINDING AffiNITIES OF SMALL ORGANIC MOLECULES TO BIOLOGICAL 16 MACROMOLECULES WOULD GREATLY ACCELERATE DRUG DISCOVERY BY REDUCING THE NUMBER OF COMPOUNDS THAT MUST 17 BE SYNTHESIZED TO REALIZE DESIRED POTENCY AND SELECTIVITY goals.
    [Show full text]