Independent-Trajectories Thermodynamic-Integration Free-Energy Changes for Biomolecular Systems: Determinants of H5N1 Avian Influenza Virus Neuraminidase Inhibition by Peramivir

Free-energy changes are essential physicochemical quantities for understanding most biochemical processes. Yet, the application of accurate thermodynamic-integration (TI) computation to biological and macromolecular systems is limited by finite-sampling artifacts. In this paper, we employ independent-trajectories thermodynamic-integration (IT-TI) computation to estimate improved free-energy changes and their uncertainties for (bio)molecular systems. IT-TI aids sampling statistics of the thermodynamic macrostates for flexible associating partners by ensemble averaging of multiple, independent simulation trajectories. We study peramivir (PVR) inhibition of the H5N1 avian influenza virus neuraminidase flexible receptor (N1). Binding site loops 150 and 119 are highly mobile, as revealed by N1-PVR 20-ns molecular dynamics. Due to such heterogeneous sampling, standard TI binding free-energy estimates span a rather large free-energy range, from a 19% underestimation to a 29% overestimation of the experimental reference value (−62.2 ± 1.8 kJ mol−1). Remarkably, our IT-TI binding free-energy estimate (−61.1 ± 5.4 kJ mol−1) agrees with a 2% relative difference. In addition, IT-TI runs provide a statistics-based free-energy uncertainty for the process of interest. Using ∼800 ns of overall sampling, we investigate N1-PVR binding determinants by IT-TI alchemical modifications of PVR moieties. These results emphasize the dominant electrostatic contribution, particularly through the N1 E277−PVR guanidinium interaction. Future drug development may be also guided by properly tuning ligand flexibility and hydrophobicity. IT-TI will allow estimation of relative free energies for systems of increasing size, with improved reliability by employing large-scale distributed computing.

[1]  D. Frenkel.  Free-energy calculations , 1991 .

[2]  Markus Christen,et al.  The GROMOS software for biomolecular simulation: GROMOS05 , 2005, J. Comput. Chem..

[3]  A. Mark,et al.  Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , 1994 .

[4]  J. Taylor.  An Introduction to Error Analysis , 1982 .

[5]  David J. Stevens,et al.  The structure of H5N1 avian influenza neuraminidase suggests new opportunities for drug design , 2006, Nature.

[6]  T. Straatsma,et al.  Separation‐shifted scaling, a new scaling method for Lennard‐Jones interactions in thermodynamic integration , 1994 .

[7]  Vincent S Stoll,et al.  Structure-based characterization and optimization of novel hydrophobic binding interactions in a series of pyrrolidine influenza neuraminidase inhibitors. , 2005, Journal of medicinal chemistry.

[8]  Ricky Chachra,et al.  Origins of Resistance Conferred by the R292K Neuraminidase Mutation via Molecular Dynamics and Free Energy Calculations. , 2008, Journal of chemical theory and computation.

[9]  R. Woods,et al.  Involvement of water in carbohydrate-protein binding. , 2001, Journal of the American Chemical Society.

[10]  Hideo Goto,et al.  Avian flu: Isolation of drug-resistant H5N1 virus , 2005, Nature.

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

[12]  A. Mclachlan.  Gene duplications in the structural evolution of chymotrypsin. , 1979, Journal of molecular biology.

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

[14]  J A McCammon,et al.  Determinants of ligand binding to cAMP-dependent protein kinase. , 1999, Biochemistry.

[15]  Jung-Hsin Lin,et al.  Remarkable loop flexibility in avian influenza N1 and its implications for antiviral drug design. , 2007, Journal of the American Chemical Society.

[16]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

[17]  Wilfred F van Gunsteren,et al.  Computational Analysis of the Mechanism and Thermodynamics of Inhibition of Phosphodiesterase 5A by Synthetic Ligands. , 2007, Journal of chemical theory and computation.

[18]  Michael R. Shirts,et al.  Direct calculation of the binding free energies of FKBP ligands. , 2005, The Journal of chemical physics.

[19]  M. Gilson,et al.  The statistical-thermodynamic basis for computation of binding affinities: a critical review. , 1997, Biophysical journal.

[20]  E. Clercq.  Emerging antiviral drugs , 2008 .

[21]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[22]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[23]  J A McCammon,et al.  Theory of biomolecular recognition. , 1998, Current opinion in structural biology.

[24]  Larisa V. Gubareva,et al.  Comparison of the Activities of Zanamivir, Oseltamivir, and RWJ-270201 against Clinical Isolates of Influenza Virus and Neuraminidase Inhibitor-Resistant Variants , 2001, Antimicrobial Agents and Chemotherapy.

[25]  Thanyada Rungrotmongkol,et al.  Understanding of known drug‐target interactions in the catalytic pocket of neuraminidase subtype N1 , 2008, Proteins.

[26]  Wilfred F van Gunsteren,et al.  Comparison of thermodynamic properties of coarse-grained and atomic-level simulation models. , 2007, Chemphyschem : a European journal of chemical physics and physical chemistry.

[27]  Wilfred F. van Gunsteren,et al.  An improved GROMOS96 force field for aliphatic hydrocarbons in the condensed phase , 2001, J. Comput. Chem..

[28]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[29]  Mark A. Miller,et al.  Why is it so difficult to simulate entropies, free energies, and their differences? , 2001, Accounts of chemical research.

[30]  R. Webster,et al.  Importance of Neuraminidase Active-Site Residues to the Neuraminidase Inhibitor Resistance of Influenza Viruses , 2006, Journal of Virology.

[31]  J Andrew McCammon,et al.  Dynamics, hydration, and motional averaging of a loop-gated artificial protein cavity: the W191G mutant of cytochrome c peroxidase in water as revealed by molecular dynamics simulations. , 2007, Biochemistry.

[32]  J A McCammon,et al.  Theoretical calculation of relative binding affinity in host-guest systems. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[33]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

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

[35]  A J Elliott,et al.  BCX-1812 (RWJ-270201): discovery of a novel, highly potent, orally active, and selective influenza neuraminidase inhibitor through structure-based drug design. , 2000, Journal of medicinal chemistry.

[36]  Y. Cheng,et al.  Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction. , 1973, Biochemical pharmacology.

[37]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[38]  David Ozonoff,et al.  Novel Druggable Hot Spots in Avian Influenza Neuraminidase H5N1 Revealed by Computational Solvent Mapping of a Reduced and Representative Receptor Ensemble , 2008, Chemical biology & drug design.

[39]  Wilfred F van Gunsteren,et al.  Principles of carbopeptoid folding: a molecular dynamics simulation study , 2005, Journal of peptide science : an official publication of the European Peptide Society.

[40]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[41]  Alan J. Hay,et al.  Crystal structures of oseltamivir-resistant influenza virus neuraminidase mutants , 2008, Nature.

[42]  I. Barr,et al.  Susceptibility of highly pathogenic A(H5N1) avian influenza viruses to the neuraminidase inhibitors and adamantanes. , 2007, Antiviral research.

[43]  Stewart A. Adcock,et al.  Molecular dynamics: survey of methods for simulating the activity of proteins. , 2006, Chemical reviews.

[44]  Christophe Chipot,et al.  Comprar Free Energy Calculations · Theory and Applications in Chemistry and Biology | Chipot, Christophe | 9783540736172 | Springer , 2007 .

[45]  Tim N. Heinz,et al.  Comparison of four methods to compute the dielectric permittivity of liquids from molecular dynamics simulations , 2001 .

[46]  Donald Hamelberg,et al.  Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method. , 2004, Journal of the American Chemical Society.

[47]  Philippe H. Hünenberger,et al.  A fast pairlist‐construction algorithm for molecular simulations under periodic boundary conditions , 2004, J. Comput. Chem..

[48]  N. Skeik,et al.  Influenza viruses and the evolution of avian influenza virus H5N1 , 2007, International Journal of Infectious Diseases.

[49]  Roberto D. Lins,et al.  A new GROMOS force field for hexopyranose‐based carbohydrates , 2005, J. Comput. Chem..

[50]  Pooran Chand,et al.  Comparison of the anti-influenza virus activity of cyclopentane derivatives with oseltamivir and zanamivir in vivo. , 2005, Bioorganic & medicinal chemistry.

[51]  E. R. Cohen.  An Introduction to Error Analysis: The Study of Uncertainties in Physical Measurements , 1998 .