Capacities of computer evaluation of hidden potential of phytochemicals of medicinal plants of the traditional Indian Ayurvedic medicine

Applicability of our computer programs PASS and PharmaExpert for prediction of biological activity spectra of rather complex and structurally diverse phytocomponents of medicinal plants, both separately and in combinations has been evaluated. For this purpose we have created the web-resource containing known information about structural formulas and biological activity of 1906 phytocomponents of 50 Ayurvedic medicinal plants used in Traditional Indian Medicine (TIM) (http://ayurveda.pharmaexpert.ru). The PASS training set was updated by addition of information about structure and biological activity of 946 natural compounds; then the training procedure and validation were performed, to estimate the quality of PASS prediction. It was shown that the differences between the average accuracy of prediction obtained in leave-5%-out cross-validation (94.467%) and in leave-one-out cross-validation (94.605%) are very small thus demonstrating high predictive ability of the program. Results of biological activity spectra prediction for all phytocomponents included in our database coincided in 83.5% of cases with known experimental data. Additional types of biological activity predicted with high probability indicate further promising directions for further studies of certain phytocomponents of some medicinal plants. The analysis of prediction results of sets of phytocomponents in each of 50 medicinal plants was made by the PharmaExpert software. Based on this analysis, we found that the combination of phytocomponents from Passiflora incarnata may exhibit nootropic, anticonvulsant, and antidepressant effects. Experiments carried out in mice models confirmed the predicted effects of P. incarnata extracts.     

Quantitative proteomic analysis of rat liver for carcinogenicity prediction in a 28-day repeated dose study

The potential of quantitative proteomic analysis to predict carcinogenicity of chemical compounds was investigated. Using 2D-DIGE, we analyzed the effects of 63 chemical compounds on protein expression in the rat liver after 28 daily doses. Types of carcinogens were categorized depending on the species and organ specificity. The carcinogen characteristic proteins for each classification were identified by Welch's t value. For evaluation of the predictive concordance we used support vector machines. The rat hepatic carcinogen-specific classification gave higher concordance than the other classification. The generalization performance was measured by leave-one-out cross-validation. For genotoxic and non-genotoxic compounds, a concordance of 79.3 and 76.5%, respectively, was obtained by the top 30 ranked proteins with Welch's t value. Furthermore, we found that the increase of the expression level of the stress response proteins as the common feature of poorly predicted chemical compounds in the leave-20%-out cross-validation. Quantitative proteomics could be promising technique for developing biomarker panels that can be used for carcinogenicity prediction. The list of proteins identified in this study and the zoomed gel images of the top ranked proteins in statistic analysis are provided in Supplementary Data.     

Genomic Models of Short-Term Exposure Accurately Predict Long-Term Chemical Carcinogenicity and Identify Putative Mechanisms of Action

Background

Despite an overall decrease in incidence of and mortality from cancer, about 40% of Americans will be diagnosed with the disease in their lifetime, and around 20% will die of it. Current approaches to test carcinogenic chemicals adopt the 2-year rodent bioassay, which is costly and time-consuming. As a result, fewer than 2% of the chemicals on the market have actually been tested. However, evidence accumulated to date suggests that gene expression profiles from model organisms exposed to chemical compounds reflect underlying mechanisms of action, and that these toxicogenomic models could be used in the prediction of chemical carcinogenicity.

Results

In this study, we used a rat-based microarray dataset from the NTP DrugMatrix Database to test the ability of toxicogenomics to model carcinogenicity. We analyzed 1,221 gene-expression profiles obtained from rats treated with 127 well-characterized compounds, including genotoxic and non-genotoxic carcinogens. We built a classifier that predicts a chemical''s carcinogenic potential with an AUC of 0.78, and validated it on an independent dataset from the Japanese Toxicogenomics Project consisting of 2,065 profiles from 72 compounds. Finally, we identified differentially expressed genes associated with chemical carcinogenesis, and developed novel data-driven approaches for the molecular characterization of the response to chemical stressors.

Conclusion

Here, we validate a toxicogenomic approach to predict carcinogenicity and provide strong evidence that, with a larger set of compounds, we should be able to improve the sensitivity and specificity of the predictions. We found that the prediction of carcinogenicity is tissue-dependent and that the results also confirm and expand upon previous studies implicating DNA damage, the peroxisome proliferator-activated receptor, the aryl hydrocarbon receptor, and regenerative pathology in the response to carcinogen exposure.     

Predicting chemical carcinogenesis in rodents: The state of the art in light of a comparative exercise

Within a recent comparative exercise, different approaches to the prediction of rodent carcinogenicity were challenged on a common set of chemicals bioassayed by the U.S. National Toxicology Program. The approaches were of very different natures. Some prediction systems looked for relationships between carcinogenicity and other, more quickly detectable biological events (activity-activity relationships, AAR). Some approaches tended to find structure-activity relationships (SAR). To give an objective evaluation of the results of the exercise, we have analyzed the rodent results and the predictions with the multivariate data analysis methods. The calculated performances varied according to the adopted carcinogenicity classification of the chemicals. When the four rodent results were summarized into a final + or − call, the Tennant approach (AAR method) showed the best performance (about 75% accuracy), whereas the best SAR systems had 60–65% accuracy. A common limitation of almost all the systems was the lack of specificity (too many false positives). Based on these results, better concordance was obtained when the input information was the very costly (and closer to the final endpoint) biological data, rather than the inexpensive (and farther from the endpoint) knowledge of the chemical structure. However, when the rodent results were summarized into a carcinogenicity classification that maintained, to some extent, the gradation intrinsic to the original experimental data, the performance of the AAR systems declined, and the SAR approaches showed a better performance. The difficulty in evaluating the various approaches was further complicated because of a fundamental difference in the approaches themselves: some approaches were ‘pure’ prediction methods (i.e. their predictions were rigorously based on information not inclusive of carcinogenicity); other approaches (e.g. Tennant, Weisburger) used ‘mixed’ information, inclusive of known carcinogenicity results from experiments performed before the NTP bioassays. As far as the SAR systems are concerned, their sets of predictions showed a fundamental similarity. This happened in spite of the extremely different procedures adopted to treat the chemical formula (initial information): very simple calculations (Benigni), intuition of the experts (Weisburger and Lijinsky), sophisticated computer programs (TOPKAT and CASE). The results of the Bakale Ke method, based on the experimental measurement of the chemical electrophilicity, and of the Salmonella typhimurium mutagenicity assay were similar to the patterns of predictions of the SAR methods.     

Markovian chemicals "in silico" design (MARCH-INSIDE), a promising approach for computer aided molecular design II: experimental and theoretical assessment of a novel method for virtual screening of fasciolicides

A novel method for in silico selection of fluckicidal drugs is introduced. Two QSARs that permit us to discriminate between fasciolicide and non-fasciolicide drugs (the first) and to outline some conclusions about the possible mechanism of action of a chemical (the second) are performed. The first model correctly classified 93.85% of compounds in the training series and 89.5% of the compounds in the predicting one. This model correctly classified 87.7, 93.8, 92.2 and 93.9% of compounds in leave- n-out cross validation procedures when n takes values from 2 to until 6. The model seems to be stable in around 92% of good classification in leave- n-out cross validation analysis when n>6. The second model correctly classified 70% of non-fasciolicide compounds, 85.71% of beta-tubulin inhibitors and 100% of proton ionophores in the training set. This model recognizes as proton ionophores 100% of any nitrosalicylanilides in the predicting series. Both models have a low p-level <0.05. Finally, the experimental assay of six organic chemicals by an in vivo test permit us to carry out an assessment of the model with a fairly good 100% agreement between experiment and theoretical prediction.     

Definitive relationships among chemical structure, carcinogenicity and mutagenicity for 301 chemicals tested by the U.S. NTP.

An analysis is presented in which are evaluated correlations among chemical structure, mutagenicity to Salmonella, and carcinogenicity to rats and mice among 301 chemicals tested by the U.S. NTP. Overall, there was a high correlation between structural alerts to DNA reactivity and mutagenicity, but the correlation of either property with carcinogenicity was low. If rodent carcinogenicity is regarded as a singular property of chemicals, then neither structural alerts nor mutagenicity to Salmonella are effective in its prediction. Given this, the database was fragmented and new correlations sought between the derived sub-groups. First, the 301 chemicals were segregated into six broad chemical groupings. Second, the rodent cancer data were partially segregated by target tissue. Using the previously assigned structural alerts to DNA reactivity (electrophilicity), the chemicals were split into 154 alerting chemicals and 147 non-alerting chemicals. The alerting chemicals were split into three chemical groups; aromatic amino/nitro-types, alkylating agents and miscellaneous structurally-alerting groups. The non-alerting chemicals were subjectively split into three broad categories; non-alerting, non-alerting containing a non-reactive halogen group, and non-alerting chemical with minor concerns about a possible structural alert. The tumor data for all 301 chemicals are re-presented according to these six chemical groupings. The most significant findings to emerge from comparisons among these six groups of chemicals were as follows: (a) Most of the rodent carcinogens, including most of the 2-species and/or multiple site carcinogens, were among the structurally alerting chemicals. (b) Most of the structurally alerting chemicals were mutagenic; 84% of the carcinogens and 66% of the non-carcinogens. 100% of the 33 aromatic amino/nitro-type 2-species carcinogens were mutagenic. Thus, for structurally alerting chemicals, the Salmonella assay showed high sensitivity and low specificity (0.84 and 0.33, respectively). (c) Among the 147 non-alerting chemicals less than 5% were mutagenic, whether they were carcinogens or non-carcinogens (sensitivity 0.04).     

The second National Toxicology Program comparative exercise on the prediction of rodent carcinogenicity: definitive results

Chemical carcinogenicity has been the target of a large array of attempts to create alternative predictive models, ranging from short-term biological assays (e.g. mutagenicity tests) to theoretical models. Among the theoretical models, the application of the science of structure-activity relationships (SAR) has earned special prominence. A crucial element is the independent evaluation of the predictive ability. In the past decade, there have been two fundamental comparative exercises on the prediction of chemical carcinogenicity, held under the aegis to the US National Toxicology Program (NTP). In both exercises, the predictions were published before the animal data were known, thus using a most stringent criterion of predictivity. We analyzed the results of the first comparative exercise in a previous paper [Mutat. Res. 387 (1997) 35]; here, we present the complete results of the second exercise, and we analyze and compare the prediction sets. The range of accuracy values was quite large: the systems that performed best in this prediction exercise were in the range 60-65% accuracy. They included various human experts approaches (e.g. Oncologic) and biologically based approaches (e.g. the experimental transformation assay in Syrian hamster embryo (SHE) cells). The main difficulty for the structure-activity relationship-based approaches was the discrimination between real carcinogens, and non-carcinogens containing structural alerts (SA) for genotoxic carcinogenicity. It is shown that the use of quantitative structure-activity relationship models, when possible, can contribute to overcome the above problem. Overall, given the uncertainty linked to the predictions, the predictions for the individual chemicals cannot be taken at face value; however, the general level of knowledge available today (especially for genotoxic carcinogens) allows qualified human experts to operate a very efficient priority setting of large sets of chemicals.     

Knowledge-based battery design of short-term tests based on dose information

A construction of batteries of short-term tests (STTs) is described which is based on a classification of 73 chemicals in regard to their carcinogenicity. The 73 chemicals were studied within the U.S. National Toxicology Program (Ashby and Tennant, 1988). The batteries are validated using the classification of 35 additional chemicals. They are defined by logically structured combinations of rules. The single rules are defined by the z-scores of the logarithmic values of the limiting doses obtained from the 4 in vitro STTs used in the study by Ashby and Tennant. The limiting dose is defined as the lowest effective dose or the highest ineffective dose (Waters et al., 1987). The batteries are constructed by minimizing the number of disagreements with the classification by Ashby and Tennant. Compared with the results obtained from single STTs, 2 batteries of 3 STTs have higher concordances with the carcinogenicity data, namely 70% for the NTP data and 74-77% for the independent test data. In addition, a theoretical result shows that the proposed battery design, for a large enough learning set of chemicals, leads to results which are replicated with high probability on a large enough validation set. Based on the first results obtained with a limited number of chemicals it is concluded that the knowledge-based battery design is worth further development.     

Simple stochastic fingerprints towards mathematical modeling in biology and medicine. 3. ocular irritability classification model

MARCH-INSIDE methodology and a statistical classification method—linear discriminant analysis (LDA)—is proposed as an alternative method to the Draize eye irritation test. This methodology has been successfully applied to a set of 46 neutral organic chemicals, which have been defined as ocular irritant or nonirritant. The model allow to categorize correctly 37 out of 46 compounds, showing an accuracy of 80.46%. Specifically, this model demonstrates the existence of a good categorization average of 91.67 and 76.47% for irritant and nonirritant compounds, respectively. Validation of the model was carried out using two cross-validation tools: Leave-one-out (LOO) and leave-group-out (LGO), showing a global predictability of the model of 71.7 and 70%, respectively. The average of coincidence of the predictions between leave-one-out/leave-group-out studies and train set were 91.3% (42 out of 46 cases)/89.1% (41 out of 46 cases) proving the robustness of the model obtained. Ocular irritancy distribution diagram is carried out in order to determine the intervals of the property where the probability of finding an irritant compound is maximal relating to the choice of find a false nonirritant one. It seems that, until today, the present model may be the first predictive linear discriminant equation able to discriminate between eye irritant and nonirritant chemicals.     

Metabonomic models of human pancreatic cancer using 1D proton NMR spectra of lipids in plasma

In this study, we hypothesized that the altered insulin and glucose levels in male pancreatic cancer patients reported in a recent JAMA article would result in an altered lipid profile in the blood of pancreatic cancer patients when compared to controls (Stolzenberg-Solomon et al., 2005). Proton nuclear magnetic resonance (NMR) spectra of human lipophilic plasma extracts were used in order to build partial least squares discriminant function (PLS-DF) models that classified samples as belonging to the pancreatic control group or to the pancreatic cancer group. The sensitivity, specificity, and overall accuracy of the PLS-DF models based on 4 bins were 96%, 88%, and 92%, respectively. The sensitivity, specificity, and overall accuracy of the PLS-DF models based on 5 bins were 98%, 94%, and 96%, respectively. The sensitivity, specificity and overall accuracy of both the 4-bin and 5-bin PLS-DF models dropped only 1–2% during leave-25%-out cross-validation testing. Mass spectrometric profiling of phospholipids in plasma found three phosphatidylinositols that were significantly lower in pancreatic cancer patients than in healthy controls. The cancer models are based upon changes in lipid profiles that may provide a more sensitive and accurate diagnosis of pancreatic cancer than current methods that are based upon a single biomarker.     

Classification according to chemical structure, mutagenicity to Salmonella and level of carcinogenicity of a further 42 chemicals tested for carcinogenicity by the U.S. National Toxicology Program

This paper is an extension and update of an earlier review published in this journal (Ashby and Tennant, 1988). A summary of the rodent carcinogenicity bioassay data on a further 42 chemicals tested by the U.S. National Toxicology Program (NTP) is presented. An evaluation of each chemical for structural alerts to DNA-reactivity is also provided, together with a summary of its mutagenicity to Salmonella. The 42 chemicals were numbered and evaluated as an extension of the earlier analysis of 222 NTP chemicals. The activity patterns and conclusions derived from the earlier study remain unchanged for the larger group of 264 chemicals. Based on the extended database of 264 NTP chemicals, the sensitivity of the Salmonella assay for rodent carcinogens is 58% and the specificity for the non-carcinogens is 73%. A total of 32 chemicals were defined as equivocal for carcinogenicity and, of these, 11 (34%) are mutagenic to Salmonella. An evaluation is made of instances where predictions of carcinogenicity, based on structural alerts, disagree with the Salmonella mutagenicity result (12% of the database). The majority of the disagreements are for structural alerts on non-mutagens, and that places these alerts as a sensitive primary screen with a specificity lower than that of the Salmonella assay. That analysis indicates some need for assays complementary to the Salmonella test when screening for potential genotoxic carcinogens. It also reveals that the correlation between structural alerts and mutagenicity to Salmonella is probably greater than 90%. Chemicals predicted to show Michael-type alkylating activity (i.e., CH2 = CHX; where X = an electron-withdrawing group, e.g. acrylamide) have been confirmed as a structural alert, and the halomethanes (624 are possible) have been classified as structurally-alerting. To this end an extended carcinogen-alert model structure is presented. Among the 138 NTP carcinogens now reviewed, 45 (33%) are non-mutagenic to Salmonella and possess a chemical structure that does not alert to DNA-reactivity. These carcinogens therefore either illustrate the need for complementary genetic screening tests to the Salmonella assay, or they represent the group of non-genotoxic carcinogens referred to most specifically by Weisburger and Williams (1981); the latter concept is favoured.     

Characteristic substructures and properties in chemical carcinogens studied by the cascade model

MOTIVATION: Chemical carcinogenicity is an important subject in health and environmental sciences, and a reliable method is expected to identify characteristic factors for carcinogenicity. The predictive toxicology challenge (PTC) 2000-2001 has provided the opportunity for various data mining methods to evaluate their performance. The cascade model, a data mining method developed by the author, has the capability to mine for local correlations in data sets with a large number of attributes. The current paper explores the effectiveness of the method on the problem of chemical carcinogenicity. RESULTS: Rodent carcinogenicity of 417 compounds examined by the National Toxicology Program (NTP) was used as the training set. The analysis by the cascade model, for example, could obtain a rule 'Highly flexible molecules are carcinogenic, if they have no hydrogen bond acceptors in halogenated alkanes and alkenes'. Resulting rules are applied to predict the activity of 185 compounds examined by the FDA. The ROC analysis performed by the PTC organizers has shown that the current method has excellent predictive power for the female rat data. AVAILABILITY: The binary program of DISCAS 2.1 and samples of input data sets on Windows PC are available at http://www.clab.kwansei.ac.jp/mining/discas/discas.html upon request from the author. SUPPLEMENTARY INFORMATION: Summary of prediction results and cross validations is accessible via http://www.clab.kwansei.ac.jp/~okada/BIJ/BIJsupple.htm. Used rules and the prediction results for each molecule are also provided.     

A survey of the Predictive Toxicology Challenge 2000-2001

MOTIVATION: The Predictive Toxicology Challenge (PTC) was initiated to stimulate the development of advanced techniques for predictive toxicology models. The goal of this challenge was to compare different approaches for the prediction of rodent carcinogenicity, based on the experimental results of the US National Toxicology Program (NTP). RESULTS: 111 sets of predictions for 185 compounds have been evaluated on quantitative and qualitative scales to select the most predictive models and those with the highest toxicological relevance. The accuracy of the submitted predictions was between 25 and 79 %. An evaluation of the most accurate models by toxicological experts showed, that it is still hard for domain experts to interpret the submitted models and to put them into relation with toxicological knowledge. AVAILABILITY: PTC details and data can be found at: http://www.predictive-toxicology.org/ptc/.     

The SOS Chromotest, a colorimetric bacterial assay for genotoxins: validation study with 83 compounds

The SOS Chromotest is a simple bacterial colorimetric assay for genotoxicity. It is based on the measure of the induction of sfiA, a gene controlled by the general repressor of the SOS system in E. coli. Expression of sfiA is monitored by means of a gene fusion with lacZ, the structural gene for beta-galactosidase. We have examined 83 compounds of various chemical classes with the SOS Chromotest using a standard procedure. Comparison of the results with those obtained in the Mutatest (the Ames test) showed that most (90%) of the mutagenic compounds were also SOS inducers. For these compounds a quantitative correlation was observed between the mutagenic potency and the SOS-inducing potency (SOSIP). The case of the 10% remaining compounds giving conflicting results in the two tests is discussed. Sensitivity, specificity and accuracy for carcinogenicity prediction have been evaluated for the SOS Chromotest and the Mutatest using 73 chemicals for which carcinogenicity data were available. In spite of some differences, similar results were obtained in the two tests. The present data indicate that the SOS Chromotest has many practical advantages and may be used as a primary screening tool or as part of a battery of short-term tests for carcinogens.