QSAR Prediction of Estrogen Activity for a Large Set of Diverse Chemicals under the Guidance of OECD Principles

Machine learning for predicting halogen radical reactivity toward aqueous organic chemicals

Rapid advances in machine learning (ML) provide fast, accurate, and widely applicable methods for predicting free radical-mediated organic pollutant reactivity. In this study, the rate constants (logk) of four halogen radicals were predicted using Morgan fingerprint (MF) and Mordred descriptor (MD) in combination with a series of ML models. The findings highlighted that making accurate predictions for various datasets depended on an effective combination of descriptors and algorithms. To further alleviate the challenge of limited sample size, we introduced a data combination strategy that improved prediction accuracy and mitigated overfitting by combining different datasets. The Light Gradient Boosting Machine (LightGBM) with MF and Random Forest (RF) with MD models based on the unified dataset were finally selected as the optimal models. The SHapley Additive exPlanations revealed insights: the MF-LightGBM model successfully captured the influence of electron-withdrawing/donating groups, while autocorrelation, walk count and information content descriptors in the MD-RF model were identified as key features. Furthermore, the important contribution of pH was emphasized. The results of the applicability domain analysis further supported that the developed model can make reliable predictions for query compounds across a broader range. Finally, a practical web application for logk calculations was built.

Estrogen Receptor Alpha Binders for Hormone-Dependent Forms of Breast Cancer: e-QSAR and Molecular Docking Supported by X-ray Resolved Structures

Cancer, a life-disturbing and lethal disease with a high global impact, causes significant economic, social, and health challenges. Breast cancer refers to the abnormal growth of cells originating from breast tissues. Hormone-dependent forms of breast cancer, such as those influenced by estrogen, prompt the exploration of estrogen receptors as targets for potential therapeutic interventions. In this study, we conducted e-QSAR molecular docking and molecular dynamics analyses on a diverse set of inhibitors targeting estrogen receptor alpha (ER-α). The e-QSAR model is based on a genetic algorithm combined with multilinear regression analysis. The newly developed model possesses a balance between predictive accuracy and mechanistic insights adhering to the OECD guidelines. The e-QSAR model pointed out that sp2-hybridized carbon and nitrogen atoms are important atoms governing binding profiles. In addition, a specific combination of H-bond donors and acceptors with carbon, nitrogen, and ring sulfur atoms also plays a crucial role. The results are supported by molecular docking, MD simulations, and X-ray-resolved structures. The novel results could be useful for future drug development for ER-α.

From Molecular Descriptors to Intrinsic Fish Toxicity of Chemicals: An Alternative Approach to Chemical Prioritization

The European and U.S. chemical agencies have listed approximately 800k chemicals about which knowledge of potential risks to human health and the environment is lacking. Filling these data gaps experimentally is impossible, so in silico approaches and prediction are essential. Many existing models are however limited by assumptions (e.g., linearity and continuity) and small training sets. In this study, we present a supervised direct classification model that connects molecular descriptors to toxicity. Categories can be driven by either data (using k-means clustering) or defined by regulation. This was tested via 907 experimentally defined 96 h LC50 values for acute fish toxicity. Our classification model explained ≈90% of the variance in our data for the training set and ≈80% for the test set. This strategy gave a 5-fold decrease in the frequency of incorrect categorization compared to a quantitative structure-activity relationship (QSAR) regression model. Our model was subsequently employed to predict the toxicity categories of ≈32k chemicals. A comparison between the model-based applicability domain (AD) and the training set AD was performed, suggesting that the training set-based AD is a more adequate way to avoid extrapolation when using such models. The better performance of our direct classification model compared to that of QSAR methods makes this approach a viable tool for assessing the hazards and risks of chemicals.

Targeting MDM2-p53 Axis through Drug Repurposing for Cancer Therapy: A Multidisciplinary Approach

Cancer remains a major cause of morbidity and mortality worldwide, and while current therapies, such as chemotherapy, immunotherapy, and cell therapy, have been effective in many patients, the development of novel therapeutic options remains an urgent priority. Mouse double minute 2 (MDM2) is a key regulator of the tumor suppressor protein p53, which plays a critical role in regulating cellular growth, apoptosis, and DNA repair. Consequently, MDM2 has been the subject of extensive research aimed at developing novel cancer therapies. In this study, we employed a machine learning-based approach to establish a quantitative structure-activity relationship model capable of predicting the potential in vitro efficacy of small molecules as MDM2 inhibitors. Our model was used to screen 5883 FDA-approved drugs, resulting in the identification of promising hits that were subsequently evaluated using molecular docking and molecular dynamics simulations. Two antihistamine drugs, cetirizine (CZ) and rupatadine (RP), exhibited particularly favorable results in the initial in silico analyses. To further assess their potential use as the activators of the p53 pathway, we investigated the antiproliferative capability of the abovementioned drugs on human glioblastoma and neuroblastoma cell lines. Both the compounds exhibited significant antiproliferative effects on the abovementioned cell lines in a dose-dependent manner. The half-maximal inhibitory concentration (IC50) of CZ was found to be 697.87 and 941.37 μM on U87 and SH-SY5Y cell lines, respectively, while the IC50 of RP was found to be 524.28 and 617.07 μM on the same cell lines, respectively. Further investigation by quantitative reverse transcriptase polymerase chain reaction analysis revealed that the CZ-treated cell lines upregulate the expression of the p53-regulated genes involved in cell cycle arrest, apoptosis, and DNA damage response compared to their respective vehicle controls. These findings suggest that CZ activates the p53 pathway by inhibiting MDM2. Our results provide compelling preclinical evidence supporting the potential use of CZ as a modulator of the MDM2-p53 axis and its plausible repurposing for cancer treatment.

New Insights on the Activity and Selectivity of MAO-B Inhibitors through In Silico Methods

To facilitate the identification of novel MAO-B inhibitors, we elaborated a consolidated computational approach, including a pharmacophoric atom-based 3D quantitative structure-activity relationship (QSAR) model, activity cliffs, fingerprint, and molecular docking analysis on a dataset of 126 molecules. An AAHR.2 hypothesis with two hydrogen bond acceptors (A), one hydrophobic (H), and one aromatic ring (R) supplied a statistically significant 3D QSAR model reflected by the parameters: R² = 0.900 (training set); Q² = 0.774 and Pearson's R = 0.884 (test set), stability s = 0.736. Hydrophobic and electron-withdrawing fields portrayed the relationships between structural characteristics and inhibitory activity. The quinolin-2-one scaffold has a key role in selectivity towards MAO-B with an AUC of 0.962, as retrieved by ECFP4 analysis. Two activity cliffs showing meaningful potency variation in the MAO-B chemical space were observed. The docking study revealed interactions with crucial residues TYR:435, TYR:326, CYS:172, and GLN:206 responsible for MAO-B activity. Molecular docking is in consensus with and complementary to pharmacophoric 3D QSAR, ECFP4, and MM-GBSA analysis. The computational scenario provided here will assist chemists in quickly designing and predicting new potent and selective candidates as MAO-B inhibitors for MAO-B-driven diseases. This approach can also be used to identify MAO-B inhibitors from other libraries or screen top molecules for other targets involved in suitable diseases.

Estrogenic flavonoids and their molecular mechanisms of action

Flavonoids are a major group of phytoestrogens associated with physiological effects, and ecological and social impacts. Although the estrogenic activity of flavonoids was reported by researchers in the fields of medical, environmental and food studies, their molecular mechanisms of action have not been comprehensively reviewed. The estrogenic activity of the respective classes of flavonoids, anthocyanidins/anthocyanins, 2-arylbenzofurans/3-arylcoumarins/α-methyldeoxybenzoins, aurones/chalcones/dihydrochalcones, coumaronochromones, coumestans, flavans/flavan-3-ols/flavan-4-ols, flavanones/dihydroflavonols, flavones/flavonols, homoisoflavonoids, isoflavans, isoflavanones, isoflavenes, isoflavones, neoflavonoids, oligoflavonoids, pterocarpans/pterocarpenes, and rotenone/rotenoids, was summarized through a comprehensive literature search, and their structure-activity relationship, biological activities, signaling pathways, and applications were discussed. Although the respective classes of flavonoids contained at least one chemical mimicking estrogen, the mechanisms varied, such as those with estrogenic, anti-estrogenic, non-estrogenic, and biphasic activities, and additional activities through crosstalk/bypassing, which exert biological activities through cell signaling pathways. Such mechanistic variations of estrogen action are not limited to flavonoids and are observed among other broad categories of chemicals, thus this group of chemicals can be termed as the "estrogenome". This review article focuses on the connection of estrogen action mainly between the outer and the inner environments, which represent variations of chemicals and biological activities/signaling pathways, respectively, and form the basis to understand their applications. The applications of chemicals will markedly progress due to emerging technologies, such as artificial intelligence for precision medicine, which is also true of the study of the estrogenome including estrogenic flavonoids.

Molecular breast cancer subtype identification using photoacoustic spectral analysis and machine learning at the biomacromolecular level

Breast cancer threatens the health of women worldwide, and its molecular subtypes largely determine the therapy and prognosis of patients. However, an uncomplicated and accurate method to identify subtypes is currently lacking. This study utilized photoacoustic spectral analysis (PASA) based on the partial least squares discriminant algorithm (PLS-DA) to identify molecular breast cancer subtypes at the biomacromolecular level in vivo. The area of power spectrum density (APSD) was extracted to semi-quantify the biomacromolecule content. The feature wavelengths were obtained via the variable importance in projection (VIP) score and the selectivity ratio (Sratio), to identify the biomarkers. The PASA achieved an accuracy of 84%. Most of the feature wavelengths fell into the collagen-dominated absorption waveband, which was consistent with the histopathological results. This paper proposes a successful method for identifying molecular breast cancer subtypes and proves that collagen can be treated as a biomarker for molecular breast cancer subtyping.

Modelling PIP4K2A inhibitory activity of 1,7-naphthyridine analogues using machine learning and molecular docking studies

PIP4K2A is a type II lipid kinase that catalyzed the rate-limiting step of the conversion of phosphatidylinositol-5-phosphate (PI5P) into phosphatidylinositol 4,5-bisphosphate (PI4,5P2). PIP4K2A has been intricately linked to the inhibition of various types of tumors via reactive oxygen species-mediated apoptosis, making it an important therapeutic target. In the quest of finding biologically active substances with efficient PIP4K2A inhibitory activity, machine learning algorithms were used to investigate the quantitative relationship between structures and inhibitory activities of 1,7-naphthyridine analogues. Three machine learning algorithms (MLR, ANN, and SVM) were used to develop QSAR models that can effectively predict the PIP4K2A inhibitory activity of a library of 1,7-naphthyridine analogues. The cascaded feature selection method was performed by sequential application of GFA and MP5 algorithms to identify a molecular descriptor subset that can best describe the PIP4K2A inhibitory activity of 1,7-naphthyridine analogues. PIP4K2A inhibitory activities predicted by the ML models were strongly correlated with the experimental values. The QSAR Modelling indicates that the best-performing ML model was SVM with the RBF kernel function. The SVM model performed very well in predicting PIP4K2A inhibitory activity of the 1,7-naphthyridine analogues with RTR and QEX values of 0.9845 and 0.8793 respectively. To further gain more structural insight into the origin of PIP4K2A inhibitory activity of 1,7-naphthyridine analogues, molecular docking studies were performed. The results indicate that five compounds; 15, 25, 13, 09, and 28 were found to have a high binding affinity with the receptor molecules. Hydrogen bonding, pi-pi interaction, and pi-cation interactions were found to modulate the binding interaction of the inhibitors. Although the SVM gives essentially a black-box model which cannot be readily interpreted, using SVM in tandem with MLR and ANN provides a unique perspective in building robust QSAR predictive models. The superior predictive performance of the ML models and the explanatory power of MLR models were combined to provide a unique insight into the structure-activity relationship of 1,7-naphthyridine inhibitors. This is relevant in that it provides information that can be invaluable as guidelines for the design of novel PIP4K2A inhibitors.

Versatile in silico modeling of XAD-air partition coefficients for POPs based on abraham descriptor and temperature

The concentration of persistent organic pollutants (POPs) makes remarkable difference to environmental fate. In the field of passive sampling, the partition coefficients between polystyrene-divinylbenzene resin (XAD) and air (i.e., K_XAD-A) are indispensable to obtain POPs concentration, and the K_XAD-A is generally thought to be governed by temperature and molecular structure of POPs. However, experimental determination of K_XAD-A is unrealistic for countless and novel chemicals. Herein, the Abraham solute descriptors of poly parameter linear free energy relationship (pp-LFER) and temperature were utilized to develop models, namely pp-LFER-T, for predicting K_XAD-A values. Two linear (MLR and LASSO) and four nonlinear (ANN, SVM, kNN and RF) machine learning algorithms were employed to develop models based on a data set of 307 sample points. For the aforementioned six models, R²_adj and Q²_ext were both beyond 0.90, indicating distinguished goodness-of-fit and robust generalization ability. By comparing the established models, the best model was observed as the RF model with R²_adj = 0.991, Q²_ext = 0.935, RMSE_tra = 0.271 and RMSE_ext = 0.868. The mechanism interpretation revealed that the temperature, size of molecules and dipole-type interactions were the predominant factors affecting K_XAD-A values. Concurrently, the developed models with the broad applicability domain provide available tools to fill the experimental data gap for untested chemicals. In addition, the developed models were helpful to preliminarily evaluate the environmental ecological risk and understand the adsorption behavior of POPs between XAD membrane and air.

Alignment-Free Antimicrobial Peptide Predictors: Improving Performance by a Thorough Analysis of the Largest Available Data Set

In the last two decades, a large number of machine-learning-based predictors for the activities of antimicrobial peptides (AMPs) have been proposed. These predictors differ from one another in the learning method and in the training and testing data sets used. Unfortunately, the training data sets present several drawbacks, such as a low representativeness regarding the experimentally validated AMP space, and duplicated peptide sequences between negative and positive data sets. These limitations give a low confidence to most of the approaches to be used in prospective studies. To address these weaknesses, we propose novel modeling and assessing data sets from the largest experimentally validated nonredundant peptide data set reported to date. From these novel data sets, alignment-free quantitative sequence-activity models (AF-QSAMs) based on Random Forest are created to identify general AMPs and their antibacterial, antifungal, antiparasitic, and antiviral functional types. An applicability domain analysis is carried out to determine the reliability of the predictions obtained, which, to the best of our knowledge, is performed for the first time for AMP recognition. A benchmarking is undertaken between the models proposed and several models from the literature that are freely available in 13 programs (ClassAMP, iAMP-2L, ADAM, MLAMP, AMPScanner v2.0, AntiFP, AMPfun, PEPred-suite, AxPEP, CAMPR3, iAMPpred, APIN, and Meta-iAVP). The models proposed are those with the best performance in all of the endpoints modeled, while most of the methods from the literature have weak-to-random predictive agreements. The models proposed are also assessed through Y-scrambling and repeated k-fold cross-validation tests, demonstrating that the outcomes obtained by them are not given by chance. Three chemometric analyses also confirmed the relevance of the peptides descriptors used in the modeling. Therefore, it can be concluded that the models built by fixing the drawbacks existing in the literature contribute to identifying antibacterial, antifungal, antiparasitic, and antiviral peptides with high effectivity and reliability. Models are freely available via the AMPDiscover tool at https://biocom-ampdiscover.cicese.mx/.

Integrating endocrine-related health effects into comparative human toxicity characterization

Endocrine-disrupting chemicals have the ability to interfere with and alter functions of the hormone system, leading to adverse effects on reproduction, growth and development. Despite growing concerns over their now ubiquitous presence in the environment, endocrine-related human health effects remain largely outside of comparative human toxicity characterization frameworks as applied for example in life cycle impact assessments. In this paper, we propose a new methodological framework to consistently integrate endocrine-related health effects into comparative human toxicity characterization. We present two quantitative and operational approaches for extrapolating towards a common point of departure from both in vivo and dosimetry-adjusted in vitro endocrine-related effect data and deriving effect factors as well as corresponding characterization factors for endocrine-active/endocrine-disrupting chemicals. Following the proposed approaches, we calculated effect factors for 323 chemicals, reflecting their endocrine potency, and related characterization factors for 157 chemicals, expressing their relative endocrine-related human toxicity potential. Developed effect and characterization factors are ready for use in the context of chemical prioritization and substitution as well as life cycle impact assessment and other comparative assessment frameworks. Endocrine-related effect factors were found comparable to existing effect factors for cancer and non-cancer effects, indicating that (1) the chemicals' endocrine potency is not necessarily higher or lower than other effect potencies and (2) using dosimetry-adjusted effect data to derive effect factors does not consistently overestimate the effect of potential endocrine disruptors. Calculated characterization factors span over 8-11 orders of magnitude for different substances and emission compartments and are dominated by the range in endocrine potencies.

Endocrine disruption: the noise in available data adversely impacts the models' performance

This paper is devoted to the analysis of available experimental data and preparation of predictive models for binding affinity of molecules with respect to two nuclear receptors involved in endocrine disruption (ED): the oestrogen (ER) and the androgen (AR) receptors. The ED-relevant data were retrieved from multiple sources, including the CERAPP, CoMPARA, and the Tox21 projects as well as ChEMBL and PubChem databases. Data analysis performed with the help of generative topographic mapping revealed the problem of low agreement between experimental values from different sources. Collected data were used to train both classification models for ER and AR binding activities and regression models for relative binding affinity (RBA) and median inhibition concentration (IC50). These models displayed relatively poor performance in classification (sensitivities ER = 0.34, AR = 0.49) and in regression (determination coefficient r ² for the RBA and IC50 models in external validation varied from 0.44 to 0.76). Our analysis demonstrates that low models' performance resulted from misinterpreted experimental endpoints or wrongly reported values, thus confirming the observations reported in CERAPP and CoMPARA studies. Developed models and collected data sets included of 6215 (ER) and 3789 (AR) unique compounds, which are freely available.

Binary and multi-class classification for androgen receptor agonists, antagonists and binders

Androgens and androgen receptor regulate a variety of biological effects in the human body. The impaired functioning of androgen receptor may have different adverse health effects from cancer to infertility. Therefore, it is important to determine whether new chemicals have any binding activity and act as androgen agonists or antagonists before commercial use. Due to the large number of chemicals that require experimental testing, the computational methods are a viable alternative. Therefore, the aim of the present study was to develop predictive QSAR models for classifying compounds according to their activity at the androgen receptor. A large data set of chemicals from the CoMPARA project was used for this purpose and random forest classification models have been developed for androgen binding, agonistic, and antagonistic activity. In addition, a unique effort has been made for multi-class approach that discriminates between inactive compounds, agonists and antagonists simultaneously. For the evaluation set, the classification models predicted agonists with 80% of accuracy and for the antagonists' and binders' the respective metrics were 72% and 78%. Combining agonists, antagonists and inactive compounds into a multi-class approach added complexity to the modelling task and resulted to 64% prediction accuracy for the evaluation set. Considering the size of the training data sets and their imbalance, the achieved evaluation accuracy is very good. The final classification models are available for exploring and predicting at QsarDB repository (https://doi.org/10.15152/QDB.236).

Principles of QSAR Modeling

At the end of her academic career, the author summarizes the main aspects of QSAR modeling, giving comments and suggestions according to her 23 years' experience in QSAR research on environmental topics. The focus is mainly on Multiple Linear Regression, particularly Ordinary Least Squares, using a Genetic Algorithm for variable selection from various theoretical molecular descriptors, but the comments can be useful also for other QSAR methods. The need for rigorous validation, also external, and for applicability domain check to guarantee predictivity and reliability of QSAR models is particularly highlighted. The commented approach is the “predictive” one, based on chemometrics, and is usefully applied to the prioritization of environmental pollutants. All the discussed points and the author's ideas are implemented in the software QSARINS, as a legacy to the QSAR community.

Application of quantitative structure-activity relationship (QSAR) model in comprehensive human health risk assessment of PAHs, and alkyl-, nitro-, carbonyl-, and hydroxyl-PAHs laden in urban road dust

The carcinogenic human health risks (CHHR) posed by the exposure to PAHs and transformed PAH products (TPPs) are currently inconclusive due to the lack of toxicity equivalency factors (TEFs) for most TPPs although some of these pollutants are more potent carcinogens. The applicability of quantitative structure-activity relationship (QSAR) model in predicting TEF of PAHs and TPPs to holistically evaluate the CHHR posed by the exposure to these pollutants in road dust from Gold Coast, Australia was examined. Statistical evaluation via ten metrics shows that partial least-squares regression (PLSR1) model has more statistical power in predicting TEF than multiple linear regression (MLR) within relevant applicability domain. For instance, the predicted residual sum of squares (PRESS) and standard deviation of error of prediction (SDEP) for PLSR is closer to zero than that of MLR. The total cancer risk estimated using the QSAR model derived TEFs and original TEFs for outliers gives a more holistic incremental lifetime cancer risk in relation to children and adults. Potential cancer risk exists for adults with this approach whereas reliance on only the originally available TEFs lead to a negligible risk diagnosis. The application of QSAR model in assessing CHHR due to PAHs and TPPs exposures is very viable.

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.

Cultivar Classification of Single Sweet Corn Seed Using Fourier Transform Near-Infrared Spectroscopy Combined with Discriminant Analysis

Seed purity is a key indicator of crop seed quality. The conventional methods for cultivar identification are time-consuming, expensive, and destructive. Fourier transform near-infrared (FT-NIR) spectroscopy combined with discriminant analyses, was studied as a rapid and nondestructive technique to classify the cultivars of sweet corn seeds. Spectra with a range of 1000–2500 nm collected from 760 seeds of two cultivars were used for the discriminant analyses. Thereafter, 126 feature wavelengths were identified from 1557 wavelengths using a genetic algorithm (GA) to build simplified classification models. Four classification algorithms, namely K-nearest neighbor (KNN), soft independent method of class analogy (SIMCA), partial least-squares discriminant analysis (PLS-DA), and support vector machine discriminant analysis (SVM-DA) were tested on full-range wavelengths and feature wavelengths, respectively. With the full-range wavelengths, all four algorithms achieved a high classification accuracy range from 97.56% to 99.59%, and the SVM-DA worked better than other models. From the feature wavelengths, no significant decline in accuracies was observed in most of the models and a high accuracy of 99.19% was still obtained by the PLS-DA model. This study demonstrated that using the FT-NIR technique with discriminant analyses could be a feasible way to classify sweet corn seed cultivars and the proper classification model could be embedded in seed sorting machinery to select high-purity seeds.

QSAR Analysis of a Series of Hydantoin-based Androgen Receptor Modulators and Corresponding Binding Affinities

Androgen receptor (AR), a member of the nuclear hormone receptor superfamily of intracellular ligand-dependent transcription factors, plays an indispensable role in normal male development through the regulation of androgen through the binding with endogenous androgens. Inappropriate amounts of androgens have a severe adverse effect on men. Excessive androgen may contribute to accelerate prostatic hypertrophy, even prostate cancer, while the absence of androgen may result in reduced muscle mass and strength, decreased bone mass, low energy, diminished sexual function and an increased risk of osteoporosis and fracture. In these cases, androgen receptor modulators are important to maintain the normal biological function of AR. So androgen receptor modulators are necessary for human being to improve their happy life index. To explore the relationships between molecular structures and corresponding binding abilities to aid the new AR modulator design, multiple linear regressions (MLR) are employed to analyze a series of hydantoin analogues, which can bind to androgen receptor acting as AR modulators. The obtained optimum model presents wonderful reliabilities and strong predictive abilities with R² =0.858, Q L O O 2 =0.822, Q L M O 2 =0.813, Q F 1 2 =0.840, Q F 2 2 =0.807, Q F 3 2 =0.814, CCC=0.893, respectively. The derived model can be used to predict the binding abilities of unknown chemicals and may help to design novel molecules with better AR affinity activity.

Single-Kernel FT-NIR Spectroscopy for Detecting Supersweet Corn (Zea mays L. Saccharata Sturt) Seed Viability with Multivariate Data Analysis

The viability and vigor of crop seeds are crucial indicators for evaluating seed quality, and high-quality seeds can increase agricultural yield. The conventional methods for assessing seed viability are time consuming, destructive, and labor intensive. Therefore, a rapid and nondestructive technique for testing seed viability has great potential benefits for agriculture. In this study, single-kernel Fourier transform near-infrared (FT-NIR) spectroscopy with a wavelength range of 1000-2500 nm was used to distinguish viable and nonviable supersweet corn seeds. Various preprocessing algorithms coupled with partial least squares discriminant analysis (PLS-DA) were implemented to test the performance of classification models. The FT-NIR spectroscopy technique successfully differentiated viable seeds from seeds that were nonviable due to overheating or artificial aging. Correct classification rates for both heat-damaged kernels and artificially aged kernels reached 98.0%. The comprehensive model could also attain an accuracy of 98.7% when combining heat-damaged samples and artificially aged samples into one category. Overall, the FT-NIR technique with multivariate data analysis methods showed great potential capacity in rapidly and nondestructively detecting seed viability in supersweet corn.

Chemical space, diversity and activity landscape analysis of estrogen receptor binders

Understanding the structure–activity relationships (SAR) of endocrine-disrupting chemicals has a major importance in toxicology. Despite the fact that classifiers and predictive models have been developed for estrogens for the past 20 years, to the best of our knowledge, there are no studies of their activity landscape or the identification of activity cliffs. Herein, we report the first SAR of a public dataset of 121 chemicals with reported estrogen receptor binding affinities using activity landscape modeling. To this end, we conducted a systematic quantitative and visual analysis of the chemical space of the 121 chemicals. The global diversity of the dataset was characterized by means of Consensus Diversity Plot, a recently developed method. Adding pairwise activity difference information to the chemical space gave rise to the activity landscape of the data set uncovering a heterogeneous SAR, in particular for some structural classes. At least eight compounds were identified with high propensity to form activity cliffs. The findings of this work further expand the current knowledge of the underlying SAR of estrogenic compounds and can be the starting point to develop novel and potentially improved predictive models.

Global diversity and activity landscape analysis of endocrine-disrupting chemicals identifies activity cliffs that are rationalized at the structure level.

Investigation and risk evaluation of the occurrence of carbamazepine, oxcarbazepine, their human metabolites and transformation products in the urban water cycle

Trace organic contaminants such as pharmaceuticals, personal care products and industrial chemicals are frequently detected in the urban water cycle, including wastewater, surface water and groundwater, as well as drinking water. These also include human metabolites (HMs), which are formed in the human body and then excreted via urine or feces, as well as transformation products (TPs) formed in engineered treatment systems and the aquatic environment. In the current study, the occurrence of HMs as well as their TPs of the anticonvulsants carbamazepine (CBZ) and oxcarbazepine (OXC) were investigated using LC tandem MS in effluents of wastewater treatment plants (WWTPs), surface water and groundwater. Highest concentrations were observed in raw wastewater for 10,11-dihydro-10,11-dihydroxycarbamazepine (DiOHCBZ), 10,11-dihydro-10-hydroxy-cabamazepine (10OHCBZ) and CBZ with concentrations ranging up to 2.7 ± 0.4, 1.7 ± 0.2 and 1.07 ± 0.06 μg L^-1, respectively. Predictions of different toxicity endpoints using a Distributed Structure-Searchable Toxicity (DSSTox) expert system query indicated that several HMs and TPs, in particular 9-carboxy-acridine (9-CA-ADIN) and acridone (ADON), may exhibit an increased genotoxicity compared to the parent compound CBZ. As 9-CA-ADIN was also detected in groundwater, a detailed investigation of the genotoxicity of 9-CA-ADIN is warranted. Investigations of an advanced wastewater treatment plant further revealed that the discharge of the investigated compounds into the aquatic environment could be substantially reduced by ozonation followed by granular activated carbon (GAC) filtration.

Rigorous 3-dimensional spectral data activity relationship approach modeling strategy for ToxCast estrogen receptor data classification, validation, and feature extraction

The estrogenic potential (expressed as a score composite of 18 high throughput screening bioassays) of 1528 compounds from the ToxCast database was modeled by a 3-dimensional spectral data activity relationship approach (3D-SDAR). Due to a lack of ¹⁷ O nuclear magnetic resonance (NMR) simulation software, the most informative carbon-carbon 3D-SDAR fingerprints were augmented with indicator variables representing oxygen atoms from carbonyl and carboxamide, ester, sulfonyl, nitro, aliphatic hydroxyl, and phenolic hydroxyl groups. To evaluate the true predictive performance of the authors' model the United States Environmental Protection Agency provided them with a blind test set consisting of 2008 compounds. Of these, 543 had available literature data-their binding affinity served to estimate the external classification accuracy of the developed model: predictive accuracy of 0.62, sensitivity of 0.71, and specificity of 0.53 were obtained. Compared with alternative modeling techniques, the authors' model displayed very little reduction in performance between the modeling and the prediction set. A 3D-SDAR mapping technique allowed identification of structural features essential for estrogenicity: 1) the presence of a phenolic OH group or cyclohexenone, 2) a second aromatic or phenolic ring at a distance of 6 Å to 8 Å from the oxygen of the first phenol ring, 3) the presence of a methyl group approximately 6 Å away from the centroid of a phenol ring, and 4) a carbonyl group in close proximity (∼4 Å measured to the centroid) to 1 of the phenol rings. Environ Toxicol Chem 2017;36:823-830. Published 2016 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

In silico investigation of agonist activity of a structurally diverse set of drugs to hPXR using HM-BSM and HM-PNN

The human pregnane X receptor (hPXR) plays a critical role in the metabolism, transport and clearance of xenobiotics in the liver and intestine. The hPXR can be activated by a structurally diverse of drugs to initiate clinically relevant drug-drug interactions. In this article, in silico investigation was performed on a structurally diverse set of drugs to identify critical structural features greatly related to their agonist activity towards hPXR. Heuristic method (HM)-Best Subset Modeling (BSM) and HM-Polynomial Neural Networks (PNN) were utilized to develop the linear and non-linear quantitative structure-activity relationship models. The applicability domain (AD) of the models was assessed by Williams plot. Statistically reliable models with good predictive power and explain were achieved (for HM-BSM, r (2)=0.881, q LOO (2) =0.797, q EXT (2) =0.674; for HM-PNN, r (2)=0.882, q LOO (2) =0.856, q EXT (2) =0.655). The developed models indicated that molecular aromatic and electric property, molecular weight and complexity may govern agonist activity of a structurally diverse set of drugs to hPXR.

QSPR prediction of the hydroxyl radical rate constant of water contaminants

In advanced oxidation processes (AOPs), the aqueous hydroxyl radical (HO) acts as a strong oxidant to react with organic contaminants. The hydroxyl radical rate constant (kHO) is important for evaluating and modelling of the AOPs. In this study, quantitative structure-property relationship (QSPR) method is applied to model the hydroxyl radical rate constant for a diverse dataset of 457 water contaminants from 27 various chemical classes. The constricted binary particle swarm optimization and multiple-linear regression (BPSO-MLR) are used to obtain the best model with eight theoretical descriptors. An optimized feed forward neural network (FFNN) is developed to investigate the complex performance of the selected molecular parameters with kHO. Although the FFNN prediction results are more accurate than those obtained using BPSO-MLR, the application of the latter is much more convenient. Various internal and external validation techniques indicate that the obtained models could predict the logarithmic hydroxyl radical rate constants of a large number of water contaminants with less than 4% absolute relative error. Finally, the above-mentioned proposed models are compared to those reported earlier and the structural factors contributing to the AOP degradation efficiency are discussed.

Conformal prediction to define applicability domain - A case study on predicting ER and AR binding

A fundamental element when deriving a robust and predictive in silico model is not only the statistical quality of the model in question but, equally important, the estimate of its predictive boundaries. This work presents a new method, conformal prediction, for applicability domain estimation in the field of endocrine disruptors. The method is applied to binders and non-binders related to the oestrogen and androgen receptors. Ensembles of decision trees are used as statistical method and three different sets (dragon, rdkit and signature fingerprints) are investigated as chemical descriptors. The conformal prediction method results in valid models where there is an excellent balance in quality between the internally validated training set and the corresponding external test set, both in terms of validity and with respect to sensitivity and specificity. With this method the level of confidence can be readily altered by the user and the consequences thereof immediately inspected. Furthermore, the predictive boundaries for the derived models are rigorously defined by using the conformal prediction framework, thus no ambiguity exists as to the level of similarity needed for new compounds to be in or out of the predictive boundaries of the derived models where reliable predictions can be expected.

3D QSAR studies of hydroxylated polychlorinated biphenyls as potential xenoestrogens

Mono-hydroxylated polychlorinated biphenyls (OH-PCBs) are found in human biological samples and lack of data on their potential estrogenic activity has been a source of concern. We have extended our previous in silico 2D QSAR study through the application of advance techniques such as docking and 3D QSAR to gain insights into their estrogen receptor (ERα) binding. The results support our earlier findings that the hydroxyl group is the most important feature on the compounds; its position, orientation and surroundings in the structure are influential for the binding of OH-PCBs to ERα. This study has also revealed the following additional interactions that influence estrogenicity of these chemicals (a) the aromatic interactions of the biphenyl moieties with the receptor, (b) hydrogen bonding interactions of the p-hydroxyl group with key amino acids ARG394 and GLU353, (c) low or no electronegative substitution at para-positions of the p-hydroxyl group, (d) enhanced electrostatic interactions at the meta position on the B ring, and (e) co-planarity of the hydroxyl group on the A ring. In combination the 2D and 3D QSAR approaches have led us to the support conclusion that the hydroxyl group is the most important feature on the OH-PCB influencing the binding to estrogen receptors, and have enhanced our understanding of the mechanistic details of estrogenicity of this class of chemicals. Such in silico computational methods could serve as useful tools in risk assessment of chemicals.

Spoilt for choice: A critical review on the chemical and biological assessment of current wastewater treatment technologies

The knowledge we have gained in recent years on the presence and effects of compounds discharged by wastewater treatment plants (WWTPs) brings us to a point where we must question the appropriateness of current water quality evaluation methodologies. An increasing number of anthropogenic chemicals is detected in treated wastewater and there is increasing evidence of adverse environmental effects related to WWTP discharges. It has thus become clear that new strategies are needed to assess overall quality of conventional and advanced treated wastewaters. There is an urgent need for multidisciplinary approaches combining expertise from engineering, analytical and environmental chemistry, (eco)toxicology, and microbiology. This review summarizes the current approaches used to assess treated wastewater quality from the chemical and ecotoxicological perspective. Discussed chemical approaches include target, non-target and suspect analysis, sum parameters, identification and monitoring of transformation products, computational modeling as well as effect directed analysis and toxicity identification evaluation. The discussed ecotoxicological methodologies encompass in vitro testing (cytotoxicity, genotoxicity, mutagenicity, endocrine disruption, adaptive stress response activation, toxicogenomics) and in vivo tests (single and multi species, biomonitoring). We critically discuss the benefits and limitations of the different methodologies reviewed. Additionally, we provide an overview of the current state of research regarding the chemical and ecotoxicological evaluation of conventional as well as the most widely used advanced wastewater treatment technologies, i.e., ozonation, advanced oxidation processes, chlorination, activated carbon, and membrane filtration. In particular, possible directions for future research activities in this area are provided.

Development and Validation of Decision Forest Model for Estrogen Receptor Binding Prediction of Chemicals Using Large Data Sets

Some chemicals in the environment possess the potential to interact with the endocrine system in the human body. Multiple receptors are involved in the endocrine system; estrogen receptor α (ERα) plays very important roles in endocrine activity and is the most studied receptor. Understanding and predicting estrogenic activity of chemicals facilitates the evaluation of their endocrine activity. Hence, we have developed a decision forest classification model to predict chemical binding to ERα using a large training data set of 3308 chemicals obtained from the U.S. Food and Drug Administration's Estrogenic Activity Database. We tested the model using cross validations and external data sets of 1641 chemicals obtained from the U.S. Environmental Protection Agency's ToxCast project. The model showed good performance in both internal (92% accuracy) and external validations (∼ 70-89% relative balanced accuracies), where the latter involved the validations of the model across different ER pathway-related assays in ToxCast. The important features that contribute to the prediction ability of the model were identified through informative descriptor analysis and were related to current knowledge of ER binding. Prediction confidence analysis revealed that the model had both high prediction confidence and accuracy for most predicted chemicals. The results demonstrated that the model constructed based on the large training data set is more accurate and robust for predicting ER binding of chemicals than the published models that have been developed using much smaller data sets. The model could be useful for the evaluation of ERα-mediated endocrine activity potential of environmental chemicals.

Reprint of "In silico methods in the discovery of endocrine disrupting chemicals"

The prevalence of sex hormone-dependent cancers, reproductive problems, obesity, and cardiovascular complications has risen especially in the Western world. It has been suggested, that the exposure to various endocrine disrupting chemicals (EDCs) contributes to the development and progression of these diseases. EDCs can interfere with various proteins: nuclear steroid hormone receptors, such as estrogen-, androgen-, glucocorticoid- and mineralocorticoid receptors (ER, AR, GR, MR), and enzymes that are involved in steroid hormone synthesis and metabolism, for example hydroxysteroid dehydrogenases (HSDs). Numerous chemicals are known as endocrine disruptors. However, the mechanism of action for most of these EDCs is still unknown. It is exhaustive and time consuming to test in vitro all chemicals - potential EDCs - used in industry, agriculture or as food preservatives against their effects on the endocrine system. Computational methods, such as virtual screening, quantitative structure activity relationships and docking, are already well recognized and used in drug development. The same methods could also aid the research on EDCs. So far, the computational methods in the search of EDCs have been retrospective. There are, however, some prospective studies reporting the use of in silico methods: five studies reporting the identification of previously unknown 17β-HSD3 inhibitors, MR agonists, and ER antagonists/agonists. This review provides an overview of case studies and in silico methods that are used in the search of EDCs. This article is part of a Special Issue entitled 'CSR 2013'.

Multicriteria approach to select polyaromatic river mutagen candidates

The identification of unknown compounds remains one of the most challenging tasks to link observed toxic effects in complex environmental mixtures to responsible toxicants in effect-directed analysis (EDA). Here, a workflow is presented based on nontarget liquid chromatography-high resolution mass spectrometry (LC-HRMS) starting with molecular formulas determined in a previous study. A compound database search (ChemSpider) was performed to retrieve candidates for each formula. Subsequently, the number of candidates was reduced by applying MS-, physical-chemical, and chromatography-based selection criteria including HRMS/MS fragmentation and plausibility, ionization efficiency with different ion sources and detection modes, acid/base behavior, octanol/water partitioning, retention time prediction and finally toxic effects (mutagenicity caused by aromatic amines). The workflow strongly decreased the number of possible candidates and resulted in the tentative identification of possible mutagens and the positive identification of the nonmutagen benzyl(diphenyl) phosphine oxide in a mutagenic fraction. The positive identification of mutagens was hampered by a lack of commercially available standards. The workflow is an innovative and promising approach and forms an excellent basis for possible further advancements.

In silico methods in the discovery of endocrine disrupting chemicals

The prevalence of sex hormone-dependent cancers, reproductive problems, obesity, and cardiovascular complications has risen especially in the Western world. It has been suggested, that the exposure to various endocrine disrupting chemicals (EDCs) contributes to the development and progression of these diseases. EDCs can interfere with various proteins: nuclear steroid hormone receptors, such as estrogen-, androgen-, glucocorticoid- and mineralocorticoid receptors (ER, AR, GR, MR), and enzymes that are involved in steroid hormone synthesis and metabolism, for example hydroxysteroid dehydrogenases (HSDs). Numerous chemicals are known as endocrine disruptors. However, the mechanism of action for most of these EDCs is still unknown. It is exhaustive and time consuming to test in vitro all chemicals - potential EDCs - used in industry, agriculture or as food preservatives against their effects on the endocrine system. Computational methods, such as virtual screening, quantitative structure activity relationships and docking, are already well recognized and used in drug development. The same methods could also aid the research on EDCs. So far, the computational methods in the search of EDCs have been retrospective. There are, however, some prospective studies reporting the use of in silico methods: five studies reporting the identification of previously unknown 17β-HSD3 inhibitors, MR agonists, and ER antagonists/agonists. This review provides an overview of case studies and in silico methods that are used in the search of EDCs. This article is part of a Special Issue entitled 'CSR 2013'.

Assessment of hydroxylated metabolites of polychlorinated biphenyls as potential xenoestrogens: a QSAR comparative analysis∗

Alternative methods, including quantitative structure-activity relationships (QSAR), are being used increasingly when appropriate data for toxicity evaluation of chemicals are not available. Approximately 40 mono-hydroxylated polychlorinated biphenyls (OH-PCBs) have been identified in humans. They represent a health and environmental concern because some of them have been shown to have agonist or antagonist interactions with human hormone receptors. This could lead to modulation of steroid hormone receptor pathways and endocrine system disruption. We performed QSAR analyses using available estrogenic activity (human estrogen receptor ER alpha) data for 71 OH-PCBs. The modelling was performed using multiple molecular descriptors including electronic, molecular, constitutional, topological, and geometrical endpoints. Multiple linear regressions and recursive partitioning were used to best fit descriptors. The results show that the position of the hydroxyl substitution, polarizability, and meta adjacent un-substituted carbon pairs at the phenolic ring contribute towards greater estrogenic activity for these chemicals. These comparative QSAR models may be used for predictive toxicity, and identification of health consequences of PCB metabolites that lack empirical data. Such information will help prioritize such molecules for additional testing, guide future basic laboratory research studies, and help the health/risk assessment community understand the complex nature of chemical mixtures.

On the development and validation of QSAR models

The fundamental and more critical steps that are necessary for the development and validation of QSAR models are presented in this chapter as best practices in the field. These procedures are discussed in the context of predictive QSAR modelling that is focused on achieving models of the highest statistical quality and with external predictive power. The most important and most used statistical parameters needed to verify the real performances of QSAR models (of both linear regression and classification) are presented. Special emphasis is placed on the validation of models, both internally and externally, as well as on the need to define model applicability domains, which should be done when models are employed for the prediction of new external compounds.

QSAR Modeling is not "Push a Button and Find a Correlation": A Case Study of Toxicity of (Benzo-)triazoles on Algae

A case study of toxicity of (benzo)triazoles ((B)TAZs) to the algae Pseudokirchneriella subcapitata is used to discuss some problems and solutions in QSAR modeling, particularly in the environmental context. The relevance of data curation (not only of experimental data, but also of chemical structures and input formats for the calculation of molecular descriptors), the crucial points of QSAR model validation and the potential application for new chemicals (internal robustness, exclusion of chance correlation, external predictivity, applicability domain) are described, while developing MLR-OLS models based on molecular descriptors, calculated by various QSAR software tools (commercial DRAGON, free PaDEL-Descriptor and QSPR-THESAURUS). Additionally, the utility of consensus models is highlighted. This work summarizes a methodology for a rigorous statistical approach to obtain reliable QSAR predictions, also for a large number of (B)TAZs in the ECHA preregistration list of REACH (even if starting from limited experimental data availability), and has evidenced some ambiguities and discrepancies related to SMILES notations from different databases; furthermore it highlighted some general problems related to QSAR model generation and was useful in the implementation of the PaDEL-Descriptor software.

Prediction of acute mammalian toxicity using QSAR methods: a case study of sulfur mustard and its breakdown products

Predicting toxicity quantitatively, using Quantitative Structure Activity Relationships (QSAR), has matured over recent years to the point that the predictions can be used to help identify missing comparison values in a substance's database. In this manuscript we investigate using the lethal dose that kills fifty percent of a test population (LD₅₀) for determining relative toxicity of a number of substances. In general, the smaller the LD₅₀ value, the more toxic the chemical, and the larger the LD₅₀ value, the lower the toxicity. When systemic toxicity and other specific toxicity data are unavailable for the chemical(s) of interest, during emergency responses, LD₅₀ values may be employed to determine the relative toxicity of a series of chemicals. In the present study, a group of chemical warfare agents and their breakdown products have been evaluated using four available rat oral QSAR LD₅₀ models. The QSAR analysis shows that the breakdown products of Sulfur Mustard (HD) are predicted to be less toxic than the parent compound as well as other known breakdown products that have known toxicities. The QSAR estimated break down products LD₅₀ values ranged from 299 mg/kg to 5,764 mg/kg. This evaluation allows for the ranking and toxicity estimation of compounds for which little toxicity information existed; thus leading to better risk decision making in the field.

Docking and QSAR study on the binding interactions between polycyclic aromatic hydrocarbons and estrogen receptor

Little is known about the estrogenic activities of polycyclic aromatic hydrocarbons (PAHs) and the underlying mechanisms on estrogenic activities are still unclear. Molecular docking and quantitative structure-activity relationship (QSAR) were used to understand the relationship between molecular structural features and estrogenic activity, and to predict the binding affinity of PAHs to estrogen receptor α (ERα). From molecular docking analysis, hydrogen bonding as well as hydrophobic and π interactions were found between PAHs and ERα. Based on the docking results, appropriate molecular structural parameters were adopted to develop a QSAR model. Five descriptors were included in the QSAR model, which indicated that the estrogenic activity was related to molecular size, van der Waals volumes, shape profiles, polarizabilities and electropological states were significant parameters explaining the estrogenicity. Comparatively, the developed QSAR model had good robustness, predictive ability and mechanistic interpretability. Moreover, the applicability domain of the model was described.

Conformation depends on 4D-QSAR analysis using EC-GA method: pharmacophore identification and bioactivity prediction of TIBOs as non-nucleoside reverse transcriptase inhibitors

The electron conformational and genetic algorithm methods (EC-GA) were integrated for the identification of the pharmacophore group and predicting the anti HIV-1 activity of tetrahydroimidazo[4,5,1-jk][1,4]benzodiazepinone (TIBO) derivatives. To reveal the pharmacophore group, each conformation of all compounds was arranged by electron conformational matrices of congruity. Multiple comparisons of these matrices, within given tolerances for high active and low active TIBO derivatives, allow the identification of the pharmacophore group that refers to the electron conformational submatrix of activity. The effects of conformations, internal and external validation were investigated by four different models based on an ensemble of conformers and a single conformer, both with and without a test set. Model 1 using an ensemble of conformers for the training (39 compounds) and test sets (13 compounds), obtained by the optimum seven parameters, gave satisfactory results (R²(training) = 0.878, R²(test)= 0.910, q² = 0.840, q²(ext1) = 0.926 and q²(ext2) = 0.900).