Interactions of omeprazole-based analogues with cytochrome P450 2C19: a computational study.

Cytochrome P450 2C19 (CYP2C19) is one of 57 drug metabolizing enzymes in humans and is responsible for the metabolism of ∼7-10% of drugs in clinical use. Recently omeprazole-based analogues were reported to be the potent inhibitors of CYP2C19 and have the potential to be used as the tool compounds for studying the substrate selectivity of CYP2C19. However, the binding modes of these compounds with CYP2C19 remain to be elucidated. In this study, a combination of molecular docking, molecular dynamics (MD), and MM/GBSA calculations was employed to systematically investigate the interactions between these compounds and CYP2C19. The binding modes of these analogues were analyzed in detail. The results indicated that the inclusion of explicit active site water molecules could improve binding energy prediction when the water molecules formed a hydrogen bonding network between the ligand and protein. We also found that the effect of active site water molecules on binding free energy prediction was dependent on the ligand binding modes. Our results unravel the interactions of these omeprazole-based analogues with CYP2C19 and might be helpful for the future design of potent CYP2C19 inhibitors with improved metabolic properties.

[1]  C. Howden,et al.  Omeprazole. Overview and opinion. , 1991, Digestive diseases and sciences.

[2]  Jeremy R. Greenwood,et al.  Epik: a software program for pKa prediction and protonation state generation for drug-like molecules , 2007, J. Comput. Aided Mol. Des..

[3]  Xiao Zhu,et al.  Fast and accurate computation schemes for evaluating vibrational entropy of proteins , 2011, J. Comput. Chem..

[4]  Alessandro Contini,et al.  Explicit Ligand Hydration Shells Improve the Correlation between MM-PB/GBSA Binding Energies and Experimental Activities. , 2013, Journal of chemical theory and computation.

[5]  N. Ford,et al.  Clopidogrel, CYP2C19, and a Black Box , 2013, Journal of clinical pharmacology.

[6]  Themis Lazaridis,et al.  Water at biomolecular binding interfaces. , 2007, Physical chemistry chemical physics : PCCP.

[7]  Woody Sherman,et al.  ConfGen: A Conformational Search Method for Efficient Generation of Bioactive Conformers , 2010, J. Chem. Inf. Model..

[8]  Charles L. Brooks,et al.  Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures , 2004, J. Comput. Chem..

[9]  Weihua Li,et al.  Reduced Catalytic Activity of P450 2A6 Mutants with Coumarin: A Computational Investigation. , 2009, Journal of chemical theory and computation.

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

[11]  Elizabeth Yuriev,et al.  Latest developments in molecular docking: 2010–2011 in review , 2013, Journal of molecular recognition : JMR.

[12]  J Andrew McCammon,et al.  Generalized Born model with a simple, robust molecular volume correction. , 2007, Journal of chemical theory and computation.

[13]  L. Wienkers,et al.  Predicting in vivo drug interactions from in vitro drug discovery data , 2005, Nature Reviews Drug Discovery.

[14]  Chris Oostenbrink,et al.  Computational prediction of drug binding and rationalisation of selectivity towards cytochromes P450 , 2008 .

[15]  Shiew-Mei Huang,et al.  Evaluation of cytochrome P450 probe substrates commonly used by the pharmaceutical industry to study in vitro drug interactions. , 2002, Drug metabolism and disposition: the biological fate of chemicals.

[16]  E. Carosati.  Modelling cytochromes P450 binding modes to predict P450 inhibition, metabolic stability and isoform selectivity. , 2013, Drug discovery today. Technologies.

[17]  L. Waskell,et al.  Theoretical study of the ligand–CYP2B4 complexes: Effect of structure on binding free energies and heme spin state , 2004, Proteins.

[18]  Chris Oostenbrink,et al.  Cooperative binding of aflatoxin B1 by cytochrome P450 3A4: a computational study. , 2014, Chemical research in toxicology.

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

[20]  Gregory D. Hawkins,et al.  Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium , 1996 .

[21]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[22]  Daniel A. Gschwend,et al.  Molecular docking towards drug discovery , 1996, Journal of molecular recognition : JMR.

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

[24]  B. Jayaram,et al.  Predicting the binding modes and sites of metabolism of xenobiotics. , 2015, Molecular bioSystems.

[25]  Jan H. Jensen,et al.  Very fast empirical prediction and rationalization of protein pKa values , 2005, Proteins.

[26]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[27]  Alessandra Magistrato,et al.  Computational Approaches Elucidate the Allosteric Mechanism of Human Aromatase Inhibition: A Novel Possible Route to Small-Molecule Regulation of CYP450s Activities? , 2014, J. Chem. Inf. Model..

[28]  R Brookmeyer,et al.  Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy. , 1997, Science.

[29]  Tingjun Hou,et al.  Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking , 2011, J. Comput. Chem..

[30]  Weihua Li,et al.  Effects of protein flexibility and active site water molecules on the prediction of sites of metabolism for cytochrome P450 2C19 substrates. , 2016, Molecular bioSystems.

[31]  M. Bodkin,et al.  Evaluation of water displacement energetics in protein binding sites with grid cell theory. , 2015, Physical chemistry chemical physics : PCCP.

[32]  K. Kim,et al.  The Effect of CYP2C19 Polymorphism on the Pharmacokinetics and Pharmacodynamics of Clopidogrel: A Possible Mechanism for Clopidogrel Resistance , 2008, Clinical pharmacology and therapeutics.

[33]  F. Guengerich,et al.  Cytochrome P450s and other enzymes in drug metabolism and toxicity , 2006, The AAPS Journal.

[34]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[35]  Eric F. Johnson,et al.  Structural Characterization of Human Cytochrome P450 2C19 , 2012, The Journal of Biological Chemistry.

[36]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[37]  W. Sippl,et al.  Docking and binding free energy calculations of sirtuin inhibitors. , 2015, European journal of medicinal chemistry.

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

[39]  U. Meyer.  Metabolic interactions of the proton-pump inhibitors lansoprazole, omeprazole and pantoprazole with other drugs , 1996, European journal of gastroenterology & hepatology.

[40]  Slobodan Rendic,et al.  Survey of Human Oxidoreductases and Cytochrome P450 Enzymes Involved in the Metabolism of Xenobiotic and Natural Chemicals , 2014, Chemical research in toxicology.

[41]  R. Skeel,et al.  Langevin stabilization of molecular dynamics , 2001 .

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

[43]  Georgios Archontis,et al.  Kinetics, in silico docking, molecular dynamics, and MM‐GBSA binding studies on prototype indirubins, KT5720, and staurosporine as phosphorylase kinase ATP‐binding site inhibitors: The role of water molecules examined , 2011, Proteins.

[44]  Jacob Kongsted,et al.  An improved method to predict the entropy term with the MM/PBSA approach , 2009, J. Comput. Aided Mol. Des..

[45]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[46]  Dieter Lang,et al.  Predicting drug metabolism: experiment and/or computation? , 2015, Nature Reviews Drug Discovery.

[47]  Colin D. Brown,et al.  LogD: lipophilicity for ionisable compounds. , 2008, Chemosphere.

[48]  L. Wienkers,et al.  Ligand-based design of a potent and selective inhibitor of cytochrome P450 2C19. , 2012, Journal of medicinal chemistry.

[49]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[50]  S. Genheden,et al.  The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities , 2015, Expert opinion on drug discovery.

[51]  Weihua Li,et al.  Computational insights into inhibitory mechanism of azole compounds against human aromatase , 2015 .

[52]  F. Guengerich.  Cytochrome p450 and chemical toxicology. , 2008, Chemical research in toxicology.

[53]  Yong-Liang Zhu,et al.  Including Explicit Water Molecules as Part of the Protein Structure in MM/PBSA Calculations , 2014, J. Chem. Inf. Model..