QUANTITATIVE STRUCTURE-METABOLISM RELATIONSHIP MODELING OF METABOLICN-DEALKYLATION REACTION RATES

In-Silico Modeling in Drug Metabolism and Interaction: Current Strategies of Lead Discovery

BACKGROUND

Drug metabolism is a complex mechanism of human body systems to detoxify foreign particles, chemicals, and drugs through bio alterations. It involves many biochemical reactions carried out by invivo enzyme systems present in the liver, kidney, intestine, lungs, and plasma. After drug administration, it crosses several biological membranes to reach into the target site for binding and produces the therapeutic response. After that, it may undergo detoxification and excretion to get rid of the biological systems. Most of the drugs and its metabolites are excreted through kidney via urination. Some drugs and their metabolites enter into intestinal mucosa and excrete through feces. Few of the drugs enter into hepatic circulation where they go into the intestinal tract. The drug leaves the liver via the bile duct and is excreted through feces. Therefore, the study of total methodology of drug biotransformation and interactions with various targets is costly.

METHODS

To minimize time and cost, in-silico algorithms have been utilized for lead-like drug discovery. Insilico modeling is the process where a computer model with a suitable algorithm is developed to perform a controlled experiment. It involves the combination of both in-vivo and in-vitro experimentation with virtual trials, eliminating the non-significant variables from a large number of variable parameters. Whereas, the major challenge for the experimenter is the selection and validation of the preferred model, as well as precise simulation in real physiological status.

RESULTS

The present review discussed the application of in-silico models to predict absorption, distribution, metabolism, and excretion (ADME) properties of drug molecules and also access the net rate of metabolism of a compound.

CONCLUSION

It helps with the identification of enzyme isoforms; which are likely to metabolize a compound, as well as the concentration dependence of metabolism and the identification of expected metabolites. In terms of drug-drug interactions (DDIs), models have been described for the inhibition of metabolism of one compound by another, and for the compound-dependent induction of drug-metabolizing enzymes.

Computationally Assessing the Bioactivation of Drugs by N-Dealkylation

Cytochromes P450 (CYPs) oxidize alkylated amines commonly found in drugs and other biologically active molecules, cleaving them into an amine and an aldehyde. Metabolic studies usually neglect to report or investigate aldehydes, even though they can be toxic. It is assumed that they are efficiently detoxified into carboxylic acids and alcohols. Nevertheless, some aldehydes are reactive and escape detoxification pathways to cause adverse events by forming DNA and protein adducts. Herein, we modeled N-dealkylations that produce both amine and aldehyde metabolites and then predicted the reactivity of the aldehyde. This model used a deep learning approach previously developed by our group to predict other types of drug metabolism. In this study, we trained the model to predict N-dealkylation by human liver microsomes (HLM), finding that including isozyme-specific metabolism data alongside HLM data significantly improved results. The final HLM model accurately predicted the site of N-dealkylation within metabolized substrates (97% top-two and 94% area under the ROC curve). Next, we combined the metabolism, metabolite structure prediction, and previously published reactivity models into a bioactivation model. This combined model predicted the structure of the most likely reactive metabolite of a small validation set of drug-like molecules known to be bioactivated by N-dealkylation. Applying this model to approved and withdrawn medicines, we found that aldehyde metabolites produced from N-dealkylation may explain the hepatotoxicity of several drugs: indinavir, piperacillin, verapamil, and ziprasidone. Our results suggest that N-dealkylation may be an under-appreciated bioactivation pathway, especially in clinical contexts where aldehyde detoxification pathways are inhibited. Moreover, this is the first report of a bioactivation model constructed by combining a metabolism and reactivity model. These results raise hope that more comprehensive models of bioactivation are possible. The model developed in this study is available at http://swami.wustl.edu/xenosite/ .

Modeling Reactivity to Biological Macromolecules with a Deep Multitask Network

Most small-molecule drug candidates fail before entering the market, frequently because of unexpected toxicity. Often, toxicity is detected only late in drug development, because many types of toxicities, especially idiosyncratic adverse drug reactions (IADRs), are particularly hard to predict and detect. Moreover, drug-induced liver injury (DILI) is the most frequent reason drugs are withdrawn from the market and causes 50% of acute liver failure cases in the United States. A common mechanism often underlies many types of drug toxicities, including both DILI and IADRs. Drugs are bioactivated by drug-metabolizing enzymes into reactive metabolites, which then conjugate to sites in proteins or DNA to form adducts. DNA adducts are often mutagenic and may alter the reading and copying of genes and their regulatory elements, causing gene dysregulation and even triggering cancer. Similarly, protein adducts can disrupt their normal biological functions and induce harmful immune responses. Unfortunately, reactive metabolites are not reliably detected by experiments, and it is also expensive to test drug candidates for potential to form DNA or protein adducts during the early stages of drug development. In contrast, computational methods have the potential to quickly screen for covalent binding potential, thereby flagging problematic molecules and reducing the total number of necessary experiments. Here, we train a deep convolution neural network-the XenoSite reactivity model-using literature data to accurately predict both sites and probability of reactivity for molecules with glutathione, cyanide, protein, and DNA. On the site level, cross-validated predictions had area under the curve (AUC) performances of 89.8% for DNA and 94.4% for protein. Furthermore, the model separated molecules electrophilically reactive with DNA and protein from nonreactive molecules with cross-validated AUC performances of 78.7% and 79.8%, respectively. On both the site- and molecule-level, the model's performances significantly outperformed reactivity indices derived from quantum simulations that are reported in the literature. Moreover, we developed and applied a selectivity score to assess preferential reactions with the macromolecules as opposed to the common screening traps. For the entire data set of 2803 molecules, this approach yielded totals of 257 (9.2%) and 227 (8.1%) molecules predicted to be reactive only with DNA and protein, respectively, and hence those that would be missed by standard reactivity screening experiments. Site of reactivity data is an underutilized resource that can be used to not only predict if molecules are reactive, but also show where they might be modified to reduce toxicity while retaining efficacy. The XenoSite reactivity model is available at http://swami.wustl.edu/xenosite/p/reactivity.

Tales from the war on error: the art and science of curating QSAR data

Curating the data underlying quantitative structure-activity relationship models is a never-ending struggle. Some curation can now be automated but much cannot, especially where data as complex as those pertaining to molecular absorption, distribution, metabolism, excretion, and toxicity are concerned (vide infra). The authors discuss some particularly challenging problem areas in terms of specific examples involving experimental context, incompleteness of data, confusion of units, problematic nomenclature, tautomerism, and misapplication of automated structure recognition tools.

Open Source Bayesian Models. 1. Application to ADME/Tox and Drug Discovery Datasets

On the order of hundreds of absorption, distribution, metabolism, excretion, and toxicity (ADME/Tox) models have been described in the literature in the past decade which are more often than not inaccessible to anyone but their authors. Public accessibility is also an issue with computational models for bioactivity, and the ability to share such models still remains a major challenge limiting drug discovery. We describe the creation of a reference implementation of a Bayesian model-building software module, which we have released as an open source component that is now included in the Chemistry Development Kit (CDK) project, as well as implemented in the CDD Vault and in several mobile apps. We use this implementation to build an array of Bayesian models for ADME/Tox, in vitro and in vivo bioactivity, and other physicochemical properties. We show that these models possess cross-validation receiver operator curve values comparable to those generated previously in prior publications using alternative tools. We have now described how the implementation of Bayesian models with FCFP6 descriptors generated in the CDD Vault enables the rapid production of robust machine learning models from public data or the user’s own datasets. The current study sets the stage for generating models in proprietary software (such as CDD) and exporting these models in a format that could be run in open source software using CDK components. This work also demonstrates that we can enable biocomputation across distributed private or public datasets to enhance drug discovery.

Predicting CYP2C19 catalytic parameters for enantioselective oxidations using artificial neural networks and a chirality code

Cytochromes P450 (CYP for isoforms) play a central role in biological processes especially metabolism of chiral molecules; thus, development of computational methods to predict parameters for chiral reactions is important for advancing this field. In this study, we identified the most optimal artificial neural networks using conformation-independent chirality codes to predict CYP2C19 catalytic parameters for enantioselective reactions. Optimization of the neural networks required identifying the most suitable representation of structure among a diverse array of training substrates, normalizing distribution of the corresponding catalytic parameters (k(cat), K(m), and k(cat)/K(m)), and determining the best topology for networks to make predictions. Among different structural descriptors, the use of partial atomic charges according to the CHelpG scheme and inclusion of hydrogens yielded the most optimal artificial neural networks. Their training also required resolution of poorly distributed output catalytic parameters using a Box-Cox transformation. End point leave-one-out cross correlations of the best neural networks revealed that predictions for individual catalytic parameters (k(cat) and K(m)) were more consistent with experimental values than those for catalytic efficiency (k(cat)/K(m)). Lastly, neural networks predicted correctly enantioselectivity and comparable catalytic parameters measured in this study for previously uncharacterized CYP2C19 substrates, R- and S-propranolol. Taken together, these seminal computational studies for CYP2C19 are the first to predict all catalytic parameters for enantioselective reactions using artificial neural networks and thus provide a foundation for expanding the prediction of cytochrome P450 reactions to chiral drugs, pollutants, and other biologically active compounds.

Quantitative Property-Property Relationship for Screening-Level Prediction of Intrinsic Clearance of Volatile Organic Chemicals in Rats and Its Integration within PBPK Models to Predict Inhalation Pharmacokinetics in Humans

The objectives of this study were (i) to develop a screening-level Quantitative property-property relationship (QPPR) for intrinsic clearance (CL_int) obtained from in vivo animal studies and (ii) to incorporate it with human physiology in a PBPK model for predicting the inhalation pharmacokinetics of VOCs. CL_int, calculated as the ratio of the in vivo V_max (μmol/h/kg bw rat) to the K_m (μM), was obtained for 26 VOCs from the literature. The QPPR model resulting from stepwise linear regression analysis passed the validation step (R² = 0.8; leave-one-out cross-validation Q² = 0.75) for CL_int normalized to the phospholipid (PL) affinity of the VOCs. The QPPR facilitated the calculation of CL_int (L PL/h/kg bw rat) from the input data on log P_ow, log blood: water PC and ionization potential. The predictions of the QPPR as lower and upper bounds of the 95% mean confidence intervals (LMCI and UMCI, resp.) were then integrated within a human PBPK model. The ratio of the maximum (using LMCI for CL_int) to minimum (using UMCI for CL_int) AUC predicted by the QPPR-PBPK model was 1.36 ± 0.4 and ranged from 1.06 (1,1-dichloroethylene) to 2.8 (isoprene). Overall, the integrated QPPR-PBPK modeling method developed in this study is a pragmatic way of characterizing the impact of the lack of knowledge of CL_int in predicting human pharmacokinetics of VOCs, as well as the impact of prediction uncertainty of CL_int on human pharmacokinetics of VOCs.

Prediction of human pharmacokinetics—evaluation of methods for prediction of hepatic metabolic clearance

Methods for prediction of hepatic clearance (CLH) in man have been evaluated. A physiologically-based in-vitro to in-vivo (PB-IVIV) method with human unbound fraction in blood (fu,bl) and hepatocyte intrinsic clearance (CLint)-data has a good rationale and appears to give the best predictions (maximum ∼2-fold errors; < 25% errors for half of CL-predictions; appropriate ranking). Inclusion of an empirical scaling factor is, however, needed, and reasons include the use of cryopreserved hepatocytes with low activity, and inappropriate CLint- and fu,bl-estimation methods. Thus, an improvement of this methodology is possible and required. Neglect of fu,bl or incorporation of incubation binding does not seem appropriate. When microsome CLint-data are used with this approach, the CLH is underpredicted by 5- to 9-fold on average, and a 106-fold underprediction (attrition potential) has been observed. The poor performance could probably be related to permeation, binding and low metabolic activity. Inclusion of scaling factors and neglect of fu,bl for basic and neutral compounds improve microsome predictions. The performance is, however, still not satisfactory. Allometry incorrectly assumes that the determinants for CLH relate to body weight and overpredicts human liver blood flow rate. Consequently, allometric methods have poor predictability. Simple allometry has an average overprediction potential, > 2-fold errors for ∼1/3 of predictions, and 140-fold underprediction to 5800-fold overprediction (potential safety risk) range. In-silico methodologies are available, but these need further development. Acceptable prediction errors for compounds with low and high CLH should be ∼50 and ∼10%, respectively. In conclusion, it is recommended that PB-IVIV with human hepatocyte CLint and fu,bl is applied and improved, limits for acceptable errors are decreased, and that animal CLH-studies and allometry are avoided.

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.

SPORCalc: A development of a database analysis that provides putative metabolic enzyme reactions for ligand-based drug design

Understanding both the enzyme reactions that contribute to intermediate metabolism and the biochemical fate of candidate therapeutic and toxic agents are essential for drug design. Traditional metabolic databases indicate whether reactions have been observed but do not provide the likelihoods of reactions occurring, for example those of mixed function oxygenases and oxidases, during phase I metabolism. The desire for more quantitative predictions motivated the development of the recently introduced Substrate Product Occurrence Ratio Calculator (SPORCalc) that identifies metabolically labile atom positions in candidate compounds. This paper describes a further development and provides a clearer explanation of SPORCalc for the computational pharmacology, medicinal chemistry and drug design communities interested in metabolic prediction of xenobiotics using chemical databases of biotransformations. Examples of reaction centre detection in Metabolite are described followed by a demonstration of almokalant, an anti-arrhythmic agent, undergoing phase I metabolism. In general, occurrence ratio (OR) values are calculated throughout a compound and its transformed metabolites to give propensity (p) values at each atom position. The OR values from substrates and products in the database are essential for addition and elimination reactions. For almokalant, the resulting p values ranged from 10(-1) to 10(-5) and their order of magnitude reflected the known and experimentally observed metabolites. SPORCalc depends entirely on the level of detail from isoform- or species-specific reaction classes in Metabolite. Labile atom positions (sites of metabolism) are identified in both the candidate compound and its metabolites. In general, the likelihood of one enzyme isoform-dependent reaction occurring relative to another and the putative metabolic routes from different isoforms can be investigated. SPORCalc can be developed further to include suitable three-dimensional, structure-activity and physiochemical information.

Molecular characterization of CYP2B6 substrates

CYP2B6 has not been as fully characterized at the molecular level as other members of the human cytochrome P450 family. As more widely used in vitro probes for characterizing the involvement of this enzyme in the metabolism of xenobiotics have become available, the number of molecules identified as CYP2B6 substrates has increased. In this study we have analyzed the available kinetic data generated by multiple laboratories with human recombinant expressed CYP2B6 and along with calculated molecular properties derived from the ChemSpider database, we have determined the molecular features that appear to be important for CYP2B6 substrates. In addition we have applied 2D and 3D QSAR methods to generate predictive pharmacophore and 2D models. For 28 molecules with K(m) data, the molecular weight (mean +/- SD) is 253.78+/-74.03, ACD/logP is 2.68+/-1.51, LogD(pH 5.5) is 1.51+/-1.43, LogD(pH 7.4) is 2.02+/-1.25, hydrogen bond donor (HBD) count is 0.57 +/-0.57, hydrogen bond acceptor (HBA) count is 2.57+/-1.37, rotatable bonds is 3.50+/-2.71 and total polar surface area (TPSA) is 27.63+/-19.42. A second set of 15 molecules without K(m) data possessed similar mean molecular property values. These properties are comparable to those of a set of 21 molecules used in a previous pharmacophore modeling study (Ekins et al., J Pharmacol Exp Ther 288 (1), 21-29, 1999). Only the LogD and HBD values were statistically significantly different between these different datasets. We have shown that CYP2B6 substrates are generally small hydrophobic molecules that are frequently central nervous system active, which may be important for drug discovery research.

In silico pharmacology for drug discovery: applications to targets and beyond

Computational (in silico) methods have been developed and widely applied to pharmacology hypothesis development and testing. These in silico methods include databases, quantitative structure-activity relationships, similarity searching, pharmacophores, homology models and other molecular modeling, machine learning, data mining, network analysis tools and data analysis tools that use a computer. Such methods have seen frequent use in the discovery and optimization of novel molecules with affinity to a target, the clarification of absorption, distribution, metabolism, excretion and toxicity properties as well as physicochemical characterization. The first part of this review discussed the methods that have been used for virtual ligand and target-based screening and profiling to predict biological activity. The aim of this second part of the review is to illustrate some of the varied applications of in silico methods for pharmacology in terms of the targets addressed. We will also discuss some of the advantages and disadvantages of in silico methods with respect to in vitro and in vivo methods for pharmacology research. Our conclusion is that the in silico pharmacology paradigm is ongoing and presents a rich array of opportunities that will assist in expediating the discovery of new targets, and ultimately lead to compounds with predicted biological activity for these novel targets.

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.

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

Inhibition of cytochrome P450 (CYP) enzymes is unwanted because of the risk of severe side effects due to drug-drug interactions. We present two in silico Gaussian kernel weighted k-nearest neighbor models based on extended connectivity fingerprints that classify CYP2D6 and CYP3A4 inhibition. Data used for modeling consisted of diverse sets of 1153 and 1382 drug candidates tested for CYP2D6 and CYP3A4 inhibition in human liver microsomes. For CYP2D6, 82% of the classified test set compounds were predicted to the correct class. For CYP3A4, 88% of the classified compounds were correctly classified. CYP2D6 and CYP3A4 inhibition were additionally classified for an external test set on 14 drugs, and multidimensional scaling plots showed that the drugs in the external test set were in the periphery of the training sets. Furthermore, fragment analyses were performed and structural fragments frequent in CYP2D6 and CYP3A4 inhibitors and noninhibitors are presented.

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.

Screening for 2-quinolinone-derived selective androgen receptor agonists in doping control analysis

Selective androgen receptor modulators (SARMs) represent a class of emerging drugs with high potential for misuse in sports, and therefore members of this group are banned as anabolic agents by the World Anti-Doping Agency. Preventive approaches to restrict their use include early implementation of target analytes into doping control screening assays and evaluation of the mass spectrometric behavior of these drugs to allow their unequivocal identification as well as the characterization of structurally related compounds and metabolic products. Four model SARMs with the 6-alkylamino-2-quinolinone structure, including the advanced drug candidate LGD-2226, were synthesized. Fragmentation pathways after positive electrospray ionization and collision-induced dissociation were studied using an LTQ Orbitrap mass analyzer, and diagnostic product ions and common dissociation pathways were employed to establish a screening procedure targeting intact quinolinone-based SARMs as well as putative metabolic products such as dealkylated analogues. Therefore, features of a triple quadrupole mass analyzer such as multiple reaction monitoring and precursor ion scanning were utilized. Sample preparation based on commonly employed liquid-liquid extraction and subsequent liquid chromatographic/tandem mass spectrometric measurement allowed for detection limits of 0.01-0.2 ng/mL, and intra- and interday precisions between 3.2 and 8.5% and between 6.3 and 16.6%, respectively. Recoveries varied from 81 to 98%, and tests for ion suppression or enhancement effects were negative for all analytes.

Quantitative structure-pharmacokinetic/pharmacodynamic relationships

Quantitative structure-activity relationships have long been considered a vital component of drug discovery and development, providing insight into the role of molecular properties in the biological activity of similar and unrelated compounds. Recognition that in vitro bioassay and/or pre-clinical activity are insufficient for anticipating which compounds are suitable leads for further development has shifted the focus toward integrated pharmacokinetic (PK) and pharmacodynamic (PD) processes. Over the last decade, considerable progress has been made in constructing empirical and mechanistic quantitative structure-PK relationships (QSPKR), as well as diverse mechanism-based pharmacodynamic models of drug effects. In this review, traditional and contemporary approaches to developing QSPKR models are discussed, along with selected examples of attempts to couple QSPKR and pharmacodynamic models to anticipate the intensity and time-course of the pharmacological effects of new or related compounds, or quantitative structure-pharmacodynamic relationships modeling. Such models are in accordance with the goals of systems biology and the ideal of designing drugs and delivery systems from first principles.

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.

A COMBINED APPROACH TO DRUG METABOLISM AND TOXICITY ASSESSMENT

The challenge of predicting the metabolism or toxicity of a drug in humans has been approached using in vivo animal models, in vitro systems, high throughput genomics and proteomics methods, and, more recently, computational approaches. Understanding the complexity of biological systems requires a broader perspective rather than focusing on just one method in isolation for prediction. Multiple methods may therefore be necessary and combined for a more accurate prediction. In the field of drug metabolism and toxicology, we have seen the growth, in recent years, of computational quantitative structure-activity relationships (QSARs), as well as empirical data from microarrays. In the current study we have further developed a novel computational approach, MetaDrug, that 1) predicts metabolites for molecules based on their chemical structure, 2) predicts the activity of the original compound and its metabolites with various absorption, distribution, metabolism, excretion, and toxicity models, 3) incorporates the predictions with human cell signaling and metabolic pathways and networks, and 4) integrates networks and metabolites, with relevant toxicogenomic or other high throughput data. We have demonstrated the utility of such an approach using recently published data from in vitro metabolism and microarray studies for aprepitant, 2(S)-((3,5-bis(trifluoromethyl)benzyl)-oxy)-3(S)phenyl-4-((3-oxo-1,2,4-triazol-5-yl)methyl)morpholine (L-742694), trovofloxacin, 4-hydroxytamoxifen, and artemisinin and other artemisinin analogs to show the predicted interactions with cytochromes P450, pregnane X receptor, and P-glycoprotein, and the metabolites and the networks of genes that are affected. As a comparison, we used a second computational approach, MetaCore, to generate statistically significant gene networks with the available expression data. These case studies demonstrate the combination of QSARs and systems biology methods.

Computational prediction of human drug metabolism

There is an urgent requirement within the pharmaceutical and biotechnology industries, regulatory authorities and academia to improve the success of molecules that are selected for clinical trials. Although absorption, distribution, metabolism, excretion and toxicity (ADME/Tox) properties are some of the many components that contribute to successful drug discovery and development, they represent factors for which we currently have in vitro and in vivo data that can be modelled computationally. Understanding the possible toxicity and the metabolic fate of xenobiotics in the human body is particularly important in early drug discovery. There is, therefore, a need for computational methodologies for uncovering the relationships between the structure and the biological activity of novel molecules. The convergence of numerous technologies, including high-throughput techniques, databases, ADME/Tox modelling and systems biology modelling, is leading to the foundation of systems-ADME/Tox. Results from experiments can be integrated with predictions to globally simulate and understand the likely complete effects of a molecule in humans. The development and early application of major components of MetaDrug (GeneGo, Inc.) software will be described, which includes rule-based metabolite prediction, quantitative structure-activity relationship models for major drug metabolising enzymes, and an extensive database of human protein-xenobiotic interactions. This represents a combined approach to predicting drug metabolism. MetaDrug can be readily used for visualising Phase I and II metabolic pathways, as well as interpreting high-throughput data derived from microarrays as networks of interacting objects. This will ultimately aid in hypothesis generation and the early triaging of molecules likely to have undesirable predicted properties or measured effects on key proteins and cellular functions.

Systems-ADME/Tox: resources and network approaches

The increasing cost of drug development is partially due to our failure to identify undesirable compounds at an early enough stage of development. The application of higher throughput screening methods have resulted in the generation of very large datasets from cells in vitro or from in vivo experiments following the treatment with drugs or known toxins. In recent years the development of systems biology, databases and pathway software has enabled the analysis of the high-throughput data in the context of the whole cell. One of the latest technology paradigms to be applied alongside the existing in vitro and computational models for absorption, distribution, metabolism, excretion and toxicology (ADME/Tox) involves the integration of complex multidimensional datasets, termed toxicogenomics. The goal is to provide a more complete understanding of the effects a molecule might have on the entire biological system. However, due to the sheer complexity of this data it may be necessary to apply one or more different types of computational approaches that have as yet not been fully utilized in this field. The present review describes the data generated currently and introduces computational approaches as a component of ADME/Tox. These methods include network algorithms and manually curated databases of interactions that have been separately classified under systems biology methods. The integration of these disparate tools will result in systems-ADME/Tox and it is important to understand exactly what data resources and technologies are available and applicable. Examples of networks derived with important drug transporters and drug metabolizing enzymes are provided to demonstrate the network technologies.