BP-Dock: A Flexible Docking Scheme for Exploring Protein-Ligand Interactions Based on Unbound Structures

Molecular docking serves as an important tool in modeling protein-ligand interactions. However, it is still challenging to incorporate overall receptor flexibility, especially backbone flexibility, in docking due to the large conformational space that needs to be sampled. To overcome this problem, we developed a novel flexible docking approach, BP-Dock (Backbone Perturbation-Dock) that can integrate both backbone and side chain conformational changes induced by ligand binding through a multi-scale approach. In the BP-Dock method, we mimic the nature of binding-induced events as a first-order approximation by perturbing the residues along the protein chain with a small Brownian kick one at a time. The response fluctuation profile of the chain upon these perturbations is computed using the perturbation response scanning method. These response fluctuation profiles are then used to generate binding-induced multiple receptor conformations for ensemble docking. To evaluate the performance of BP-Dock, we applied our approach on a large and diverse data set using unbound structures as receptors. We also compared the BP-Dock results with bound and unbound docking, where overall receptor flexibility was not taken into account. Our results highlight the importance of modeling backbone flexibility in docking for recapitulating the experimental binding affinities, especially when an unbound structure is used. With BP-Dock, we can generate a wide range of binding site conformations realized in nature even in the absence of a ligand that can help us to improve the accuracy of unbound docking. We expect that our fast and efficient flexible docking approach may further aid in our understanding of protein-ligand interactions as well as virtual screening of novel targets for rational drug design.

[1]  Martin Zacharias,et al.  Rapid protein–ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: Binding of FK506 to FKBP , 2004, Proteins.

[2]  R. Friesner,et al.  Novel procedure for modeling ligand/receptor induced fit effects. , 2006, Journal of medicinal chemistry.

[3]  Ruth Nussinov,et al.  Geometry‐based flexible and symmetric protein docking , 2005, Proteins.

[4]  M. Noble,et al.  The adaptability of the active site of trypanosomal triosephosphate isomerase as observed in the crystal structures of three different complexes , 1991, Proteins.

[5]  Jens Meiler,et al.  A Correspondence Between Solution-State Dynamics of an Individual Protein and the Sequence and Conformational Diversity of its Family , 2009, PLoS Comput. Biol..

[6]  Martin Zacharias,et al.  Accounting for conformational changes during protein-protein docking. , 2010, Current opinion in structural biology.

[7]  D. Case,et al.  Molecular Dynamics Simulations of Nucleic Acids with a Generalized Born Solvation Model , 2000 .

[8]  J. Bajorath,et al.  Docking and scoring in virtual screening for drug discovery: methods and applications , 2004, Nature Reviews Drug Discovery.

[9]  Ian W. Davis,et al.  The backrub motion: how protein backbone shrugs when a sidechain dances. , 2006, Structure.

[10]  Katrina W Lexa,et al.  Protein flexibility in docking and surface mapping , 2012, Quarterly Reviews of Biophysics.

[11]  Ian W. Davis,et al.  RosettaLigand docking with full ligand and receptor flexibility. , 2009, Journal of molecular biology.

[12]  Tirion,et al.  Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis. , 1996, Physical review letters.

[13]  Tanja Kortemme,et al.  Structure-based prediction of the peptide sequence space recognized by natural and synthetic PDZ domains. , 2010, Journal of molecular biology.

[14]  L. Cantley,et al.  Recognition of Unique Carboxyl-Terminal Motifs by Distinct PDZ Domains , 1997, Science.

[15]  Michel F. Sanner,et al.  Protein–ligand docking with multiple flexible side chains , 2008, J. Comput. Aided Mol. Des..

[16]  Thomas Lengauer,et al.  On the Applicability of Elastic Network Normal Modes in Small-Molecule Docking , 2012, J. Chem. Inf. Model..

[17]  Christopher R. Corbeil,et al.  Docking Ligands into Flexible and Solvated Macromolecules. 2. Development and Application of Fitted 1.5 to the Virtual Screening of Potential HCV Polymerase Inhibitors , 2008, J. Chem. Inf. Model..

[18]  Guang Song,et al.  Protein elastic network models and the ranges of cooperativity , 2009, Proceedings of the National Academy of Sciences.

[19]  Ali Rana Atilgan,et al.  Perturbation-Response Scanning Reveals Ligand Entry-Exit Mechanisms of Ferric Binding Protein , 2009, PLoS Comput. Biol..

[20]  Gerhard Klebe,et al.  Probing flexibility and “induced‐fit” phenomena in aldose reductase by comparative crystal structure analysis and molecular dynamics simulations , 2004, Proteins.

[21]  Rachelle J Bienstock,et al.  Computational drug design targeting protein-protein interactions. , 2012, Current pharmaceutical design.

[22]  Christopher R. Corbeil,et al.  Docking Ligands into Flexible and Solvated Macromolecules. 3. Impact of Input Ligand Conformation, Protein Flexibility, and Water Molecules on the Accuracy of Docking Programs , 2009, J. Chem. Inf. Model..

[23]  Oliver Korb,et al.  Potential and Limitations of Ensemble Docking , 2012, J. Chem. Inf. Model..

[24]  Martin Zacharias,et al.  Combining Elastic Network Analysis and Molecular Dynamics Simulations by Hamiltonian Replica Exchange. , 2008, Journal of chemical theory and computation.

[25]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[26]  Z. Nevin Gerek,et al.  Change in Allosteric Network Affects Binding Affinities of PDZ Domains: Analysis through Perturbation Response Scanning , 2011, PLoS Comput. Biol..

[27]  I. Bahar,et al.  Global dynamics of proteins: bridging between structure and function. , 2010, Annual review of biophysics.

[28]  R. Jernigan,et al.  Anisotropy of fluctuation dynamics of proteins with an elastic network model. , 2001, Biophysical journal.

[29]  Ruth Nussinov,et al.  A Method for Biomolecular Structural Recognition and Docking Allowing Conformational Flexibility , 1998, J. Comput. Biol..

[30]  R Nussinov,et al.  Flexible docking allowing induced fit in proteins: Insights from an open to closed conformational isomers , 1998, Proteins.

[31]  Jens Meiler,et al.  ROSETTALIGAND: Protein–small molecule docking with full side‐chain flexibility , 2006, Proteins.

[32]  W. Lipscomb,et al.  Comparison of the structures of three carboxypeptidase A-phosphonate complexes determined by X-ray crystallography. , 1994, Biochemistry.

[33]  A. Kidera,et al.  Protein structural change upon ligand binding: linear response theory. , 2005, Physical review letters.

[34]  A. Atilgan,et al.  Manipulation of conformational change in proteins by single-residue perturbations. , 2010, Biophysical journal.

[35]  Markus Wagener,et al.  A flexible approach to induced fit docking. , 2007, Journal of medicinal chemistry.

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

[37]  S. Ozkan,et al.  A flexible docking scheme to explore the binding selectivity of PDZ domains , 2010, Protein science : a publication of the Protein Society.

[38]  Tanja Kortemme,et al.  RosettaBackrub—a web server for flexible backbone protein structure modeling and design , 2010, Nucleic Acids Res..

[39]  K. Hinsen.  Analysis of domain motions by approximate normal mode calculations , 1998, Proteins.

[40]  Jacek Otlewski,et al.  PDZ tandem of human syntenin: crystal structure and functional properties. , 2003, Structure.

[41]  B. Shoichet,et al.  Information decay in molecular docking screens against holo, apo, and modeled conformations of enzymes. , 2003, Journal of medicinal chemistry.

[42]  O. Keskin,et al.  Identification of specificity and promiscuity of PDZ domain interactions through their dynamic behavior , 2009, Proteins.

[43]  J A McCammon,et al.  Accommodating protein flexibility in computational drug design. , 2000, Molecular pharmacology.

[44]  David R Cooper,et al.  Molecular roots of degenerate specificity in syntenin's PDZ2 domain: reassessment of the PDZ recognition paradigm. , 2003, Structure.

[45]  Ruth Nussinov,et al.  An integrated suite of fast docking algorithms , 2010, Proteins.

[46]  Claudio N. Cavasotto,et al.  Representing receptor flexibility in ligand docking through relevant normal modes. , 2005, Journal of the American Chemical Society.

[47]  G. Klebe,et al.  Expect the unexpected or caveat for drug designers: multiple structure determinations using aldose reductase crystals treated under varying soaking and co-crystallisation conditions. , 2006, Journal of molecular biology.

[48]  W. Delano.  The PyMOL Molecular Graphics System , 2002 .

[49]  Christopher R. Corbeil,et al.  Docking Ligands into Flexible and Solvated Macromolecules, 1. Development and Validation of FITTED 1.0 , 2007, J. Chem. Inf. Model..

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

[51]  S Mangani,et al.  High-resolution structure of the complex between carboxypeptidase A and L-phenyl lactate. , 1993, Acta crystallographica. Section D, Biological crystallography.

[52]  X. Zou,et al.  Ensemble docking of multiple protein structures: Considering protein structural variations in molecular docking , 2006, Proteins.

[53]  Paul N. Mortenson,et al.  Diverse, high-quality test set for the validation of protein-ligand docking performance. , 2007, Journal of medicinal chemistry.

[54]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[55]  Thomas Lengauer,et al.  FlexE: efficient molecular docking considering protein structure variations. , 2001, Journal of molecular biology.

[56]  H. Wolfson,et al.  Principles of flexible protein–protein docking , 2008, Proteins.

[57]  J. MacQueen.  Some methods for classification and analysis of multivariate observations , 1967 .

[58]  Ruth Nussinov,et al.  Automatic prediction of protein interactions with large scale motion , 2007, Proteins.

[59]  Luhua Lai,et al.  Further development and validation of empirical scoring functions for structure-based binding affinity prediction , 2002, J. Comput. Aided Mol. Des..

[60]  Renxiao Wang,et al.  Comparative evaluation of 11 scoring functions for molecular docking. , 2003, Journal of medicinal chemistry.

[61]  Michel F Sanner,et al.  FLIPDock: Docking flexible ligands into flexible receptors , 2007, Proteins.

[62]  M. Philippopoulos,et al.  Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI , 1999, Proteins.

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

[64]  David S. Goodsell,et al.  AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility , 2009, J. Comput. Chem..

[65]  Chaok Seok,et al.  GalaxyDock: Protein-Ligand Docking with Flexible Protein Side-chains , 2012, J. Chem. Inf. Model..

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

[67]  Feng Ding,et al.  Incorporating Backbone Flexibility in MedusaDock Improves Ligand-Binding Pose Prediction in the CSAR2011 Docking Benchmark , 2013, J. Chem. Inf. Model..

[68]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[69]  O. Keskin,et al.  The binding affinities of proteins interacting with the PDZ domain of PICK1 , 2012, Proteins.

[70]  Woody Sherman,et al.  Use of an Induced Fit Receptor Structure in Virtual Screening , 2006, Chemical biology & drug design.

[71]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[72]  Heather A Carlson,et al.  Exploring experimental sources of multiple protein conformations in structure-based drug design. , 2007, Journal of the American Chemical Society.

[73]  J M Blaney,et al.  A geometric approach to macromolecule-ligand interactions. , 1982, Journal of molecular biology.

[74]  C L Brooks,et al.  Ligand-protein database: linking protein-ligand complex structures to binding data. , 2001, Journal of medicinal chemistry.

[75]  R. Nussinov,et al.  Exploiting conformational ensembles in modeling protein-protein interactions on the proteome scale. , 2013, Journal of proteome research.