Optimizing QSAR Models for Predicting Ligand Binding to the Drug-Metabolizing Cytochrome P450 Isoenzyme CYP2D6

Computational prediction of cytochrome P450 inhibition and induction

Cytochrome P450 (CYP) enzymes play an important role in the phase I metabolism of many xenobiotics. Most drug-drug interactions (DDIs) associated with CYP are caused by either CYP inhibition or induction. The early detection of potential DDIs is highly desirable in the pharmaceutical industry because DDIs can cause serious adverse events, which can lead to poor patient health and drug development failures. Recently, many computational studies predicting CYP inhibition and induction have been reported. The current computational modeling approaches for CYP metabolism are classified as ligand- and structure-based; various techniques, such as quantitative structure-activity relationships, machine learning, docking, and molecular dynamic simulation, are involved in both the approaches. Recently, combining these two approaches have resulted in improvements in the prediction accuracy of DDIs. In this review, we present important, recent developments in the computational prediction of the inhibition of four clinically crucial CYP isoforms (CYP1A2, 2C9, 2D6, and 3A4) and three nuclear receptors (aryl hydrocarbon receptor, constitutive androstane receptor, and pregnane X receptor) involved in the induction of CYP1A2, 2B6, and 3A4, respectively.

Exposition and reactivity optimization to predict sites of metabolism in chemicals

Chemical modifications of drugs induced by phase I biotransformations significantly affect their pharmacokinetic properties. Because the metabolites produced can themselves have a pharmacological effect and an intrinsic toxicity, medicinal chemists need to accurately predict the sites of metabolism (SoM) of drugs as early as possible. However, site of metabolism prediction is rarely accompanied by a prediction of the relative abundance of the various metabolites. Such a prediction would be a great help in the study of drug– drug interactions and in the process of reducing the toxicity of potential drug candidates. The aim of this paper is to present recent developments in the prediction of xenobiotic metabolism and to use concrete examples to explain the computational mechanism employed.

Development of a computational tool to rival experts in the prediction of sites of metabolism of xenobiotics by p450s

The metabolism of xenobiotics--and more specifically drugs--in the liver is a critical process controlling their half-life. Although there exist experimental methods, which measure the metabolic stability of xenobiotics and identify their metabolites, developing higher throughput predictive methods is an avenue of research. It is expected that predicting the chemical nature of the metabolites would be an asset for designing safer drugs and/or drugs with modulated half-lives. We have developed IMPACTS (In-silico Metabolism Prediction by Activated Cytochromes and Transition States), a computational tool combining docking to metabolic enzymes, transition state modeling, and rule-based substrate reactivity prediction to predict the site of metabolism (SoM) of xenobiotics. Its application to sets of CYP1A2, CYP2C9, CYP2D6, and CYP3A4 substrates and comparison to experts' predictions demonstrates its accuracy and significance. IMPACTS identified an experimentally observed SoM in the top 2 predicted sites for 77% of the substrates, while the accuracy of biotransformation experts' prediction was 65%. Application of IMPACTS to external sets and comparison of its accuracy to those of eleven other methods further validated the method implemented in IMPACTS.

Identification of cytochrome P450 2D6 and 2C9 substrates and inhibitors by QSAR analysis

This paper presents four new QSAR models for CYP2C9 and CYP2D6 substrate recognition and inhibitor identification based on human clinical data. The models were used to screen a large data set of environmental chemicals for CYP activity, and to analyze the frequency of CYP activity among these compounds. A large fraction of these chemicals were found to be CYP active, and thus potentially capable of affecting human physiology. 20% of the compounds within applicability domain of the models were predicted to be CYP2C9 substrates, and 17% to be inhibitors. The corresponding numbers for CYP2D6 were 9% and 21%. Where the majority of CYP2C9 active compounds were predicted to be both a substrate and an inhibitor at the same time, the CYP2D6 active compounds were primarily predicted to be only inhibitors. It was demonstrated that the models could identify compound classes with a high occurrence of specific CYP activity. An overrepresentation was seen for poly-aromatic hydrocarbons (group of procarcinogens) among CYP2C9 active and mutagenic compounds compared to CYP2C9 inactive and mutagenic compounds. The mutagenicity was predicted with a QSAR model based on Ames in vitro test data.