HBonanza: a computer algorithm for molecular-dynamics-trajectory hydrogen-bond analysis.

In the current work, we present a hydrogen-bond analysis of 2673 ligand-receptor complexes that suggests the total number of hydrogen bonds formed between a ligand and its receptor is a poor predictor of ligand potency; furthermore, even that poor prediction does not suggest a statistically significant correlation between hydrogen-bond formation and potency. While we are not the first to suggest that hydrogen bonds on average do not generally contribute to ligand binding affinities, this additional evidence is nevertheless interesting. The primary role of hydrogen bonds may instead be to ensure specificity, to correctly position the ligand within the active site, and to hold the protein active site in a ligand-friendly conformation. We also present a new computer program called HBonanza (hydrogen-bond analyzer) that aids the analysis and visualization of hydrogen-bond networks. HBonanza, which can be used to analyze single structures or the many structures of a molecular dynamics trajectory, is open source and python implemented, making it easily editable, customizable, and platform independent. Unlike many other freely available hydrogen-bond analysis tools, HBonanza provides not only a text-based table describing the hydrogen-bond network, but also a Tcl script to facilitate visualization in VMD, a popular molecular visualization program. Visualization in other programs is also possible. A copy of HBonanza can be obtained free of charge from http://www.nbcr.net/hbonanza.

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

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

[3]  G. Desiraju,et al.  Defining the hydrogen bond: an account , 2011 .

[4]  T. S. Moore,et al.  CLXXVII.—The state of amines in aqueous solution , 1912 .

[5]  K N Houk,et al.  The origins of femtomolar protein-ligand binding: hydrogen-bond cooperativity and desolvation energetics in the biotin-(strept)avidin binding site. , 2007, Journal of the American Chemical Society.

[6]  J Andrew McCammon,et al.  BINANA: a novel algorithm for ligand-binding characterization. , 2011, Journal of molecular graphics & modelling.

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

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

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

[10]  Nathan A. Baker,et al.  PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations , 2004, Nucleic Acids Res..

[11]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

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

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

[14]  F. Lombardo,et al.  Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings , 1997 .

[15]  A. Davis,et al.  Hydrogen Bonding, Hydrophobic Interactions, and Failure of the Rigid Receptor Hypothesis. , 1999, Angewandte Chemie.

[16]  A. Schnaufer,et al.  High resolution crystal structure of a key editosome enzyme from Trypanosoma brucei: RNA editing ligase 1. , 2004, Journal of molecular biology.

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

[18]  Bruce Tidor,et al.  Optimal charges in lead progression: a structure-based neuraminidase case study. , 2006, Journal of medicinal chemistry.

[19]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[20]  David S. Wishart,et al.  DrugBank: a comprehensive resource for in silico drug discovery and exploration , 2005, Nucleic Acids Res..

[21]  Robert V. Swift,et al.  Novel Naphthalene-Based Inhibitors of Trypanosoma brucei RNA Editing Ligase 1 , 2010, PLoS neglected tropical diseases.

[22]  Michael G. Lerner,et al.  Binding MOAD (Mother Of All Databases) , 2005, Proteins.

[23]  Gerhard Klebe,et al.  PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations , 2007, Nucleic Acids Res..

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

[25]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[26]  Jacob D. Durrant,et al.  NNScore: A Neural-Network-Based Scoring Function for the Characterization of Protein−Ligand Complexes , 2010, J. Chem. Inf. Model..