Modern drug design: the implication of using artificial neuronal networks and multiple molecular dynamic simulations

We report the implementation of molecular modeling approaches developed as a part of the 2016 Grand Challenge 2, the blinded competition of computer aided drug design technologies held by the D3R Drug Design Data Resource (https://drugdesigndata.org/). The challenge was focused on the ligands of the farnesoid X receptor (FXR), a highly flexible nuclear receptor of the cholesterol derivative chenodeoxycholic acid. FXR is considered an important therapeutic target for metabolic, inflammatory, bowel and obesity related diseases (Expert Opin Drug Metab Toxicol 4:523-532, 2015), but in the context of this competition it is also interesting due to the significant ligand-induced conformational changes displayed by the protein. To deal with these conformational changes we employed multiple simulations of molecular dynamics (MD). Our MD-based protocols were top-ranked in estimating the free energy of binding of the ligands and FXR protein. Our approach was ranked second in the prediction of the binding poses where we also combined MD with molecular docking and artificial neural networks. Our approach showed mediocre results for high-throughput scoring of interactions.

[1]  Bert L. de Groot,et al.  g_wham—A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates , 2010 .

[2]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[3]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[4]  J. Kästner.  Umbrella sampling , 2011 .

[5]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[6]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[7]  Guigang Zhang,et al.  Deep Learning , 2016, Int. J. Semantic Comput..

[8]  M. Parrinello,et al.  Funnel metadynamics as accurate binding free-energy method , 2013, Proceedings of the National Academy of Sciences.

[9]  Geoffrey E. Hinton,et al.  Rectified Linear Units Improve Restricted Boltzmann Machines , 2010, ICML.

[10]  Jens Carlsson,et al.  Combining docking, molecular dynamics and the linear interaction energy method to predict binding modes and affinities for non-nucleoside inhibitors to HIV-1 reverse transcriptase. , 2008, Journal of medicinal chemistry.

[11]  Andreas Bender,et al.  Recognizing Pitfalls in Virtual Screening: A Critical Review , 2012, J. Chem. Inf. Model..

[12]  Michael K. Gilson,et al.  Overcoming dissipation in the calculation of standard binding free energies by ligand extraction , 2013, J. Comput. Chem..

[13]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[14]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[15]  Ryan D. Morin,et al.  Theoretical Investigation of the D83V Mutation within the Myocyte-Specific Enhancer Factor-2 Beta and Its Role in Cancer , 2013 .

[16]  E. Patsouris,et al.  Farnesoid x receptor in human metabolism and disease: the interplay between gene polymorphisms, clinical phenotypes and disease susceptibility , 2015, Expert opinion on drug metabolism & toxicology.

[17]  Scott Kirkpatrick,et al.  Optimization by Simmulated Annealing , 1983, Sci..

[18]  B. Roux,et al.  Calculation of absolute protein-ligand binding free energy from computer simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[19]  B. Kuhn,et al.  Optimization of a novel class of benzimidazole-based farnesoid X receptor (FXR) agonists to improve physicochemical and ADME properties. , 2011, Bioorganic & medicinal chemistry letters.

[20]  M. Gilson,et al.  The statistical-thermodynamic basis for computation of binding affinities: a critical review. , 1997, Biophysical journal.

[21]  B. Roux,et al.  Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials. , 2006, Biophysical journal.

[22]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

[23]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[24]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[25]  Geoffrey E. Hinton,et al.  ImageNet classification with deep convolutional neural networks , 2012, Commun. ACM.

[26]  KumarShankar,et al.  The weighted histogram analysis method for free-energy calculations on biomolecules. I , 1992 .

[27]  Yoshua Bengio,et al.  Gradient-based learning applied to document recognition , 1998, Proc. IEEE.

[28]  Christophe Chipot,et al.  The Adaptive Biasing Force Method: Everything You Always Wanted To Know but Were Afraid To Ask , 2014, The journal of physical chemistry. B.

[29]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[30]  Andrew T. Fenley,et al.  Computational Calorimetry: High-Precision Calculation of Host–Guest Binding Thermodynamics , 2015, Journal of chemical theory and computation.

[31]  Dariusz Plewczynski,et al.  Can we trust docking results? Evaluation of seven commonly used programs on PDBbind database , 2011, J. Comput. Chem..

[32]  B. MacLennan,et al.  Combining Genetic Algorithms and Neural Networks : The Encoding Problem , 1994 .

[33]  A. Cavalli,et al.  The ligand binding mechanism to purine nucleoside phosphorylase elucidated via molecular dynamics and machine learning , 2015, Nature Communications.

[34]  Volodymyr G. Bdzhola,et al.  Kirchhoff atomic charges fitted to multipole moments: Implementation for a virtual screening system , 2008, J. Comput. Chem..

[35]  Michael K Gilson,et al.  On the theory of noncovalent binding. , 2004, Biophysical journal.

[36]  Luca Mollica,et al.  Kinetics of protein-ligand unbinding via smoothed potential molecular dynamics simulations , 2015, Scientific Reports.

[37]  Izhar Wallach,et al.  AtomNet: A Deep Convolutional Neural Network for Bioactivity Prediction in Structure-based Drug Discovery , 2015, ArXiv.

[38]  Richard H. Henchman,et al.  Standard Free Energy of Binding from a One-Dimensional Potential of Mean Force. , 2009, Journal of chemical theory and computation.

[39]  Hans Richter,et al.  Identification of an N-oxide pyridine GW4064 analog as a potent FXR agonist. , 2009, Bioorganic & medicinal chemistry letters.

[40]  Niel M. Henriksen,et al.  Attach-Pull-Release Calculations of Ligand Binding and Conformational Changes on the First BRD4 Bromodomain , 2017, Journal of chemical theory and computation.

[41]  Richard P. Brent,et al.  An Algorithm with Guaranteed Convergence for Finding a Zero of a Function , 1971, Comput. J..

[42]  Christophe Chipot,et al.  Standard binding free energies from computer simulations: What is the best strategy? , 2013, Journal of chemical theory and computation.

[43]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[44]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[45]  Jens Carlsson,et al.  Improving the Accuracy of the Linear Interaction Energy Method for Solvation Free Energies. , 2007, Journal of chemical theory and computation.

[46]  Steven J. M. Jones,et al.  Ab initio parameterization of YFF1, a universal force field for drug-design applications , 2012, Journal of Molecular Modeling.

[47]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[48]  Melanie Keller,et al.  Essentials Of Computational Chemistry Theories And Models , 2016 .

[49]  Bradley M Dickson,et al.  A fast, open source implementation of adaptive biasing potentials uncovers a ligand design strategy for the chromatin regulator BRD4. , 2016, The Journal of chemical physics.

[50]  A. Cavalli,et al.  Molecular Dynamics Simulations and Kinetic Measurements to Estimate and Predict Protein-Ligand Residence Times. , 2016, Journal of medicinal chemistry.