A novel approach to predicting P450 mediated drug metabolism. CYP2D6 catalyzed N-dealkylation reactions and qualitative metabolite predictions using a combined protein and pharmacophore model for CYP2D6

Strategies for Assessing Acceptable Intakes for Novel N-Nitrosamines Derived from Active Pharmaceutical Ingredients

The detection of N-nitrosamines, derived from solvents and reagents and, on occasion, the active pharmaceutical ingredient (API) at higher than acceptable levels in drug products, has led regulators to request a detailed review for their presence in all medicinal products. In the absence of rodent carcinogenicity data for novel N-nitrosamines derived from amine-containing APIs, a conservative class limit of 18 ng/day (based on the most carcinogenic N-nitrosamines) or the derivation of acceptable intakes (AIs) using structurally related surrogates with robust rodent carcinogenicity data is recommended. The guidance has implications for the pharmaceutical industry given the vast number of marketed amine-containing drugs. In this perspective, the rate-limiting step in N-nitrosamine carcinogenicity, involving cytochrome P450-mediated α-carbon hydroxylation to yield DNA-reactive diazonium or carbonium ion intermediates, is discussed with reference to the selection of read-across analogs to derive AIs. Risk-mitigation strategies for managing putative N-nitrosamines in the preclinical discovery setting are also presented.

Metabolic stability studies of lead compounds supported by separation techniques and chemometrics analysis

With metabolism being one of the main routes of drug elimination from the body (accounting for removal of around 75% of known drugs), it is crucial to understand and study metabolic stability of drug candidates. Metabolically unstable compounds are uncomfortable to administer (requiring repetitive dosage during therapy), while overly stable drugs increase risk of adverse drug reactions. Additionally, biotransformation reactions can lead to formation of toxic or pharmacologically active metabolites (either less-active than parent drug, or even with different action). There were numerous approaches in estimating metabolic stability, including in vitro, in vivo, in silico, and high-throughput screening to name a few. This review aims at describing separation techniques used in in vitro metabolic stability estimation, as well as chemometric techniques allowing for creation of predictive models which enable high-throughput screening approach for estimation of metabolic stability. With a very low rate of drug approval, it is important to understand in silico methods that aim at supporting classical in vitro approach. Predictive models that allow assessment of certain biological properties of drug candidates allow for cutting not only cost, but also time required to synthesize compounds predicted to be unstable or inactive by in silico models.

Scaffold-hopping as a strategy to address metabolic liabilities of aromatic compounds

Understanding and minimizing oxidative metabolism of aromatic compounds is a key hurdle in lead optimization. Metabolic processes not only clear compounds from the body, but they can also transform parent compounds into reactive metabolites. One particularly useful strategy when addressing metabolically labile or oxidation-prone structures is scaffold-hopping. Replacement of an aromatic system with a more electron-deficient ring system can often increase robustness towards cytochrome P450-mediated oxidation while conserving the structural requirements of the pharmacophore. The most common example of this substitution strategy, replacement of a phenyl ring with a pyridyl substituent, is prevalent throughout the literature; however scaffold-hopping encompasses a much wider scope of heterocycle replacement. This review will showcase recent examples where different scaffold-hopping approaches were used to reduce metabolic clearance or block the formation of reactive metabolites. Additionally, we will highlight considerations that should be made to garner the most benefit from a scaffold-hopping strategy for lead optimization.

Linker Variation and Structure-Activity Relationship Analyses of Carboxylic Acid-based Small Molecule STAT3 Inhibitors

The molecular determinants for the activities of the reported benzoic acid (SH4-54), salicylic acid (BP-1-102), and benzohydroxamic acid (SH5-07)-based STAT3 inhibitors were investigated to design optimized analogues. All three leads are based on an N-methylglycinamide scaffold, with its two amine groups condensed with three different functionalities. The three functionalities and the CH₂ group of the glycinamide scaffold were separately modified. The replacement of the pentafluorobenzene or cyclohexylbenzene, or replacing the benzene ring of the aromatic carboxylic or hydroxamic acid motif with heterocyclic components (containing nitrogen and oxygen elements) all decreased potency. Notably, the Ala-linker analogues, 1a and 2v, and the Pro-based derivative 5d, all with (R)-configuration at the chiral center, had improved inhibitory activity and selectivity against STAT3 DNA-binding activity in vitro, with IC₅₀ of 3.0 ± 0.9, 1.80 ± 0.94, and 2.4 ± 0.2 μM, respectively. Compounds 1a, 2v, 5d, and other analogues inhibited constitutive STAT3 phosphorylation and activation in human breast cancer and melanoma lines, and blocked tumor cell viability, growth, colony formation, and migration in vitro. Pro-based analogue, 5h, with a relatively polar tetrahydropyranyl (THP) ring, instead of the cyclohexyl, showed improved permeability. In general, the (R)-configuration Pro-based analogs showed the overall best profile, including physicochemical properties (e.g., microsomal metabolic stability, Caco-2 permeability), and in particular, 5d showed improved tumor-cell specificity.

Drug Metabolism in Preclinical Drug Development: A Survey of the Discovery Process, Toxicology, and Computational Tools

BACKGROUND

While establishing efficacy in translational models and humans through clinically-relevant endpoints for disease is of great interest, assessing the potential toxicity of a putative therapeutic drug is critical. Toxicological assessments in the pre-clinical discovery phase help to avoid future failure in the clinical phases of drug development. Many in vitro assays exist to aid in modular toxicological assessment, such as hepatotoxicity and genotoxicity. While these methods have provided tremendous insight into human toxicity by investigational new drugs, they are expensive, require substantial resources, and do not account for pharmacogenomics as well as critical ADME properties. Computational tools can fill this niche in toxicology if in silico models are accurate in relating drug molecular properties to toxicological endpoints as well as reliable in predicting important drug-target interactions that mediate known adverse events or adverse outcome pathways (AOPs).

METHODS

We undertook an unstructured search of multiple bibliographic databases for peer-reviewed literature regarding computational methods in predictive toxicology for in silico drug discovery. As this review paper is meant to serve as a survey of available methods for the interested reader, no focused criteria were applied. Literature chosen was based on the writers' expertise and intent in communicating important aspects of in silico toxicology to the interested reader.

CONCLUSION

This review provides a purview of computational methods of pre-clinical toxicologic assessments for novel small molecule drugs that may be of use for novice and experienced investigators as well as academic and commercial drug discovery entities.

Construction of Metabolism Prediction Models for CYP450 3A4, 2D6, and 2C9 Based on Microsomal Metabolic Reaction System

During the past decades, there have been continuous attempts in the prediction of metabolism mediated by cytochrome P450s (CYP450s) 3A4, 2D6, and 2C9. However, it has indeed remained a huge challenge to accurately predict the metabolism of xenobiotics mediated by these enzymes. To address this issue, microsomal metabolic reaction system (MMRS)—a novel concept, which integrates information about site of metabolism (SOM) and enzyme—was introduced. By incorporating the use of multiple feature selection (FS) techniques (ChiSquared (CHI), InfoGain (IG), GainRatio (GR), Relief) and hybrid classification procedures (Kstar, Bayes (BN), K-nearest neighbours (IBK), C4.5 decision tree (J48), RandomForest (RF), Support vector machines (SVM), AdaBoostM1, Bagging), metabolism prediction models were established based on metabolism data released by Sheridan et al. Four major biotransformations, including aliphatic C-hydroxylation, aromatic C-hydroxylation, N-dealkylation and O-dealkylation, were involved. For validation, the overall accuracies of all four biotransformations exceeded 0.95. For receiver operating characteristic (ROC) analysis, each of these models gave a significant area under curve (AUC) value >0.98. In addition, an external test was performed based on dataset published previously. As a result, 87.7% of the potential SOMs were correctly identified by our four models. In summary, four MMRS-based models were established, which can be used to predict the metabolism mediated by CYP3A4, 2D6, and 2C9 with high accuracy.

Predicting the binding modes and sites of metabolism of xenobiotics

Metabolism studies are an essential integral part of ADMET profiling of drug candidates to evaluate their safety and efficacy. Cytochrome P-450 (CYP) metabolizes a wide variety of xenobiotics/drugs. The binding modes of these compounds with CYP and their intrinsic reactivities decide the metabolic products. We report here a novel computational protocol, which comprises docking of ligands to heme-containing CYPs and prediction of binding energies through a newly developed scoring function, followed by analyses of the docked structures and molecular orbitals of the ligand molecules, for predicting the sites of metabolism (SOM) of ligands. The calculated binding free energies of 121 heme-containing protein-ligand docked complexes yielded a correlation coefficient of 0.84 against experiment. Molecular orbital analyses of the resultant top three unique poses of the docked complexes provided a success rate of 87% in identifying the experimentally known sites of metabolism of the xenobiotics. The SOM prediction methodology is freely accessible at .

Further Optimization and Evaluation of Bioavailable, Mixed-Efficacy μ-Opioid Receptor (MOR) Agonists/δ-Opioid Receptor (DOR) Antagonists: Balancing MOR and DOR Affinities

In a previously described peptidomimetic series, we reported the development of bifunctional μ-opioid receptor (MOR) agonist and δ-opioid receptor (DOR) antagonist ligands with a lead compound that produced antinociception for 1 h after intraperitoneal administration in mice. In this paper, we expand on our original series by presenting two modifications, both of which were designed with the following objectives: (1) probing bioavailability and improving metabolic stability, (2) balancing affinities between MOR and DOR while reducing affinity and efficacy at the κ-opioid receptor (KOR), and (3) improving in vivo efficacy. Here, we establish that, through N-acetylation of our original peptidomimetic series, we are able to improve DOR affinity and increase selectivity relative to KOR while maintaining the desired MOR agonist/DOR antagonist profile. From initial in vivo studies, one compound (14a) was found to produce dose-dependent antinociception after peripheral administration with an improved duration of action of longer than 3 h.

Identification and in silico prediction of metabolites of the model compound, tebufenozide by human CYP3A4 and CYP2C19

The metabolites of tebufenozide, a model compound, formed by the yeast-expressed human CYP3A4 and CYP2C19 were identified to clarify the substrate recognition mechanism of the human cytochrome P450 (CYP) isozymes. We then determined whether tebufenozide metabolites may be predicted in silico. Hydrogen abstraction energies were calculated with the density functional theory method B3LYP/6-31G(∗). A docking simulation was performed using FRED software. Several alkyl sites of tebufenozide were hydroxylated by CYP3A4 whereas only one site was modified by CYP2C19. The accessibility of each site of tebufenozide to the reaction center of CYP enzymes and the susceptibility of each hydrogen atom for metabolism by CYP enzymes were evaluated by a docking simulation and hydrogen abstraction energy estimation, respectively.

Antihypertensive drugs metabolism: an update to pharmacokinetic profiles and computational approaches

Drug discovery and development is a high-risk enterprise that requires significant investments in capital, time and scientific expertise. The studies of xenobiotic metabolism remain as one of the main topics in the research and development of drugs, cosmetics and nutritional supplements.

Antihypertensive drugs are used for the treatment of high blood pressure, which is one the most frequent symptoms of the patients that undergo cardiovascular diseases such as myocardial infraction and strokes. In current cardiovascular disease pharmacology, four drug clusters - Angiotensin Converting Enzyme Inhibitors, Beta-Blockers, Calcium Channel Blockers and Diuretics - cover the major therapeutic characteristics of the most antihypertensive drugs. The pharmacokinetic and specifically the metabolic profile of the antihypertensive agents are intensively studied because of the broad inter-individual variability on plasma concentrations and the diversity on the efficacy response especially due to the P450 dependent metabolic status they present.

Several computational methods have been developed with the aim to: (i) model and better understand the human drug metabolism; and (ii) enhance the experimental investigation of the metabolism of small xenobiotic molecules. The main predictive tools these methods employ are rule-based approaches, quantitative structure metabolism/activity relationships and docking approaches.

This review paper provides detailed metabolic profiles of the major clusters of antihypertensive agents, including their metabolites and their metabolizing enzymes, and it also provides specific information concerning the computational approaches that have been used to predict the metabolic profile of several antihypertensive drugs.

Modeling of interactions between xenobiotics and cytochrome P450 (CYP) enzymes

The adverse effects to humans and environment of only few chemicals are well known. Absorption, distribution, metabolism, and excretion (ADME) are the steps of pharmaco/toxicokinetics that determine the internal dose of chemicals to which the organism is exposed. Of all the xenobiotic-metabolizing enzymes, the cytochrome P450 (CYP) enzymes are the most important due to their abundance and versatility. Reactions catalyzed by CYPs usually turn xenobiotics to harmless and excretable metabolites, but sometimes an innocuous xenobiotic is transformed into a toxic metabolite. Data on ADME and toxicity properties of compounds are increasingly generated using in vitro and modeling (in silico) tools. Both physics-based and empirical modeling approaches are used. Numerous ligand-based and target-based as well as combined modeling methods have been employed to evaluate determinants of CYP ligand binding as well as predicting sites of metabolism and inhibition characteristics of test molecules. In silico prediction of CYP–ligand interactions have made crucial contributions in understanding (1) determinants of CYP ligand binding recognition and affinity; (2) prediction of likely metabolites from substrates; (3) prediction of inhibitors and their inhibition potency. Truly predictive models of toxic outcomes cannot be created without incorporating metabolic characteristics; in silico methods help producing such information and filling gaps in experimentally derived data. Currently modeling methods are not mature enough to replace standard in vitro and in vivo approaches, but they are already used as an important component in risk assessment of drugs and other chemicals.

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.

Selective and Dual Targeting of CCR2 and CCR5 Receptors: A Current Overview

The chemokine receptor 2 (CCR2) and chemokine receptor 5 (CCR5) are important mediators of leukocyte trafficking in inflammatory processes. The emerging evidence for a role of CCR2 and CCR5 receptors in human inflammatory diseases led to a growing interest in CCR2- and CCR5-selective antagonists. In this review, we focus on the recent development of selective CCR2/CCR5 receptor ligands and dual antagonists. Several compounds targeting CCR2, e.g., INCB8761 and MK0812, were developed as promising candidates for clinical trials, but failed to show clinical efficacy as presumed from preclinical models. The role of CCR5 receptors as the second co-receptor for the HIV-host cell fusion led to the development of various CCR5-selective ligands. Maraviroc is the first CCR5-targeting drug for the treatment of HIV-1 infections on the market. The role of CCR5 receptors in the progression of inflammatory processes fueled the use of CCR5 antagonists for the treatment of rheumatoid arthritis. Unfortunately, the use of maraviroc for the treatment of rheumatoid arthritis failed due to its inefficacy. Some of the ligands, e.g., TAK-779 and TAK-652, were also found to be dual antagonists of CCR2 and CCR5 receptors. The fact that CCR2 and CCR5 receptor antagonists contribute to the treatment of inflammatory diseases renders the development of dual antagonists as promising novel therapeutic strategy.

Metabolism site prediction based on xenobiotic structural formulas and PASS prediction algorithm

A new ligand-based method for the prediction of sites of metabolism (SOMs) for xenobiotics has been developed on the basis of the LMNA (labeled multilevel neighborhoods of atom) descriptors and the PASS (prediction of activity spectra for substances) algorithm and applied to predict the SOMs of the 1A2, 2C9, 2C19, 2D6, and 3A4 isoforms of cytochrome P450. An average IAP (invariant accuracy of prediction) of SOMs calculated by the leave-one-out cross-validation procedure was 0.89 for the developed method. The external validation was made with evaluation sets containing data on biotransformations for 57 cardiovascular drugs. An average IAP of regioselectivity for evaluation sets was 0.83. It was shown that the proposed method exceeds accuracy of SOM prediction by RS-Predictor for CYP 1A2, 2D6, 2C9, 2C19, and 3A4 and is comparable to or better than SMARTCyp for CYP 2C9 and 2D6.

Characterization of Cupriavidus metallidurans CYP116B1--a thiocarbamate herbicide oxygenating P450-phthalate dioxygenase reductase fusion protein

The novel cytochrome P450/redox partner fusion enzyme CYP116B1 from Cupriavidus metallidurans was expressed in and purified from Escherichia coli. Isolated CYP116B1 exhibited a characteristic Fe(II)CO complex with Soret maximum at 449 nm. EPR and resonance Raman analyses indicated low-spin, cysteinate-coordinated ferric haem iron at both 10 K and ambient temperature, respectively, for oxidized CYP116B1. The EPR of reduced CYP116B1 demonstrated stoichiometric binding of a 2Fe-2S cluster in the reductase domain. FMN binding in the reductase domain was confirmed by flavin fluorescence studies. Steady-state reduction of cytochrome c and ferricyanide were supported by both NADPH/NADH, with NADPH used more efficiently (K(m[NADPH]) = 0.9 ± 0.5 μM and K(m[NADH]) = 399.1 ± 52.1 μM). Stopped-flow studies of NAD(P)H-dependent electron transfer to the reductase confirmed the preference for NADPH. The reduction potential of the P450 haem iron was -301 ± 7 mV, with retention of haem thiolate ligation in the ferrous enzyme. Redox potentials for the 2Fe-2S and FMN cofactors were more positive than that of the haem iron. Multi-angle laser light scattering demonstrated CYP116B1 to be monomeric. Type I (substrate-like) binding of selected unsaturated fatty acids (myristoleic, palmitoleic and arachidonic acids) was shown, but these substrates were not oxidized by CYP116B1. However, CYP116B1 catalysed hydroxylation (on propyl chains) of the herbicides S-ethyl dipropylthiocarbamate (EPTC) and S-propyl dipropylthiocarbamate (vernolate), and the subsequent N-dealkylation of vernolate. CYP116B1 thus has similar thiocarbamate-oxidizing catalytic properties to Rhodoccocus erythropolis CYP116A1, a P450 involved in the oxidative degradation of EPTC.

Ligand-Based Site of Metabolism Prediction for Cytochrome P450 2D6

A ligand-based method based on the SMARTCyp approach that predicts the sites of cytochrome P450 2D6-mediated metabolism of druglike molecules has been developed. The method uses only two descriptors besides the reactivity from SMARTCyp: the distance to a protonated nitrogen atom and the distance to the end of the molecule. Hence, the site of metabolism is predicted directly from the 2D structure of a molecule, without requiring calculation of electronic properties or generation of 3D structures. Testing on an independent test set gives an area under the curve value of 0.94, and a site of metabolism is found among the top two ranked atoms for 91% of the compounds.

Selecting Relevant Descriptors for Classification by Bayesian Estimates: A Comparison with Decision Trees and Support Vector Machines Approaches for Disparate Data Sets

Classification algorithms suffer from the curse of dimensionality, which leads to overfitting, particularly if the problem is over-determined. Therefore it is of particular interest to identify the most relevant descriptors to reduce the complexity. We applied Bayesian estimates to model the probability distribution of descriptors values used for binary classification using n-fold cross-validation. As a measure for the discriminative power of the classifiers, the symmetric form of the Kullback-Leibler divergence of their probability distributions was computed. We found that the most relevant descriptors possess a Gaussian-like distribution of their values, show the largest divergences, and therefore appear most often in the cross-validation scenario. The results were compared to those of the LASSO feature selection method applied to multiple decision trees and support vector machine approaches for data sets of substrates and nonsubstrates of three Cytochrome P450 isoenzymes, which comprise strongly unbalanced compound distributions. In contrast to decision trees and support vector machines, the performance of Bayesian estimates is less affected by unbalanced data sets. This strategy reveals those descriptors that allow a simple linear separation of the classes, whereas the superior accuracy of decision trees and support vector machines can be attributed to nonlinear separation, which are in turn more prone to overfitting.

RS-predictor: a new tool for predicting sites of cytochrome P450-mediated metabolism applied to CYP 3A4

This article describes RegioSelectivity-Predictor (RS-Predictor), a new in silico method for generating predictive models of P450-mediated metabolism for drug-like compounds. Within this method, potential sites of metabolism (SOMs) are represented as "metabolophores": A concept that describes the hierarchical combination of topological and quantum chemical descriptors needed to represent the reactivity of potential metabolic reaction sites. RS-Predictor modeling involves the use of metabolophore descriptors together with multiple-instance ranking (MIRank) to generate an optimized descriptor weight vector that encodes regioselectivity trends across all cases in a training set. The resulting pathway-independent (O-dealkylation vs N-oxidation vs Csp(3) hydroxylation, etc.), isozyme-specific regioselectivity model may be used to predict potential metabolic liabilities. In the present work, cross-validated RS-Predictor models were generated for a set of 394 substrates of CYP 3A4 as a proof-of-principle for the method. Rank aggregation was then employed to merge independently generated predictions for each substrate into a single consensus prediction. The resulting consensus RS-Predictor models were shown to reliably identify at least one observed site of metabolism in the top two rank-positions on 78% of the substrates. Comparisons between RS-Predictor and previously described regioselectivity prediction methods reveal new insights into how in silico metabolite prediction methods should be compared.

Probing small-molecule binding to cytochrome P450 2D6 and 2C9: An in silico protocol for generating toxicity alerts

Drug metabolism, toxicity, and their interaction profiles are major issues in the drug-discovery and lead-optimization processes. The cytochromes P450 (CYPs) 2D6 and 2C9 are enzymes involved in the oxidative metabolism of a majority of marketed drugs. Therefore, the prediction of the binding affinity towards CYP2D6 and CYP2C9 would be beneficial for identifying cytochrome-mediated adverse effects triggered by drugs or chemicals (e.g., toxic reactions, drug-drug, and food-drug interactions). By identifying the binding mode by using pharmacophore prealignment, automated flexible docking, and by quantifying the binding affinity by multidimensional QSAR (mQSAR), we validated a model family of 56 compounds (46 training, 10 test) and 85 compounds (68 training, 17 test) for CYP2D6 and CYP2C9, respectively. The correlation with the experimental data (cross-validated r²=0.811 for CYP2D6 and 0.687 for CYP2C9) suggests that our approach is suited for predicting the binding affinity of compounds towards CYP2D6 and CYP2C9. The models were challenged by Y-scrambling and by testing an external dataset of binding compounds (15 compounds for CYP2D6 and 40 for CYP2C9). To assess the probability of false-positive predictions, datasets of nonbinders (64 compounds for CYP2D6 and 56 for CYP2C9) were tested by using the same protocol. The two validated mQSAR models were subsequently added to the VirtualToxLab (VTL, http://www.virtualtoxlab.org).

The challenges of in silico contributions to drug metabolism in lead optimization

IMPORTANCE OF THE FIELD

The site of metabolism (SOM) predictions by CYP 3A4 are extremely important during the drug discovery process especially during the lead discovery or library design phases. With the ability to rapidly characterize metabolites from these enzymes, the challenges facing in silico contribution change during the drug optimization phase. Some of the challenges are addressed in this article. Some aspects of the SOM prediction software and methodology are discussed in this opinion article and examples of software utility in overcoming metabolic instability in drug optimization are shown.

AREAS COVERED IN THIS REVIEW

SOM prediction by various approaches is discussed. Two ways of overcoming metabolic instability, blocking the metabolic softspots and rational modification of the instable molecule to avoid interaction with the CYP pocket, are discussed. The contribution plot in MetaSite and its use are discussed.

WHAT THE READER WILL GAIN

The reader will gain an understanding of possible approaches to either blocking the metabolic softspot or rationally modifying the molecule using MetaSite software or docking approaches. Blocking metabolism using fluorination has risks especially introducing multifluorinated benzene rings in the molecule.

TAKE HOME MESSAGE

During the lead optimization phase of drug discovery, when metabolic instability is an issue in a series, in silico approaches can be used to modify the molecule in order to decrease clearance due to metabolism, even that due to CYP3A4.

Using a homology model of cytochrome P450 2D6 to predict substrate site of metabolism

CYP2D6 is an important enzyme that is involved in first pass metabolism and is responsible for metabolizing ~25% of currently marketed drugs. A homology model of CYP2D6 was built using X-ray structures of ligand-bound CYP2C5 complexes as templates. This homology model was used in docking studies to rationalize and predict the site of metabolism of known CYP2D6 substrates. While the homology model was generally found to be in good agreement with the recently solved apo (ligand-free) X-ray structure of CYP2D6, significant differences between the structures were observed in the B' and F-G helical region. These structural differences are similar to those observed between ligand-free and ligand-bound structures of other CYPs and suggest that these conformational changes result from induced-fit adaptations upon ligand binding. By docking to the homology model using Glide, it was possible to identify the correct site of metabolism for a set of 16 CYP2D6 substrates 85% of the time when the 5 top scoring poses were examined. On the other hand, docking to the apo CYP2D6 X-ray structure led to a loss in accuracy in predicting the sites of metabolism for many of the CYP2D6 substrates considered in this study. These results demonstrate the importance of describing substrate-induced conformational changes that occur upon binding. The best results were obtained using Glide SP with van der Waals scaling set to 0.8 for both the receptor and ligand atoms. A discussion of putative binding modes that explain the distribution of metabolic sites for substrates, as well as a relationship between the number of metabolic sites and substrate size, are also presented. In addition, analysis of these binding modes enabled us to rationalize the typical hydroxylation and O-demethylation reactions catalyzed by CYP2D6 as well as the less common N-dealkylation.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

Understanding CYP2D6 interactions

Owing to the polymorphic nature of CYP2D6, clinically significant issues can arise when drugs rely on that enzyme either for clearance, or metabolism to an active metabolite. Available screening methods to determine if the compound is likely to cause drug-drug interactions, or is likely to be a victim of inhibition of CYP2D6 by other compounds will be described. Computational models and examples will be given on strategies to design out the CYP2D6 liabilities for both heme-binding compounds and non-heme-binding compounds.

QSAR models for P450 (2D6) substrate activity

Human Cytochrome P450 (CYP) is a large group of enzymes that possess an essential function in metabolising different exogenous and endogenous compounds. Humans have more than 50 different genes encoding CYP enzymes, among these a gene encoding for the CYP isoenzyme 2D6, a CYP able to metabolise drugs and other chemicals. A training set of 747 chemicals primarily based on in vivo human data for the CYP isoenzyme 2D6 was collected from the literature. QSAR models focusing on substrate/non-substrate activity were constructed by the use of MultiCASE, Leadscope and MDL quantitative structure-activity relationship (QSAR) modelling systems. They cross validated (leave-groups-out) with concordances of 71%, 81% and 82%, respectively. Discrete organic European Inventory of Existing Commercial Chemical Substances (EINECS) chemicals were screened to predict an approximate percentage of CYP 2D6 substrates. These chemicals are potentially present in the environment. The biological importance of the CYP 2D6 and the use of the software mentioned above were discussed.

AmineDB: Large scale docking of amines with CYP2D6 and scoring for druglike properties—towards defining the scope of the chemical defense against foreign amines in humans

Organic amines are prevalent in nature and in drugs, especially the psychotherapeutic agents, and a major defense against potentially toxic amines is metabolism by CYP2D6. In order to understand better the constraints on the broad specificity of CYP2D6, 4207 amines were docked into the binding site of this enzyme. Docking poses were found predominantly with the positively charged amino groups closest to Asp301, with aromatic rings close to Phe120 and sometimes extending as far as Phe483. Organic amines that bind best to CYP2D6 tend to have larger molecular weights and logP values. Organic amines that score highly as being druglike, based on a Bayesian model constructed using a 5223-drug training set, are least likely to bind to CYP2D6. This correlation suggests that the set of known drugs, which have been largely designed or selected to avoid high affinity CYP binding, partially encodes the binding site preferences (or rather anti-preferences) of CYP2D6. Finally, in order to benchmark our docking and druglike scoring procedures, an analysis of psychotherapeutic agents is presented. All of these data, including the 4207 AM1-optimized ligand structures in proper ionization states, docking poses and scores, Druglike Scores and Lipinski properties, can be viewed from an online database, the AmineDB.

Designing rapid onset selective serotonin re-uptake inhibitors. Part 3: Site-directed metabolism as a strategy to avoid active circulating metabolites: structure-activity relationships of (thioalkyl)phenoxy benzylamines

A series of thio-alkyl containing diphenylethers were designed and evaluated, as a strategy to competitively direct metabolism away from unwanted amine N-demethylation and deliver a pharmacologically inactive S-oxide metabolite. Overall, sulfonamide 20 was found to possess the best balance of target pharmacology, pharmacokinetics and metabolism profile.

Regioselectivity prediction of CYP1A2-mediated phase I metabolism

A kinetic, reactivity-binding model has been proposed to predict the regioselectivity of substrates meditated by the CYP1A2 enzyme, which is responsible for the metabolism of planar-conjugated compounds such as caffeine. This model consists of a docking simulation for binding energy and a semiempirical molecular orbital calculation for activation energy. Possible binding modes of CYP1A2 substrates were first examined using automated docking based on the crystal structure of CYP1A2, and binding energy was calculated. Then, activation energies for CYP1A2-mediated metabolism reactions were calculated using the semiempirical molecular orbital calculation, AM1. Finally, the metabolic probability obtained from two energy terms, binding and activation energies, was used for predicting the most probable metabolic site. This model predicted 8 out of 12 substrates accurately as the primary preferred site among all possible metabolic sites, and the other four substrates were predicted into the secondary preferred site. This method can be applied for qualitative prediction of drug metabolism mediated by CYP1A2 and other CYP450 family enzymes, helping to develop drugs efficiently.

Insights into drug metabolism by cytochromes P450 from modelling studies of CYP2D6-drug interactions

The cytochromes P450 (CYPs) comprise a vast superfamily of enzymes found in virtually all life forms. In mammals, xenobiotic metabolizing CYPs provide crucial protection from the effects of exposure to a wide variety of chemicals, including environmental toxins and therapeutic drugs. Ideally, the information on the possible metabolism by CYPs required during drug development would be obtained from crystal structures of all the CYPs of interest. For some years only crystal structures of distantly related bacterial CYPs were available and homology modelling techniques were used to bridge the gap and produce structural models of human CYPs, and thereby obtain useful functional information. A significant step forward in the reliability of these models came seven years ago with the first crystal structure of a mammalian CYP, rabbit CYP2C5, followed by the structures of six human enzymes, CYP1A2, CYP2A6, CYP2C8, CYP2C9, CYP2D6 and CYP3A4, and a second rabbit enzyme, CYP2B4. In this review we describe as a case study the evolution of a CYP2D6 model, leading to the validation of the model as an in silico tool for predicting binding and metabolism. This work has led directly to the successful design of CYP2D6 mutants with novel activity-including creating a testosterone hydroxylase, converting quinidine from inhibitor to substrate, creating a diclofenac hydroxylase and creating a dextromethorphan O-demethylase. Our modelling-derived hypothesis-driven integrated interdisciplinary studies have given key insight into the molecular determinants of CYP2D6 and other important drug metabolizing enzymes.

Biotransformation of metoprolol by the fungus Cunninghamella blakesleeana

AIM

To investigate the biotransformation of metoprolol, a beta1-cardioselective adrenoceptor antagonist, by filamentous fungus, and to compare the parallels between microbial transformation and mammalian metabolism.

METHODS

Five strains of Cunninghamella (C elegans AS 3.156, C elegans AS 3.2028, C echinulata AS 3.2004, C blakesleeana AS 3.153 and AS 3.910) were screened for the ability to transform metoprolol. The metabolites of metoprolol produced by C blakesleeana AS 3.153 were separated and assayed by liquid chromatography-tandem mass spectrometry (LC/MS(n)). The major metabolites were isolated by semipreparative HPLC and the structures were identified by a combination of LC/MS(n) and nuclear magnetic resonance analysis.

RESULTS

Metoprolol was transformed to 7 metabolites; 2 were identified as new metabolites and 5 were known metabolites in mammals.

CONCLUSION

The microbial transformation of metoprolol was similar to the metabolism in mammals. The fungi belonging to Cunninghamella species could be used as complementary models for predicting in vivo metabolism and producing quantities of metabolite references for drugs like metoprolol.

State-of-the-art tools for computational site of metabolism predictions: comparative analysis, mechanistical insights, and future applications

In drug design, it is crucial to have reliable information on how a chemical entity behaves in the presence of metabolizing enzymes. This requires substantial experimental efforts. Consequently, being able to predict the likely site/s of metabolism in any compound, synthesized or virtual, would be highly beneficial and time efficient. In this work, six different methodologies for predictions of the site of metabolism (SOM) have been compared and validated using structurally diverse data sets of drug-like molecules with well-established metabolic pattern in CYP3A4, CYP2C9, or both. Three of the methods predict the SOM based on the ligand's chemical structure, two additional methods use structural information of the enzymes, and the sixth method combines structure and ligand similarity and reactivity. The SOM is correctly predicted in 50 to 90% of the cases, depending on method and enzyme, which is an encouraging rate. We also discuss the underlying mechanisms of cytochrome P450 metabolism in the light of the results from this comparison.

Computational Approaches That Predict Metabolic Intermediate Complex Formation with CYP3A4 (+b5)

Some mechanism-based inhibitors cause irreversible inhibition by forming a metabolic intermediate complex (MIC) with cytochrome P450. In the present study, 54 molecules (substrates of CYP3A and amine-containing compounds that are not known substrates of CYP3A) were spectrophotometrically assessed for their propensity to cause MIC formation with recombinant CYP3A4 (+b(5)). Comparisons of common physicochemical properties showed that mean (+/-S.D.) mol. wt. of MIC-forming compounds was significantly greater than mean mol. wt. of non-MIC-forming compounds, 472 (+/-173) versus 307 (+/-137), respectively. Computational pharmacophores, logistic regression, and recursive partitioning (RP) approaches were applied to predict MIC formation from molecular structure and to generate a quantitative structure activity relationship. A pharmacophore built with SKF-525A (2-diethylaminoethyl 2:2-diphenylvalerate hydrochloride), erythromycin, amprenavir, and norverapamil indicated that four hydrophobic features and a hydrogen bond acceptor were important for these MIC-forming compounds. Two different RP methods using either simple descriptors or 2D augmented atom descriptors indicated that hydro-phobic and hydrogen bond acceptor features were required for MIC formation. Both of these RP methods correctly predicted the MIC formation status with CYP3A4 for 10 of 12 literature molecules in an independent test set. Logistic multiple regression and a third classification tree model predicted 11 of 12 molecules correctly. Both models possessed a hydrogen bond acceptor and represent an approach for predicting CYP3A4 MIC formation that can be improved using more data and molecular descriptors. The preliminary pharmacophores provide structural insights that complement those for CYP3A4 inhibitors and substrates.

Predicting drug metabolism--an evaluation of the expert system METEOR

The paper begins with a discussion of the goals of metabolic predictions in early drug research, and some difficulties toward this objective, mainly the various substrate and product selectivities characteristic of drug metabolism. The major in silico approaches to predict drug metabolism are then classified and summarized. A discrimination is, thus, made between 'local' and 'global' systems. In its second part, an evaluation of METEOR, a rule-based expert system used to predict the metabolism of drugs and other xenobiotics, is reported. The published metabolic data of ten substrates were used in this evaluation, the overall results being discussed in terms of correct vs. disputable (i.e., false-positive and false-negative) predictions. The predictions for four representative substrates are presented in detail (Figs. 1-4), illustrating the interest of such an evaluation in identifying where and how predictive rules can be improved.

Methods for Predicting Human Drug Metabolism

Drug metabolism information is a necessary component of drug discovery and development. The key issues in drug metabolism include identifying: the enzyme(s) involved, the site(s) of metabolism, the resulting metabolite(s), and the rate of metabolism. Methods for predicting human drug metabolism from in vitro and computational methodologies and determining relationships between the structure and metabolic activity of molecules are also critically important for understanding potential drug interactions and toxicity. There are numerous experimental and computational approaches that have been developed in order to predict human metabolism which have their own limitations. It is apparent that few of the computational tools for metabolism prediction alone provide the major integrated functions needed to assist in drug discovery. Similarly the different in vitro methods for human drug metabolism themselves have implicit limitations. The utilization of these methods for pharmaceutical and other applications as well as their integration is discussed as it is likely that hybrid methods will provide the most success.

Pharmacophore-based discovery of ligands for drug transporters

The ability to identify ligands for drug transporters is an important step in drug discovery and development. It can both improve accurate profiling of lead pharmacokinetic properties and assist in the discovery of new chemical entities targeting transporters. In silico approaches, especially pharmacophore-based database screening methods have great potential in improving the throughput of current transporter ligand identification assays, leading to a higher hit rate by focusing in vitro testing to the most promising hits. In this review, the potential of different in silico methods in transporter ligand identification studies are compared and summarized with an emphasis on pharmacophore modeling. Various implementations of pharmacophore model generation, database compilation and flexible screening algorithms are also introduced. Recent successful utilization of database searching with pharmacophores to identify novel ligands for the pharmaceutically significant transporters hPepT1, P-gp, BCRP, MRP1 and DAT are reviewed and the challenges encountered with current approaches are discussed.

Designing better drugs: predicting cytochrome P450 metabolism

Many 3D ligand-based and structure-based computational approaches have been used to predict, and thus help explain, the metabolism catalyzed by the enzymes of the cytochrome P450 superfamily (P450s). P450s are responsible for >90% of the metabolism of all drugs, so the computational prediction of metabolism can help to design out drug-drug interactions in the early phases of the drug discovery process. Computational methodologies have focused on a few P450s that are directly involved in drug metabolism. The recently derived crystal structures for human P450s enable better 3D modelling of these important metabolizing enzymes. Models derived for P450s have evolved from simple comparisons of known substrates to more-elaborate experiments that require considerable computer power involving 3D overlaps and docking experiments. These models help to explain and, more importantly, predict the involvement of P450s in the metabolism of specific compounds and guide the drug-design process.

Computational approaches for predicting CYP-related metabolism properties in the screening of new drugs

The site of biotransformation, the extent and rate of metabolism and the number of active metabolic pathways are among the most important characteristics of the pharmacokinetics of a drug. The catalytic activity of drug metabolizing enzymes is likely the most influential determinant of the pharmacokinetic variability. Metabolic stability is the prerequisite for sustaining the therapeutically relevant concentrations. Metabolic inhibition and induction can give rise to clinically important drug-drug interactions. A variety of computational approaches are currently available for predicting different cytochrome P450 (CYP)-related metabolism endpoints. The present review will describe these approaches and their impact on drug development process. Indications on the available software for the implementation will also be given.

Functional interaction of nitrogenous organic bases with cytochrome P450: A critical assessment and update of substrate features and predicted key active-site elements steering the access, binding, and orientation of amines

The widespread use of nitrogenous organic bases as environmental chemicals, food additives, and clinically important drugs necessitates precise knowledge about the molecular principles governing biotransformation of this category of substrates. In this regard, analysis of the topological background of complex formation between amines and P450s, acting as major catalysts in C- and N-oxidative attack, is of paramount importance. Thus, progress in collaborative investigations, combining physico-chemical techniques with chemical-modification as well as genetic engineering experiments, enables substantiation of hypothetical work resulting from the design of pharmacophores or homology modelling of P450s. Based on a general, CYP2D6-related construct, the majority of prospective amine-docking residues was found to cluster near the distal heme face in the six known SRSs, made up by the highly variant helices B', F and G as well as the N-terminal portion of helix C and certain beta-structures. Most of the contact sites examined show a frequency of conservation < 20%, hinting at the requirement of some degree of conformational versatility, while a limited number of amino acids exhibiting a higher level of conservation reside close to the heme core. Some key determinants may have a dual role in amine binding and/or maintenance of protein integrity. Importantly, a series of non-SRS elements are likely to be operative via long-range effects. While hydrophobic mechanisms appear to dominate orientation of the nitrogenous compounds toward the iron-oxene species, polar residues seem to foster binding events through H-bonding or salt-bridge formation. Careful uncovering of structure-function relationships in amine-enzyme association together with recently developed unsupervised machine learning approaches will be helpful in both tailoring of novel amine-type drugs and early elimination of potentially toxic or mutagenic candidates. Also, chimeragenesis might serve in the construction of more efficient P450s for activation of amine drugs and/or bioremediation.