Graph convolutional and message-passing networks can be a powerful tool for predicting physical properties of small molecules when coupled to a simple physical model that encodes the relevant invariances. Here, we show the ability of graph nets to predict partial atomic charges for use in molecular dynamics simulations and physical docking.
Abraham MJ, Apostolov R, Barnoud J, Bauer P, Blau C, Bonvin AMJJ, Chavent M, Chodera JD, Condic-Jurkic K, Delemotte L, Grubmüller H, Howard RJ, Lindahl E, Ollila S, Salent J, Smith D, Stansfeld PJ, Tiemann J, Trellet M, Woods C, and Zhmurov A.
Preprint ahead of publication. [chemRxiv]
There is a dire need to establish standards for sharing data in the molecular sciences. Here, we review the findings of a workshop held in Stockholm in Nov 2018 to discuss this need.
We describe the design and data collection (and associated challenges) for the SAMPL6 part II logP octanol-water blind prediction challenge, where the goal was to benchmark the accuracy of force fields for druglike molecules (here, molecules resembling kinase inhibitors).
Dajun Sang, Sudarshan Pinglay, Rafal P Wiewiora, Myvizhi E Selvan, Hua Jane Lou, John D Chodera, Benjamin E Turk, Zeynep H Gümüş, and Liam J Holt. eLife 2019;8:e38805 [DOI] [eLife] [PDF] [Folding@home data]
To understand how kinase regulation by phosphorylation emerged, we reconstruct the common ancestor of CDKs and MAPKs, using biochemical experiments and massively parallel molecular simulations to study how a few mutations were sufficient to switch ERK-family kinases from high- to low-autophosphorylation.
We assess the accuracy of the SMIRNOFF99Frosst 1.0.5 force field in reproducing host-guest binding thermodynamics in comparison with the GAFF force field, demonstrating how the SMIRNOFF format for compactly specifying force fields provide comparable accuracy with 20x fewer parameters.
Gerard Minuesa, Steven K Albanese, Arthur Chow, Alexandra Schurer, Sun-Mi Park, Christina Z. Rotsides, James Taggart, Andrea Rizzi, Levi N. Naden, Timothy Chou, Saroj Gourkanti, Daniel Cappel, Maria C Passarelli, Lauren Fairchild, Carolina Adura, Fraser J Glickman, Jessica Schulman, Christopher Famulare, Minal Patel, Joseph K Eibl, Gregory M Ross, Derek S Tan, Christina S Leslie, Thijs Beuming, Yehuda Goldgur, John D Chodera, Michael G Kharas
Nature Communications 10:2691, 2019. [DOI] [bioRxiv] [GitHub] [MSKCC blog post]
We use absolute alchemical free energy calculations to identify the likely interaction site for a small hydrophobic ligand that shows activity against MUSASHI in AML.
Rafal P. Wiewiora*, Shi Chen*, Fanwang Meng, Nicolas Babault, Anqi Ma, Wenyu Yu, Kun Qian, Hao Hu, Hua Zou, Junyi Wang, Shijie Fan, Gil Blum, Fabio Pittella-Silva, Kyle A. Beauchamp, Wolfram Tempel, Hualing Jiang, Kaixian Chen, Robert Skene, Y. George Zheng, Peter J. Brown, Jian Jin, John D. Chodera+, and Minkui Luo+.
eLife 8:e45403, 2019. [DOI] [bioRxiv] [GitHub] [OSF] [movies] [MSKCC blog post]
* These authors contributed equally to this work
+ Co-corresponding authors
In this work, we show how targeted X-ray crystallography using covalent inhibitors and depletion of native ligands to reveal structures of low-population hidden conformations can be combined with massively distributed molecular simulation to resolve the functional dynamic landscape of the protein methyltransferase SETD8 in unprecedented atomistic detail. Using an aggregate of six milliseconds of fully atomistic simulation from Folding@home, we use Markov state models to illuminate the conformational dynamics of this important epigenetic protein.
The trajectories generated for this project were used as the source for a unique musical composition 'Metastable' by George Holloway, performed by the Ligeti String Quartet with visual accompaniment from Robert Arbon.
A class of kinases are particularly promiscuous binders of small molecule inhibitors. Using combined biomolecular simulations and biochemical studies, we show that the promiscuity of DDR1, one of the major members of this class, is likely due to an unusually stable DFG-out conformation.
To make powerful path sampling techniques broadly accessible and efficient, we have produced a new Python framework for easily implementing path sampling strategies (such as transition path and interface sampling) in Python. This second publication describes advanced aspects of the theory and details of how to customize path ensembles.
Paulina M. Wojnarowicz, Raquel Lima e Silva, Masayuki Ohanka, Sang Bae Lee, Yvette Chin, Anita Kulukian, Sung-Hee Chang, Bina Desai, Marta Garcia Escolano, Riddhi Shah, Marta Garcia-Cao, Sijia Xu, Rashmi Thakar, Yehuda Goldgur, Meredith A. Miller, Ouathek Ouerfelli, Guangli Yang, Tsutomu Arakawa, Steven K. Albanese, William A. Garland, Glenn Stoller, Jaideep Chaudhary, Rajesh Soni, John Philip, Ronald C. Hendrickson, Antonio Iavarone, Andrew J. Dannenberg, John D. Chodera, Nikola Pavletich, Anna Lasorella, Peter A. Campochiaro, and Robert Benezra
We report the discovery and characterization of a small molecule, AGX51, with the surprising ability to inhibit the interaction of Id1 with E47, which leads to ubiquitin-mediated degradation of Ids.
To make powerful path sampling techniques broadly accessible and efficient, we have produced a new Python framework for easily implementing path sampling strategies (such as transition path and interface sampling) in Python. This first publication describes some of the theory and capabilities behind the approach.
Camila Zanette, Caitlin C. Bannan, Christopher I. Bayly, Josh Fass, Michael K. Gilson, Michael R. Shirts, John Chodera, and David L. Mobley
Journal of Chemical Theory and Computation, 15:402, 2019. [DOI] [ChemRxiv] [GitHub]
We show how machine learning can learn typing rules for molecular mechanics force fields within a Bayesian statistical framework.
Andrea Rizzi, Steven Murkli, John N. McNeill, Wei Yao, Matthew Sullivan, Michael K. Gilson, Michael W. Chiu, Lyle Isaacs, Bruce C. Gibb, David L. Mobley*, John D. Chodera*
* denotes co-corresponding authors
Journal of Computer-Aided Molecular Design special issue on SAMPL6, 32:937, 2018. [DOI] [bioRxiv] [GitHub]
We present an overview of the host-guest systems and participant performance for the SAMPL6 host-guest blind affinity prediction challenges, assessing how well various physical modeling approaches were able to predict ligand binding affinities for simple ligand recognition problems where receptor sampling and protonation state effects are eliminated due to the simplicity of supramolecular hosts. We find that progress is now stagnated likely due to force field limitations.
Steven K. Albanese*, Daniel L. Parton*, Mehtap Isik**, Lucelenie Rodríguez-Laureano**, Sonya M. Hanson, Julie M. Behr, Scott Gradia, Chris Jeans, Nicholas M. Levinson, Markus A. Seeliger, and John D. Chodera.
* co-first author; ** co-second author
Biochemistry 57:4675, 2018. [DOI] [PDF] [bioRxiv] [GitHub]
Interactive data browser: [github.io]
Plasmids available via AddGene
Human kinase catalytic domains---the therapeutic target of selective kinase inhibitors used in the treatment of cancer and other diseases---are notoriously difficult and expensive to express in insect or human cells. Here, we utilize the phosphatase co-expression technology developed by Markus Seeliger (now at Stony Brook) to develop a library of human kinase catalytic domains for facile and inexpensive expression in bacteria.
Mehtap Işık, Dorothy Levorse, Ariën S. Rustenburg, Ikenna E. Ndukwe, Heather Wang , Xiao Wang , Mikhail Reibarkh , Gary E. Martin , Alexey A. Makarov , David L. Mobley, Timothy Rhodes*, John D. Chodera*.
* co-corresponding authors
Journal of Computer-Aided Molecular Design special issue on SAMPL6 32:1117, 2018.
[DOI] [PDF] [bioRxiv] [Supplementary Tables and Figures] [Supplementary Data (includes Sirius T3 reports on all measurements)]
The SAMPL5 blind challenge exercises identified neglect of protonation state effects as a major accuracy-limiting factor in physical modeling of biomolecular interactions. In this study, we report the experimental measurements behind a SAMPL6 blind challenges in which we assess the ability of community codes to predict small molecule pKas for small molecule resembling fragments of selective kinase inhibitors.
Andrew M. Intlekofer, Alan H. Shih, Bo Wang, Abbas Nazir, Ariën S. Rustenburg, Steven K. Albanese, Minal Patel, Christopher Famulare, Fabian M. Correa, Naofumi Takemoto, Vidushi Durani, Hui Liu, Justin Taylor, Noushin Farnoud, Elli Papaemmanuil, Justin R. Cross, Martin S. Tallman, Maria E. Arcila, Mikhail Roshal, Gregory A. Petsko, Bin Wu, Sung Choe, Zenon D. Konteatis, Scott A. Biller, John D. Chodera, Craig B. Thompson, Ross L. Levine, and Eytan M. Stein
Nature 559:125, 2018. [DOI] [PDF]
We probe the biochemical and biophysical mechanisms underlying an unusual set of clinical resistance mutations that appeared in trans in the IDH2 dimer interface in cancer patients treated with IDH2 inhibitors.
Kevin Hauser, Christopher Negron, Steven K. Albanese, Soumya Ray, Thomas Steinbrecher, Robert Abel, John D. Chodera, and Lingle Wang.
Communications Biology 1:70, 2018 [DOI] [PDF] [input files and analysis scripts]
In our first collaborative paper with Schrödinger, we present the first comprehensive benchmark assessing the ability for alchemical free energy calculations to predict clinical mutational resistance or susceptibility to targeted kinase inhibitors using the well-studied kinase Abl, the target of therapy for chronic myelogenous leukemia (CML).
Trung Hai Nguyen, Arien S. Rustenburg, Stefan G. Krimmer, Hexi Zhang, John D. Clark, Paul A. Novick, Kim Branson, Vijay S. Pande, John D Chodera, David D. L. Minh.
PLoS One 13:e0203224, 2018. [DOI] [bioRxiv] [GitHub]
We show how Bayesian inference can produce greatly improved estimates of statistical uncertainty from isothermal titration calorimetry (ITC) experiments, allowing the joint distribution of thermodynamic parameter uncertainties to be inferred.
Molecular dynamics simulations necessarily use a finite timestep, which introduces error or bias in the sampled configuration space density that grows rapidly with increasing timestep. For the first time, we show how to compute a natural measure of this error---the KL divergence---in both phase and configuration space for a widely used family of Langevin integrators, and show that VRORV is generally superior for simulation of molecular systems.
David L. Mobley, Caitlin C. Bannan, Andrea Rizzi, Christopher I. Bayly, John D. Chodera, Victoria T Lim, Nathan M. Lim, Kyle A. Beauchamp, Michael R. Shirts, Michael K. Gilson, Peter K. Eastman.
Journal of Chemical Theory and Computation 14:6076, 2018 [DOI] [bioRxiv]
We describe the philosophy behind a modern approach to molecular mechanics forcefield parameterization, and present initial results for the first SMIRNOFF-encoded forcefield: SMIRNOFF99Frosst.