REACH and in silico methods: an attractive opportunity for medicinal chemists

Where developmental toxicity meets explainable artificial intelligence: state-of-the-art and perspectives

INTRODUCTION

The application of Artificial Intelligence (AI) to predictive toxicology is rapidly increasing, particularly aiming to develop non-testing methods that effectively address ethical concerns and reduce economic costs. In this context, Developmental Toxicity (Dev Tox) stands as a key human health endpoint, especially significant for safeguarding maternal and child well-being.

AREAS COVERED

This review outlines the existing methods employed in Dev Tox predictions and underscores the benefits of utilizing New Approach Methodologies (NAMs), specifically focusing on eXplainable Artificial Intelligence (XAI), which proves highly efficient in constructing reliable and transparent models aligned with recommendations from international regulatory bodies.

EXPERT OPINION

The limited availability of high-quality data and the absence of dependable Dev Tox methodologies render XAI an appealing avenue for systematically developing interpretable and transparent models, which hold immense potential for both scientific evaluations and regulatory decision-making.

Chemometrics-driven prediction and prioritization of diverse pesticides on chickens for addressing hazardous effects on public health

The extensive use of various pesticides in the agriculture field badly affects both chickens and humans, primarily through residues in food products and environmental exposure. This study offers the first quantitative structure-toxicity relationship (QSTR) and quantitative read-across-structure toxicity relationship (q-RASTR) models encompassing the LOEL and NOEL endpoints for acute toxicity in chicken, a widely consumed protein. The study's significance lies in the direct link between chemical toxicity in chicken, human intake, and environmental damage. Both the QSTR and the similarity-based read-across algorithms are applied concurrently to improve the predictability of the models. The q-RASTR models were generated by combining read-across derived similarity and error-based parameters, alongside structural and physicochemical descriptors. Machine Learning approaches (SVM and RR) were also employed with the optimization of relevant hyperparameters based on the cross-validation approach, and the final test set prediction results were compared. The PLS-based q-RASTR models for NOEL and LOEL endpoints showed good statistical performance, as traced from the external validation metrics Q2F1: 0.762-0.844; Q2F2: 0.759-0.831 and MAEtest: 0.195-0.214. The developed models were further used to screen the Pesticide Properties DataBase (PPDB) for potential toxicants in chickens. Thus, established models can address eco-toxicological data gaps and development of novel and safe eco-friendly pesticides.

Recent advances and current challenges of new approach methodologies in developmental and adult neurotoxicity testing

Adult neurotoxicity (ANT) and developmental neurotoxicity (DNT) assessments aim to understand the adverse effects and underlying mechanisms of toxicants on the human nervous system. In recent years, there has been an increasing focus on the so-called new approach methodologies (NAMs). The Organization for Economic Co-operation and Development (OECD), together with European and American regulatory agencies, promote the use of validated alternative test systems, but to date, guidelines for regulatory DNT and ANT assessment rely primarily on classical animal testing. Alternative methods include both non-animal approaches and test systems on non-vertebrates (e.g., nematodes) or non-mammals (e.g., fish). Therefore, this review summarizes the recent advances of NAMs focusing on ANT and DNT and highlights the potential and current critical issues for the full implementation of these methods in the future. The status of the DNT in vitro battery (DNT IVB) is also reviewed as a first step of NAMs for the assessment of neurotoxicity in the regulatory context. Critical issues such as (i) the need for test batteries and method integration (from in silico and in vitro to in vivo alternatives, e.g., zebrafish, C. elegans) requiring interdisciplinarity to manage complexity, (ii) interlaboratory transferability, and (iii) the urgent need for method validation are discussed.

NeuTox: A weighted ensemble model for screening potential neuronal cytotoxicity of chemicals based on various types of molecular representations

Chemical-induced neurotoxicity has been widely brought into focus in the risk assessment of chemical safety. However, the traditional in vivo animal models to evaluate neurotoxicity are time-consuming and expensive, which cannot completely represent the pathophysiology of neurotoxicity in humans. Cytotoxicity to human neuroblastoma cell line (SH-SY5Y) is commonly used as an alternative to animal testing for the assessment of neurotoxicity, yet it is still not appropriate for high throughput screening of potential neuronal cytotoxicity of chemicals. In this study, we constructed an ensemble prediction model, termed NeuTox, by combining multiple machine learning algorithms with molecular representations based on the weighted score of Particle Swarm Optimization. For the test set, NeuTox shows excellent performance with an accuracy of 0.9064, which are superior to the top-performing individual models. The subsequent experimental verifications reveal that 5,5'-isopropylidenedi-2-biphenylol and 4,4'-cyclo-hexylidenebisphenol exhibited stronger SH-SY5Y-based cytotoxicity compared to bisphenol A, suggesting that NeuTox has good generalization ability in the first-tier assessment of neuronal cytotoxicity of BPA analogs. For ease of use, NeuTox is presented as an online web server that can be freely accessed via http://www.iehneutox-predictor.cn/NeuToxPredict/Predict.

TISBE: A Public Web Platform for the Consensus-Based Explainable Prediction of Developmental Toxicity

Despite being extremely relevant for the protection of prenatal and neonatal health, the developmental toxicity (Dev Tox) is a highly complex endpoint whose molecular rationale is still largely unknown. The lack of availability of high-quality data as well as robust nontesting methods makes its understanding even more difficult. Thus, the application of new explainable alternative methods is of utmost importance, with Dev Tox being one of the most animal-intensive research themes of regulatory toxicology. Descending from TIRESIA (Toxicology Intelligence and Regulatory Evaluations for Scientific and Industry Applications), the present work describes TISBE (TIRESIA Improved on Structure-Based Explainability), a new public web platform implementing four fundamental advancements for in silico analyses: a three times larger dataset, a transparent XAI (explainable artificial intelligence) framework employing a fragment-based fingerprint coding, a novel consensus classifier based on five independent machine learning models, and a new applicability domain (AD) method based on a double top-down approach for better estimating the prediction reliability. The training set (TS) includes as many as 1008 chemicals annotated with experimental toxicity values. Based on a 5-fold cross-validation, a median value of 0.410 for the Matthews correlation coefficient was calculated; TISBE was very effective, with a median value of sensitivity and specificity equal to 0.984 and 0.274, respectively. TISBE was applied on two external pools made of 1484 bioactive compounds and 85 pediatric drugs taken from ChEMBL (Chemical European Molecular Biology Laboratory) and TEDDY (Task-Force in Europe for Drug Development in the Young) repositories, respectively. Notably, TISBE gives users the option to clearly spot the molecular fragments responsible for the toxicity or the safety of a given chemical query and is available for free at https://prometheus.farmacia.uniba.it/tisbe.

Harmonization Risks and Rewards: Nano-QSAR for Agricultural Nanomaterials

This comprehensive review explores the emerging landscape of Nano-QSAR (quantitative structure-activity relationship) for assessing the risk and potency of nanomaterials in agricultural settings. The paper begins with an introduction to Nano-QSAR, providing background and rationale, and explicitly states the hypotheses guiding the review. The study navigates through various dimensions of nanomaterial applications in agriculture, encompassing their diverse properties, types, and associated challenges. Delving into the principles of QSAR in nanotoxicology, this article elucidates its application in evaluating the safety of nanomaterials, while addressing the unique limitations posed by these materials. The narrative then transitions to the progression of Nano-QSAR in the context of agricultural nanomaterials, exemplified by insightful case studies that highlight both the strengths and the limitations inherent in this methodology. Emerging prospects and hurdles tied to Nano-QSAR in agriculture are rigorously examined, casting light on important pathways forward, existing constraints, and avenues for research enhancement. Culminating in a synthesis of key insights, the review underscores the significance of Nano-QSAR in shaping the future of nanoenabled agriculture. It provides strategic guidance to steer forthcoming research endeavors in this dynamic field.

Predictive binary mixture toxicity modeling of fluoroquinolones (FQs) and the projection of toxicity of hypothetical binary FQ mixtures: a combination of 2D-QSAR and machine-learning approaches

All sorts of chemicals get degraded under various environmental stresses, and the degradates coexist with the parent compounds as mixtures in the environment. Antibiotics emerge as an additional concern due to the bioactive nature of both the parent compound and degradation products and their combined exposure to the environment. Therefore, environmental risk assessment of antibiotics and their degradation products is very much necessary. In this direction, we made use of in silico new approach methodologies (NAMs) and machine-learning algorithms. In this study, we have developed a robust and predictive mixture-quantitative structure-activity relationship (QSAR) model with promising quality and predictability (internal: MAETrain = 0.085, QLOO2 = 0.849, external: MAETest = 0.090, and QF12 = 0.859) for predicting the toxicity of the mixtures of a class of antibiotics and their degradation products. To obtain the predictive model, toxicity data of 78 binary fluoroquinolone mixtures in E. coli (endpoint: log 1/IC50 in molar) have been utilized. We have used only 0D-2D descriptors to efficiently encode the structural features of mixture components without any additional complexities. The optimization of the class of mixture descriptors has been performed in this study by using three different mixing rules (linear combination of molecular contributions, the squared molecular contributions, and the norm of molecular contributions). Different machine-learning approaches namely, random forest (RF), ada boost, gradient boost (GB), extreme gradient boost (XGB), support vector machine (SVM), linear support vector machine (LSVM), and ridge regression (RR) have been employed here apart from the conventional partial least squares (PLS) regression to optimize the modeling approach. A rigorous validation protocol has been used for assessing the goodness-of-fit, robustness, and external predictability of the models. Finally, the toxicity of possible untested mixtures of different photodegradation products of fluoroquinolones has been predicted using the best model reported in this study.

QSAR assessment of aquatic toxicity potential of diverse agrochemicals

The fast-increasing number of commercially produced chemicals challenges the experimental ecotoxicity assessment methods, which are costly, time-consuming, and dependent on the sacrifice of animals. In this regard, Quantitative Structure-Property/Activity Relationships (QSPR/QSAR) have led the way in developing ecotoxicity assessment models. In this study, QSAR models have been developed using the pEC50 values of 82 diverse agrochemicals or agro-molecules against a planktonic crustacean Daphnia magna with easily interpretable 2D descriptors. Moreover, a link among octanol-water partition coefficient (KOW), bio-concentration factor (BCF), and critical body residue (CBR) has been addressed, and their imputation for the prediction of the toxicity endpoint (EC50) has been done with an objective of the advanced exploration of several ecotoxicological parameters for toxic chemicals. The developed partial least squares (PLS) models were validated rigorously and proved to be robust, sound, and immensely well-predictive. The final Daphnia toxicity model derived from experimental derived properties along with computed descriptors emerged better in statistical quality and predictivity than those obtained solely from computed descriptors. Additionally, the pEC50 and other important properties (log KOW, log BCF, and log CBR) for a set of external agro-molecules, not employed in model development, were predicted to show the predictive ability of the models.

"Data fusion" quantitative read-across structure-activity-activity relationships (q-RASAARs) for the prediction of toxicities of binary and ternary antibiotic mixtures toward three bacterial species

Antibiotics are often found in the environment as pollutants. They are usually found as mixtures in the environment and may produce toxicity against different ecological species due to joint exposure in the sub-optimal range. Sometimes the degradation products of parent chemicals also interact with it and cause mixture toxicity. In this study, we have developed three different mixture-Quantitative Structure-Activity Relationship (mixture-QSAR) models for three different bacterial species (Vibrio fischeri, Escherichia coli, and Bacillus subtilis). The toxicity data were collected from a previous experimental report in the literature, which comprised binary and ternary mixtures of sulfonamides (SAs), sulfonamide potentiators (SAPs), and tetracyclines (TCs). We have also explored the interspecies modeling to find inter-correlation among the toxicity of these studied organisms and have developed quantitative structure activity-activity relationship (QSAAR) models by employing the "data fusion" quantitative read-across structure-activity-activity relationship (q-RASAAR) and partial least squares (PLS) regression algorithms. All the models are strictly validated using both internal and external validation tests as suggested in the OECD guidelines. Three different mixing rules have been used in this study for descriptor computations to incorporate the additive and interaction effects among the mixture components. To the best of our knowledge, this is the first report of interspecies mixture toxicity models which can predict the cellular toxicity of binary and ternary mixtures against any of the three above-mentioned organisms.

In silico prediction of the mutagenicity of nitroaromatic compounds using correlation weights of fragments of local symmetry

Most quantitative structure-property/activity relationships (QSPRs/QSARs) techniques involve using different programs separately for generating molecular descriptors and separately for building models based on available descriptors. Here, the capabilities of the CORAL program are evaluated. A user of the program should apply as the basis for models the representation of the molecular structure by means of the simplified molecular input-line entry system (SMILES) as well as experimental data on the endpoint of interest. The local symmetry of SMILES is a novel composition of symmetrically represented symbols, which are three 'xyx', four 'xyyx', or five symbols 'xyzyx'. We updated our CORAL software using this optimal, new flexible descriptor, sensitive to the symmetric composition of a specific part of the molecule. Computational experiments have shown that taking account of these attributes of SMILES can improve the predictive potential of models for the mutagenicity of nitroaromatic compounds. In addition, the above computational experiments have confirmed the advantage of using the index of ideality of correlation (IIC) and the correlation intensity index (CII) for Monte Carlo optimization of the correlation weights for various attributes of SMILES, including the local symmetry. The average value of the coefficient of determination for the validation set (five different models) without fragments of local symmetry is 0.8589 ± 0.025, whereas using fragments of local symmetry improves this criterion of the predictive potential up to 0.9055 ± 0.010.

Aquatic toxicity prediction of diverse pesticides on two algal species using QSTR modeling approach

With the aim of identification of toxic nature of the diverse pesticides on the aquatic compartment, a large dataset of pesticides (n = 325) with experimental toxicity data on two algal test species (Pseudokirchneriella subcapitata (PS) (synonym: Raphidocelis subcapitata, Selenastrum capricornutum) and Scenedemus subspicatus (SS)) was gathered and subjected to quantitative structure toxicity relationship (QSTR) analysis to predict aquatic toxicity of pesticides. The QSTR models were developed by multiple linear regressions (MLRs), and the genetic algorithm (GA) was used for the variable selection. The developed GA-MLR models were statistically robust enough internally (Q²_LOO = 0.620-0.663) and externally (Q²F_n = 0.693-0.868, CCC_ext = 0.843-0.877). The leverage approach of applicability domain (AD) and prediction reliability indicator assured the reliability of the developed models. The mechanistic interpretation highlighted that the presence of SO₂, F and aromatic rings influenced the toxicity of pesticides towards PS species while the presence of alkyl, alkyl halide, aromatic rings and carbonyl was responsible for the toxicity of pesticides towards SS species. Additionally, we have reported the application of developed models to pesticides without experimental value and the cumulative toxicity of pesticides on the aquatic environment by using principal component analysis (PCA). The reliable prediction and prioritization of toxic compounds from the developed models will be useful in the aquatic toxicity assessment of pesticides.

Experimental and Computational Nanotoxicology-Complementary Approaches for Nanomaterial Hazard Assessment

The growing development and applications of nanomaterials lead to an increasing release of these materials in the environment. The adverse effects they may elicit on ecosystems or human health are not always fully characterized. Such potential toxicity must be carefully assessed with the underlying mechanisms elucidated. To that purpose, different approaches can be used. First, experimental toxicology consisting of conducting in vitro or in vivo experiments (including clinical studies) can be used to evaluate the nanomaterial hazard. It can rely on variable models (more or less complex), allowing the investigation of different biological endpoints. The respective advantages and limitations of in vitro and in vivo models are discussed as well as some issues associated with experimental nanotoxicology. Perspectives of future developments in the field are also proposed. Second, computational nanotoxicology, i.e., in silico approaches, can be used to predict nanomaterial toxicity. In this context, we describe the general principles, advantages, and limitations especially of quantitative structure–activity relationship (QSAR) models and grouping/read-across approaches. The aim of this review is to provide an overview of these different approaches based on examples and highlight their complementarity.

Ecotoxicological QSTR and QSTTR Modeling for the Prediction of Acute Oral Toxicity of Pesticides against Multiple Avian Species

The ever-increasing use of pesticides in response to the rising agricultural demand has threatened the existence of nontarget organisms like avian species, disrupting the global ecological integrity. Therefore, it is critical to protect and restore different endangered bird species from the perspective of ecosystem safety. In the present work, we have developed regression-based two-dimensional quantitative structure toxicity relationship (2D QSTR) and quantitative structure toxicity-toxicity relationship (QSTTR) models to estimate the toxicity of pesticides on five different avian species following the Organization for Economic Co-operation and Development (OECD) guidelines. Rigorous validation has been performed using different statistical internal and external validation parameters to ensure the robustness and interpretability of the developed models. From the developed models, it can be stated that the presence of electronegative and lipophilic features greatly enhance pesticide toxicity, whereas the hydrophilic characters are shown to have a detrimental impact on the toxicity of pesticides. Moreover, the developed QSTTR models have been employed to the in silico toxicity prediction of 124, 154, and 250 pesticides against bobwhite quail, ring-necked pheasant, and mallard duck species, respectively, extracted from the Office of Pesticides Program (OPP) Pesticide Ecotoxicity Database. The information obtained from the modeled descriptors might be used for pesticide risk assessment in the future, with the added benefit of providing an early caution of their possible negative impact on birds for regulatory purposes.

Trimethoxylated Halogenated Chalcones as Dual Inhibitors of MAO-B and BACE-1 for the Treatment of Neurodegenerative Disorders

Six halogenated trimethoxy chalcone derivatives (CH1-CH6) were synthesized and spectrally characterized. The compounds were further evaluated for their inhibitory potential against monoamine oxidases (MAOs) and β-secretase (BACE-1). Six compounds inhibited MAO-B more effectively than MAO-A, and the 2',3',4'-methoxy moiety in CH4-CH6 was more effective for MAO-B inhibition than the 2',4',6'-methoxy moiety in CH1-CH3. Compound CH5 most potently inhibited MAO-B, with an IC₅₀ value of 0.46 µM, followed by CH4 (IC₅₀ = 0.84 µM). In 2',3',4'-methoxy derivatives (CH4-CH6), the order of inhibition was -Br in CH5 > -Cl in CH4 > -F in CH6 at the para-position in ring B of chalcone. CH4 and CH5 were selective for MAO-B, with selectivity index (SI) values of 15.1 and 31.3, respectively, over MAO-A. CH4 and CH5 moderately inhibited BACE-1 with IC₅₀ values of 13.6 and 19.8 µM, respectively. When CH4 and CH5 were assessed for their cell viability studies on the normal African Green Monkey kidney cell line (VERO) using MTT assays, it was noted that both compounds were found to be safe, and only a slightly toxic effect was observed in concentrations above 200 µg/mL. CH4 and CH5 decreased reactive oxygen species (ROS) levels of VERO cells treated with H₂O₂, indicating both compounds retained protective effects on the cells by antioxidant activities. All compounds showed high blood brain barrier permeabilities analyzed by a parallel artificial membrane permeability assay (PAMPA). Molecular docking and ADME prediction of the lead compounds provided more insights into the rationale behind the binding and the CNS drug likeness. From non-test mutagenicity and cardiotoxicity studies, CH4 and CH5 were non-mutagenic and non-/weak-cardiotoxic. These results suggest that CH4 and CH5 could be considered candidates for the cure of neurological dysfunctions.

Predictive classification-based QSTR models for toxicity study of diverse pesticides on multiple avian species

Protection and restoration of different endangered bird species from pesticide exposure is crucial from the point of safety assessment of ecosystem. Toxicity predictions or risk assessment of pesticides by chemometric tools is one of the challenging fields in recent era. In the present study, classification-based quantitative structure toxicity relationship (QSTR) models were developed for a large dataset (516) of diverse pesticides on multiple avian species mallard duck, bobwhite quail, and zebra finch according to the Organization for Economic Co-operation and Development guidelines. The QSTR models were developed by linear discriminant analysis method with genetic algorithm for feature selection from 2D descriptors using QSAR-Co software. Different statistical metrics assured the reliability and robustness of the developed models. External compound prediction highlighted predictive nature of the models. The mechanistic interpretation suggested that presence of phosphate, halogens (Cl, Br), ether linkage, and NCOO influence the avian toxicity. Furthermore, model reliability was checked by the application of the standardization approach of the applicability domain (AD). Finally, the developed models provided a priori toxic and non-toxic classification for unknown pesticides (inside AD), with particular emphasis on organophosphate pesticides. The interspecies toxicity correlation and predictions encouraged for their further applicability for the fulfilment of data gaps in vital missing species.

Virtual screening of potentially endocrine-disrupting chemicals against nuclear receptors and its application to identify PPARγ-bound fatty acids

Nuclear receptors (NRs) are key regulators of energy homeostasis, body development, and sexual reproduction. Xenobiotics binding to NRs may disrupt natural hormonal systems and induce undesired adverse effects in the body. However, many chemicals of concerns have limited or no experimental data on their potential or lack-of-potential endocrine-disrupting effects. Here, we propose a virtual screening method based on molecular docking for predicting potential endocrine-disrupting chemicals (EDCs) that bind to NRs. For 12 NRs, we systematically analyzed how multiple crystal structures can be used to distinguish actives and inactives found in previous high-throughput experiments. Our method is based on (i) consensus docking scores from multiple structures at a single functional state (agonist-bound or antagonist-bound), (ii) multiple functional states (agonist-bound and antagonist-bound), and (iii) multiple pockets (orthosteric site and alternative sites) of these NRs. We found that the consensus enrichment from multiple structures is better than or comparable to the best enrichment from a single structure. The discriminating power of this consensus strategy was further enhanced by a chemical similarity-weighted scoring scheme, yielding better or comparable enrichment for all studied NRs. Applying this optimized method, we screened 252 fatty acids against peroxisome proliferator-activated receptor gamma (PPARγ) and successfully identified 3 previously unknown fatty acids with Kd = 100-250 μM including two furan fatty acids: furannonanoic acid (FNA) and furanundecanoic acid (FUA), and one cyclopropane fatty acid: phytomonic acid (PTA). These results suggested that the proposed method can be used to rapidly screen and prioritize potential EDCs for further experimental evaluations.

In silico prediction of chemical neurotoxicity using machine learning

Neurotoxicity is one of the main causes of drug withdrawal, and the biological experimental methods of detecting neurotoxic toxicity are time-consuming and laborious. In addition, the existing computational prediction models of neurotoxicity still have some shortcomings. In response to these shortcomings, we collected a large number of data set of neurotoxicity and used PyBioMed molecular descriptors and eight machine learning algorithms to construct regression prediction models of chemical neurotoxicity. Through the cross-validation and test set validation of the models, it was found that the extra-trees regressor model had the best predictive effect on neurotoxicity ([Formula: see text] = 0.784). In addition, we get the applicability domain of the models by calculating the standard deviation distance and the lever distance of the training set. We also found that some molecular descriptors are closely related to neurotoxicity by calculating the contribution of the molecular descriptors to the models. Considering the accuracy of the regression models, we recommend using the extra-trees regressor model to predict the chemical autonomic neurotoxicity.

Prediction of Acute Oral Systemic Toxicity Using a Multifingerprint Similarity Approach

The implementation of nonanimal approaches is of particular importance to regulatory agencies for the prediction of potential hazards associated with acute exposures to chemicals. This work was carried out in the framework of an international modeling initiative organized by the Acute Toxicity Workgroup (ATWG) of the Interagency Coordinating Committee on the Validation of Alternative Methods (ICCVAM) with the participation of 32 international groups across government, industry, and academia. Our contribution was to develop a multifingerprints similarity approach for predicting five relevant toxicology endpoints related to the acute oral systemic toxicity that are: the median lethal dose (LD50) point prediction, the "nontoxic" (LD50 > 2000 mg/kg) and "very toxic" (LD50<50 mg/kg) binary classification, and the multiclass categorization of chemicals based on the United States Environmental Protection Agency and Globally Harmonized System of Classification and Labeling of Chemicals schemes. Provided by the ICCVAM's ATWG, the training set used to develop the models consisted of 8944 chemicals having high-quality rat acute oral lethality data. The proposed approach integrates the results coming from a similarity search based on 19 different fingerprint definitions to return a consensus prediction value. Moreover, the herein described algorithm is tailored to properly tackling the so-called toxicity cliffs alerting that a large gap in LD50 values exists despite a high structural similarity for a given molecular pair. An external validation set made available by ICCVAM and consisting in 2896 chemicals was employed to further evaluate the selected models. This work returned high-accuracy predictions based on the evaluations conducted by ICCVAM's ATWG.

Human ether-à-go-go-related potassium channel: exploring SAR to improve drug design

hERG is best known as a primary anti-target, the inhibition of which is responsible for serious side effects. A renewed interest in hERG as a desired target, especially in oncology, was sparked because of its role in cellular proliferation and apoptosis. In this study, we survey the most recent advances regarding hERG by focusing on SAR in the attempt to elucidate, at a molecular level, off-target and on-target actions of potential hERG binders, which are highly promiscuous and largely varying in structure. Understanding the rationale behind hERG interactions and the molecular determinants of hERG activity is a real challenge and comprehension of this is of the utmost importance to prioritize compounds in early stages of drug discovery and to minimize cardiotoxicity attrition in preclinical and clinical studies.

Best Practices for QSAR Model Reporting: Physical and Chemical Properties, Ecotoxicity, Environmental Fate, Human Health, and Toxicokinetics Endpoints

BACKGROUND

Quantitative and qualitative structure–activity relationships (QSARs) have been used to understand chemical behavior for almost a century. The main source of QSAR models is the scientific literature, but the open question is how well these models are documented.

OBJECTIVES

The main aim of this study was to critically analyze the publication practices of QSARs with regard to transparency, potential reproducibility, and independent verification. The focus was on the level of technical completeness of the published QSARs.

METHODS

A total of 1,533 QSAR articles reporting 79 individual endpoints, mostly in environmental and health science, were reviewed. The QSAR parameters required for technical completeness were grouped into five categories: chemical structures, experimental endpoint values, descriptor values, mathematical representation of the model, and predicted endpoint values. The data were summarized and discussed using Circos plots.

RESULTS

Altogether, 42.5% of the reviewed articles were found to be potentially reproducible. The potential reproducibility for different endpoint groups varied; the respective rates were 39% for physical and chemical properties, 52% for ecotoxicity, 56% for environmental fate, 30% for human health, and 32% for toxicokinetics. The reproducibility of QSARs is discussed and placed in the context of the reproducibility of the experimental methods. Included are 65 references to open QSAR datasets as examples of models restored from scientific articles.

DISCUSSION

Strikingly poor documentation of QSARs was observed, which reduces the transparency, availability, and consequently, the application of research results in scientific, industrial, and regulatory areas. A list of the components needed to ensure the best practices for QSAR reporting is provided, allowing long-term use and preservation of the models. This list also allows an assessment of the reproducibility of models by interested parties such as journal editors, reviewers, regulators, evaluators, and potential users. https://doi.org/10.1289/EHP3264.

A systems biology approach to predictive developmental neurotoxicity of a larvicide used in the prevention of Zika virus transmission

The need to prevent developmental brain disorders has led to an increased interest in efficient neurotoxicity testing. When an epidemic of microcephaly occurred in Brazil, Zika virus infection was soon identified as the likely culprit. However, the pathogenesis appeared to be complex, and a larvicide used to control mosquitoes responsible for transmission of the virus was soon suggested as an important causative factor. Yet, it is challenging to identify relevant and efficient tests that are also in line with ethical research defined by the 3Rs rule (Replacement, Reduction and Refinement). Especially in an acute situation like the microcephaly epidemic, where little toxicity documentation is available, new and innovative alternative methods, whether in vitro or in silico, must be considered. We have developed a network-based model using an integrative systems biology approach to explore the potential developmental neurotoxicity, and we applied this method to examine the larvicide pyriproxyfen widely used in the prevention of Zika virus transmission. Our computational model covered a wide range of possible pathways providing mechanistic hypotheses between pyriproxyfen and neurological disorders via protein complexes, thus adding to the plausibility of pyriproxyfen neurotoxicity. Although providing only tentative evidence and comparisons with retinoic acid, our computational systems biology approach is rapid and inexpensive. The case study of pyriproxyfen illustrates its usefulness as an initial or screening step in the assessment of toxicity potentials of chemicals with incompletely known toxic properties.

Considerations of nano-QSAR/QSPR models for nanopesticide risk assessment within the European legislative framework

The European market for pesticides is currently legislated through the well-developed Regulation (EC) No. 1107/2009. This regulation promotes the competitiveness of European agriculture, recognizing the necessity of safe pesticides for human and animal health and the environment to protect crops against pests, diseases and weeds. In this sense, nanotechnology can provide a tremendous opportunity to achieve a more rational use of pesticides. However, the lack of information regarding nanopesticides and their fate and behavior in the environment and their effects on human and animal health is inhibiting rapid nanopesticide incorporation into European Union agriculture. This review analyzes the recent state of knowledge on nanopesticide risk assessment, highlighting the challenges that need to be overcame to accelerate the arrival of these new tools for plant protection to European agricultural professionals. Novel nano-Quantitative Structure-Activity/Structure-Property Relationship (nano-QSAR/QSPR) tools for risk assessment are analyzed, including modeling methods and validation procedures towards the potential of these computational instruments to meet the current requirements for authorization of nanoformulations. Future trends on these issues, of pressing importance within the context of the current European pesticide legislative framework, are also discussed. Standard protocols to make high-quality and well-described datasets for the series of related but differently sized nanoparticles/nanopesticides are required.

Automated identification of structurally heterogeneous and patentable antiproliferative hits as potential tubulin inhibitors

By employing a recently developed hierarchical computational platform, we identified 37 novel and structurally diverse tubulin targeting compounds. In particular, hierarchical molecular filters, based on molecular shape similarity, structure-based pharmacophore, and molecular docking, were applied on a large chemical collection of commercial compounds to identify unexplored and patentable microtubule-destabilizing candidates. The herein proposed 37 novel hits, showing new molecular scaffolds (such as 1,3,3a,4-tetraaza-1,2,3,4,5,6,7,7a-octahydroindene or dihydropyrrolidin-2-one fused to a chromen-4-one), are provided with antiproliferative activity in the μm range toward MCF-7 (human breast cancer lines). Importantly, there is a likely causative relationship between cytotoxicity and the inhibition of tubulin polymerization at the colchicine binding site, assessed through fluorescence polymerization assays.

Strategies of Virtual Screening in Medicinal Chemistry

Virtual screening represents an effective computational strategy to rise-up the chances of finding new bioactive compounds by accelerating the time needed to move from an initial intuition to market. Classically, the most pursued approaches rely on ligand- and structure-based studies, the former employed when structural data information about the target is missing while the latter employed when X-ray/NMR solved or homology models are instead available for the target. The authors will focus on the most advanced techniques applied in this area. In particular, they will survey the key concepts of virtual screening by discussing how to properly select chemical libraries, how to make database curation, how to applying and- and structure-based techniques, how to wisely use post-processing methods. Emphasis will be also given to the most meaningful databases used in VS protocols. For the ease of discussion several examples will be presented.

Nontest Methods to Predict Acute Toxicity: State of the Art for Applications of In Silico Methods

The assessment of acute toxicity of chemicals by in silico methods is actually done by two methodologies, read-across and QSAR. The two approaches are strongly based on the similarity between the chemical for which a risk assessment is required and the reference chemical(s) for which the experimental data are known. Here, we describe the two methodologies with some main publications as illustrations and the in silico data associated with acute toxicity endpoints (ECHA, REACH) accessible via eChemPortal.

Applications of Chemoinformatics in Predictive Toxicology for Regulatory Purposes, Especially in the Context of the EU REACH Legislation

Chemoinformatics methodologies such as QSAR/QSPR have been used for decades in drug discovery projects, especially for the finding of new compounds with therapeutic properties and the optimization of ADME properties on chemical series. The application of computational techniques in predictive toxicology is much more recent, and they are experiencing an increasingly interest because of the new legal requirements imposed by national and international regulations. In the pharmaceutical field, the US Food and Drug Administration (FDA) support the use of predictive models for regulatory decision-making when assessing the genotoxic and carcinogenic potential of drug impurities. In Europe, the REACH legislation promotes the use of QSAR in order to reduce the huge amount of animal testing needed to demonstrate the safety of new chemical entities subjected to registration, provided they meet specific conditions to ensure their quality and predictive power. In this review, the authors summarize the state of art of in silico methods for regulatory purposes, with especial emphasis on QSAR models.

Molecular Docking for Predictive Toxicology

Molecular docking is an in silico method widely applied in drug discovery programs to predict the binding mode of a given molecule interacting with a specific biological target. This computational technique is today emerging also in the field of predictive toxicology for regulatory purposes, being for instance successfully applied to develop classification models for the prediction of the endocrine disruptor potential of chemicals. Herein, we describe the protocol for adapting molecular docking to the purposes of predictive toxicology.

(Q)SAR Methods for Predicting Genotoxicity and Carcinogenicity: Scientific Rationale and Regulatory Frameworks

Knowledge of the genotoxicity and carcinogenicity potential of chemical substances is one of the key scientific elements able to better protect human health. Genotoxicity assessment is also considered as prescreening of carcinogenicity. The assessment of both endpoints is a fundamental component of national and international legislations, for all types of substances, and has stimulated the development of alternative, nontesting methods. Over the recent decades, much attention has been given to the use and further development of structure-activity relationships-based approaches, to be used in isolation or in combination with in vitro assays for predictive purposes. In this chapter, we briefly introduce the rationale for the main (Q)SAR approaches, and detail the most important regulatory initiatives and frameworks. It appears that the existence and needs of regulatory frameworks stimulate the development of better predictive tools; in turn, this allows the regulators to fine-tune their requirements for an improved defense of human health.

A Round Trip from Medicinal Chemistry to Predictive Toxicology

Predictive toxicology is a new emerging multifaceted research field aimed at protecting human health and environment from risks posed by chemicals. Such issue is of extreme public relevance and requires a multidisciplinary approach where the experience in medicinal chemistry is of utmost importance. Herein, we will survey some basic recommendations to gather good data and then will review three recent case studies to show how strategies of ligand- and structure-based molecular design, widely applied in medicinal chemistry, can be adapted to meet the more restrictive scientific and regulatory goals of predictive toxicology. In particular, we will report: Docking-based classification models to predict the estrogenic potentials of chemicals. Predicting the bioconcentration factor using biokinetics descriptors. Modeling oral sub-chronic toxicity using a customized k-nearest neighbors (k-NN) approach.

Predictive Structure-Based Toxicology Approaches To Assess the Androgenic Potential of Chemicals

We present a practical and easy-to-run in silico workflow exploiting a structure-based strategy making use of docking simulations to derive highly predictive classification models of the androgenic potential of chemicals. Models were trained on a high-quality chemical collection comprising 1689 curated compounds made available within the CoMPARA consortium from the US Environmental Protection Agency and were integrated with a two-step applicability domain whose implementation had the effect of improving both the confidence in prediction and statistics by reducing the number of false negatives. Among the nine androgen receptor X-ray solved structures, the crystal 2PNU (entry code from the Protein Data Bank) was associated with the best performing structure-based classification model. Three validation sets comprising each 2590 compounds extracted by the DUD-E collection were used to challenge model performance and the effectiveness of Applicability Domain implementation. Next, the 2PNU model was applied to screen and prioritize two collections of chemicals. The first is a small pool of 12 representative androgenic compounds that were accurately classified based on outstanding rationale at the molecular level. The second is a large external blind set of 55450 chemicals with potential for human exposure. We show how the use of molecular docking provides highly interpretable models and can represent a real-life option as an alternative nontesting method for predictive toxicology.

Ligand efficiency metrics in drug discovery: the pros and cons from a practical perspective

INTRODUCTION

Ligand efficiency metrics are almost universally accepted as a valuable indicator of compound quality and an aid to reduce attrition. Areas covered: In this review, the authors describe ligand efficiency metrics giving a balanced overview on their merits and points of weakness in order to enable the readers to gain an informed opinion. Relevant theoretical breakthroughs and drug-like properties are also illustrated. Several recent exemplary case studies are discussed in order to illustrate the main fields of application of ligand efficiency metrics. Expert opinion: As a medicinal chemist guide, ligand efficiency metrics perform in a context- and chemotype-dependent manner; thus, they should not be used as a magic box. Since the 'big bang' of efficiency metrics occurred more or less ten years ago and the average time to develop a new drug is over the same period, the next few years will give a clearer outlook on the increased rate of success, if any, gained by means of these new intriguing tools.

Chemical toxicity prediction for major classes of industrial chemicals: Is it possible to develop universal models covering cosmetics, drugs, and pesticides?

Computational models have earned broad acceptance for assessing chemical toxicity during early stages of drug discovery or environmental safety assessment. The majority of publicly available QSAR toxicity models have been developed for datasets including mostly drugs or drug-like compounds. We have evaluated and compared chemical spaces occupied by cosmetics, drugs, and pesticides, and explored whether current computational models of toxicity endpoints can be universally applied to all these chemicals. Our analysis of the chemical space overlap and applicability domain (AD) of models built previously for twenty different toxicity endpoints showed that most of these models afforded high coverage (>90%) for all three classes of compounds analyzed herein. Only T. pyriformis models demonstrated lower coverage for drugs and pesticides (38% and 54%, respectively). These results show that, for the most part, historical QSAR models built with data available for different toxicity endpoints can be used for toxicity assessment of novel chemicals irrespective of the intended commercial use; however, the AD restriction is necessary to assure the expected prediction accuracy. Local models may need to be developed to capture chemicals that appear as outliers with respect to global models.

Active learning for computational chemogenomics

Aim: Computational chemogenomics models the compound–protein interaction space, typically for drug discovery, where existing methods predominantly either incorporate increasing numbers of bioactivity samples or focus on specific subfamilies of proteins and ligands. As an alternative to modeling entire large datasets at once, active learning adaptively incorporates a minimum of informative examples for modeling, yielding compact but high quality models. Results/methodology: We assessed active learning for protein/target family-wide chemogenomic modeling by replicate experiment. Results demonstrate that small yet highly predictive models can be extracted from only 10–25% of large bioactivity datasets, irrespective of molecule descriptors used. Conclusion: Chemogenomic active learning identifies small subsets of ligand–target interactions in a large screening database that lead to knowledge discovery and highly predictive models.

Public (Q)SAR Services, Integrated Modeling Environments, and Model Repositories on the Web: State of the Art and Perspectives for Future Development

Thousands of (Quantitative) Structure-Activity Relationships (Q)SAR models have been described in peer-reviewed publications; however, this way of sharing seldom makes models available for the use by the research community outside of the developer's laboratory. Conversely, on-line models allow broad dissemination and application representing the most effective way of sharing the scientific knowledge. Approaches for sharing and providing on-line access to models range from web services created by individual users and laboratories to integrated modeling environments and model repositories. This emerging transition from the descriptive and informative, but "static", and for the most part, non-executable print format to interactive, transparent and functional delivery of "living" models is expected to have a transformative effect on modern experimental research in areas of scientific and regulatory use of (Q)SAR models.

Preselection of A- and B- modified d-homo lactone and d-seco androstane derivatives as potent compounds with antiproliferative activity against breast and prostate cancer cells - QSAR approach and molecular docking analysis

The problem with trial-and-error approach in organic synthesis of targeted anticancer compounds can be successfully avoided by computational modeling of molecules, docking studies and chemometric tools. It has been proven that A- and B- modified d-homo lactone and d-seco androstane derivatives are compounds with significant antiproliferative activity against estrogen-independent breast adenocarcinoma (ER-, MDA-MB-231) and androgen-independent prostate cancer cells (AR-, PC-3). This paper presents the quantitative structure-activity relationship (QSAR) models based on artificial neural networks (ANNs) which are able to predict whether d-homo lactone and/or d-seco androstane-based compounds will express antiproliferative activity against breast cancer cells (MDA-MB-231) or not. Also, the present paper describes the molecular docking study of 3β-acetoxy-5α,6α-epoxy- (3) and 6α,7α-epoxy-1,4-dien-3-one (24) d-homo lactone androstane derivatives, as well as 4-en-3-one (15) d-seco androstane derivative, which are compounds with strong or moderate antiproliferative activity against prostate cancer cells (PC-3), and compares them with commercially available medicament for prostate cancer - abiraterone. The obtained promising results can be used as guidelines in further syntheses of novel d-homo lactone and d-seco androstane derivatives with antiproliferative activity against breast and prostate cancer cells.

Mind the Gap! A Journey towards Computational Toxicology

Computational methods have advanced toxicology towards the development of target-specific models based on a clear cause-effect rationale. However, the predictive potential of these models presents strengths and weaknesses. On the good side, in silico models are valuable cheap alternatives to in vitro and in vivo experiments. On the other, the unconscious use of in silico methods can mislead end-users with elusive results. The focus of this review is on the basic scientific and regulatory recommendations in the derivation and application of computational models. Attention is paid to examine the interplay between computational toxicology and drug discovery and development. Avoiding the easy temptation of an overoptimistic future, we report our view on what can, or cannot, realistically be done. Indeed, studies of safety/toxicity represent a key element of chemical prioritization programs carried out by chemical industries, and primarily by pharmaceutical companies.

ADMET evaluation in drug discovery: 15. Accurate prediction of rat oral acute toxicity using relevance vector machine and consensus modeling

Background

Determination of acute toxicity, expressed as median lethal dose (LD₅₀), is one of the most important steps in drug discovery pipeline. Because in vivo assays for oral acute toxicity in mammals are time-consuming and costly, there is thus an urgent need to develop in silico prediction models of oral acute toxicity.

Results

In this study, based on a comprehensive data set containing 7314 diverse chemicals with rat oral LD₅₀ values, relevance vector machine (RVM) technique was employed to build the regression models for the prediction of oral acute toxicity in rate, which were compared with those built using other six machine learning approaches, including k-nearest-neighbor regression, random forest (RF), support vector machine, local approximate Gaussian process, multilayer perceptron ensemble, and eXtreme gradient boosting. A subset of the original molecular descriptors and structural fingerprints (PubChem or SubFP) was chosen by the Chi squared statistics. The prediction capabilities of individual QSAR models, measured by q_ext² for the test set containing 2376 molecules, ranged from 0.572 to 0.659.

Conclusion

Considering the overall prediction accuracy for the test set, RVM with Laplacian kernel and RF were recommended to build in silico models with better predictivity for rat oral acute toxicity. By combining the predictions from individual models, four consensus models were developed, yielding better prediction capabilities for the test set (q_ext² = 0.669–0.689). Finally, some essential descriptors and substructures relevant to oral acute toxicity were identified and analyzed, and they may be served as property or substructure alerts to avoid toxicity. We believe that the best consensus model with high prediction accuracy can be used as a reliable virtual screening tool to filter out compounds with high rat oral acute toxicity. Graphical abstract

Identification of structural alerts for liver and kidney toxicity using repeated dose toxicity data

BACKGROUND

The potential for a compound to cause hepatotoxicity and nephrotoxicity is a matter of extreme interest for human health risk assessment. To assess liver and kidney toxicity, repeated-dose toxicity (RDT) studies are conducted mainly on rodents. However, these tests are expensive, time-consuming and require large numbers of animals. For early toxicity screening, in silico models can be applied, reducing the costs, time and animals used. Among in silico approaches, structure-activity relationship (SAR) methods, based on the identification of chemical substructures (structural alerts, SAs) related to a particular activity (toxicity), are widely employed.

RESULTS

We identified and evaluated some SAs related to liver and kidney toxicity, using RDT data on rats taken from the hazard evaluation support system (HESS) database. We considered only SAs that gave the best percentages of true positives (TP).

CONCLUSIONS

It was not possible to assign an unambiguous mode of action for all the SAs, but a mechanistic explanation is provided for some of them. Such achievements may help in the early identification of liver and renal toxicity of substances.

Docking-based classification models for exploratory toxicology studies on high-quality estrogenic experimental data

BACKGROUND

The ethical and practical limitation of animal testing has recently promoted computational methods for the fast screening of huge collections of chemicals.

RESULTS

The authors derived 24 reliable docking-based classification models able to predict the estrogenic potential of a large collection of chemicals provided by the US Environmental Protection Agency. Model performances were challenged by considering AUC, EF1% (EFmax = 7.1), -LR (at sensitivity = 0.75); +LR (at sensitivity = 0.25) and 37 reference compounds comprised within the training set. Moreover, external predictions were made successfully on ten representative known estrogenic chemicals and on a set consisting of >32,000 chemicals.

CONCLUSION

The authors demonstrate that structure-based methods, widely applied to drug discovery programs, can be fairly adapted to exploratory toxicology studies.

Studies on the antiplatelet and antithrombotic profile of anti-inflammatory coumarin derivatives

The interest towards coumarin-based structures stems from their polypharmacological profile. Herein, we present a series of Mannich bases and 7-azomethine-linked coumarin derivatives exhibiting antiplatelet and antithrombotic activities, in addition to the already known anti-inflammatory and antioxidant activities. Among others, compounds 15 and 16 were found to be the most potent and selective inhibitors of platelet aggregation whereas compound 3 also proved to be the most potent in the clot retraction assay. Structure-activity relationship studies were conducted to elucidate the molecular determinants responsible for the herein observed activities. The chance of inhibiting cyclooxygenase-1 was also investigated for evaluating the platelet aggregation induced by arachidonic acid. Taken together, these results suggest that the investigation of other targets connected to the antiplatelet activity, such as phosphodiesterase-3 (PDE3), could be a viable strategy to shed light on the polypharmacological profile of coumarin-based compounds. Docking simulations towards PDE3 were also carried out.

Modeling of interactions between xenobiotics and cytochrome P450 (CYP) enzymes

The adverse effects to humans and environment of only few chemicals are well known. Absorption, distribution, metabolism, and excretion (ADME) are the steps of pharmaco/toxicokinetics that determine the internal dose of chemicals to which the organism is exposed. Of all the xenobiotic-metabolizing enzymes, the cytochrome P450 (CYP) enzymes are the most important due to their abundance and versatility. Reactions catalyzed by CYPs usually turn xenobiotics to harmless and excretable metabolites, but sometimes an innocuous xenobiotic is transformed into a toxic metabolite. Data on ADME and toxicity properties of compounds are increasingly generated using in vitro and modeling (in silico) tools. Both physics-based and empirical modeling approaches are used. Numerous ligand-based and target-based as well as combined modeling methods have been employed to evaluate determinants of CYP ligand binding as well as predicting sites of metabolism and inhibition characteristics of test molecules. In silico prediction of CYP–ligand interactions have made crucial contributions in understanding (1) determinants of CYP ligand binding recognition and affinity; (2) prediction of likely metabolites from substrates; (3) prediction of inhibitors and their inhibition potency. Truly predictive models of toxic outcomes cannot be created without incorporating metabolic characteristics; in silico methods help producing such information and filling gaps in experimentally derived data. Currently modeling methods are not mature enough to replace standard in vitro and in vivo approaches, but they are already used as an important component in risk assessment of drugs and other chemicals.

Biomolecular recognition of antagonists by α7 nicotinic acetylcholine receptor: Antagonistic mechanism and structure-activity relationships studies

As the key constituent of ligand-gated ion channels in the central nervous system, nicotinic acetylcholine receptors (nAChRs) and neurodegenerative diseases are strongly coupled in the human species. In recently years the developments of selective agonists by using nAChRs as the drug target have made a large progress, but the studies of selective antagonists are severely lacked. Currently these antagonists rest mainly on the extraction of partly natural products from some animals and plants; however, the production of these crude substances is quite restricted, and artificial synthesis of nAChR antagonists is still one of the completely new research fields. In the context of this manuscript, our primary objective was to comprehensively analyze the recognition patterns and the critical interaction descriptors between target α7 nAChR and a series of the novel compounds with potentially antagonistic activity by means of virtual screening, molecular docking and molecular dynamics simulation, and meanwhile these recognition reactions were also compared with the biointeraction of α7 nAChR with a commercially natural antagonist - methyllycaconitine. The results suggested clearly that there are relatively obvious differences of molecular structures between synthetic antagonists and methyllycaconitine, while the two systems have similar recognition modes on the whole. The interaction energy and the crucially noncovalent forces of the α7 nAChR-antagonists are ascertained according to the method of Molecular Mechanics/Generalized Born Surface Area. Several amino acid residues, such as B/Tyr-93, B/Lys-143, B/Trp-147, B/Tyr-188, B/Tyr-195, A/Trp-55 and A/Leu-118 played a major role in the α7 nAChR-antagonist recognition processes, in particular, residues B/Tyr-93, B/Trp-147 and B/Tyr-188 are the most important. These outcomes tally satisfactorily with the discussions of amino acid mutations. Based on the explorations of three-dimensional quantitative structure-activity relationships, the structure-antagonistic activity relationships of antagonists and the characteristics of α7 nAChR-ligand recognitions were received a reasonable summary as well. These attempts emerged herein would not only provide helpful guidance for the design of α7 nAChR antagonists, but shed new light on the subsequent researches in antagonistic mechanism.