Search for β2 Adrenergic Receptor Ligands by Virtual Screening via Grid Computing and Investigation of Binding Modes by Docking and Molecular Dynamics Simulations

We designed a program called MolGridCal that can be used to screen small molecule database in grid computing on basis of JPPF grid environment. Based on MolGridCal program, we proposed an integrated strategy for virtual screening and binding mode investigation by combining molecular docking, molecular dynamics (MD) simulations and free energy calculations. To test the effectiveness of MolGridCal, we screened potential ligands for β2 adrenergic receptor (β2AR) from a database containing 50,000 small molecules. MolGridCal can not only send tasks to the grid server automatically, but also can distribute tasks using the screensaver function. As for the results of virtual screening, the known agonist BI-167107 of β2AR is ranked among the top 2% of the screened candidates, indicating MolGridCal program can give reasonable results. To further study the binding mode and refine the results of MolGridCal, more accurate docking and scoring methods are used to estimate the binding affinity for the top three molecules (agonist BI-167107, neutral antagonist alprenolol and inverse agonist ICI 118,551). The results indicate agonist BI-167107 has the best binding affinity. MD simulation and free energy calculation are employed to investigate the dynamic interaction mechanism between the ligands and β2AR. The results show that the agonist BI-167107 also has the lowest binding free energy. This study can provide a new way to perform virtual screening effectively through integrating molecular docking based on grid computing, MD simulations and free energy calculations. The source codes of MolGridCal are freely available at http://molgridcal.codeplex.com.

[1]  I. Hall,et al.  Influence of agonist efficacy and receptor phosphorylation on antagonist affinity measurements: differences between second messenger and reporter gene responses. , 2003, Molecular pharmacology.

[2]  L. Tan,et al.  Polymorphisms at the β2-adrenergic receptor gene influence Alzheimer's disease susceptibility , 2008, Brain Research.

[3]  David S. Goodsell,et al.  Protein Flexibility in Virtual Screening: The BACE-1 Case Study , 2012, J. Chem. Inf. Model..

[4]  Kaihsu Tai,et al.  Grid computing and biomolecular simulation , 2005, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[5]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[6]  Claude Lecomte,et al.  Hydrogen bond strengths revealed by topological analyses of experimentally observed electron densities , 1998 .

[7]  S. Liggett,et al.  Molecular Mechanisms of β2-Adrenergic Receptor Function and Regulation , 2005 .

[8]  Vijay S. Pande,et al.  Screen Savers of the World Unite! , 2000, Science.

[9]  Thomas Sander,et al.  Integration of distributed computing into the drug discovery process , 2011, Expert opinion on drug discovery.

[10]  Neil Geddes,et al.  The Large Hadron Collider and Grid computing , 2012, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[11]  Markus H. J. Seifert,et al.  Virtual high-throughput screening of molecular databases. , 2007, Current opinion in drug discovery & development.

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

[13]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[14]  Malcolm W Johnson.  Molecular mechanisms of beta(2)-adrenergic receptor function, response, and regulation. , 2006, The Journal of allergy and clinical immunology.

[15]  Ryan G. Coleman,et al.  ZINC: A Free Tool to Discover Chemistry for Biology , 2012, J. Chem. Inf. Model..

[16]  Klaus Kayser,et al.  Grid computing in image analysis , 2011, Diagnostic pathology.

[17]  Steven W. Muchmore,et al.  High-Throughput Calculation of Protein-Ligand Binding Affinities: Modification and Adaptation of the MM-PBSA Protocol to Enterprise Grid Computing , 2006, J. Chem. Inf. Model..

[18]  Rommie E. Amaro,et al.  AutoGrow: A Novel Algorithm for Protein Inhibitor Design , 2009, Chemical biology & drug design.

[19]  J. Baker,et al.  The selectivity of β‐adrenoceptor agonists at human β1‐, β2‐ and β3‐adrenoceptors , 2010, British journal of pharmacology.

[20]  Qifeng Bai,et al.  Molecular modeling study on the dynamical structural features of human smoothened receptor and binding mechanism of antagonist LY2940680 by metadynamics simulation and free energy calculation. , 2014, Biochimica et biophysica acta.

[21]  Brian K. Shoichet,et al.  ZINC - A Free Database of Commercially Available Compounds for Virtual Screening , 2005, J. Chem. Inf. Model..

[22]  Qifeng Bai,et al.  Computational Study on the Different Ligands Induced Conformation Change of β2 Adrenergic Receptor-Gs Protein Complex , 2013, PloS one.

[23]  Davide Provasi,et al.  Ligand-Induced Modulation of the Free-Energy Landscape of G Protein-Coupled Receptors Explored by Adaptive Biasing Techniques , 2011, PLoS Comput. Biol..

[24]  Stefano Vanni,et al.  Predicting Novel Binding Modes of Agonists to β Adrenergic Receptors Using All-Atom Molecular Dynamics Simulations , 2011, PLoS Comput. Biol..

[25]  Charles L. Brooks,et al.  Detailed analysis of grid‐based molecular docking: A case study of CDOCKER—A CHARMm‐based MD docking algorithm , 2003, J. Comput. Chem..

[26]  Zhihui Wang,et al.  Accelerating cancer systems biology research through Semantic Web technology , 2013, Wiley interdisciplinary reviews. Systems biology and medicine.

[27]  N. Jacq,et al.  Grid-Enabled High-Throughput In Silico Screening Against Influenza A Neuraminidase , 2006, IEEE Transactions on NanoBioscience.

[28]  Adam Hunter,et al.  Yabi: An online research environment for grid, high performance and cloud computing , 2012, Source Code for Biology and Medicine.

[29]  Malcolm W Johnson.  Molecular mechanisms of β2-adrenergic receptor function, response, and regulation , 2006 .

[30]  Dimitrios Zikos,et al.  Distributed and Grid Computing Projects with Research Focus in Human Health , 2012, EFMI-STC.

[31]  T. Nero,et al.  LK 204-545, a highly selective beta1-adrenoceptor antagonist at human beta-adrenoceptors. , 1999, European journal of pharmacology.

[32]  Tingjun Hou,et al.  Studies on the Interactions between β2 Adrenergic Receptor and Gs Protein by Molecular Dynamics Simulations , 2012, J. Chem. Inf. Model..

[33]  W. Graham Richards,et al.  Virtual screening using grid computing: the screensaver project , 2002, Nature Reviews Drug Discovery.

[34]  E. Bleecker,et al.  Effect of β2-adrenergic receptor gene (ADRB2) 3′ untranslated region polymorphisms on inhaled corticosteroid/long-acting β2-adrenergic agonist response , 2012, Respiratory Research.

[35]  Jillian G. Baker,et al.  The selectivity of β‐adrenoceptor antagonists at the human β1, β2 and β3 adrenoceptors , 2005 .

[36]  Stephen R. Johnson,et al.  Grid computing in large pharmaceutical molecular modeling. , 2008, Drug discovery today.

[37]  Hermann-Georg Holzhütter,et al.  Molecular dioxygen enters the active site of 12/15-lipoxygenase via dynamic oxygen access channels , 2007, Proceedings of the National Academy of Sciences.

[38]  Dariusz Leniowski,et al.  Acorn: A grid computing system for constraint based modeling and visualization of the genome scale metabolic reaction networks via a web interface , 2011, BMC Bioinformatics.

[39]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[40]  H Binder,et al.  Grid and Cloud Computing Methods in Biomedical Research , 2013, Methods of Information in Medicine.

[41]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[42]  S. Liggett,et al.  Molecular mechanisms of beta2-adrenergic receptor function and regulation. , 2005, Proceedings of the American Thoracic Society.

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

[44]  K. Garcia,et al.  Adrenaline-activated structure of the β2-adrenoceptor stabilized by an engineered nanobody , 2013, Nature.

[45]  C. Janson,et al.  Pneumonia and pneumonia related mortality in patients with COPD treated with fixed combinations of inhaled corticosteroid and long acting β2 agonist: observational matched cohort study (PATHOS) , 2013, BMJ.

[46]  Albert C. Pan,et al.  Pathway and mechanism of drug binding to G-protein-coupled receptors , 2011, Proceedings of the National Academy of Sciences.

[47]  Antony J. Williams,et al.  A perspective of publicly accessible/open-access chemistry databases. , 2008, Drug discovery today.

[48]  Jonathan W. Essex,et al.  A review of protein-small molecule docking methods , 2002, J. Comput. Aided Mol. Des..

[49]  Leonardo Pardo,et al.  Molecular Basis of Ligand Dissociation in β-Adrenergic Receptors , 2011, PloS one.

[50]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[51]  M. Mihăşan,et al.  What in silico molecular docking can do for the ‘bench-working biologists’ , 2012, Journal of Biosciences.

[52]  J. Shuster,et al.  Methacholine Challenge as a Clinical Bioassay of Pulmonary Delivery of a Long‐Acting β2‐Adrenergic Agonist , 2011, Pharmacotherapy.

[53]  Rafael Mayo,et al.  Characterization of antigenetic serotypes from the dengue virus in Venezuela by means of Grid Computing , 2010, HealthGrid.

[54]  Eric Darve,et al.  Adaptive biasing force method for scalar and vector free energy calculations. , 2008, The Journal of chemical physics.

[55]  Christophe Chipot,et al.  Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables. , 2010, Journal of chemical theory and computation.

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

[57]  T. Patterson,et al.  Vulnerable caregivers of Alzheimer disease patients have a deficit in beta 2-adrenergic receptor sensitivity and density. , 2004, The American journal of geriatric psychiatry : official journal of the American Association for Geriatric Psychiatry.

[58]  Francesca Fanelli,et al.  Update 1 of: computational modeling approaches to structure-function analysis of G protein-coupled receptors. , 2011, Chemical reviews.

[59]  Martha S. Head,et al.  Validation Studies of the Site-Directed Docking Program LibDock , 2007, J. Chem. Inf. Model..

[60]  Frank Dickmann,et al.  Parallel High-Performance Grid Computing: Capabilities and Opportunities of a Novel Demanding Service and Business Class Allowing Highest Resource Efficiency , 2010, HealthGrid.

[61]  Michael Nilges,et al.  Grid computing for improving conformational sampling in NMR structure calculation , 2011, Bioinform..

[62]  Z. Oltvai,et al.  Identifying Ligand Binding Conformations of the β2-Adrenergic Receptor by Using Its Agonists as Computational Probes , 2012, PloS one.

[63]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[64]  Qifeng Bai,et al.  Modeling a New Water Channel That Allows SET9 to Dimethylate p53 , 2011, PLoS ONE.

[65]  Charles Perkins,et al.  Hydra: A Self Regenerating High Performance Computing Grid for Drug Discovery , 2008, J. Chem. Inf. Model..

[66]  A. Woodcock,et al.  Safety and tolerability of the novel inhaled corticosteroid fluticasone furoate in combination with the β2 agonist vilanterol administered once daily for 52 weeks in patients ≥12 years old with asthma: a randomised trial , 2013, Thorax.

[67]  Javier Bajo,et al.  Automatic knowledge extraction in sequencing analysis with multiagent system and grid computing , 2012, J. Integr. Bioinform..

[68]  Miguel Ángel Guevara-López,et al.  A Software Framework for Building Biomedical Machine Learning Classifiers through Grid Computing Resources , 2012, Journal of Medical Systems.

[69]  Ian Foster,et al.  Grid technologies empowering drug discovery. , 2002, Drug discovery today.

[70]  Michael R. Shirts,et al.  COMPUTING: Screen Savers of the World Unite! , 2000, Science.

[71]  José M. García,et al.  High-Throughput parallel blind Virtual Screening using BINDSURF , 2012, BMC Bioinformatics.

[72]  Toni Giorgino,et al.  MEMBPLUGIN: studying membrane complexity in VMD , 2014, Bioinform..

[73]  Bryan L. Roth,et al.  Structure of the human smoothened receptor bound to an antitumour agent , 2013, Nature.

[74]  Danfeng Shi,et al.  Ligand induced change of β2 adrenergic receptor from active to inactive conformation and its implication for the closed/open state of the water channel: insight from molecular dynamics simulation, free energy calculation and Markov state model analysis. , 2014, Physical chemistry chemical physics : PCCP.

[75]  R. Abagyan,et al.  Conserved binding mode of human beta2 adrenergic receptor inverse agonists and antagonist revealed by X-ray crystallography. , 2010, Journal of the American Chemical Society.

[76]  Y. Iwaki,et al.  Evaluation of Bedoradrine Sulfate (MN-221), a Novel, Highly Selective Beta2-Adrenergic Receptor Agonist for the Treatment of Asthma via Intravenous Infusion , 2012, The Journal of asthma : official journal of the Association for the Care of Asthma.

[77]  Samuel V. Angiuoli,et al.  Resources and Costs for Microbial Sequence Analysis Evaluated Using Virtual Machines and Cloud Computing , 2011, PloS one.

[78]  Kai Huang,et al.  PharmMapper server: a web server for potential drug target identification using pharmacophore mapping approach , 2010, Nucleic Acids Res..

[79]  S. Rasmussen,et al.  Crystal Structure of the β2Adrenergic Receptor-Gs protein complex , 2011, Nature.

[80]  M F Sanner,et al.  Python: a programming language for software integration and development. , 1999, Journal of molecular graphics & modelling.

[81]  S. Rasmussen,et al.  Structure of a nanobody-stabilized active state of the β2 adrenoceptor , 2010, Nature.

[82]  Tingjun Hou,et al.  Feasibility of Using Molecular Docking-Based Virtual Screening for Searching Dual Target Kinase Inhibitors , 2013, J. Chem. Inf. Model..

[83]  J. Baker.  The selectivity of beta-adrenoceptor antagonists at the human beta1, beta2 and beta3 adrenoceptors. , 2005, British journal of pharmacology.

[84]  L. Pernomian,et al.  Cross‐talk with β2‐adrenoceptors enhances ligand affinity properties from endothelial alpha1D‐adrenoceptors that mediates carotid relaxation , 2013, The Journal of pharmacy and pharmacology.

[85]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[86]  Susumu Date,et al.  Grid Heterogeneity in In-silico Experiments: An Exploration of Drug Screening Using DOCK on Cloud Environments , 2010, HealthGrid.