Chemicals In Training Sets Research Articles

Study on the molecular mechanism of anti-liver cancer effect of Evodiae fructus by network pharmacology and QSAR model.

Introduction: Evodiae Fructus (EF) is the dried, near ripe fruit of Euodia rutaecarpa (Juss.) Benth in Rutaceae. Numerous studies have demonstrated its anti-liver cancer properties. However, the molecular mechanism of Evodiae fructus against liver cancer and its structure-activity connection still require clarification. Methods: We utilized network pharmacology and a QSAR (2- and 3-dimensional) model to study the anti-liver cancer effect of Evodiae fructus. First, by using network pharmacology to screen the active substances and targets of Evodiae fructus, we investigated the signaling pathways involved in the anti-liver cancer actions of Evodiae fructus. The 2D-QSAR pharmacophore model was then used to predict the pIC50 values of compounds. The hiphop method was used to create an ideal 3D-QSAR pharmacophore model for the prediction of Evodiae fructus compounds. Finally, molecular docking was used to validate the rationality of the pharmacophore, and molecular dynamics was used to disclose the stability of the compounds by assessing the trajectories in 10ns using RMSD, RMSF, Rg, and hydrogen bonding metrics. Results: In total, 27 compounds were acquired from the TCMSP and TCM-ID databases, and 45 intersection targets were compiled using Venn diagrams. Network integration analysis was used in this study to identify SRC as a primary target. Key pathways were discovered by KEGG pathway analysis, including PD-L1 expression and PD-1 checkpoint pathway, EGFR tyrosine kinase inhibitor resistance, and ErbB signaling pathway. Using a 2D-QSAR pharmacophore model and the MLR approach to predict chemical activity, ten highly active compounds were found. Two hydrophobic features and one hydrogen bond acceptor feature in the 3D-QSAR pharmacophore model were validated by training set chemicals. The results of molecular docking revealed that 10 active compounds had better docking scores with SRC and were linked to residues via hydrogen and hydrophobic bonds. Molecular dynamics was used to show the structural stability of obacunone, beta-sitosterol, and sitosterol. Conclusion:Pharmacophore 01 has high selectivity and the ability to distinguish active and inactive compounds, which is the optimal model for this study. Obacunone has the optimal binding ability with SRC. The pharmacophore model proposed in this study provides theoretical support for further screening effective anti-cancer Chinese herbal compounds and optimizing the compound structure.

Aneugen Molecular Mechanism Assay: Proof-of-Concept With 27 Reference Chemicals.

A tiered bioassay and data analysis scheme is described for elucidating the most common molecular targets responsible for chemical-induced in vitro aneugenicity: tubulin destabilization, tubulin stabilization, and inhibition of mitotic kinase(s). To evaluate this strategy, TK6 cells were first exposed to each of 27 presumed aneugens over a range of concentrations. After 4 and 24 h of treatment, γH2AX, p53, phospho-histone H3 (p-H3), and polyploidization biomarkers were evaluated using the MultiFlow DNA Damage Assay Kit. The assay identified 27 of 27 chemicals as genotoxic, with 25 exhibiting aneugenic signatures, 1 aneugenic and clastogenic, and 1 clastogenic. Subsequently, a newly described follow-up assay was employed to investigate the aneugenic agents' molecular targets. For these experiments, TK6 cells were exposed to each of 26 chemicals in the presence of 488 Taxol. After 4 h, cells were lysed and the liberated nuclei and mitotic chromosomes were stained with a nucleic acid dye and labeled with fluorescent antibodies against p-H3 and Ki-67. Flow cytometric analyses revealed that alterations to 488 Taxol-associated fluorescence were only observed with tubulin binders-increases in the case of tubulin stabilizers, decreases with destabilizers. Mitotic kinase inhibitors with known Aurora kinase B inhibiting activity were the only aneugens that dramatically decreased the ratio of p-H3-positive to Ki-67-positive nuclei. Unsupervised hierarchical clustering based on 488 Taxol fluorescence and p-H3: Ki-67 ratios clearly distinguished compounds with these disparate molecular mechanisms. Furthermore, a classification algorithm based on an artificial neural network was found to effectively predict molecular target, as leave-one-out cross-validation resulted in 25/26 agreement with a priori expectations. These results are encouraging, as they suggest that an adequate number of training set chemicals, in conjunction with a machine learning algorithm based on 488 Taxol, p-H3, and Ki-67 responses, can reliably elucidate the most commonly encountered aneugenic molecular targets.

Genotoxic mode of action predictions from a multiplexed flow cytometric assay and a machine learning approach.

Several endpoints associated with cellular responses to DNA damage as well as overt cytotoxicity were multiplexed into a miniaturized, "add and read" type flow cytometric assay. Reagents included a detergent to liberate nuclei, RNase and propidium iodide to serve as a pan-DNA dye, fluorescent antibodies against γH2AX, phospho-histone H3, and p53, and fluorescent microspheres for absolute nuclei counts. The assay was applied to TK6 cells and 67 diverse reference chemicals that served as a training set. Exposure was for 24 hrs in 96-well plates, and unless precipitation or foreknowledge about cytotoxicity suggested otherwise, the highest concentration was 1 mM. At 4- and 24-hrs aliquots were removed and added to microtiter plates containing the reagent mix. Following a brief incubation period robotic sampling facilitated walk-away data acquisition. Univariate analyses identified biomarkers and time points that were valuable for classifying agents into one of three groups: clastogenic, aneugenic, or non-genotoxic. These mode of action predictions were optimized with a forward-stepping process that considered Wald test p-values, receiver operator characteristic curves, and pseudo R(2) values, among others. A particularly high performing multinomial logistic regression model was comprised of four factors: 4 hr γH2AX and phospho-histone H3 values, and 24 hr p53 and polyploidy values. For the training set chemicals, the four-factor model resulted in 94% concordance with our a priori classifications. Cross validation occurred via a leave-one-out approach, and in this case 91% concordance was observed. A test set of 17 chemicals that were not used to construct the model were evaluated, some of which utilized a short-term treatment in the presence of a metabolic activation system, and in 16 cases mode of action was correctly predicted. These initial results are encouraging as they suggest a machine learning strategy can be used to rapidly and reliably predict new chemicals' genotoxic mode of action based on data from an efficient and highly scalable multiplexed assay.

Acute Toxicity-Supported Chronic Toxicity Prediction: A k-Nearest Neighbor Coupled Read-Across Strategy.

A k-nearest neighbor (k-NN) classification model was constructed for 118 RDT NEDO (Repeated Dose Toxicity New Energy and industrial technology Development Organization; currently known as the Hazard Evaluation Support System (HESS)) database chemicals, employing two acute toxicity (LD50)-based classes as a response and using a series of eight PaDEL software-derived fingerprints as predictor variables. A model developed using Estate type fingerprints correctly predicted the LD50 classes for 70 of 94 training set chemicals and 19 of 24 test set chemicals. An individual category was formed for each of the chemicals by extracting its corresponding k-analogs that were identified by k-NN classification. These categories were used to perform the read-across study for prediction of the chronic toxicity, i.e., Lowest Observed Effect Levels (LOEL). We have successfully predicted the LOELs of 54 of 70 training set chemicals (77%) and 14 of 19 test set chemicals (74%) to within an order of magnitude from their experimental LOEL values. Given the success thus far, we conclude that if the k-NN model predicts LD50 classes correctly for a certain chemical, then the k-analogs of such a chemical can be successfully used for data gap filling for the LOEL. This model should support the in silico prediction of repeated dose toxicity.

Characterization and validation of an in silico toxicology model to predict the mutagenic potential of drug impurities*

Control and minimization of human exposure to potential genotoxic impurities found in drug substances and products is an important part of preclinical safety assessments of new drug products. The FDA's 2008 draft guidance on genotoxic and carcinogenic impurities in drug substances and products allows use of computational quantitative structure–activity relationships (QSAR) to identify structural alerts for known and expected impurities present at levels below qualified thresholds. This study provides the information necessary to establish the practical use of a new in silico toxicology model for predicting Salmonella t. mutagenicity (Ames assay outcome) of drug impurities and other chemicals. We describe the model's chemical content and toxicity fingerprint in terms of compound space, molecular and structural toxicophores, and have rigorously tested its predictive power using both cross-validation and external validation experiments, as well as case studies. Consistent with desired regulatory use, the model performs with high sensitivity (81%) and high negative predictivity (81%) based on external validation with 2368 compounds foreign to the model and having known mutagenicity. A database of drug impurities was created from proprietary FDA submissions and the public literature which found significant overlap between the structural features of drug impurities and training set chemicals in the QSAR model. Overall, the model's predictive performance was found to be acceptable for screening drug impurities for Salmonella mutagenicity.

Identification of the Structural Requirements for Mutagencitiy, by Incorporating Molecular Flexibility and Metabolic Activation of Chemicals. II. General Ames Mutagenicity Model

The tissue metabolic simulator (TIMES) modeling approach is a hybrid expert system that couples a metabolic simulator together with structure toxicity rules, underpinned by structural alerts, to predict interaction of chemicals or their metabolites with target macromolecules. Some of the structural alerts representing the reactivity pattern-causing effect could interact directly with the target whereas others necessitated a combination with two- or three-dimensional quantitative structure-activity relationship models describing the firing of the alerts from the rest of the molecules. Recently, TIMES has been used to model bacterial mutagenicity [Mekenyan, O., Dimitrov, S., Serafimova, R., Thompson, E., Kotov, S., Dimitrova, N., and Walker, J. (2004) Identification of the structural requirements for mutagenicity by incorporating molecular flexibility and metabolic activation of chemicals I: TA100 model. Chem. Res. Toxicol. 17 (6), 753-766]. The original model was derived for a single tester strain, Salmonella typhimurium (TA100), using the Ames test by the National Toxicology Program (NTP). The model correctly identified 82% of the primary acting mutagens, 94% of the nonmutagens, and 77% of the metabolically activated chemicals in a training set. The identified high correlation between activities across different strains changed the initial strategic direction to look at the other strains in the next modeling developments. In this respect, the focus of the present work was to build a general mutagenicity model predicting mutagenicity with respect to any of the Ames tester strains. The use of all reactivity alerts in the model was justified by their interaction mechanisms with DNA, found in the literature. The alerts identified for the current model were analyzed by comparison with other established alerts derived from human experts. In the new model, the original NTP training set with 1341 structures was expanded by 1626 proprietary chemicals provided by BASF AG. Eventually, the training set consisted of 435 chemicals, which are mutagenic as parents, 397 chemicals that are mutagenic after S9 metabolic activation, and 2012 nonmutagenic chemicals. The general mutagenicity model was found to have 82% sensitivity, 89% specificity, and 88% concordance for training set chemicals. The model applicability domain was introduced accounting for similarity (structural, mechanistic, etc.) between predicted chemicals and training set chemicals for which the model performs correctly.

Non-stochastic and stochastic linear indices of the molecular pseudograph’s atom-adjacency matrix: a novel approach for computational in silico screening and “rational” selection of new lead antibacterial agents

A novel approach (TOMOCOMD-CARDD) to computer-aided "rational" drug design is illustrated. This approach is based on the calculation of the non-stochastic and stochastic linear indices of the molecular pseudograph's atom-adjacency matrix representing molecular structures. These TOMOCOMD-CARDD descriptors are introduced for the computational (virtual) screening and "rational" selection of new lead antibacterial agents using linear discrimination analysis. The two structure-based antibacterial-activity classification models, including non-stochastic and stochastic indices, classify correctly 91.61% and 90.75%, respectively, of 1525 chemicals in training sets. These models show high Matthews correlation coefficients (MCC=0.84 and 0.82). An external validation process was carried out to assess the robustness and predictive power of the model obtained. These QSAR models permit the correct classification of 91.49% and 89.31% of 505 compounds in an external test set, yielding MCCs of 0.84 and 0.79, respectively. The TOMOCOMD-CARDD approach compares satisfactorily with respect to nine of the most useful models for antimicrobial selection reported to date. Finally, an in silico screening of 87 new chemicals reported in the anti-infective field with antibacterial activities is developed showing the ability of the TOMOCOMD-CARDD models to identify new lead antibacterial compounds.

QSAR Applicability Domain Estimation by Projection of the Training Set in Descriptor Space: A Review

As the use of Quantitative Structure Activity Relationship (QSAR) models for chemical management increases, the reliability of the predictions from such models is a matter of growing concern. The OECD QSAR Validation Principles recommend that a model should be used within its applicability domain (AD). The Setubal Workshop report provided conceptual guidance on defining a (Q)SAR AD, but it is difficult to use directly. The practical application of the AD concept requires an operational definition that permits the design of an automatic (computerised), quantitative procedure to determine a models AD. An attempt is made to address this need, and methods and criteria for estimating AD through training set interpolation in descriptor space are reviewed. It is proposed that response space should be included in the training set representation. Thus, training set chemicals are points in n-dimensional descriptor space and m-dimensional model response space. Four major approaches for estimating interpolation regions in a multivariate space are reviewed and compared: range, distance, geometrical, and probability density distribution.

QSAR prioritization of chemical inventories for endocrine disruptor testing

Binding affinity between chemicals and the estrogen receptor (ER) serves as an indicator of the potential to cause endocrine disruption through this receptor-mediated endocrine pathway. Estimating ER-binding affinity is, therefore, one strategic approach to reducing the costs of screening chemicals for potential risks of endocrine disruption. While measuring ER binding with in vitro assays may be the first choice in prioritizing chemicals for additional in vitro or in vivo estrogenicity testing, the time and costs associated with screening thousands of chemicals is prohibitive. Recent advances in 3D modeling of the reactivity of flexible structures make in silico methods for estimating ER binding possible. One technique, the common reactivity pattern (COREPA) approach, was applied to development of reactivity patterns for ER relative binding affinity based on global nucleophilicity, interatomic distances between nucleophilic sites, and local electron donor capability of the nucleophilic sites. The reactivity patterns provided descriptor profiles for order-of-magnitude RBA ranges of training set chemicals. An exploratory expert system was subsequently developed to predict RBA and rank chemicals with respect to potential estrogenicity. A strategy is presented for extending initial exploratory 3D QSAR models beyond current training sets to increase applicability to more diverse structures in large chemical inventories.