A framework for virtual screening

Recent advances in Bioinformatics and in Computer simulation and modelling have positively impacted the drug discovery process by turning viable the rational drug design (RDD). One of the major challenges in RDD is the understanding about protein-ligand interaction simulated at the atomic level by molecular docking algorithms. Virtual screening (VS) is defined as a computational approach applied to the analyses of large libraries of chemical structures in order to identify possible drug candidates to a target. The major challenge of VS based on molecular docking is the time required to run each experiment and the countless parameters and characteristics that should be defined by the researcher such as: the target(s) receptor, one or a set of ligands, the receptor binding site and so on. In order to perform more realistic docking simulations it is also necessary to account for the receptor and ligand flexibility. Therefore, this paper presents a framework for VS, where the user configure an experiment in a Web based platform informing the path of input and output files as well as the size, center and variation of the binding site(s). Then, the proposed framework generates a Python script that performs the VS experiment on the users personal computer. We expect that researchers from diverse backgrounds as Biology, Physics, Pharmacy, etc. can easily prepare VS experiments without the necessity of learning how to write scripts. To validate our proposed framework we performed five different case studies considering the AcrB protein as target receptor. All the case studies were easily realized using the proposed framework. The results show that the framework is effective to configure the VS experiments with different characteristics. Besides, the experiments can help on the search for new drug candidates for this important target.

[1]  Yanli Wang,et al.  Structure-Based Virtual Screening for Drug Discovery: a Problem-Centric Review , 2012, The AAPS Journal.

[2]  Holger Gohlke,et al.  Target flexibility: an emerging consideration in drug discovery and design. , 2008, Journal of medicinal chemistry.

[3]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[4]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

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

[6]  I. Kuntz Structure-Based Strategies for Drug Design and Discovery , 1992, Science.

[7]  Satoshi Murakami,et al.  Crystal structure of bacterial multidrug efflux transporter AcrB , 2002, Nature.

[8]  Karina S Machado,et al.  Effect of the explicit flexibility of the InhA enzyme from Mycobacterium tuberculosis in molecular docking simulations , 2011, BMC Genomics.

[9]  Janet Dodd DIGITAL BRIEFS: NEW SOFTWARE AND WEBSITES FOR THE CHEMICAL ENTERPRISE , 2005 .

[10]  Brian K. Shoichet,et al.  Virtual screening of chemical libraries , 2004, Nature.

[11]  Karina S Machado,et al.  FReDoWS: a method to automate molecular docking simulations with explicit receptor flexibility and snapshots selection , 2011, BMC Genomics.

[12]  M. Mezei,et al.  Molecular docking: a powerful approach for structure-based drug discovery. , 2011, Current computer-aided drug design.

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

[14]  Scott P. Brown,et al.  A unified, probabilistic framework for structure- and ligand-based virtual screening. , 2011, Journal of medicinal chemistry.

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

[16]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[17]  J M Thornton,et al.  LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. , 1995, Protein engineering.