Flexible receptor docking for drug discovery

Introduction: Molecular docking has become a popular method for virtual screening. Docking small molecules to a rigid biological receptor is fast but could produce many false negatives and identify less diverse compounds. Flexible receptor docking has alleviated this problem. Areas covered: This article focuses on reviewing ensemble docking as an approximate but inexpensive method to incorporate receptor flexibility in molecular docking. It outlines key features and recent advances of this method and points out problem areas that need to be addressed to make it even more useful in drug discovery. Expert opinion: Among the different methods introduced for flexible receptor docking, ensemble docking represents one of the most popular approaches, especially for high-throughput virtual screening. One can generate structural ensembles by using experimental structures, by structural modeling and by various types of molecular simulations. In building a structural ensemble, a judicious choice of the structures to be included can improve performance. Furthermore, reducing the size of the structural ensemble can cut computational costs, and removing the structures that can bind few ligands well could enrich the number of true actives identified by ensemble docking. The ability of ensemble docking to identify more true positives at the top of a rank-ordered list also depends on the choice of the methods to score and rank compounds, an area that needs further research.

[1]  J. P. Grossman,et al.  Simulation and embedded software development for Anton, a parallel machine with heterogeneous multicore ASICs , 2008, CODES+ISSS '08.

[2]  Kristine Steen Jensen,et al.  Dynamics and mechanisms of coupled protein folding and binding reactions. , 2012, Current opinion in structural biology.

[3]  F. Rao,et al.  Replica exchange molecular dynamics simulations of amyloid peptide aggregation. , 2004, The Journal of chemical physics.

[4]  D. Chipman,et al.  Solution of the linearized Poisson-Boltzmann equation. , 2004, The Journal of chemical physics.

[5]  J. P. Grossman,et al.  A 32x32x32, spatially distributed 3D FFT in four microseconds on Anton , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[6]  Robert V. Swift,et al.  Molecular docking to flexible targets. , 2015, Methods in molecular biology.

[7]  Wilfred F van Gunsteren,et al.  Practical Aspects of Free-Energy Calculations: A Review. , 2014, Journal of chemical theory and computation.

[8]  R. Copeland,et al.  Residence time of receptor-ligand complexes and its effect on biological function. , 2008, Biochemistry.

[9]  J. P. Grossman,et al.  Incorporating flexibility in Anton, a specialized machine for molecular dynamics simulation , 2008, 2008 IEEE 14th International Symposium on High Performance Computer Architecture.

[10]  Martin Zacharias,et al.  Efficient inclusion of receptor flexibility in grid‐based protein–ligand docking* , 2011, J. Comput. Chem..

[11]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[12]  J. P. Grossman,et al.  Unifying on-chip and inter-node switching within the Anton 2 network , 2014, 2014 ACM/IEEE 41st International Symposium on Computer Architecture (ISCA).

[13]  D. Swinney The role of binding kinetics in therapeutically useful drug action. , 2009, Current opinion in drug discovery & development.

[14]  P. A. Bash,et al.  Free energy calculations by computer simulation. , 1987, Science.

[15]  García,et al.  Large-amplitude nonlinear motions in proteins. , 1992, Physical review letters.

[16]  K. Fichthorn,et al.  Accelerated molecular dynamics of infrequent events , 1999 .

[17]  E. Di Cera,et al.  Essential role of conformational selection in ligand binding. , 2014, Biophysical chemistry.

[18]  Edmond Chow,et al.  Extending the Generality of Molecular Dynamics Simulations on a Special-Purpose Machine , 2013, 2013 IEEE 27th International Symposium on Parallel and Distributed Processing.

[19]  Erratum: The role of dynamic conformational ensembles in biomolecular recognition (Nature Chemical Biology (2009) 5 (789-796) , 2009 .

[20]  J Andrew McCammon,et al.  Molecular docking of balanol to dynamics snapshots of protein kinase A , 2005, Proteins.

[21]  G. Clore,et al.  Interplay between conformational selection and induced fit in multidomain protein-ligand binding probed by paramagnetic relaxation enhancement. , 2014, Biophysical chemistry.

[22]  R. Nussinov,et al.  Induced Fit, Conformational Selection and Independent Dynamic Segments: an Extended View of Binding Events Opinion , 2022 .

[23]  W. L. Jorgensen Free energy calculations: a breakthrough for modeling organic chemistry in solution , 1989 .

[24]  J Andrew McCammon,et al.  Conformational selection in G-proteins: lessons from Ras and Rho. , 2010, Biophysical journal.

[25]  György M. Keserü,et al.  Ensemble Docking into Flexible Active Sites. Critical Evaluation of FlexE against JNK-3 and beta-Secretase , 2006, J. Chem. Inf. Model..

[26]  Chung F. Wong Molecular simulation of drug-binding kinetics , 2014 .

[27]  Huan‐Xiang Zhou,et al.  Both protein dynamics and ligand concentration can shift the binding mechanism between conformational selection and induced fit , 2014, Proceedings of the National Academy of Sciences.

[28]  A Mitsutake,et al.  Generalized-ensemble algorithms for molecular simulations of biopolymers. , 2000, Biopolymers.

[29]  Nathan A. Baker,et al.  Poisson-Boltzmann Methods for Biomolecular Electrostatics , 2004, Numerical Computer Methods, Part D.

[30]  Heinz Sklenar,et al.  Harmonic modes as variables to approximately account for receptor flexibility in ligand-receptor docking simulations: Application to DNA minor groove ligand complex , 1999, J. Comput. Chem..

[31]  Chung F Wong,et al.  A mining minima approach to exploring the docking pathways of p-nitrocatechol sulfate to YopH. , 2007, Biophysical journal.

[32]  Arieh Warshel,et al.  Modeling electrostatic effects in proteins. , 2006, Biochimica et biophysica acta.

[33]  R. Elber,et al.  Reaction path study of conformational transitions in flexible systems: Applications to peptides , 1990 .

[34]  Chung F. Wong,et al.  Conformational selection of protein kinase A revealed by flexible‐ligand flexible‐protein docking , 2009, J. Comput. Chem..

[35]  C. MacKintosh,et al.  Inhibitors of protein kinases and phosphatases. , 1994, Trends in biochemical sciences.

[36]  M R Lee,et al.  Use of MM‐PB/SA in estimating the free energies of proteins: Application to native, intermediates, and unfolded villin headpiece , 2000, Proteins.

[37]  Jung-Hsin Lin,et al.  The relaxed complex method: Accommodating receptor flexibility for drug design with an improved scoring scheme. , 2003, Biopolymers.

[38]  D. Swinney,et al.  Applications of Binding Kinetics to Drug Discovery , 2012, Pharmaceutical Medicine.

[39]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[40]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[41]  John L. Klepeis,et al.  High-throughput pairwise point interactions in Anton, a specialized machine for molecular dynamics simulation , 2008, 2008 IEEE 14th International Symposium on High Performance Computer Architecture.

[42]  In-Hee Park,et al.  Dynamic ligand-induced-fit simulation via enhanced conformational samplings and ensemble dockings: a survivin example. , 2010, The journal of physical chemistry. B.

[43]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[44]  R. Abagyan,et al.  Systematic Exploitation of Multiple Receptor Conformations for Virtual Ligand Screening , 2011, PloS one.

[45]  Edmond Chow,et al.  Exploiting 162-Nanosecond End-to-End Communication Latency on Anton , 2010, 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis.

[46]  Chung F Wong,et al.  Flexible protein–flexible ligand docking with disrupted velocity simulated annealing , 2008, Proteins.

[47]  Y. Mu,et al.  Reduction of False Positives in Structure-Based Virtual Screening When Receptor Plasticity Is Considered , 2015, Molecules.

[48]  Ruth Nussinov,et al.  Protein dynamics and conformational selection in bidirectional signal transduction , 2012, BMC Biology.

[49]  K. Sanbonmatsu,et al.  Structure of Met‐enkephalin in explicit aqueous solution using replica exchange molecular dynamics , 2002, Proteins.

[50]  Dan Li,et al.  Importance of protein flexibility in ranking inhibitor affinities: modeling the binding mechanisms of piperidine carboxamides as Type I1/2 ALK inhibitors. , 2015, Physical chemistry chemical physics : PCCP.

[51]  A. Cavalli,et al.  Investigating drug-target association and dissociation mechanisms using metadynamics-based algorithms. , 2015, Accounts of chemical research.

[52]  L. R. Scott,et al.  Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program , 1995 .

[53]  John L. Klepeis,et al.  Anton, a special-purpose machine for molecular dynamics simulation , 2007, ISCA '07.

[54]  Chung F. Wong,et al.  A case study of scoring and rescoring in peptide docking. , 2012, Methods in molecular biology.

[55]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998, J. Comput. Chem..

[56]  Sally R. Ellingson,et al.  Multi-conformer ensemble docking to difficult protein targets. , 2015, The journal of physical chemistry. B.

[57]  Alessandro Pandini,et al.  Predicting the accuracy of protein–ligand docking on homology models , 2011, J. Comput. Chem..

[58]  J. P. Grossman,et al.  Hardware support for fine-grained event-driven computation in Anton 2 , 2013, ASPLOS '13.

[59]  Daniel-Adriano Silva,et al.  A Role for Both Conformational Selection and Induced Fit in Ligand Binding by the LAO Protein , 2011, PLoS Comput. Biol..

[60]  B. Shoichet,et al.  Flexible ligand docking using conformational ensembles , 1998, Protein science : a publication of the Protein Society.

[61]  M. Zacharias,et al.  Protein-ligand docking accounting for receptor side chain and global flexibility in normal modes: evaluation on kinase inhibitor cross docking. , 2008, Journal of medicinal chemistry.

[62]  Jung-Hsin Lin,et al.  Bridging implicit and explicit solvent approaches for membrane electrostatics. , 2002, Biophysical journal.

[63]  M Mezei,et al.  Free Energy Simulations a , 1986, Annals of the New York Academy of Sciences.

[64]  R. Nussinov,et al.  The role of dynamic conformational ensembles in biomolecular recognition. , 2009, Nature chemical biology.

[65]  J. Irwin,et al.  Benchmarking sets for molecular docking. , 2006, Journal of medicinal chemistry.

[66]  J. Mccammon,et al.  Computational drug design accommodating receptor flexibility: the relaxed complex scheme. , 2002, Journal of the American Chemical Society.

[67]  John L. Klepeis,et al.  Millisecond-scale molecular dynamics simulations on Anton , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[68]  Donald Hamelberg,et al.  Towards fast, rigorous and efficient conformational sampling of biomolecules: Advances in accelerated molecular dynamics. , 2015, Biochimica et biophysica acta.

[69]  Rommie E. Amaro,et al.  Emerging methods for ensemble-based virtual screening. , 2010, Current topics in medicinal chemistry.

[70]  Ron Elber,et al.  A method for determining reaction paths in large molecules: application to myoglobin , 1987 .

[71]  J. Andrew McCammon,et al.  Dynamics and design of enzymes and inhibitors , 1986 .

[72]  J. Mccammon,et al.  Induced Fit or Conformational Selection? The Role of the Semi-closed State in the Maltose Binding Protein , 2011, Biochemistry.

[73]  M. Karplus,et al.  Evaluation of the configurational entropy for proteins: application to molecular dynamics simulations of an α-helix , 1984 .

[74]  J. P. Grossman,et al.  The role of Cascade, a cycle-based simulation infrastructure, in designing the Anton special-purpose supercomputers , 2013, 2013 50th ACM/EDAC/IEEE Design Automation Conference (DAC).

[75]  P. Kollman,et al.  Computational Alanine Scanning To Probe Protein−Protein Interactions: A Novel Approach To Evaluate Binding Free Energies , 1999 .

[76]  J. Mongan,et al.  Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules. , 2004, The Journal of chemical physics.

[77]  Yuji Sugita,et al.  Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape , 2000, cond-mat/0009119.

[78]  Ruben Abagyan,et al.  Recipes for the Selection of Experimental Protein Conformations for Virtual Screening , 2010, J. Chem. Inf. Model..

[79]  A. Voter Hyperdynamics: Accelerated Molecular Dynamics of Infrequent Events , 1997 .

[80]  A. Voter,et al.  Accelerating the dynamics of infrequent events: Combining hyperdynamics and parallel replica dynamics to treat epitaxial layer growth , 1998 .

[81]  R. Elber,et al.  Reaction Path Study of Ligand Diffusion in Proteins: Application of the Self Penalty Walk (SPW) Method to Calculate Reaction Coordinates for the Motion of CO through Leghemoglobin , 1991 .

[82]  Edmond Chow,et al.  Overcoming Communication Latency Barriers in Massively Parallel Scientific Computation , 2011, IEEE Micro.

[83]  Chung F Wong,et al.  Docking flexible peptide to flexible protein by molecular dynamics using two implicit-solvent models: an evaluation in protein kinase and phosphatase systems. , 2009, The journal of physical chemistry. B.

[84]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[85]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[86]  Chung F. Wong,et al.  Beyond thermodynamics: drug binding kinetics could influence epidermal growth factor signaling. , 2009, Journal of medicinal chemistry.

[87]  Peter A. Kollman,et al.  A Ligand That Is Predicted to Bind Better to Avidin than Biotin: Insights from Computational Fluorine Scanning , 2000 .

[88]  P. Kollman,et al.  Binding of a diverse set of ligands to avidin and streptavidin: an accurate quantitative prediction of their relative affinities by a combination of molecular mechanics and continuum solvent models. , 2000, Journal of medicinal chemistry.

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

[90]  Claudio N. Cavasotto,et al.  Protein flexibility in ligand docking and virtual screening to protein kinases. , 2004, Journal of molecular biology.

[91]  Ron Elber,et al.  Reaction path study of helix formation in tetrapeptides: Effect of side chains , 1991 .

[92]  L. Chong,et al.  Reaching biological timescales with all-atom molecular dynamics simulations. , 2010, Current opinion in pharmacology.

[93]  I. Kuntz,et al.  Molecular docking to ensembles of protein structures. , 1997, Journal of molecular biology.

[94]  Irwin Oppenheim,et al.  Statistical Mechanical Theory of Transport Processes. VII. The Coefficient of Thermal Conductivity of Monatomic Liquids , 1954 .

[95]  Heinz Sklenar,et al.  Harmonic modes as variables to approximately account for receptor flexibility in ligand–receptor docking simulations: Application to DNA minor groove ligand complex , 1999 .

[96]  E. Di Cera,et al.  Conformational selection is a dominant mechanism of ligand binding. , 2013, Biochemistry.

[97]  J. Feeney,et al.  The use of saturation transfer NMR experiments to monitor the conformational selection accompanying ligand—protein interactions , 1980, FEBS letters.

[98]  Woody Sherman,et al.  Exploring protein flexibility: incorporating structural ensembles from crystal structures and simulation into virtual screening protocols. , 2012, The journal of physical chemistry. B.

[99]  Chung F. Wong,et al.  Cytochrome c: A Molecular Proving Ground for Computer Simulations , 1993 .

[100]  X. Barril,et al.  Unveiling the full potential of flexible receptor docking using multiple crystallographic structures. , 2005, Journal of medicinal chemistry.

[101]  B. Shoichet,et al.  Soft docking and multiple receptor conformations in virtual screening. , 2004, Journal of medicinal chemistry.

[102]  R. Abagyan,et al.  Representing receptor flexibility in ligand docking through relevant normal modes. , 2005, Journal of the American Chemical Society.

[103]  J. Changeux,et al.  Conformational selection or induced fit? 50 years of debate resolved , 2011, F1000 biology reports.

[104]  Michael D. Daily,et al.  Many local motions cooperate to produce the adenylate kinase conformational transition. , 2010, Journal of molecular biology.

[105]  Chung F. Wong Conformational transition paths harbor structures useful for aiding drug discovery and understanding enzymatic mechanisms in protein kinases , 2016, Protein science : a publication of the Protein Society.

[106]  Rafael C. Bernardi,et al.  Enhanced sampling techniques in molecular dynamics simulations of biological systems. , 2015, Biochimica et biophysica acta.

[107]  Eric T. Kim,et al.  How does a drug molecule find its target binding site? , 2011, Journal of the American Chemical Society.

[108]  Shantenu Jha,et al.  Computing Clinically Relevant Binding Free Energies of HIV-1 Protease Inhibitors , 2014, Journal of chemical theory and computation.

[109]  M. Zacharias,et al.  Accounting for global protein deformability during protein-protein and protein-ligand docking. , 2005, Biochimica et biophysica acta.

[110]  S. M. Sullivan,et al.  Enzymes with lid-gated active sites must operate by an induced fit mechanism instead of conformational selection , 2008, Proceedings of the National Academy of Sciences.

[111]  J. Andrew McCammon,et al.  Computer Simulation and the Design of New Biological Molecules , 1986 .

[112]  R. Elber,et al.  Molecular Dynamics at Extended Timescales , 2014 .

[113]  A. Voter A method for accelerating the molecular dynamics simulation of infrequent events , 1997 .

[114]  Brian K Shoichet,et al.  Testing a flexible-receptor docking algorithm in a model binding site. , 2004, Journal of molecular biology.

[115]  D. Swinney,et al.  A study of the molecular mechanism of binding kinetics and long residence times of human CCR5 receptor small molecule allosteric ligands , 2014, British journal of pharmacology.

[116]  D. Koshland Application of a Theory of Enzyme Specificity to Protein Synthesis. , 1958, Proceedings of the National Academy of Sciences of the United States of America.

[117]  Johan Åqvist,et al.  Ligand binding affinity prediction by linear interaction energy methods , 1998, J. Comput. Aided Mol. Des..

[118]  C. Brooks,et al.  Large-scale allosteric conformational transitions of adenylate kinase appear to involve a population-shift mechanism , 2007, Proceedings of the National Academy of Sciences.

[119]  Simulation reveals two major docking pathways between the hexapeptide GDYMNM and the catalytic domain of the insulin receptor protein kinase , 2012, Proteins.

[120]  R. Copeland Conformational adaptation in drug-target interactions and residence time. , 2011, Future medicinal chemistry.

[121]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[122]  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..

[123]  Haruki Nakamura,et al.  Virtual‐system‐coupled adaptive umbrella sampling to compute free‐energy landscape for flexible molecular docking , 2015, J. Comput. Chem..

[124]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[125]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[126]  A. Warshel,et al.  Electrostatic energy and macromolecular function. , 1991, Annual review of biophysics and biophysical chemistry.

[127]  Martin Zacharias,et al.  ReFlexIn: A Flexible Receptor Protein-Ligand Docking Scheme Evaluated on HIV-1 Protease , 2012, PloS one.

[128]  Christopher I. Bayly,et al.  Evaluating Virtual Screening Methods: Good and Bad Metrics for the "Early Recognition" Problem , 2007, J. Chem. Inf. Model..

[129]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[130]  Hugo Gutiérrez-de-Terán,et al.  Linear interaction energy: method and applications in drug design. , 2012, Methods in molecular biology.

[131]  Adrian E. Roitberg,et al.  MOIL: A program for simulations of macromolecules , 1995 .

[132]  D. Goodsell,et al.  Automated docking to multiple target structures: Incorporation of protein mobility and structural water heterogeneity in AutoDock , 2002, Proteins.

[133]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[134]  R. Abagyan,et al.  Flexible ligand docking to multiple receptor conformations: a practical alternative. , 2008, Current opinion in structural biology.

[135]  P. Kollman,et al.  Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding , 2000 .