Development of an information-intensive structure-activity relationship model and its application to human respiratory chemical sensitizers

Structure-activity relationship (SAR) models are recognized as powerful tools to predict the toxicologic potential of new or untested chemicals and also provide insight into possible mechanisms of toxicity. Models have been based on physicochemical attributes and structural features of chemicals. We describe herein the development of a new SAR modeling algorithm called cat-SAR that is capable of analyzing and predicting chemical activity from divergent biological response data. The cat-SAR program develops chemical fragment-based SAR models from categorical biological response data (e.g. toxicologically active and inactive compounds). The database selected for model development was a published set of chemicals documented to cause respiratory hypersensitivity in humans. Two models were generated that differed only in that one model included explicate hydrogen containing fragments. The predictive abilities of the models were tested using leave-one-out cross-validation tests. One model had a sensitivity of 0.94 and specificity of 0.87 yielding an overall correct prediction of 91%. The second model had a sensitivity of 0.89, specificity of 0.95 and overall correct prediction of 92%. The demonstrated predictive capabilities of the cat-SAR approach, together with its modeling flexibility and design transparency, suggest the potential for its widespread applicability to toxicity prediction and for deriving mechanistic insight into toxicologic effects.     

Comparison of metabolic soft spot predictions of CYP3A4, CYP2C9 and CYP2D6 substrates using MetaSite and StarDrop

Metabolite identification study plays an important role in determining the sites of metabolic liability of new chemical entities (NCEs) in drug discovery for lead optimization. Here we compare the two predictive software, MetaSite and StarDrop, available for this purpose. They work very differently but are used to predict the site of oxidation by major human cytochrome P450 (CYP) isoforms. Neither software can predict non-CYP catalyzed metabolism nor the rates of metabolism. For the purpose of comparing the two software packages, we tested known probe substrate for these enzymes, which included 12 substrates of CYP3A4 and 18 substrates of CYP2C9 and CYP2D6 were analyzed by each software and the results were compared. It is possible that these known substrates were part of the training set but we are not aware of it. To assess the performance of each software we assigned a point system for each correct prediction. The total points assigned for each CYP isoform experimentally were compared as a percentage of the total points assigned theoretically for the first choice prediction for all substrates for each isoform. Our results show that MetaSite and StarDrop are similar in predicting the correct site of metabolism by CYP3A4 (78% vs 83%, respectively). StarDrop appears to do slightly better in predicting the correct site of metabolism by CYP2C9 and CYP2D6 metabolism (89% and 93%, respectively) compared to MetaSite (63% and 70%, respectively). The sites of metabolism (SOM) from 34 in-house NCEs incubated in human liver microsomes or human hepatocytes were also evaluated using two prediction software packages and the results showed comparable SOM predictions. What makes this comparison challenging is that the contribution of each isoform to the intrinsic clearance (Clint) is not known. Overall the software were comparable except for MetaSite performing better for CYP2D6 and that MetaSite has a liver model that is absent in StarDrop that predicted with 82% accuracy.     

Development of an information-intensive structure–activity relationship model and its application to human respiratory chemical sensitizers

Structure–activity relationship (SAR) models are recognized as powerful tools to predict the toxicologic potential of new or untested chemicals and also provide insight into possible mechanisms of toxicity. Models have been based on physicochemical attributes and structural features of chemicals. We describe herein the development of a new SAR modeling algorithm called cat-SAR that is capable of analyzing and predicting chemical activity from divergent biological response data. The cat-SAR program develops chemical fragment-based SAR models from categorical biological response data (e.g. toxicologically active and inactive compounds). The database selected for model development was a published set of chemicals documented to cause respiratory hypersensitivity in humans. Two models were generated that differed only in that one model included explicate hydrogen containing fragments. The predictive abilities of the models were tested using leave-one-out cross-validation tests. One model had a sensitivity of 0.94 and specificity of 0.87 yielding an overall correct prediction of 91%. The second model had a sensitivity of 0.89, specificity of 0.95 and overall correct prediction of 92%. The demonstrated predictive capabilities of the cat-SAR approach, together with its modeling flexibility and design transparency, suggest the potential for its widespread applicability to toxicity prediction and for deriving mechanistic insight into toxicologic effects.     

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.     

Use of 13C NMR spectrometric data to produce a predictive model of estrogen receptor binding activity

We have developed a spectroscopic data-activity relationship (SDAR) model based on 13C NMR spectral data for 30 estrogenic chemicals whose relative binding affinities (RBA) are available for the alpha (ERalpha) and beta (ERbeta) estrogen receptors. The SDAR models segregated the 30 compounds into strong and medium binding affinities. The SDAR model gave a leave-one-out (LOO) cross-validation of 90%. Two compounds that were classified incorrectly in the SDAR model were in the transition zone between classifications. Real and predicted 13C NMR chemical shifts were used with test compounds to evaluate the predictive behavior of the SDAR model. The 13C NMR SDAR model using predicted 13C NMR data for the test compounds provides a rapid, reliable, and simple way to screen whether a compound binds to the estrogen receptors.     

Computational toxicology and the generation of mechanistic hypotheses: gamma-butyrolactone

In this study, we use SAR approaches in an attempt to elucidate the action of gamma-butyrolactone (GBL), an illicit drug and a dietary supplement, that can cause coma and deaths in humans while exhibiting low systemic toxicity towards rodents. The lack of systemic toxicity of GBL and of its metabolite(s) was also predicted by validated SAR models. In fact using diverse SAR models, the only significant SAR prediction was that GBL had the potential for inhibiting human cytochrome P4502D6 (CYP2D6). However, inhibition of that isozyme is not necessarily associated with toxicity. It is suggested that GBL users also abuse other substances. When GBL inhibits CYP2D6 this may prevent the CYP2D6-mediated detoxification of other toxicants simultaneously consumed by the GBL user.     

Automated information extraction and structure-activity relationship analysis of cytochrome P450 substrates

Information on CYP-chemical interactions was comprehensively explored by a text-mining technique, to confirm our previous structure-activity relationship model for CYP substrates (Yamashita et al. J. Chem. Inf. Model. 2008, 48, 364-369). The text-mining technique is based on natural language processing and can extract chemical names and their interaction patterns according to sentence context. After chemicals were automatically extracted and classified into CYP substrates, inhibitors, and inducers, 709 substrates were retrieved from the PubChem database and categorized as 216, 145, 136, 217, 156, and 379 substrates for CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4, respectively. Although the previous classification model was developed using data from only 161 compounds, the model classified the substrates found by text-mining analysis with reasonable accuracy. This confirmed the validity of both the multi-objective classification model for CYP substrates and the text-mining procedure.     

Computational Toxicology and the Generation of Mechanistic Hypotheses: γ-Butyrolactone

Abstract

In this study, we use SAR approaches in an attempt to elucidate the action of γ-butyrolactone (GBL), an illicit drug and a dietary supplement, that can cause coma and deaths in humans while exhibiting low systemic toxicity towards rodents.

The lack of systemic toxicity of GBL and of its metabolite(s) was also predicted by validated SAR models. In fact using diverse SAR models, the only significant SAR prediction was that GBL had the potential for inhibiting human cytochrome P4502D6 (CYP2D6). However, inhibition of that isozyme is not necessarily associated with toxicity. It is suggested that GBL users also abuse other substances. When GBL inhibits CYP2D6 this may prevent the CYP2D6-mediated detoxification of other toxicants simultaneously consumed by the GBL user.     

Predictive models for cytochrome p450 isozymes based on quantitative high throughput screening data

The human cytochrome P450 (CYP450) isozymes are the most important enzymes in the body to metabolize many endogenous and exogenous substances including environmental toxins and therapeutic drugs. Any unnecessary interactions between a small molecule and CYP450 isozymes may raise a potential to disarm the integrity of the protection. Accurately predicting the potential interactions between a small molecule and CYP450 isozymes is highly desirable for assessing the metabolic stability and toxicity of the molecule. The National Institutes of Health Chemical Genomics Center (NCGC) has screened a collection of over 17,000 compounds against the five major isozymes of CYP450 (1A2, 2C9, 2C19, 2D6, and 3A4) in a quantitative high throughput screening (qHTS) format. In this study, we developed support vector classification (SVC) models for these five isozymes using a set of customized generic atom types. The CYP450 data sets were randomly split into equal-sized training and test sets. The optimized SVC models exhibited high predictive power against the test sets for all five CYP450 isozymes with accuracies of 0.93, 0.89, 0.89, 0.85, and 0.87 for 1A2, 2C9, 2C19, 2D6, and 3A4, respectively, as measured by the area under the receiver operating characteristic (ROC) curves. The important atom types and features extracted from the five models are consistent with the structural preferences for different CYP450 substrates reported in the literature. We also identified novel features with significant discerning power to separate CYP450 actives from inactives. These models can be useful in prioritizing compounds in a drug discovery pipeline or recognizing the toxic potential of environmental chemicals.