We propose the development of a novel theoretical framework for the creation of next-generation quantum mechanical force fields (QMFFs) designed to meet the challenges of biomolecular and drug discovery applications. The proposed QMFFs are orders of magnitude faster than comparable linear-scaling electronic structure methods because the explicit interfragment orbital coupling is replaced by an empirical term that depends, at most, only on the electron density. The empirical parameters are tuned to reproduce experimental condensed phase properties such as ligand binding affinities, pKa shifts, solvation free energies, and lattice energies. In these ways, the proposed QMFFs will offer improved accuracy and predictive capability for multiscale modeling applications of complex biological and pharmaceutical problems. Specifically, we propose to 1. Create a general theoretical framework for development of new QMFFs for molecular simulations. This framework will allow arbitrary quantum mechanical Hamiltonian models to be integrated to form new QMFFs with accurate, balanced intermolecular interactions. Focus will be placed on the parameterization of intermolecular interactions using highly efficient low-level approximate density-functional tight-binding models, and high-level ab initio density-functional and coupled cluster methods that will be accessible through a specific interface with the PSI4 quantum chemistry package. The methods will take advantage of new linear-scaling methods for electrostatic interactions and generalized solvation effects, and will be implemented into a flexible set of software libraries that have interfaces with molecular simulation programs. 2. Apply QMFFs to thermodynamic properties of crystalline pharmaceuticals. The bioavailability of drugs are affected by crystal polymorph and co-crystal properties, such as solubility, and their thermodynamic characterization is thus an essential part of the drug discovery and chemical process pipeline. Density-functional QMFFs will be developed to study crystalline polymorphs by parameterizing them to reproduce known lattice energies. We will exploit the decomposition of the QMFF energy to compute the relative free energy between polymorphs, and absolute sublimation and solvation free energies of crystals using novel simulation methods that scale the intermolecular interactions. Our QMFF's ability to describe bonding events will be used to explore pKa shifts of drug compounds in co-crystals and make predictions about protonation states that are of pharmaceutical importance. 3. Apply QMFFs to aid in the interpretation of spectroscopic data. We will apply density-functional QMFFs to aid in the interpretation of IR, Raman, and NMR molecular spectroscopy of proteins and nucleic acids. First, small biologically relevant molecules will be used to benchmark the ability of density-functional QMFFs to predict IR, Raman and NMR spectra from molecular simulations. Second, we will apply the QMFFs to probe electric field heterogeneity on ketosteroid isomerase (KSI) with IR and 13C NMR chemical shifts. Third, we will examine the NMR chemical shift distributions in DNA and RNA systems derived from representative ensembles generated from MD simulation, and compare with recent and ongoing experiments.

Public Health Relevance

We propose a novel strategy to develop a class of integrated quantum mechanical force fields that can be used efficiently in molecular simulations of complex biological problems that involve reactivity and catalysis. The theoretical framework developed through this research will overcome a critical barrier to progress in the field by providing a foundation from which accurate, robust quantum mechanical methods can be applied in large scale biomolecular simulations. These methods will be utilized to interpret and predict various spectroscopic data, including IR, and NMR, of important biological systems, as well as to explore different polymorphs of crystalline pharmaceuticals through their thermodynamic properties.

Agency
National Institute of Health (NIH)
Institute
National Institute of General Medical Sciences (NIGMS)
Type
Research Project (R01)
Project #
5R01GM107485-03
Application #
9251846
Study Section
Special Emphasis Panel (ZRG1-BCMB-W (02)M)
Program Officer
Wehrle, Janna P
Project Start
2015-08-01
Project End
2019-03-31
Budget Start
2017-04-01
Budget End
2018-03-31
Support Year
3
Fiscal Year
2017
Total Cost
$268,649
Indirect Cost
$90,899
Name
Rutgers University
Department
Chemistry
Type
Schools of Arts and Sciences
DUNS #
001912864
City
Piscataway
State
NJ
Country
United States
Zip Code
08854
König, Gerhard; Pickard, Frank C; Huang, Jing et al. (2018) A Comparison of QM/MM Simulations with and without the Drude Oscillator Model Based on Hydration Free Energies of Simple Solutes. Molecules 23:
König, Gerhard; Brooks, Bernard R; Thiel, Walter et al. (2018) On the convergence of multi-scale free energy simulations. Mol Simul 44:1062-1081
Lee, Tai-Sung; Cerutti, David S; Mermelstein, Dan et al. (2018) GPU-Accelerated Molecular Dynamics and Free Energy Methods in Amber18: Performance Enhancements and New Features. J Chem Inf Model 58:2043-2050
Giese, Timothy J; York, Darrin M (2018) A GPU-Accelerated Parameter Interpolation Thermodynamic Integration Free Energy Method. J Chem Theory Comput 14:1564-1582
Lee, Tai-Sung; Hu, Yuan; Sherborne, Brad et al. (2017) Toward Fast and Accurate Binding Affinity Prediction with pmemdGTI: An Efficient Implementation of GPU-Accelerated Thermodynamic Integration. J Chem Theory Comput 13:3077-3084
Huang, Ming; Dissanayake, Thakshila; Kuechler, Erich et al. (2017) A Multidimensional B-Spline Correction for Accurate Modeling Sugar Puckering in QM/MM Simulations. J Chem Theory Comput 13:3975-3984
Giese, Timothy J; York, Darrin M (2017) Quantum mechanical force fields for condensed phase molecular simulations. J Phys Condens Matter 29:383002
Kuechler, Erich R; Giese, Timothy J; York, Darrin M (2016) VR-SCOSMO: A smooth conductor-like screening model with charge-dependent radii for modeling chemical reactions. J Chem Phys 144:164115
Giese, Timothy J; York, Darrin M (2016) Ambient-Potential Composite Ewald Method for ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulation. J Chem Theory Comput 12:2611-32
Hu, Yuan; Sherborne, Brad; Lee, Tai-Sung et al. (2016) The importance of protonation and tautomerization in relative binding affinity prediction: a comparison of AMBER TI and Schrödinger FEP. J Comput Aided Mol Des 30:533-9

Showing the most recent 10 out of 11 publications