Role of Water in Molecular Docking Simulations of Cytochrome P450 2D6

Active-site water molecules form an important component in biological systems, facilitating promiscuous binding or an increase in specificity and affinity. Taking water molecules into account in computational approaches to drug design or site-of-metabolism predictions is currently far from straightforward. In this study, the effects of including water molecules in molecular docking simulations of the important metabolic enzyme cytochrome P450 2D6 are investigated. The structure and dynamics of water molecules that are present in the active site simultaneously with a selected substrate are described, and based on this description, water molecules are selected to be included in docking experiments into multiple protein conformations. Apart from the parent substrate, 11 similar and 53 dissimilar substrates are included to investigate the transferability of active-site hydration sites between substrates. The role of water molecules appears to be highly dependent on the protein conformation and the substrate.

[1]  B. McConkey,et al.  The performance of current methods in ligand-protein docking , 2002 .

[2]  J. Ladbury.  Just add water! The effect of water on the specificity of protein-ligand binding sites and its potential application to drug design. , 1996, Chemistry & biology.

[3]  Ricardo L. Mancera,et al.  Ligand-Protein Docking with Water Molecules , 2008, J. Chem. Inf. Model..

[4]  Philip M. Dean,et al.  Hydration in drug design. 1. Multiple hydrogen-bonding features of water molecules in mediating protein-ligand interactions , 1995, J. Comput. Aided Mol. Des..

[5]  Philip M. Dean,et al.  Hydration in drug design. 2. Influence of local site surface shape on water binding , 1995, J. Comput. Aided Mol. Des..

[6]  Pietro Cozzini,et al.  Water: How to evaluate its contribution in protein–ligand interactions , 2006 .

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

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

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

[10]  Julien Michel,et al.  Prediction of the water content in protein binding sites. , 2009, The journal of physical chemistry. B.

[11]  Marcel L. Verdonk,et al.  Sensitivity of molecular docking to induced fit effects in influenza virus neuraminidase , 2002, J. Comput. Aided Mol. Des..

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

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

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

[15]  J. Dawson,et al.  Heme-Containing Oxygenases. , 1996, Chemical reviews.

[16]  Chris de Graaf,et al.  Binding mode prediction of cytochrome p450 and thymidine kinase protein-ligand complexes by consideration of water and rescoring in automated docking. , 2005, Journal of medicinal chemistry.

[17]  Stewart B Kirton,et al.  Prediction of binding modes for ligands in the cytochromes P450 and other heme‐containing proteins , 2005, Proteins.

[18]  Christopher R. Corbeil,et al.  Docking Ligands into Flexible and Solvated Macromolecules. 3. Impact of Input Ligand Conformation, Protein Flexibility, and Water Molecules on the Accuracy of Docking Programs , 2009, J. Chem. Inf. Model..

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

[20]  Richard D. Taylor,et al.  Modeling water molecules in protein-ligand docking using GOLD. , 2005, Journal of medicinal chemistry.

[21]  G. Cruciani,et al.  MetaSite: understanding metabolism in human cytochromes from the perspective of the chemist. , 2005, Journal of medicinal chemistry.

[22]  Volkhard Helms,et al.  Protein dynamics tightly connected to the dynamics of surrounding and internal water molecules. , 2007, Chemphyschem : a European journal of chemical physics and physical chemistry.

[23]  F. Guengerich,et al.  A malleable catalyst dominates the metabolism of drugs , 2006, Proceedings of the National Academy of Sciences.

[24]  Christopher R. Corbeil,et al.  Towards the development of universal, fast and highly accurate docking/scoring methods: a long way to go , 2008, British journal of pharmacology.

[25]  Pedro Alexandrino Fernandes,et al.  Protein–ligand docking: Current status and future challenges , 2006, Proteins.

[26]  Chris Oostenbrink,et al.  Catalytic site prediction and virtual screening of cytochrome P450 2D6 substrates by consideration of water and rescoring in automated docking. , 2006, Journal of medicinal chemistry.

[27]  Niu Huang,et al.  Exploiting ordered waters in molecular docking. , 2008, Journal of medicinal chemistry.

[28]  Jozef Hritz,et al.  Impact of plasticity and flexibility on docking results for cytochrome P450 2D6: a combined approach of molecular dynamics and ligand docking. , 2008, Journal of medicinal chemistry.

[29]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[30]  Kenneth R. Harris,et al.  Pressure and temperature dependence of the self diffusion coefficient of water and oxygen-18 water , 1980 .

[31]  B. Berne,et al.  Role of the active-site solvent in the thermodynamics of factor Xa ligand binding. , 2008, Journal of the American Chemical Society.

[32]  K. Korttila,et al.  Diazepam effects and kinetics in Caucasians and Orientals , 1981, Clinical pharmacology and therapeutics.

[33]  Chris Oostenbrink,et al.  Computational prediction of drug binding and rationalisation of selectivity towards cytochromes P450. , 2008, Expert opinion on drug metabolism & toxicology.

[34]  Ilme Schlichting,et al.  Structure and chemistry of cytochrome P450. , 2005, Chemical reviews.

[35]  G. V. Paolini,et al.  Empirical scoring functions: I. The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes , 1997, J. Comput. Aided Mol. Des..

[36]  Barry C. Jones,et al.  DRUG-DRUG INTERACTIONS FOR UDP-GLUCURONOSYLTRANSFERASE SUBSTRATES: A PHARMACOKINETIC EXPLANATION FOR TYPICALLY OBSERVED LOW EXPOSURE (AUCI/AUC) RATIOS , 2004, Drug Metabolism and Disposition.

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

[38]  Lars Olsen,et al.  Dynamics of water molecules in the active-site cavity of human cytochromes P450. , 2007, The journal of physical chemistry. B.