Development of classification and regression based QSAR models to predict rodent carcinogenic potency using oral slope factor

ARKA: a framework of dimensionality reduction for machine-learning classification modeling, risk assessment, and data gap-filling of sparse environmental toxicity data

Due to the lack of experimental toxicity data for environmental chemicals, there arises a need to fill data gaps by in silico approaches. One of the most commonly used in silico approaches for toxicity assessment of small datasets is the Quantitative Structure-Activity Relationship (QSAR), which generates predictive models for the efficient prediction of query compounds. However, the reliability of the predictions from QSARs derived from small datasets is often questionable from a statistical point of view. This is due to the presence of a larger number of descriptors as compared to the number of training compounds, which reduces the degree of freedom of the developed model. To reduce the overall prediction error for a particular QSAR model, we have proposed here the computation of the novel Arithmetic Residuals in K-groups Analysis (ARKA) descriptors. We have reduced the number of modeling descriptors in a supervised manner by partitioning them into K classes (K = 2 here) depending on the higher mean normalized values of the descriptors to a particular response class, thus preventing the loss of chemical information. A scatter plot of the data points using the values of two ARKA descriptors (ARKA_2 vs. ARKA_1) can potentially identify activity cliffs, less confident data points, and less modelable data points. We have used here five representative environmentally relevant endpoints (skin sensitization, earthworm toxicity, milk/plasma partitioning, algal toxicity, and rodent carcinogenicity of hazardous chemicals) with graded responses to which the ARKA framework was applied for classification modeling. On comparing the performance of the models generated using conventional QSAR descriptors and the ARKA descriptors, the prediction quality of the models derived from ARKA descriptors was found, based on multiple graded-data validation metrics-derived decision criteria, much better than the models derived from QSAR descriptors signifying the potential of ARKA descriptors in ecotoxicological classification modeling of small data sets. Additionally, this holds true for the Read-Across approach as well, since the Read-Across predictions using ARKA descriptors supersede the predictions generated from QSAR descriptors. For the ease of users, a Java-based expert system has been developed that computes the ARKA descriptors from the input of QSAR descriptors.

In Silico Methods for Carcinogenicity Assessment

Screening compounds for potential carcinogenicity is of major importance for prevention of environmentally induced cancers. A large sequence of predictive models, ranging from short-term biological assays (e.g., mutagenicity tests) to theoretical models, has been attempted in this field. Theoretical approaches such as (Q)SAR are highly desirable for identifying carcinogens, since they actively promote the replacement, reduction, and refinement of animal tests. This chapter reports and describes some of the most noted (Q)SAR models based on human expert knowledge and statistical approaches, aiming at predicting the carcinogenicity of chemicals. Additionally, the performance of the selected models has been evaluated, and the results are interpreted in details by applying these predictive models to some pharmaceutical molecules.

First structure-activity relationship analysis of SARS-CoV-2 virus main protease (Mpro) inhibitors: an endeavor on COVID-19 drug discovery

Main protease (Mpro) of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) intervenes in the replication and transcription processes of the virus. Hence, it is a lucrative target for anti-viral drug development. In this study, molecular modeling analyses were performed on the structure activity data of recently reported diverse SARS-CoV-2 Mpro inhibitors to understand the structural requirements for higher inhibitory activity. The classification-based quantitative structure-activity relationship (QSAR) models were generated between SARS-CoV-2 Mpro inhibitory activities and different descriptors. Identification of structural fingerprints to increase or decrease in the inhibitory activity was mapped for possible inclusion/exclusion of these fingerprints in the lead optimization process. Challenges in ADME properties of protease inhibitors were also discussed to overcome the problems of oral bioavailability. Further, depending on the modeling results, we have proposed novel as well as potent SARS-CoV-2 Mpro inhibitors.

QSAR Models for Human Carcinogenicity: An Assessment Based on Oral and Inhalation Slope Factors

Carcinogenicity is a crucial endpoint for the safety assessment of chemicals and products. During the last few decades, the development of quantitative structure-activity relationship ((Q)SAR) models has gained importance for regulatory use, in combination with in vitro testing or expert-based reasoning. Several classification models can now predict both human and rat carcinogenicity, but there are few models to quantitatively assess carcinogenicity in humans. To our knowledge, slope factor (SF), a parameter describing carcinogenicity potential used especially for human risk assessment of contaminated sites, has never been modeled for both inhalation and oral exposures. In this study, we developed classification and regression models for inhalation and oral SFs using data from the Risk Assessment Information System (RAIS) and different machine learning approaches. The models performed well in classification, with accuracies for the external set of 0.76 and 0.74 for oral and inhalation exposure, respectively, and r2 values of 0.57 and 0.65 in the regression models for oral and inhalation SFs in external validation. These models might therefore support regulators in (de)prioritizing substances for regulatory action and in weighing evidence in the context of chemical safety assessments. Moreover, these models are implemented on the VEGA platform and are now freely downloadable online.

An in vitro and in silico investigation of human pregnane X receptor agonistic activity of poly- and perfluorinated compounds using the heuristic method-best subset and comparative similarity indices analysis

Poly- and perfluorinated compounds (PFCs) may induce potential endocrine-disrupting hormonal effects. However, the molecular mechanism of the toxicology of PFCs remains unclear, and the insufficient information is available on the biological activities of PFCs at present. In this study, the cell-based reporter gene assays were used to determine the agonistic activity of PFCs on the human pregnane X receptor (hPXR). The heuristic method combined with best subset modeling (HM-BSM) based on Dragon descriptors and comparative similarity indices analysis (CoMSIA) were employed to build classical quantitative structure-activity relationship (QSAR) and three-dimensional QSAR models, respectively. The applicability domain (AD) of the classical QSAR model was assessed. Both the HM-BSM and CoMSIA approaches demonstrated good robustness, predictive ability, and mechanistic interpretability. The r² and leave-one-out cross-validation squared correlated coefficient (q²_LOO) values were 0.872 and 0.759 for the HM-BSM, and 0.976 and 0.751 for the CoMSIA model, respectively. The hPXR agonistic activity of the PFCs predicted by the built HM-BSM and CoMSIA agreed well with experimental activity, with root mean square error (RMSE) values of 0.0803 and 0.117, respectively, and external validation squared correlated coefficients (q²_EXT) of 0.972 and 0.932, respectively. The hPXR agonistic activity of PFCs was related to their molecular polarizability, charge and atomic mass. Hydrogen bonding and hydrophobic interactions constituted the primary intermolecular forces between PFCs and the hPXR. The developed models were used to screen the PFCs with high hPXR agonistic activity.

Toward comprehension of multiple human cells uptake of engineered nano metal oxides: quantitative inter cell line uptake specificity (QICLUS) modeling

To address the nanomaterial exposure threat, it is imperative to understand how nanomaterials are recognized, internalized, and distributed within diverse cell systems. Targeting of nanomaterials to a specific cell type is generally attained through the modification of the nanoparticle (NP) surface leading to required cellular uptake. The enhanced cellular uptake to normal cells can direct to the higher interaction of NPs with subcellular organelles resulting the provocation of various signaling pathways. The successes of NPs rely on the prospect for the synthesis of functionalized NPs with necessary properties and their enhanced potential for cellular uptake for specific targeting. In the present study, we have modeled the cellular uptake of 109 surface modifiers of metal oxide nanoparticles (MNPs) for three different cell lines: HUVEC (Human endothelial cells), U937 (human macrophage cells), and PaCa2 (cancer cell lines). Along with the quantitative structure-activity relationship (QSAR) models, for the very first time we have developed and performed quantitative inter cell line uptake specificity (QICLUS) modeling to identify the physicochemical properties, as well as majorly structural fragments responsible for cellular uptake differences between two specific cell lines. The present work provides a comprehensive understanding of the cellular uptake of MNPs and the underlying structural parameters controlling the nano-cellular interactions. This phenomenon has also been analyzed from the QSAR and QICLUS models that concluded the functional groups of surface modifiers like amine, anhydride, halogen atoms, nitro group, acids have the dominating roles for the uptake of MNPs into the cell lines. Thus, the developed models may be used for designing of novel surface modifiers of MNPs of desired characteristics for proper cell-NPs interactions, as well as in the context of virtual screening aspect. Moreover, the MNP-cell interactions can give some idea about the toxicity for target-specific drug delivery treatment as higher cellular uptake is required for specific cells to treat the disease and lower uptake to the neighboring cells for lower toxicity.

In Silico Estimation of Chemical Carcinogenicity with Binary and Ternary Classification Methods

Carcinogenicity is one of the most concerned properties of chemicals to human health, thus it is important to identify chemical carcinogenicity as early as possible. In this study, 829 diverse compounds with rat carcinogenicity were collected from Carcinogenic Potency Database (CPDB). Using six types of fingerprints to represent the molecules, 30 binary and ternary classification models were generated to predict chemical carcinogenicity by five machine learning methods. The models were evaluated by an external validation set containing 87 chemicals from ISSCAN database. The best binary model was developed by MACCS keys and kNN algorithm with predictive accuracy at 83.91 %, while the best ternary model was also generated by MACCS keys and kNN algorithm with overall accuracy at 80.46 %. Furthermore, the best binary and ternary classification models were used to estimate carcinogenicity of tobacco smoke components containing 2251 compounds. 981 ones were predicted as carcinogens by binary classification model, while 110 compounds were predicted as strong carcinogens and 807 ones as weak carcinogens by ternary classification model. The results indicated that our models would be helpful for prediction of chemical carcinogenicity.

Computational toxicology, friend or foe?

Computational approaches offer the attraction of being both fast and cheap to run being able to process thousands of chemical structures in a few minutes. As with all new technology, there is a tendency for these approaches to be hyped up and claims of reliability and performance may be exaggerated. So just how good are these computational methods?

Predictive modeling of chemical toxicity towards Pseudokirchneriella subcapitata using regression and classification based approaches

Biodiversity nurturing may be a valuable pathway in controlling chemical stress on the ecosystem. In the present work, in silico studies have been performed to develop regression based quantitative structure toxicity relationship (QSTR) models using a data set containing 105 organic chemicals for the prediction of 48-h chemical toxicity towards Pseudokirchneriella subcapitata. Classification based linear discriminant analysis (LDA) was also performed to distinguish chemicals into toxic and nontoxic groups using the same data set. The developed models were found to possess good predictive quality in terms of internal, external and overall validation parameters. The regression based QSTR model suggests that second order molecular connectivity index (molecular size and lipophilicity), density (aromaticity), relative shape of molecules (cyclicity/aromaticity), and specific molecular fragments of the chemicals are important properties of chemicals to exert their toxicity on P. subcapitata. The classification based LDA QSTR model suggested that fused ring aromatic systems, secondary carbon atom fragments, second order valence molecular connectivity indices (molecular size and branching) and molecular weight are the distinguishing features to differentiate chemicals into toxic and nontoxic groups.

Predicting the DPP-IV inhibitory activity pIC₅₀ based on their physicochemical properties

The second development program developed in this work was introduced to obtain physicochemical properties of DPP-IV inhibitors. Based on the computation of molecular descriptors, a two-stage feature selection method called mRMR-BFS (minimum redundancy maximum relevance-backward feature selection) was adopted. Then, the support vector regression (SVR) was used in the establishment of the model to map DPP-IV inhibitors to their corresponding inhibitory activity possible. The squared correlation coefficient for the training set of LOOCV and the test set are 0.815 and 0.884, respectively. An online server for predicting inhibitory activity pIC50 of the DPP-IV inhibitors as described in this paper has been given in the introduction.