In silico prediction of potential chemical reactions mediated by human enzymes

BackgroundAdministered drugs are often converted into an ineffective or activated form by enzymes in our body. Conventional in silico prediction approaches focused on therapeutically important enzymes such as CYP450. However, there are more than thousands of different cellular enzymes that potentially convert administered drug into other forms.ResultWe developed an in silico model to predict which of human enzymes including metabolic enzymes as well as CYP450 family can catalyze a given chemical compound. The prediction is based on the chemical and physical similarity between known enzyme substrates and a query chemical compound. Our in silico model was developed using multiple linear regression and the model showed high performance (AUC = 0.896) despite of the large number of enzymes. When evaluated on a test dataset, it also showed significantly high performance (AUC = 0.746). Interestingly, evaluation with literature data showed that our model can be used to predict not only enzymatic reactions but also drug conversion and enzyme inhibition.ConclusionOur model was able to predict enzymatic reactions of a query molecule with a high accuracy. This may foster to discover new metabolic routes and to accelerate the computational development of drug candidates by enabling the prediction of the potential conversion of administered drugs into active or inactive forms.

[1]  David E. Williams,et al.  The Role of Flavin-Containing Monooxygenase (FMO) in the Metabolism of Tamoxifen and Other Tertiary Amines , 2006, Drug metabolism reviews.

[2]  C. Volk OCTs, OATs, and OCTNs: structure and function of the polyspecific organic ion transporters of the SLC22 family , 2014 .

[3]  Berith F. Jensen,et al.  In silico prediction of cytochrome P450 2D6 and 3A4 inhibition using Gaussian kernel weighted k-nearest neighbor and extended connectivity fingerprints, including structural fragment analysis of inhibitors versus noninhibitors. , 2007, Journal of medicinal chemistry.

[4]  N. Sucher,et al.  Searching for synergy in silico, in vitro and in vivo , 2014 .

[5]  Antony J. Williams,et al.  ChemSpider:: An Online Chemical Information Resource , 2010 .

[6]  Chris Oostenbrink,et al.  Prediction of cytochrome P450 mediated metabolism. , 2015, Advanced drug delivery reviews.

[7]  David S. Wishart,et al.  HMDB 3.0—The Human Metabolome Database in 2013 , 2012, Nucleic Acids Res..

[8]  J. Gsponer,et al.  Systems-wide Identification of cis-Regulatory Elements in Proteins. , 2016, Cell systems.

[9]  Christian von Mering,et al.  STRING: known and predicted protein–protein associations, integrated and transferred across organisms , 2004, Nucleic Acids Res..

[10]  Jie Shen,et al.  admetSAR: A Comprehensive Source and Free Tool for Assessment of Chemical ADMET Properties , 2012, J. Chem. Inf. Model..

[11]  A. Cederbaum,et al.  CYP2E1 and oxidative liver injury by alcohol. , 2008, Free radical biology & medicine.

[12]  H. Wulff,et al.  Therapeutic potential of KCa3.1 blockers: recent advances and promising trends , 2010, Expert review of clinical pharmacology.

[13]  R. Tukey,et al.  Caffeine metabolism by human hepatic cytochromes P450: contributions of 1A2, 2E1 and 3A isoforms. , 1994, Biochemical pharmacology.

[14]  CHUN WEI YAP,et al.  PaDEL‐descriptor: An open source software to calculate molecular descriptors and fingerprints , 2011, J. Comput. Chem..

[15]  G. Tucker,et al.  Variable contribution of cytochromes P450 2D6, 2C9 and 3A4 to the 4-hydroxylation of tamoxifen by human liver microsomes. , 1997, Biochemical pharmacology.

[16]  Jean-Loup Faulon,et al.  Genome scale enzyme–metabolite and drug–target interaction predictions using the signature molecular descriptor , 2008 .

[17]  Stephen R. Heller,et al.  InChI - the worldwide chemical structure identifier standard , 2013, Journal of Cheminformatics.

[18]  Y. Funae,et al.  Role of human cytochrome P4502A6 in C-oxidation of nicotine. , 1996, Drug metabolism and disposition: the biological fate of chemicals.

[19]  Bing Niu,et al.  Prediction of Substrate-Enzyme-Product Interaction Based on Molecular Descriptors and Physicochemical Properties , 2013, BioMed research international.

[20]  Q. Feng,et al.  The Role of Genipin and Geniposide in Liver Diseases: A Review , 2013 .

[21]  David S. Wishart,et al.  DrugBank 4.0: shedding new light on drug metabolism , 2013, Nucleic Acids Res..

[22]  G. Pérez,et al.  Dual Effect of Tamoxifen on Arterial KCa Channels Does Not Depend on the Presence of the β1 Subunit* , 2005, Journal of Biological Chemistry.

[23]  R. Blantz,et al.  Polyamine transport system mediates agmatine transport in mammalian cells. , 2001, American journal of physiology. Cell physiology.

[24]  Ming Wen,et al.  Deep-Learning-Based Drug-Target Interaction Prediction. , 2017, Journal of proteome research.

[25]  Martin Mozina,et al.  Orange: data mining toolbox in python , 2013, J. Mach. Learn. Res..

[26]  Antje Chang,et al.  BRENDA in 2017: new perspectives and new tools in BRENDA , 2016, Nucleic Acids Res..

[27]  Robert F. Tate,et al.  Correlation Between a Discrete and a Continuous Variable. Point-Biserial Correlation , 1954 .

[28]  Teruko Imai,et al.  Substrate Specificity of Carboxylesterase Isozymes and Their Contribution to Hydrolase Activity in Human Liver and Small Intestine , 2006, Drug Metabolism and Disposition.

[29]  Z R Li,et al.  Application of support vector machines to in silico prediction of cytochrome p450 enzyme substrates and inhibitors. , 2006, Current topics in medicinal chemistry.

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

[31]  David Weininger,et al.  SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules , 1988, J. Chem. Inf. Comput. Sci..