Two new parameters based on distances in a receiver operating characteristic chart for the selection of classification models

Exploring different classification-dependent QSAR modelling strategies for HDAC3 inhibitors in search of meaningful structural contributors

Histone deacetylase 3 (HDAC3), a Zn2+-dependent class I HDACs, contributes to numerous disorders such as neurodegenerative disorders, diabetes, cardiovascular disease, kidney disease and several types of cancers. Therefore, the development of novel and selective HDAC3 inhibitors might be promising to combat such diseases. Here, different classification-based molecular modelling studies such as Bayesian classification, recursive partitioning (RP), SARpy and linear discriminant analysis (LDA) were conducted on a set of HDAC3 inhibitors to pinpoint essential structural requirements contributing to HDAC3 inhibition followed by molecular docking study and molecular dynamics (MD) simulation analyses. The current study revealed the importance of hydroxamate function for Zn2+ chelation as well as hydrogen bonding interaction with Tyr298 residue. The importance of hydroxamate function for higher HDAC3 inhibition was noticed in the case of Bayesian classification, recursive partitioning and SARpy models. Also, the importance of substituted thiazole ring was revealed, whereas the presence of linear alkyl groups with carboxylic acid function, any type of ester function, benzodiazepine moiety and methoxy group in the molecular structure can be detrimental to HDAC3 inhibition. Therefore, this study can aid in the design and discovery of effective novel HDAC3 inhibitors in the future.

A fragment-based exploration of diverse MMP-9 inhibitors through classification-dependent structural assessment

Matrix metalloproteinases (MMPs) are belonging to the Zn2+-dependent metalloenzymes. These can degenerate the extracellular matrix (ECM) that is entailed with various biological processes. Among the MMP family members, MMP-9 is associated with several pathophysiological circumstances. Apart from wound healing, remodeling of bone, inflammatory mechanisms, and rheumatoid arthritis, MMP-9 has also significant roles in tumor invasion and metastasis. Therefore, MMP-9 has been in the spotlight of anticancer drug discovery programs for more than a decade. In this present study, classification-based QSAR techniques along with fragment-based data mining have been carried out on divergent MMP-9 inhibitors to point out the important structural attributes. This current study may be able to elucidate the importance of several pivotal molecular fragments such as sulfonamide, hydroxamate, i-butyl, and ethoxy functions for imparting potential MMP-9 inhibition. These observations are in correlation with the ligand-bound co-crystal structures of MMP-9. Therefore, these findings are beneficial for the design and discovery of effective MMP-9 inhibitors in the future.

QSAR classification model for diverse series of antifungal agents based on binary coyote optimization algorithm

One of the recently developed metaheuristic algorithms, the coyote optimization algorithm (COA), has shown to perform better in a number of difficult optimization tasks. The binary form, BCOA, is used in this study as a solution to the descriptor selection issue in classifying diverse antifungal series. Z-shape transfer functions (ZTF) are evaluated to verify their efficiency in improving BCOA performance in QSAR classification based on classification accuracy (CA), the geometric mean of sensitivity and specificity (G-mean), and the area under the curve (AUC). The Kruskal-Wallis test is also applied to show the statistical differences between the functions. The efficacy of the best suggested transfer function, ZTF4, is further assessed by comparing it to the most recent binary algorithms. The results prove that ZTF, especially ZTF4, significantly improves the performance of the original BCOA. The ZTF4 function yields the best CA and G-mean of 99.03% and 0.992%, respectively. It shows the fastest convergence behaviour compared to other binary algorithms. It takes the fewest iterations to reach high classification performance and selects the fewest descriptors. In conclusion, the obtained results indicate the ability of the ZTF4-based BCOA to find the smallest subset of descriptors while maintaining the best classification accuracy performance.

Identification of Pharmacophoric Fragments of DYRK1A Inhibitors Using Machine Learning Classification Models

Dual-specific tyrosine phosphorylation regulated kinase 1 (DYRK1A) has been regarded as a potential therapeutic target of neurodegenerative diseases, and considerable progress has been made in the discovery of DYRK1A inhibitors. Identification of pharmacophoric fragments provides valuable information for structure- and fragment-based design of potent and selective DYRK1A inhibitors. In this study, seven machine learning methods along with five molecular fingerprints were employed to develop qualitative classification models of DYRK1A inhibitors, which were evaluated by cross-validation, test set, and external validation set with four performance indicators of predictive classification accuracy (CA), the area under receiver operating characteristic (AUC), Matthews correlation coefficient (MCC), and balanced accuracy (BA). The PubChem fingerprint-support vector machine model (CA = 0.909, AUC = 0.933, MCC = 0.717, BA = 0.855) and PubChem fingerprint along with the artificial neural model (CA = 0.862, AUC = 0.911, MCC = 0.705, BA = 0.870) were considered as the optimal modes for training set and test set, respectively. A hybrid data balancing method SMOTETL, a combination of synthetic minority over-sampling technique (SMOTE) and Tomek link (TL) algorithms, was applied to explore the impact of balanced learning on the performance of models. Based on the frequency analysis and information gain, pharmacophoric fragments related to DYRK1A inhibition were also identified. All the results will provide theoretical supports and clues for the screening and design of novel DYRK1A inhibitors.

Identification of Kukoamine A, Zeaxanthin, and Clexane as New Furin Inhibitors

The endogenous protease furin is a key protein in many different diseases, such as cancer and infections. For this reason, a wide range of studies has focused on targeting furin from a therapeutic point of view. Our main objective consisted of identifying new compounds that could enlarge the furin inhibitor arsenal; secondarily, we assayed their adjuvant effect in combination with a known furin inhibitor, CMK, which avoids the SARS-CoV-2 S protein cleavage by means of that inhibition. Virtual screening was carried out to identify potential furin inhibitors. The inhibition of physiological and purified recombinant furin by screening selected compounds, Clexane, and these drugs in combination with CMK was assayed in fluorogenic tests by using a specific furin substrate. The effects of the selected inhibitors from virtual screening on cell viability (293T HEK cell line) were assayed by means of flow cytometry. Through virtual screening, Zeaxanthin and Kukoamine A were selected as the main potential furin inhibitors. In fluorogenic assays, these two compounds and Clexane inhibited both physiological and recombinant furin in a dose-dependent way. In addition, these compounds increased physiological furin inhibition by CMK, showing an adjuvant effect. In conclusion, we identified Kukoamine A, Zeaxanthin, and Clexane as new furin inhibitors. In addition, these drugs were able to increase furin inhibition by CMK, so they could also increase its efficiency when avoiding S protein proteolysis, which is essential for SARS-CoV-2 cell infection.

Applying comparative molecular modelling techniques on diverse hydroxamate-based HDAC2 inhibitors: an attempt to identify promising structural features for potent HDAC2 inhibition

Histone deacetylase 2 (HDAC2) has been implicated in a variety of cardiovascular and neurodegenerative disorders as well as in cancers. Thus, HDAC2 has become an exclusive target for anticancer drug development. Therefore, the development of newer HDAC2 inhibitors in disease conditions is a prime goal to restrain such a scenario. Although a handful of HDAC inhibitors was accepted for the treatment of HDAC-related disease conditions, the non-selective nature of these entities is one of the major setbacks in the treatment of specific HDAC isoform-related pathophysiology. In this framework, the analyses of pre-existing molecules are essential to identify the important structural features that can fulfil the requirements for the cap and linker moieties to obtain potent and effective HDAC2 inhibition. Thus, in this study, the implementation of a combined comparative 2D and 3D molecular modelling techniques was done on a group of 92 diverse hydroxamate derivatives having a wide range of HDAC2 inhibitory potency. Besides other crucial features, this study upheld the importance of groups like triazole and benzyl moieties along with the molecular fields that are crucial for regulating HDAC2 inhibition. The outcomes of this study may be employed for the designing of HDAC2 inhibitors in future.

Machine Learning Models for the Classification of CK2 Natural Products Inhibitors with Molecular Fingerprint Descriptors

Casein kinase 2 (CK2) is considered an important target for anti-cancer drugs. Given the structural diversity and broad spectrum of pharmaceutical activities of natural products, numerous studies have been performed to prove them as valuable sources of drugs. However, there has been little study relevant to identifying structural factors responsible for their inhibitory activity against CK2 with machine learning methods. In this study, classification studies were conducted on 115 natural products as CK2 inhibitors. Seven machine learning methods along with six molecular fingerprints were employed to develop qualitative classification models. The performances of all models were evaluated by cross-validation and test set. By taking predictive accuracy(CA), the area under receiver operating characteristic (AUC), and (MCC)as three performance indicators, the optimal models with high reliability and predictive ability were obtained, including the Extended Fingerprint-Logistic Regression model (CA = 0.859, AUC = 0.826, MCC = 0.520) for training test andPubChem fingerprint along with the artificial neural model (CA = 0.826, AUC = 0.933, MCC = 0.628) for test set. Meanwhile, the privileged substructures responsible for their inhibitory activity against CK2 were also identified through a combination of frequency analysis and information gain. The results are expected to provide useful information for the further utilization of natural products and the discovery of novel CK2 inhibitors.

Robust classification-based molecular modelling of diverse chemical entities as potential SARS-CoV-2 3CLpro inhibitors: theoretical justification in light of experimental evidences

COVID-19 is the most unanticipated incidence of 2020 affecting the human population worldwide. Currently, it is utmost important to produce novel small molecule anti-SARS-CoV-2 drugs urgently that can save human lives globally. Based on the earlier SARS-CoV and MERS-CoV infection along with the general characters of coronaviral replication, a number of drug molecules have been proposed. However, one of the major limitations is the lack of experimental observations with different drug molecules. In this article, 70 diverse chemicals having experimental SARS-CoV-2 3CL^proinhibitory activity were accounted for robust classification-based QSAR analysis statistically validated with 4 different methodologies to recognize the crucial structural features responsible for imparting the activity. Results obtained from all these methodologies supported and validated each other. Important observations obtained from these analyses were also justified with the ligand-bound crystal structure of SARS-CoV-2 3CL^pro enzyme. Our results suggest that molecules should contain a 2-oxopyrrolidine scaffold as well as a methylene (hydroxy) sulphonic acid warhead in proper orientation to achieve higher inhibitory potency against SARS-CoV-2 3CL^pro. Outcomes of our study may be able to design and discover highly effective SARS-CoV-2 3CL^pro inhibitors as potential anticoronaviral therapy to crusade against COVID-19.

Evaluating the cytotoxicity of a large pool of metal oxide nanoparticles to Escherichia coli: Mechanistic understanding through In Vitro and In Silico studies

The toxic effect of eight metal oxide nanoparticles (MONPs) on Escherichia coli was experimentally evaluated following standard bioassay protocols. The obtained cytotoxicity ranking of these studied MONPs is Er₂O₃, Gd₂O₃, CeO₂, Co₂O₃, Mn₂O₃, Co₃O₄, Fe₃O₄/WO₃ (in descending order). The computed EC₅₀ values from experimental data suggested that Er₂O₃ and Gd₂O₃ were the most acutely toxic MONPs to E. coli. To identify the mechanism of toxicity of these 8 MONPs along with 17 other MONPs from our previous study, we employed seven classifications and machine learning (ML) algorithms including linear discriminant analysis (LDA), naïve bayes (NB), multinomial logistic regression (MLogitR), sequential minimal optimization (SMO), AdaBoost, J48, and random forest (RF). We also employed 1^st and 2^nd generation periodic table descriptors developed by us (without any sophisticated computing facilities) along with experimentally analyzed Zeta-potential, to model the cytotoxicity of these MONPs. Based on qualitative validation metrics, the LDA model appeared to be the best among the 7 tested models. The core environment of metal defined by the ratio of the number of core electrons to the number of valence electrons and the electronegativity count of oxygen showed a positive impact on toxicity. The identified properties were important for understanding the mechanisms of nanotoxicity and for predicting the potential environmental risk associated with MONPs exposure. The developed models can be utilized for environmental risk assessment of any untested MONP to E. coli, thereby providing a scientific basis for the design and preparation of safe nanomaterials.

Sugars and Sweeteners: Structure, Properties and In Silico Modeling

Several studies report the effects of excessive use of sugars and sweeteners in the diet. These include obesity, cardiac diseases, diabetes, and even lymphomas, leukemias, cancers of the bladder and brain, chronic fatigue syndrome, Parkinson's disease, Alzheimer's disease, multiple sclerosis, autism, and systemic lupus. On the other hand, each sugar and sweetener has a distinct metabolic assimilation process, and its chemical structure plays an important role in this process. Several scientific papers present the biological effects of the sugars and sweeteners in relation to their chemical structure. One important issue dealing with the sugars is the degree of similarity in their structures, focusing mostly on optical isomerism. Finding and developing new sugars and sweeteners with desired properties is an emerging research area, in which in silico approaches play an important role.

Prediction on the mutagenicity of nitroaromatic compounds using quantum chemistry descriptors based QSAR and machine learning derived classification methods

Nitroaromatic compounds (NACs) are an important type of environmental organic pollutants. However, it is lack of sufficient information relating to their potential adverse effects on human health and the environment due to the limited resources. Thus, using in silico technologies to assess their potential hazardous effects is urgent and promising. In this study, quantitative structure activity relationship (QSAR) and classification models were constructed using a set of NACs based on their mutagenicity against Salmonella typhimurium TA100 strain. For QSAR studies, DRAGON descriptors together with quantum chemistry descriptors were calculated for characterizing the detailed molecular information. Based on genetic algorithm (GA) and multiple linear regression (MLR) analyses, we screened descriptors and developed QSAR models. For classification studies, seven machine learning methods along with six molecular fingerprints were applied to develop qualitative classification models. The goodness of fitting, reliability, robustness and predictive performance of all developed models were measured by rigorous statistical validation criteria, then the best QSAR and classification models were chosen. Moreover, the QSAR models with quantum chemistry descriptors were compared to that without quantum chemistry descriptors and previously reported models. Notably, we also obtained some specific molecular properties or privileged substructures responsible for the high mutagenicity of NACs. Overall, the developed QSAR and classification models can be utilized as potential tools for rapidly predicting the mutagenicity of new or untested NACs for environmental hazard assessment and regulatory purposes, and may provide insights into the in vivo toxicity mechanisms of NACs and related compounds.

QSAR classification model for diverse series of antifungal agents based on improved binary differential search algorithm

An improved binary differential search (improved BDS) algorithm is proposed for QSAR classification of diverse series of antimicrobial compounds against Candida albicans inhibitors. The transfer functions is the most important component of the BDS algorithm, and converts continuous values of the donor into discrete values. In this paper, the eight types of transfer functions are investigated to verify their efficiency in improving BDS algorithm performance in QSAR classification. The performance was evaluated using three metrics: classification accuracy (CA), geometric mean of sensitivity and specificity (G-mean), and area under the curve. The Kruskal-Wallis test was also applied to show the statistical differences between the functions. Two functions, S1 and V4, show the best classification achievement, with a slightly better performance of V4 than S1. The V4 function takes the lowest iterations and selects the fewest descriptors. In addition, the V4 function yields the best CA and G-mean of 98.07% and 0.977%, respectively. The results prove that the V4 transfer function significantly improves the performance of the original BDS.

In Silico Prediction of O⁶-Methylguanine-DNA Methyltransferase Inhibitory Potency of Base Analogs with QSAR and Machine Learning Methods

O⁶-methylguanine-DNA methyltransferase (MGMT), a unique DNA repair enzyme, can confer resistance to DNA anticancer alkylating agents that modify the O⁶-position of guanine. Thus, inhibition of MGMT activity in tumors has a great interest for cancer researchers because it can significantly improve the anticancer efficacy of such alkylating agents. In this study, we performed a quantitative structure activity relationship (QSAR) and classification study based on a total of 134 base analogs related to their ED50 values (50% inhibitory concentration) against MGMT. Molecular information of all compounds were described by quantum chemical descriptors and Dragon descriptors. Genetic algorithm (GA) and multiple linear regression (MLR) analysis were combined to develop QSAR models. Classification models were generated by seven machine-learning methods based on six types of molecular fingerprints. Performances of all developed models were assessed by internal and external validation techniques. The best QSAR model was obtained with Q²Loo = 0.83, R² = 0.87, Q²ext = 0.67, and R²ext = 0.69 based on 84 compounds. The results from QSAR studies indicated topological charge indices, polarizability, ionization potential (IP), and number of primary aromatic amines are main contributors for MGMT inhibition of base analogs. For classification studies, the accuracies of 10-fold cross-validation ranged from 0.750 to 0.885 for top ten models. The range of accuracy for the external test set ranged from 0.800 to 0.880 except for PubChem-Tree model, suggesting a satisfactory predictive ability. Three models (Ext-SVM, Ext-Tree and Graph-RF) showed high and reliable predictive accuracy for both training and external test sets. In addition, several representative substructures for characterizing MGMT inhibitors were identified by information gain and substructure frequency analysis method. Our studies might be useful for further study to design and rapidly identify potential MGMT inhibitors.

QSAR and Classification Study on Prediction of Acute Oral Toxicity of N-Nitroso Compounds

To better understand the mechanism of in vivo toxicity of N-nitroso compounds (NNCs), the toxicity data of 80 NNCs related to their rat acute oral toxicity data (50% lethal dose concentration, LD50) were used to establish quantitative structure-activity relationship (QSAR) and classification models. Quantum chemistry methods calculated descriptors and Dragon descriptors were combined to describe the molecular information of all compounds. Genetic algorithm (GA) and multiple linear regression (MLR) analyses were combined to develop QSAR models. Fingerprints and machine learning methods were used to establish classification models. The quality and predictive performance of all established models were evaluated by internal and external validation techniques. The best GA-MLR-based QSAR model containing eight molecular descriptors was obtained with Q²loo = 0.7533, R² = 0.8071, Q²ext = 0.7041 and R²ext = 0.7195. The results derived from QSAR studies showed that the acute oral toxicity of NNCs mainly depends on three factors, namely, the polarizability, the ionization potential (IP) and the presence/absence and frequency of C⁻O bond. For classification studies, the best model was obtained using the MACCS keys fingerprint combined with artificial neural network (ANN) algorithm. The classification models suggested that several representative substructures, including nitrile, hetero N nonbasic, alkylchloride and amine-containing fragments are main contributors for the high toxicity of NNCs. Overall, the developed QSAR and classification models of the rat acute oral toxicity of NNCs showed satisfying predictive abilities. The results provide an insight into the understanding of the toxicity mechanism of NNCs in vivo, which might be used for a preliminary assessment of NNCs toxicity to mammals.

First Report on the Validated Classification-Based Chemometric Modeling of Human Rhinovirus 3C Protease (HRV 3Cpro) Inhibitors

Human rhinoviruses (HRVs), a major cause of common cold and upper respiratory infections, may trigger severe respiratory complications like asthma and COPD. To date, no drugs are available in the market which are designed as novel HRV inhibitors despite the involvement of some pharmaceutical companies' due to economical and clinical constraints. HRV 3C protease may be a potential target for drug design as it plays crucial role in viral RNA replication and virion assembly process. Therefore, designing novel HRV 3Cpro inhibitors is necessary and demanding in the field of antiviral drug design. In this article, statistically significant and validated classification-based QSARs of a series of HRV 3Cpro inhibitors were performed for the first time as per the authors' knowledge. Results suggest that oxopyrrolidine and piperidinone rings are favored whereas carboxybenzyl and unsubstituted benzyl functions may be unfavorable. Moreover, this group, along with cyclic alkyl or aryl ring structures may favor HRV 3Cpro inhibition. These observations may be utilized for the design of a higher active anti-HRV agent in future.

Longitudinal change in the serology of antibodies to Chlamydia trachomatis pgp3 in children residing in a trachoma area

A serologic test for antibodies to chlamydial antigen pgp3 may be a useful tool for trachoma surveillance. However, little is known about the stability of antibody status over time, or factors associated with seroreversion/conversion. A cohort of 2,111 children ages 1-9 years in Tanzania were followed for one year in the absence of mass azithromycin. At baseline and follow-up, they were evaluated for trachoma, chlamydial infection, and antibodies to chlamydial antigen pgp3. At baseline, 31% of children were seropositive for pgp3 antibodies and 6.4% seroreverted to negative over one year. Of those seronegative, 9.8% seroconverted over the year. The seroreverters had lower baseline mean fluorescence intensity (MFI-BG) values compared to the seropositives who remained positive (Odds Ratio = 0.04 for every unit increase in log₁₀MFI-BG, 95% CI = 0.02-0.09), and were more likely to live in communities with trachoma <5% (p < 0.008). While seroconversion was expected, seroreversion was unexpected. The low seroprevalence rate reported from low endemic areas may be due to seroreversion as well as lack of exposure.

Multiple molecular modelling studies on some derivatives and analogues of glutamic acid as matrix metalloproteinase-2 inhibitors

Matrix metalloproteinase-2 (MMP-2) is a potential target in anticancer drug discovery due to its association with angiogenesis, metastasis and tumour progression. In this study, 67 glutamic acid derivatives, synthesized and evaluated as MMP-2 inhibitors, were taken into account for multi-QSAR modelling study (regression-based 2D-QSAR, classification-based LDA-QSAR, Bayesian classification QSAR, HQSAR, 3D-QSAR CoMFA and CoMSIA as well as Open3DQSAR). All these QSAR studies were statistically validated individually. Regarding the 3D-QSAR analysis, the Open3DQSAR results were better than CoMFA and CoMSIA, although all these 3D-QSAR models supported each other. The importance of biphenylsulphonyl moiety over phenylacetyl/naphthylacetyl moieties was established due to its association with favourable steric and hydrophobic characters. HQSAR, LDA-QSAR and Bayesian classification QSAR studies also suggested that the biphenylsulphonamido group was better than the phenylacetylcarboxamido function. Additionally, glutamines were proven to be far better inhibitors than isoglutamines. Observations obtained from the current study were revalidated and supported by the earlier reported molecular modelling studies. Depending on these observations, newer glutamic acid-based compounds may be designed further in future for potent MMP-2 inhibitory activity.

Three-Dimensional Biologically Relevant Spectrum (BRS-3D): Shape Similarity Profile Based on PDB Ligands as Molecular Descriptors

The crystallized ligands in the Protein Data Bank (PDB) can be treated as the inverse shapes of the active sites of corresponding proteins. Therefore, the shape similarity between a molecule and PDB ligands indicated the possibility of the molecule to bind with the targets. In this paper, we proposed a shape similarity profile that can be used as a molecular descriptor for ligand-based virtual screening. First, through three-dimensional (3D) structural clustering, 300 diverse ligands were extracted from the druggable protein-ligand database, sc-PDB. Then, each of the molecules under scrutiny was flexibly superimposed onto the 300 ligands. Superimpositions were scored by shape overlap and property similarity, producing a 300 dimensional similarity array termed the "Three-Dimensional Biologically Relevant Spectrum (BRS-3D)". Finally, quantitative or discriminant models were developed with the 300 dimensional descriptor using machine learning methods (support vector machine). The effectiveness of this approach was evaluated using 42 benchmark data sets from the G protein-coupled receptor (GPCR) ligand library and the GPCR decoy database (GLL/GDD). We compared the performance of BRS-3D with other 2D and 3D state-of-the-art molecular descriptors. The results showed that models built with BRS-3D performed best for most GLL/GDD data sets. We also applied BRS-3D in histone deacetylase 1 inhibitors screening and GPCR subtype selectivity prediction. The advantages and disadvantages of this approach are discussed.

In silico Prediction of Drug Induced Liver Toxicity Using Substructure Pattern Recognition Method

Drug-induced liver injury (DILI) is a leading cause of acute liver failure in the US and less severe liver injury worldwide. It is also one of the major reasons of drug withdrawal from the market. Thus, DILI has become one of the most important concerns of drugs, and should be predicted in very early stage of drug discovery process. In this study, a comprehensive data set containing 1317 diverse compounds was collected from publications. Then, high accuracy classification models were built using five machine learning methods based on MACCS and FP4 fingerprints after evaluating by substructure pattern recognition method. The best model was built using SVM method together with FP4 fingerprint at the IG value threshold of 0.0005. Its overall predictive accuracies were 79.7 % and 64.5 % for the training and test sets, separately, which yielded overall accuracy of 75.0 % for the external validation dataset, consisting of 88 compounds collected from a benchmark DILI database - the Liver Toxicity Knowledge Base. This model could be used for drug-induced liver toxicity prediction. Moreover, some key substructure patterns correlated with drug-induced liver toxicity were also identified as structural alerts.

Exploring simple, transparent, interpretable and predictive QSAR models for classification and quantitative prediction of rat toxicity of ionic liquids using OECD recommended guidelines

The present study explores the chemical attributes of diverse ionic liquids responsible for their cytotoxicity in a rat leukemia cell line (IPC-81) by developing predictive classification as well as regression-based mathematical models. Simple and interpretable descriptors derived from a two-dimensional representation of the chemical structures along with quantum topological molecular similarity indices have been used for model development, employing unambiguous modeling strategies that strictly obey the guidelines of the Organization for Economic Co-operation and Development (OECD) for quantitative structure-activity relationship (QSAR) analysis. The structure-toxicity relationships that emerged from both classification and regression-based models were in accordance with the findings of some previous studies. The models suggested that the cytotoxicity of ionic liquids is dependent on the cationic surfactant action, long alkyl side chains, cationic lipophilicity as well as aromaticity, the presence of a dialkylamino substituent at the 4-position of the pyridinium nucleus and a bulky anionic moiety. The models have been transparently presented in the form of equations, thus allowing their easy transferability in accordance with the OECD guidelines. The models have also been subjected to rigorous validation tests proving their predictive potential and can hence be used for designing novel and "greener" ionic liquids. The major strength of the present study lies in the use of a diverse and large dataset, use of simple reproducible descriptors and compliance with the OECD norms.

In silico prediction of chemical toxicity on avian species using chemical category approaches

Avian species are sensitive to pesticides and industrial chemicals, and hence used as model species in evaluation of chemical toxicity. In present study, we assessed the toxicity of more than 663 diverse chemicals on 17 avian species. All the chemicals were classified into three categories, i.e. highly toxic, slightly toxic and non-toxic, based on the toxicity classification criteria of the United States Environmental Protection Agency (EPA). To evaluate these chemicals, the toxicity prediction models were built using chemical category approaches with molecular descriptors and five commonly used fingerprints, in which five machine learning methods were performed on two standard test species: aquatic bird mallard duck and terrestrial bird northern bobwhite quail. The support vector machine (SVM) method with Pubchem fingerprint performed best as revealed by 5-fold cross-validation and the external validation set on Japanese quail. No species difference existed in our database despite several chemicals with different toxicity on some avian species. The best model had an overall accuracy at 0.851 for the prediction of toxicity on avian species, which outperformed the work of Mazzatorta et al. Furthermore, several representative substructures for characterizing avian toxicity were identified via information gain (IG) method. This study would provide a new tool for chemical safety assessment.

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.

Application of Receiver Operating Characteristic (ROC) Curves for Explosives Detection Using Different Sampling and Detection Techniques

Reported for the first time are receiver operating characteristic (ROC) curves constructed to describe the performance of a sorbent-coated disk, planar solid phase microextraction (PSPME) unit for non-contact sampling of a variety of volatiles. The PSPME is coupled to ion mobility spectrometers (IMSs) for the detection of volatile chemical markers associated with the presence of smokeless powders, model systems of explosives containing diphenylamine (DPA), 2,4-dinitrotoluene (2,4-DNT) and nitroglycerin (NG) as the target analytes. The performance of the PSPME-IMS was compared with the widely accepted solid-phase microextraction (SPME), coupled to a GC-MS. A set of optimized sampling conditions for different volume containers (1–45 L) with various sample amounts of explosives, were studied in replicates (n = 30) to determine the true positive rates (TPR) and false positive detection rates (FPR) for the different scenarios. These studies were obtained in order to construct the ROC curves for two IMS instruments (a bench-top and field-portable system) and a bench top GC-MS system in low and high clutter environments. Both static and dynamic PSPME sampling were studied in which 10–500 mg quantities of smokeless powders were detected within 10 min of static sampling and 1 min of dynamic sampling.

Topological structural alerts modulations of mammalian cell mutagenicity for halogenated derivatives

Genotoxicity is a key toxicity endpoint for current regulatory requirements regarding new and existing chemicals. However, genotoxicity testing is time-consuming and costly, and involves the use of laboratory animals. This has motivated the development of computational approaches, designed to predict genotoxicity without the need to conduct laboratory tests. Currently, many existing computational methods, like quantitative structure-activity relationship (QSAR) models, provide limited information about the possible mechanisms involved in mutagenicity or predictions based on structural alerts (SAs) do not take statistical models into account. This paper describes an attempt to address this problem by using the TOPological Substructural MOlecular Design (TOPS-MODE) approach to develop and validate improved QSAR models for predicting the mutagenicity of a range of halogenated derivatives. Our most predictive model has an accuracy of 94.12%, exhibits excellent cross-validation and external set statistics. A reasonable interpretation of the model in term of SAs was achieved by means of bond contributions to activity. The results obtained led to the following conclusions: primary halogenated derivatives are more mutagenic than secondary ones; and substitution of chlorine by bromine increases mutagenicity while polyhalogenation decreases activity. The paper demonstrates the potential of the TOPS-MODE approach in developing QSAR models for identifying structural alerts for mutagenicity, combining high predictivity with relevant mechanistic interpretation.

The use of ROC analysis for the qualitative prediction of human oral bioavailability from animal data

Purpose

To develop and evaluate a tool for the qualitative prediction of human oral bioavailability (F_human) from animal oral bioavailability (F_animal) data employing ROC analysis and to identify the optimal thresholds for such predictions.

Methods

A dataset of 184 compounds with known F_human and F_animal in at least one species (mouse, rat, dog and non-human primates (NHP)) was employed. A binary classification model for F_human was built by setting a threshold for high/low F_human at 50%. The thresholds for high/low F_animal were varied from 0 to 100 to generate the ROC curves. Optimal thresholds were derived from ‘cost analysis’ and the outcomes with respect to false negative and false positive predictions were analyzed against the BDDCS class distributions.

Results

We successfully built ROC curves for the combined dataset and per individual species. Optimal F_animal thresholds were found to be 67% (mouse), 22% (rat), 58% (dog), 35% (NHP) and 47% (combined dataset). No significant trends were observed when sub-categorizing the outcomes by the BDDCS.

Conclusions

F_animal can predict high/low F_human with adequate sensitivity and specificity. This methodology and associated thresholds can be employed as part of decisions related to planning necessary studies during development of new drug candidates and lead selection.

Electronic supplementary material

The online version of this article (doi:10.1007/s11095-013-1193-2) contains supplementary material, which is available to authorized users.

Development and validation of regression-based QSAR models for quantification of contributions of molecular fragments to skin sensitization potency of diverse organic chemicals

In our present work, we have developed regression-based QSAR models for skin sensitization potential of 51 diverse organic chemicals. The objective behind the present work is to determine the influence of different molecular features on the skin sensitizing potential of chemicals. Among several models developed, the two best ones are discussed to unveil specific information regarding the contribution of different structural and physicochemical features towards the property of skin sensitization. The QSAR models suggested that aromatic compounds are more skin sensitizing than aliphatic ones, but in the case of carbonyl compounds, aliphatic ones are more skin sensitizing than aromatic ones. Other descriptors such as LUMO and <2-Atype_H-47> signify the importance of the electrophilic and hydrophilic character respectively of the molecules for showing skin sensitizing property. Another two descriptors, <Dipole-mag-2.72> and (3)χC also exert significant contributions towards the skin sensitization potential of the chemicals. Further, it is observed that the nitrogen atoms (nN), triple bonds (nTB) and also the fragment Al-C(=X)-Al (Atype_C38) are responsible for skin sensitizing property. All the above information provides additional support for further research involving identification of the skin sensitization potential of new chemicals.

QSTR with extended topochemical atom (ETA) indices. 16. Development of predictive classification and regression models for toxicity of ionic liquids towards Daphnia magna

Ionic liquids have been judged much with respect to their wide applicability than their considerable harmful effects towards the living ecosystem which has been observed in many instances. Hence, toxicological introspection of these chemicals by the development of predictive mathematical models can be of good help. This study presents an attempt to develop predictive classification and regression models correlating the structurally derived chemical information of a group of 62 diverse ionic liquids with their toxicity towards Daphnia magna and their interpretation. We have principally used the extended topochemical atom (ETA) indices along with various topological non-ETA and thermodynamic parameters as independent variables. The developed quantitative models have been subjected to extensive statistical tests employing multiple validation strategies from which acceptable results have been reported. The best models obtained from classification and regression studies captured necessary structural information on lipophilicity, branching pattern, electronegativity and chain length of the cationic substituents for explaining ecotoxicity of ionic liquids towards D. magna. The derived information can be successfully used to design better ionic liquid analogues acquiring the qualities of a true eco-friendly green chemical.

Combining QSAR classification models for predictive modeling of human monoamine oxidase inhibitors

Due to their role in the metabolism of monoamine neurotransmitters, MAO-A and MAO-B present a significant pharmacological interest. For instance the inhibitors of human MAO-B are considered useful tools for the treatment of Parkinson Disease. Therefore, the rational design and synthesis of new MAOs inhibitors is considered of great importance for the development of new and more effective treatments of Parkinson Disease. In this work, Quantitative Structure Activity Relationships (QSAR) has been developed to predict the human MAO inhibitory activity and selectivity. The first step was the selection of a suitable dataset of heterocyclic compounds that include chromones, coumarins, chalcones, thiazolylhydrazones, etc. These compounds were previously synthesized in one of our laboratories, or elsewhere, and their activities measured by the same assays and for the same laboratory staff. Applying linear discriminant analysis to data derived from a variety of molecular representations and feature selection algorithms, reliable QSAR models were built which could be used to predict for test compounds the inhibitory activity and selectivity toward human MAO. This work also showed how several QSAR models can be combined to make better predictions. The final models exhibit significant statistics, interpretability, as well as displaying predictive power on an external validation set made up of chromone derivatives with unknown activity (that are being reported here for first time) synthesized by our group, and coumarins recently reported in the literature.

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

Carcinogenicity is among the toxicological endpoints posing the highest concern for human health. Oral slope factors (OSFs) are used to estimate quantitatively the carcinogenic potency or the risk associated with exposure to the chemical by oral route. Regulatory agencies in food and drug administration and environmental protection are employing quantitative structure-activity relationship (QSAR) models to fill the data gaps related with properties of chemicals affecting the environment and human health. In this background, we have developed quantitative structure-carcinogenicity regression models for rodents based on the carcinogenic potential of 70 chemicals with wide diversity of molecular structures, spanning a large number of chemical classes and biological mechanisms. All the developed models have been assessed according to the Organization for Economic Cooperation and Development (OECD) principles for the validation of QSAR models. We have also attempted to develop a carcinogenicity classification model based on Linear Discriminant Analysis (LDA). Developed regression and LDA models are rigorously validated internally as well as externally. Our in silico studies make it possible to obtain a quantitative interpretation of the structural information of carcinogenicity along with identification of the discriminant functions between lower and higher carcinogenic compounds by LDA. Pharmacological distribution diagrams (PDDs) are used as a visualizing technique for the identification and selection of chemicals with lower carcinogenicity. Constructive, informative and comparable interpretations have been observed in both cases of classification and regression based modeling.