Pharmaceutical drug discovery is time-consuming and expensive, with each new drug brought to market now costing roughly $1 billion on average. This cost is driven by the difficulty of drug discovery, and in part by the amount of trial and error involved in the process of finding initial """"""""hits"""""""" which modulate the function of a biomolecule, and then refining these into """"""""leads"""""""" which have adequate affinity for the biomolecular target and other desirable properties. Computational methods ideally could guide this process, reducing the amount of trial and error involved by suggesting hits in advance of experiment and predicting chemical modifications which will improve these into leads, enhancing affinity while maintaining drug-like properties. But current computational methods are not adequate to change the discovery process in this way. Recent innovations in alchemical free energy calculations based on molecular simulations show considerable promise at reaching the level of accuracy needed to help drug discovery, but these simulations require considerable expertise to set up and conduct, and a great deal of computer power. This proposal focuses on lowering these barriers, providing a new approach to automatically plan and set up these calculations, and improved computational efficiency. Alchemical free energy calculations are one of the most physically realistic computational approaches available, and one of the most promising in terms of accuracy. This project's aims are to (1) develop a new tool to automate setup of relative binding free energy calculations for drug lead optimization;(2) efficiently calculate ligand binding mode occupancies, dramatically reducing the computational expense of binding free energy predictions;and (3) use these techniques to guide experimental drug discovery of histone methyltransferase inhibitors, which show considerable promise as potential anti-cancer drugs. While considerable effort has gone into alchemical free energy calculations, one innovative aspect of this work is the focus on predicting binding mode as well as binding affinity. This is handled by using fast docking methods, in combination with exploratory simulations, to identify a variety of stable ligand binding modes, then including all of these in binding free energy calculations, so that bound structures of individual inhibitors need not be known in advance. This work will speed up promising tools for affinity calculation, and improve automation so that they can more easily be applied to problems in drug discovery. The long-term goal of these techniques is to change the early stage drug discovery process by providing robust computational affinity predictions, and this work provides an important step in that direction.

Public Health Relevance

Pharmaceutical drug discovery produces dramatic public health benefits through new and improved treatments for common diseases and disorders, but it is an expensive, time-consuming process involving much trial and error, and failure is common. This work focuses on dramatically improving computer tools to predict interactions between small molecules and molecular machines. The proposed approaches will help guide the development of new drugs, resulting potential public health rewards.

National Institute of Health (NIH)
National Institute of General Medical Sciences (NIGMS)
Research Project (R01)
Project #
Application #
Study Section
Macromolecular Structure and Function D Study Section (MSFD)
Program Officer
Preusch, Peter
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
University of California Irvine
Schools of Pharmacy
United States
Zip Code
Riquelme, Maximiliano; Lara, Alejandro; Mobley, David L et al. (2018) Hydration Free Energies in the FreeSolv Database Calculated with Polarized Iterative Hirshfeld Charges. J Chem Inf Model 58:1779-1797
I??k, Mehtap; Levorse, Dorothy; Rustenburg, Ariën S et al. (2018) pKa measurements for the SAMPL6 prediction challenge for a set of kinase inhibitor-like fragments. J Comput Aided Mol Des 32:1117-1138
Gill, Samuel C; Lim, Nathan M; Grinaway, Patrick B et al. (2018) Binding Modes of Ligands Using Enhanced Sampling (BLUES): Rapid Decorrelation of Ligand Binding Modes via Nonequilibrium Candidate Monte Carlo. J Phys Chem B 122:5579-5598
Kyrychenko, Alexander; Lim, Nathan M; Vasquez-Montes, Victor et al. (2018) Refining Protein Penetration into the Lipid Bilayer Using Fluorescence Quenching and Molecular Dynamics Simulations: The Case of Diphtheria Toxin Translocation Domain. J Membr Biol 251:379-391
Mobley, David L; Bannan, Caitlin C; Rizzi, Andrea et al. (2018) Escaping Atom Types in Force Fields Using Direct Chemical Perception. J Chem Theory Comput 14:6076-6092
Mobley, David L; Gilson, Michael K (2017) Predicting Binding Free Energies: Frontiers and Benchmarks. Annu Rev Biophys 46:531-558
Könst, Zef A; Szklarski, Anne R; Pellegrino, Simone et al. (2017) Synthesis facilitates an understanding of the structural basis for translation inhibition by the lissoclimides. Nat Chem 9:1140-1149
Yin, Jian; Henriksen, Niel M; Slochower, David R et al. (2017) Overview of the SAMPL5 host-guest challenge: Are we doing better? J Comput Aided Mol Des 31:1-19
Shirts, Michael R; Klein, Christoph; Swails, Jason M et al. (2017) Lessons learned from comparing molecular dynamics engines on the SAMPL5 dataset. J Comput Aided Mol Des 31:147-161
Bannan, Caitlin C; Burley, Kalistyn H; Chiu, Michael et al. (2016) Blind prediction of cyclohexane-water distribution coefficients from the SAMPL5 challenge. J Comput Aided Mol Des 30:927-944

Showing the most recent 10 out of 17 publications