Strategies to Calculate Water Binding Free Energies in Protein–Ligand Complexes

Water molecules are commonplace in protein binding pockets, where they can typically form a complex between the protein and a ligand or become displaced upon ligand binding. As a result, it is often of great interest to establish both the binding free energy and location of such molecules. Several approaches to predicting the location and affinity of water molecules to proteins have been proposed and utilized in the literature, although it is often unclear which method should be used under what circumstances. We report here a comparison between three such methodologies, Just Add Water Molecules (JAWS), Grand Canonical Monte Carlo (GCMC), and double-decoupling, in the hope of understanding the advantages and limitations of each method when applied to enclosed binding sites. As a result, we have adapted the JAWS scoring procedure, allowing the binding free energies of strongly bound water molecules to be calculated to a high degree of accuracy, requiring significantly less computational effort than more rigorous approaches. The combination of JAWS and GCMC offers a route to a rapid scheme capable of both locating and scoring water molecules for rational drug design

Drugs against avian influenza A virus : design of novel sulfonate inhibitors of neuraminidase N1

The outbreak of avian influenza A (H5N1) virus has raised a global concern for both the animal as well as human health. Besides vaccination, that may not achieve full protection in certain groups of patients, inhibiting neuraminidase or the transmembrane protein M2 represents the main measure of controlling the disease. Due to alarming emergence of influenza virus strains resistant to the currently available drugs, development of new neuraminidase N1 inhibitors is of utmost importance. The present paper provides an overview of the recent advances in the design of new antiviral drugs against avian influenza. It also reports findings in binding free energy calculations for nine neuraminidase N1 inhibitors (oseltamivir, zanamivir, and peramivir -carboxylate, -phosphonate, and -sulfonate) using the Linear Interaction Energy method. Molecular dynamics simulations of these inhibitors were performed in a free and two bound states the so called open and closed conformations of neuraminidase N1. Obtained results successfully reproduce the experimental binding affinities of the already known neuraminidase N1 inhibitors, i.e. peramivir being a stronger binder than zanamivir that is in turn stronger binder than oseltamivir, or phosphonate inhibitors being stronger binders than their carboxylate analogues. In addition, the newly proposed sulfonate inhibitors are predicted to be the strongest binders - a fact to be confirmed by their chemical synthesis and a subsequent test of their biological activity. Finally, contributions of individual inhibitor moieties to the overall binding affinity are explicitly evaluated to assist further drug development towards inhibition of the H5N1 avian influenza A virus