AutoDock VinaXB: implementation of XBSF, new empirical halogen bond scoring function, into AutoDock Vina

BackgroundHalogen bonding has recently come to play as a target for lead optimization in rational drug design. However, most docking program don’t account for halogen bonding in their scoring functions and are not able to utilize this new approach. In this study a new and improved halogen bonding scoring function (XBSF) is presented along with its implementation in the AutoDock Vina molecular docking software. This new improved program is termed as AutoDock VinaXB, where XB stands for the halogen bonding parameters that were added.ResultsXBSF scoring function is derived based on the X···A distance and C–X···A angle of interacting atoms. The distance term was further corrected to account for the polar flattening effect of halogens. A total of 106 protein-halogenated ligand complexes were tested and compared in terms of binding affinity and docking poses using Vina and VinaXB. VinaXB performed superior to Vina in the majority of instances. VinaXB was closer to native pose both above and below 2 Å deviation categories almost twice as frequently as Vina.ConclusionsImplementation of XBSF into AutoDock Vina has been shown to improve the accuracy of the docking result with regards to halogenated ligands. AutoDock VinaXB addresses the issues of halogen bonds that were previously being scored unfavorably due to repulsion factors, thus effectively lowering the output RMSD values.

[1]  Renxiao Wang,et al.  The PDBbind database: methodologies and updates. , 2005, Journal of medicinal chemistry.

[2]  Mahmoud A A Ibrahim,et al.  AMBER empirical potential describes the geometry and energy of noncovalent halogen interactions better than advanced semiempirical quantum mechanical method PM6-DH2X. , 2012, The journal of physical chemistry. B.

[3]  Timothy Clark,et al.  Halogen bonding: the σ-hole , 2007 .

[4]  Patric Schyman,et al.  Treatment of Halogen Bonding in the OPLS-AA Force Field; Application to Potent Anti-HIV Agents. , 2012, Journal of chemical theory and computation.

[5]  Thomas Lengauer,et al.  Computational methods for biomolecular docking. , 1996, Current opinion in structural biology.

[6]  A. Rappé,et al.  Scalable Anisotropic Shape and Electrostatic Models for Biological Bromine Halogen Bonds. , 2012, Journal of chemical theory and computation.

[7]  Pavel Hobza,et al.  Plugging the explicit σ-holes in molecular docking. , 2013, Chemical communications.

[8]  Ajay N. Jain.  Scoring noncovalent protein-ligand interactions: A continuous differentiable function tuned to compute binding affinities , 1996, J. Comput. Aided Mol. Des..

[9]  Pavel Hobza,et al.  On Extension of the Current Biomolecular Empirical Force Field for the Description of Halogen Bonds. , 2012, Journal of chemical theory and computation.

[10]  A. Rappé,et al.  Force Field Model of Periodic Trends in Biomolecular Halogen Bonds. , 2015, The journal of physical chemistry. B.

[11]  Mahmoud A. A. Ibrahim,et al.  Molecular mechanical study of halogen bonding in drug discovery , 2011, J. Comput. Chem..

[12]  Mahmoud A A Ibrahim,et al.  Molecular mechanical perspective on halogen bonding , 2012, Journal of Molecular Modeling.

[13]  Eric Westhof,et al.  Halogen bonds in biological molecules. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

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

[15]  W. L. Jorgensen.  The Many Roles of Computation in Drug Discovery , 2004, Science.

[16]  Markus O. Zimmermann,et al.  Evaluating the Potential of Halogen Bonding in Molecular Design: Automated Scaffold Decoration Using the New Scoring Function XBScore , 2015, J. Chem. Inf. Model..

[17]  Weiliang Zhu,et al.  A knowledge-based halogen bonding scoring function for predicting protein-ligand interactions , 2013, Journal of Molecular Modeling.

[18]  S. Kazmirski,et al.  Beyond the MEK-pocket: can current MEK kinase inhibitors be utilized to synthesize novel type III NCKIs? Does the MEK-pocket exist in kinases other than MEK? , 2009, Bioorganic & medicinal chemistry letters.

[19]  Renxiao Wang,et al.  The PDBbind database: collection of binding affinities for protein-ligand complexes with known three-dimensional structures. , 2004, Journal of medicinal chemistry.

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

[21]  Lili Wang,et al.  A polarizable ellipsoidal force field for halogen bonds , 2013, J. Comput. Chem..

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

[23]  Suman Sirimulla,et al.  Halogen Interactions in Protein-Ligand Complexes: Implications of Halogen Bonding for Rational Drug Design , 2013, J. Chem. Inf. Model..