A stepwise approach for defining the applicability domain of SAR and QSAR models

In silico assessment of nanoparticle toxicity powered by the Enalos Cloud Platform: Integrating automated machine learning and synthetic data for enhanced nanosafety evaluation

The rapid advance of nanotechnology has led to the development and widespread application of nanomaterials, raising concerns regarding their potential adverse effects on human health and the environment. Traditional (experimental) methods for assessing the nanoparticles (NPs) safety are time-consuming, expensive, and resource-intensive, and raise ethical concerns due to their reliance on animals. To address these challenges, we propose an in silico workflow that serves as an alternative or complementary approach to conventional hazard and risk assessment strategies, which incorporates state-of-the-art computational methodologies. In this study we present an automated machine learning (autoML) scheme that employs dose-response toxicity data for silver (Ag), titanium dioxide (TiO2), and copper oxide (CuO) NPs. This model is further enriched with atomistic descriptors to capture the NPs' underlying structural properties. To overcome the issue of limited data availability, synthetic data generation techniques are used. These techniques help in broadening the dataset, thus improving the representation of different NP classes. A key aspect of this approach is a novel three-step applicability domain method (which includes the development of a local similarity approach) that enhances user confidence in the results by evaluating the prediction's reliability. We anticipate that this approach will significantly expedite the nanosafety assessment process enabling regulation to keep pace with innovation, and will provide valuable insights for the design and development of safe and sustainable NPs. The ML model developed in this study is made available to the scientific community as an easy-to-use web-service through the Enalos Cloud Platform (www.enaloscloud.novamechanics.com/sabydoma/safenanoscope/), facilitating broader access and collaborative advancements in nanosafety.

Mechanistic QSAR modeling derived virtual screening, drug repurposing, ADMET and in-vitro evaluation to identify anticancer lead as lysine-specific demethylase 5a inhibitor

A lysine-specific demethylase is an enzyme that selectively eliminates methyl groups from lysine residues. KDM5A, also known as JARID1A or RBP2, belongs to the KDM5 Jumonji histone demethylase subfamily. To identify novel molecules that interact with the LSD5A receptor, we created a quantitative structure-activity relationship (QSAR) model. A group of 435 compounds was used in a study of the quantitative relationship between structure and activity to guess the IC50 values for blocking LASD5A. We used a genetic algorithm-multilinear regression-based quantitative structure-activity connection model to forecast the bioactivity (PIC50) of 1615 food and drug administration pharmaceuticals from the zinc database with the goal of repurposing clinically used medications. We used molecular docking, molecular dynamic simulation modelling, and molecular mechanics generalised surface area analysis to investigate the molecule's binding mechanism. A genetic algorithm and multi-linear regression method were used to make six variable-based quantitative structure-activity relationship models that worked well (R2 = 0.8521, Q2LOO = 0.8438, and Q2LMO = 0.8414). ZINC000000538621 was found to be a new hit against LSD5A after a quantitative structure-activity relationship-based virtual screening of 1615 zinc food and drug administration compounds. The docking analysis revealed that the hit molecule 11 in the KDM5A binding pocket adopted a conformation similar to the pdb-6bh1 ligand (docking score: -8.61 kcal/mol). The results from molecular docking and the quantitative structure-activity relationship were complementary and consistent. The most active lead molecule 11, which has shown encouraging results, has good absorption, distribution, metabolism, and excretion (ADME) properties, and its toxicity has been shown to be minimal. In addition, the MTT assay of ZINC000000538621 with MCF-7 cell lines backs up the in silico studies. We used molecular mechanics generalise borne surface area analysis and a 200-ns molecular dynamics simulation to find structural motifs for KDM5A enzyme interactions. Thus, our strategy will likely expand food and drug administration molecule repurposing research to find better anticancer drugs and therapies.Communicated by Ramaswamy H. Sarma.

Investigation of the Binding Fraction of PFAS in Human Plasma and Underlying Mechanisms Based on Machine Learning and Molecular Dynamics Simulation

More than 7000 per- and polyfluorinated alkyl substances (PFAS) have been documented in the U.S. Environmental Protection Agency's CompTox Chemicals database. These PFAS can be used in a broad range of industrial and consumer applications but may pose potential environmental issues and health risks. However, little is known about emerging PFAS bioaccumulation to assess their chemical safety. This study focuses specifically on the large and high-quality data set of fluorochemicals from the related environmental and pharmaceutical chemicals databases, and machine learning (ML) models were developed for the classification prediction of the unbound fraction of compounds in plasma. A comprehensive evaluation of the ML models shows that the best blending model yields an accuracy of 0.901 for the test set. The predictions suggest that most PFAS (∼92%) have a high binding fraction in plasma. Introduction of alkaline amino groups is likely to reduce the binding affinities of PFAS with plasma proteins. Molecular dynamics simulations indicate a clear distinction between the high and low binding fractions of PFAS. These computational workflows can be used to predict the bioaccumulation of emerging PFAS and are also helpful for the molecular design of PFAS to prevent the release of high-bioaccumulation compounds into the environment.

Machine learning models based on residue interaction network for ABCG2 transportable compounds recognition

As the one of the most important protein of placental transport of environmental substances, the identification of ABCG2 transport molecules is the key step for assessing the risk of placental exposure to environmental chemicals. Here, residue interaction network (RIN) was used to explore the difference of ABCG2 binding conformations between transportable and non-transportable compounds. The RIN were treated as a kind of special quantitative data of protein conformation, which not only reflected the changes of single amino acid conformation in protein, but also indicated the changes of distance and action type between amino acids. Based on the quantitative RIN, four machine learning algorithms were applied to establish the classification and recognition model for 1100 compounds with transported by ABCG2 potential. The random forest (RF) models constructed with RIN presented the best and satisfied predictive ability with an accuracy of training set of 0.97 and the test set of 0.96 respectively. In conclusion, the construction of residue interaction network provided a new perspective for the quantitative characterization of protein conformation and the establishment of prediction models for transporter molecular recognition. The ABCG2 transport molecular recognition model based on residue interaction network provides a possible way for screening environmental chemistry transported through placenta.

Merging bioactivity predictions from cell morphology and chemical fingerprint models using similarity to training data

The applicability domain of machine learning models trained on structural fingerprints for the prediction of biological endpoints is often limited by the lack of diversity of chemical space of the training data. In this work, we developed similarity-based merger models which combined the outputs of individual models trained on cell morphology (based on Cell Painting) and chemical structure (based on chemical fingerprints) and the structural and morphological similarities of the compounds in the test dataset to compounds in the training dataset. We applied these similarity-based merger models using logistic regression models on the predictions and similarities as features and predicted assay hit calls of 177 assays from ChEMBL, PubChem and the Broad Institute (where the required Cell Painting annotations were available). We found that the similarity-based merger models outperformed other models with an additional 20% assays (79 out of 177 assays) with an AUC > 0.70 compared with 65 out of 177 assays using structural models and 50 out of 177 assays using Cell Painting models. Our results demonstrated that similarity-based merger models combining structure and cell morphology models can more accurately predict a wide range of biological assay outcomes and further expanded the applicability domain by better extrapolating to new structural and morphology spaces.

Design of Potent Inhibitors Targeting the Main Protease of SARS-CoV-2 Using QSAR Modeling, Molecular Docking, and Molecular Dynamics Simulations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is a serious global public health threat. The evolving strains of SARS-CoV-2 have reduced the effectiveness of vaccines. Therefore, antiviral drugs against SARS-CoV-2 are urgently needed. The main protease (Mpro) of SARS-CoV-2 is an extremely potent target due to its pivotal role in virus replication and low susceptibility to mutation. In the present study, a quantitative structure-activity relationship (QSAR) study was performed to design new molecules that might have higher inhibitory activity against SARS-CoV-2 Mpro. In this context, a set of 55 dihydrophenanthrene derivatives was used to build two 2D-QSAR models using the Monte Carlo optimization method and the Genetic Algorithm Multi-Linear Regression (GA-MLR) method. From the CORAL QSAR model outputs, the promoters responsible for the increase/decrease in inhibitory activity were extracted and interpreted. The promoters responsible for an increase in activity were added to the lead compound to design new molecules. The GA-MLR QSAR model was used to ensure the inhibitory activity of the designed molecules. For further validation, the designed molecules were subjected to molecular docking analysis and molecular dynamics simulations along with an absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis. The results of this study suggest that the newly designed molecules have the potential to be developed as effective drugs against SARS-CoV-2.

Predicted health and environmental hazards of liquid crystal materials via quantitative structure-property relationship modeling

Liquid crystal materials (LCMs) are considered as emerging contaminants with high persistent and bioaccumulative potentials, but their toxicological effects are not well understood. To address this issue, a list of 1431 LCMs commercially available in the market was established through literature reviews and surveys of LCM suppliers. Toxicological properties of 221 target LCMs were derived from the Classification and Labeling Inventory by the European Chemicals Agency. More than 80 % of target LCMs likely pose adverse effects on human health or aquatic ecosystems. Two quantitative structure-property relationship (QSPR) models developed from the toxicological properties of LCMs achieved approximately 90 % accuracy in external data sets. The probability-based approach was more efficient in defining the applicability domain for the QSPR models than a range- or distance-based approach. The highest accuracy was achieved for chemicals within the probability-based applicability domain. The QSPR models were applied to predict health and environmental hazards of 1210 LCMs that had not been notified to the Classification and Labeling Inventory, and 301 and 94 LCMs were recognized as posing potential hazards to human health and the environment, respectively. The present study highlights the potential detrimental effects of LCMs and offers a specific in silico technique for screening hazardous LCMs.

Monte Carlo Method and GA-MLR-Based QSAR Modeling of NS5A Inhibitors against the Hepatitis C Virus

Hepatitis C virus (HCV) is a serious disease that threatens human health. Despite consistent efforts to inhibit the virus, it has infected more than 58 million people, with 300,000 deaths per year. The HCV nonstructural protein NS5A plays a critical role in the viral life cycle, as it is a major contributor to the viral replication and assembly processes. Therefore, its importance is evident in all currently approved HCV combination treatments. The present study identifies new potential compounds for possible medical use against HCV using the quantitative structure–activity relationship (QSAR). In this context, a set of 36 NS5A inhibitors was used to build QSAR models using genetic algorithm multiple linear regression (GA-MLR) and Monte Carlo optimization and were implemented in the software CORAL. The Monte Carlo method was used to build QSAR models using SMILES-based optimal descriptors. Four splits were performed and 24 QSAR models were developed and verified through internal and external validation. The model created for split 3 produced a higher value of the determination coefficients using the validation set (R2 = 0.991 and Q2 = 0.943). In addition, this model provides interesting information about the structural features responsible for the increase and decrease of inhibitory activity, which were used to develop eight novel NS5A inhibitors. The constructed GA-MLR model with satisfactory statistical parameters (R2 = 0.915 and Q2 = 0.941) confirmed the predicted inhibitory activity for these compounds. The Absorption, Distribution, Metabolism, Elimination, and Toxicity (ADMET) predictions showed that the newly designed compounds were nontoxic and exhibited acceptable pharmacological properties. These results could accelerate the process of discovering new drugs against HCV.

High-Throughput Assessment of the Abiotic Stability of Test Chemicals in In Vitro Bioassays

Abiotic stability of chemicals is not routinely tested prior to performing in vitro bioassays, although abiotic degradation can reduce the concentration of test chemicals leading to the formation of active or inactive transformation products, which may lead to misinterpretation of bioassay results. A high-throughput workflow was developed to measure the abiotic stability of 22 test chemicals in protein-rich aqueous media under typical bioassay conditions at 37 °C for 48 h. These test chemicals were degradable in the environment according to a literature review. The chemicals were extracted from the exposure media at different time points using a novel 96-pin solid-phase microextraction. The conditions were varied to differentiate between various reaction mechanisms. For most hydrolyzable chemicals, pH-dependent degradation in phosphate-buffered saline indicated that acid-catalyzed hydrolysis was less important than reactions with hydroxide ions. Reactions with proteins were mainly responsible for the depletion of the test chemicals in the media, which was simulated by bovine serum albumin (BSA) and glutathione (GSH). 1,2-Benzisothiazol-3(2H)-one, 2-methyl-4-isothiazolinone, and l-sulforaphane reacted almost instantaneously with GSH but not with BSA, indicating that GSH is a good proxy for reactivity with electrophilic amino acids but may overestimate the actual reaction with three-dimensional proteins. Chemicals such as hydroquinones or polyunsaturated chemicals are prone to autoxidation, but this reaction is difficult to differentiate from hydrolysis and could not be simulated by the oxidant N-bromosuccinimide. Photodegradation played a minor role because cells are exposed in incubators in the dark and simulations with high light intensities did not yield realistic degradation. Stability predictions from various in silico prediction models for environmental conditions can give initial indications of the stability but were not always consistent with the experimental stability in bioassays. As the presented workflow can be performed in high throughput under realistic bioassay conditions, it can be used to provide an experimental database for developing bioassay-specific stability prediction models.

Automated read-across workflow for predicting acute oral toxicity: I. The decision scheme in the QSAR toolbox

A decision-scheme outlining the steps for identifying the appropriate chemical category and subsequently appropriate tested source analog(s) for data gap filling of a target chemical by read-across is described. The primary features used in the grouping of the target chemical with source analogues within a database of 10,039 discrete organic substances include reactivity mechanisms associated with protein interactions and specific-acute-oral-toxicity-related mechanisms (e.g., mitochondrial uncoupling). Additionally, the grouping of chemicals making use of the in vivo rat metabolic simulator and neutral hydrolysis. Subsequently, a series of structure-based profilers are used to narrow the group to the most similar analogues. The scheme is implemented in the OECD QSAR Toolbox, so it automatically predicts acute oral toxicity as the rat oral LD50 value in log [1/mol/kg]. It was demonstrated that due to the inherent variability in experimental data, classification distribution should be employed as more adequate in comparison to the exact classification. It was proved that the predictions falling in the adjacent GSH categories to the experimentally-stated ones are acceptable given the variation in experimental data. The model performance estimated by adjacent accuracy was found to be 0.89 and 0.54 while based on R². The mechanistic and predictive coverages were >0.85.

Predictive Models to Identify Small Molecule Activators and Inhibitors of Opioid Receptors

Opioid receptors (OPRs) are the main targets for the treatment of pain and related disorders. The opiate compounds that activate these receptors are effective analgesics but their use leads to adverse effects, and they often are highly addictive drugs of abuse. There is an urgent need for alternative chemicals that are analgesics and to reduce/avoid the unwanted effects in order to relieve the public health crisis of opioid addiction. Here, we aim to develop computational models to predict the OPR activity of small molecule compounds based on chemical structures and apply these models to identify novel OPR active compounds. We used four different machine learning algorithms to build models based on quantitative high throughput screening (qHTS) data sets of three OPRs in both agonist and antagonist modes. The best performing models were applied to virtually screen a large collection of compounds. The model predicted active compounds were experimentally validated using the same qHTS assays that generated the training data. Random forest was the best classifier with the highest performance metrics, and the mu OPR (OPRM)-agonist model achieved the best performance measured by AUC-ROC (0.88) and MCC (0.7) values. The model predicted actives resulted in hit rates ranging from 2.3% (delta OPR-agonist) to 15.8% (OPRM-agonist) after experimental confirmation. Compared to the original assay hit rate, all models enriched the hit rate by ≥2-fold. Our approach produced robust OPR prediction models that can be applied to prioritize compounds from large libraries for further experimental validation. The models identified several novel potent compounds as activators/inhibitors of OPRs that were confirmed experimentally. The potent hits were further investigated using molecular docking to find the interactions of the novel ligands in the active site of the corresponding OPR.

Trends in predictive biodegradation for sustainable mitigation of environmental pollutants: Recent progress and future outlook

The feasibility of in-silico techniques, together with the computational framework, has been applied to predictive bioremediation aiming to clean-up contaminants, toxicity evaluation, and possibilities for the degradation of complex recalcitrant compounds. Emerging contaminants from different industries have posed a significant hazard to the environment and public health. Given current bioremediation strategies, it is often a failure or inadequate for sustainable mitigation of hazardous pollutants. However, clear-cut vital information about biodegradation is quite incomplete from a conventional remediation techniques perspective. Lacking complete information on bio-transformed compounds leads to seeking alternative methods. Only scarce information about the transformed products and toxicity profile is available in the published literature. To fulfill this literature gap, various computational or in-silico technologies have emerged as alternating techniques, which are being recognized as in-silico approaches for bioremediation. Molecular docking, molecular dynamics simulation, and biodegradation pathways predictions are the vital part of predictive biodegradation, including the Quantitative Structure-Activity Relationship (QSAR), Quantitative structure-biodegradation relationship (QSBR) model system. Furthermore, machine learning (ML), artificial neural network (ANN), genetic algorithm (GA) based programs offer simultaneous biodegradation prediction along with toxicity and environmental fate prediction. Herein, we spotlight the feasibility of in-silico remediation approaches for various persistent, recalcitrant contaminants while traditional bioremediation fails to mitigate such pollutants. Such could be addressed by exploiting described model systems and algorithm-based programs. Furthermore, recent advances in QSAR modeling, algorithm, and dedicated biodegradation prediction system have been summarized with unique attributes.

HuskinDB, a database for skin permeation of xenobiotics

Skin permeation is an essential biological property of small organic compounds our body is exposed to, such as drugs in topic formulations, cosmetics, and environmental toxins. Despite the limited availability of experimental data, there is a lack of systematic analysis and structure. We present a novel resource on skin permeation data that collects all measurements available in the literature and systematically structures experimental conditions. Besides the skin permeation value k_p, it includes experimental protocols such as skin source site, skin layer used, preparation technique, storage conditions, as well as test conditions such as temperature, pH as well as the type of donor and acceptor solution. It is important to include these parameters in the assessment of the skin permeation data. In addition, we provide an analysis of physicochemical properties and chemical space coverage, laying the basis for applicability domain determination of insights drawn from the collected data points. The database is freely accessible under https://huskindb.drug-design.de or https://doi.org/10.7303/syn21998881 .

Comparison of the predictive nature of the Genomic Allergen Rapid Detection (GARD) assay with mammalian assays in determining the skin sensitisation potential of agrochemical active ingredients

Alternatives to mammalian testing are highly desirable to predict the skin sensitisation potential of agrochemical active ingredients (AI). The GARD assay, a stimulated, dendritic cell-like, cell line measuring genomic signatures, was evaluated using twelve AIs (seven sensitisers and five non-sensitisers) and the results compared with historical results from guinea pig or local lymph node assay (LLNA) studies. Initial GARD results suggested 11/12 AIs were sensitisers and six concurred with mammalian data. Conformal predictions changed one AI to a non-sensitiser. An AI identified as non-sensitising in the GARD assay was considered a potent sensitiser in the LLNA. In total 7/12 GARD results corresponded with mammalian data. AI chemistries might not be comparable to the GARD training set in terms of applicability domains. Whilst the GARD assay can replace mammalian tests for skin sensitisation evaluation for compounds including cosmetic ingredients, further work in agrochemical chemistries is needed for this assay to be a viable replacement to animal testing. The work conducted here is, however, considered exploratory research and the methodology needs further development to be validated for agrochemicals. Mammalian and other alternative assays for regulatory safety assessments of AIs must provide confidence to assign the appropriate classification for human health protection.

Successive Statistical and Structure-Based Modeling to Identify Chemically Novel Kinase Inhibitors

Kinases are frequently studied in the context of anticancer drugs. Their involvement in cell responses, such as proliferation, differentiation, and apoptosis, makes them interesting subjects in multitarget drug design. In this study, a workflow is presented that models the bioactivity spectra for two panels of kinases: (1) inhibition of RET, BRAF, SRC, and S6K, while avoiding inhibition of MKNK1, TTK, ERK8, PDK1, and PAK3, and (2) inhibition of AURKA, PAK1, FGFR1, and LKB1, while avoiding inhibition of PAK3, TAK1, and PIK3CA. Both statistical and structure-based models were included, which were thoroughly benchmarked and optimized. A virtual screening was performed to test the workflow for one of the main targets, RET kinase. This resulted in 5 novel and chemically dissimilar RET inhibitors with remaining RET activity of <60% (at a concentration of 10 μM) and similarities with known RET inhibitors from 0.18 to 0.29 (Tanimoto, ECFP6). The four more potent inhibitors were assessed in a concentration range and proved to be modestly active with a pIC₅₀ value of 5.1 for the most active compound. The experimental validation of inhibitors for RET strongly indicates that the multitarget workflow is able to detect novel inhibitors for kinases, and hence, this workflow can potentially be applied in polypharmacology modeling. We conclude that this approach can identify new chemical matter for existing targets. Moreover, this workflow can easily be applied to other targets as well.

Novel and Predictive QSAR Model for Steroidal and Nonsteroidal 5α- Reductase Type II Inhibitors

AIMS AND OBJECTIVE

In this study, a novel quantitative structure activity relationship (QSAR) model has been developed for inhibitors of human 5-alpha reductase type II, which are used to treat benign prostate hypertrophy (BPH).

METHODS

The dataset consisted of 113 compounds-mainly nonsteroidal-with known inhibitory concentration. Then 3D structures of compounds were optimized and molecular structure descriptors were calculated. The stepwise multiple linear regression was used to select descriptors encoding the inhibitory activity of the compounds. Multiple linear regression (MLR) was used to build up the linear QSAR model.

RESULTS

The results obtained revealed that the descriptors which best describe the activity were atom type electropological state, carbon type, radial distribution function (RDF), barysz matrix and molecular linear free energy relation. The suggested model could achieve satisfied square correlation coefficient of R² = 0.72, higher than of many previous studies, indicating its superiority. Rigid validation criteria were met using external data with Q² ˃ 0.5 and R² = 0.75, reflecting the predictive power of the model.

CONCLUSION

The QSAR model was applied for screening botanical components of herbal preparations used to treat BPH, and could predict the activity of some, among others, making reasonable attribution to the proposed effect of these preparations. Gamma tocopherol was found to be an active inhibitor, in consistence with many previous studies, anticipating the power of this model in the prediction of new candidate molecules and suggesting further investigations.

Whether the Validation of the Predictive Potential of Toxicity Models is a Solved Task?

Different kinds of biological activities are defined by complex biochemical interactions, which are termed as a "mathematical function" not only of the molecular structure but also for some additional circumstances, such as physicochemical conditions, interactions via energy and information effects between a substance and organisms, organs, cells. These circumstances lead to the great complexity of prediction for biochemical endpoints, since all "details" of corresponding phenomena are practically unavailable for the accurate registration and analysis. Researchers have not a possibility to carry out and analyse all possible ways of the biochemical interactions, which define toxicological or therapeutically attractive effects via direct experiment. Consequently, a compromise, i.e. the development of predictive models of the above phenomena, becomes necessary. However, the estimation of the predictive potential of these models remains a task that is solved only partially. This mini-review presents a collection of attempts to be used for the above-mentioned task, two special statistical indices are proposed, which may be a measure of the predictive potential of models. These indices are (i) Index of Ideality of Correlation; and (ii) Correlation Contradiction Index.

Prediction and Experimental Confirmation of Novel Peripheral Cannabinoid-1 Receptor Antagonists

Small molecules targeting peripheral CB1 receptors have therapeutic potential in a variety of disorders including obesity-related, hormonal, and metabolic abnormalities, while avoiding the psychoactive effects in the central nervous system. We applied our in-house algorithm, iterative stochastic elimination, to produce a ligand-based model that distinguishes between CB1R antagonists and random molecules by physicochemical properties only. We screened ∼2 million commercially available molecules and found that about 500 of them are potential candidates to antagonize the CB1R. We applied a few criteria for peripheral activity and narrowed that set down to 30 molecules, out of which 15 could be purchased. Ten out of those 15 showed good affinity to the CB1R and two of them with nanomolar affinities (K_i of ∼400 nM). The eight molecules with top affinities were tested for activity: two compounds were pure antagonists, and five others were inverse agonists. These molecules are now being examined in vivo for their peripheral versus central distribution and subsequently will be tested for their effects on obesity in small animals.

Quantitative structure-activity relationship to predict the anti-malarial activity in a set of new imidazolopiperazines based on artificial neural networks

Background

After years of efforts on the control of malaria, it remains as a most deadly infectious disease. A major problem for the available anti-malarial drugs is the occurrence of drug resistance in Plasmodium. Developing of new compounds or modification of existing anti-malarial drugs is an effective approach to face this challenge. Quantitative structure activity relationship (QSAR) modelling plays an important role in design and modification of anti-malarial compounds by estimation of the activity of the compounds.

Methods

In this research, the QSAR study was done on anti-malarial activity of 33 imidazolopiperazine compounds based on artificial neural networks (ANN). The structural descriptors of imidazolopiperazine molecules was used as the independents variables and their activity against 3D7 and W2 strains was used as the dependent variables. During modelling process, 70% of compound was used as the training and two 15% of imidazolopiperazines were used as the validation and external test sets. In this work, stepwise multiple linear regression was applied as the valuable selection and ANN with Levenberg–Marquardt algorithm was utilized as an efficient non-linear approach to correlate between structural information of molecules and their anti-malarial activity.

Results

The sufficiency of the suggested method to estimate the anti-malarial activity of imidazolopiperazine compounds at two 3D7 and W2 strains was demonstrated using statistical parameters, such as correlation coefficient (R²), mean square error (MSE). For instance R²_train = 0.947, R²_val = 0.959, R²_test = 0.920 shows the potential of the suggested model for the prediction of 3D7 activity. Different statistical approaches such as and applicability domain (AD) and y-scrambling was also showed the validity of models.

Conclusion

QSAR can be an efficient way to virtual screening the molecules to design more efficient compounds with activity against malaria (3D7 and W2 strains). Imidazolopiperazines can be good candidates and change in the structure and functional groups can be done intelligently using QSAR approach to rich more efficient compounds with decreasing trial–error runs during synthesis.

Mechanistic relationship between biodegradation and bioaccumulation. Practical outcomes

According to the REACH Regulation, for all substances manufactured or imported in amounts of 10 or more tons per year, that are not exempted from the registration requirement, a Chemical Safety Assessment (CSA) must be conducted. According to CSA criteria, for these substances persistent, bioaccumulative and toxic (PBT), and very persistent and very bioaccumulative (vPvB) assessment is requested. In order to reduce the number of applications of the expensive bioaccumulation test it seems useful to search thresholds for other related parameters above which no bioaccumulation is observed. Given the known relationship between ready biodegradability and bioaccumulation, one such parameter is biodegradation. This article addresses this relationship in searching for BOD threshold above which no vB and B chemicals could be observed. It was found that the regulatory criteria for persistency could be used for identification of not vB and B chemicals. In addition, fish liver metabolism is determined as the most significant factor in reducing of maximum bioaccumulation potential of the chemicals. It was found that parameters associated with the models simulating fish metabolism could be also used for identification of not vB and B chemicals.

Compartment-Specific Screening Tools for Persistence: Potential Role and Application in the Regulatory Context

The persistence assessment under the European Union regulation Registration, Evaluation, Authorisation and Restriction of Chemicals (REACH) relies on compartment-specific degradation half-lives derived from laboratory simulation studies with surface water, aquatic sediment, or soil. Although these data are given priority, they are not available for most of the compounds. Therefore, according to the Integrated Assessment and Testing Strategy (ITS) for persistence assessment, results from ready biodegradability tests (RBTs) are used within a persistence screening to decide whether a substance is considered as "not persistent" or "potentially persistent." However, ready biodegradability is currently tested only in water. Consequently, there is a lack of approaches that include the soil and sediment compartments for persistence assessment at the screening level. In previous studies, compartment-specific screening tools for water-sediment (Water-Sediment Screening Tool [WSST]) and soil (Soil Screening Tool [SST]) were developed based on the existing test guideline Organisation for Economic Development and Co-operation (OECD TG 301C [MITI (Ministry of International Trade and Industry, Japan) test]). The test systems MITI, WSST, and SST were successfully applied to determine sound and reliable biodegradation data for 15 test compounds. In the present study, these results are used within the scope of a new alternative persistence screening approach, the Compartment-Specific Persistence Screening (CSPS). Compared to the persistence screening under REACH, the CSPS is a more conservative approach that provides additional reasonable results, particularly for compounds that sorb to sediment and soil, and for which the current standard persistence screening might be insufficient. Thus, the CSPS can be used to identify potentially persistent and nonpersistent compounds in the regulatory context by a comprehensive assessment that includes water, soil, and sediment. Moreover, experimentally determined half-lives from the compartment-specific screening tools can be used as input for multimedia models that estimate, for example, overall persistence (P_ov ). The application of fixed half-life factors to extrapolate from water to soil and sediment, which is here demonstrated to be inappropriate, can thereby be avoided. Integr Environ Assess Manag 2019;00:000-000. © 2019 SETAC.

Validation of the performance of TIMES genotoxicity models with EFSA pesticide data

This study validates the performance of the TIssue MEtabolism Simulator (TIMES) genotoxicity models with data on pesticide chemicals included in a recently released European Food Safety Authority (EFSA) genotoxicity database. The EFSA database is biased towards negative chemicals. A comparison of substances included in the EFSA database and TIMES genotoxicity databases showed that the majority of the EFSA pesticides is not included in the TIMES genotoxicity databases and, thus, out of the applicability domains of the current TIMES models. However, the EFSA genotoxicity database provides an opportunity to expand the TIMES models. Where there is overlap of substances, consistency between EFSA and TIMES databases for the chemicals with documented data is found to be high (>80%) with respect to the Ames data and lower than the Ames data with respect to chromosomal aberration (CA) and mouse lymphoma assay (MLA) data. No conclusion for consistency with respect to micronucleus test and comet genotoxicity data can be provided due to the limited number of overlapping substances. Specificity of the models is important, given the prevalence of negative genotoxicity data in the EFSA database. High specificity (>80%) is obtained for prediction of the EFSA pesticides with Ames data. Moreover, this high specificity of the TIMES Ames models is not dependant on pesticides being within the domains. Specificity of the TIMES CA and MLA models is lower (>40%) to pesticides for out of domain. Sensitivity of TIMES in vitro and in vivo models cannot be properly estimated due to the small number of positive chemicals in the EFSA database.

A review on machine learning methods for in silico toxicity prediction

In silico toxicity prediction plays an important role in the regulatory decision making and selection of leads in drug design as in vitro/vivo methods are often limited by ethics, time, budget, and other resources. Many computational methods have been employed in predicting the toxicity profile of chemicals. This review provides a detailed end-to-end overview of the application of machine learning algorithms to Structure-Activity Relationship (SAR)-based predictive toxicology. From raw data to model validation, the importance of data quality is stressed as it greatly affects the predictive power of derived models. Commonly overlooked challenges such as data imbalance, activity cliff, model evaluation, and definition of applicability domain are highlighted, and plausible solutions for alleviating these challenges are discussed.

QSAR Development for Plasma Protein Binding: Influence of the Ionization State

Purpose

This study explored several strategies to improve the performance of literature QSAR models for plasma protein binding (PPB), such as a suitable endpoint transformation, a correct representation of chemicals, more consistency in the dataset, and a reliable definition of the applicability domain.

Methods

We retrieved human fraction unbound (Fu) data for 670 compounds from the literature and carefully checked them for consistency. Descriptors were calculated taking account of the ionization state of molecules at physiological pH (7.4), in order to better estimate the affinity of molecules to blood proteins. We used different algorithms and chemical descriptors to explore the most suitable strategy for modeling the endpoint. SMILES (simplified molecular input line entry system)-based string descriptors were also tested with the CORAL software (CORelation And Logic). We did an outlier analysis to establish the models to use (or not to use) in case of well recognized families.

Results

Internal validation of the selected models returned Q² values close to 0.60. External validation also gave r² values always greater than 0.60. The CORAL descriptor based model for √fu was the best, with r² 0.74 in external validation.

Conclusions

Performance in prediction confirmed the robustness of all the derived models and their suitability for real-life purposes, i.e. screening chemicals for their ADMET profiling. Optimization of descriptors can be useful in order to obtain the correct results with a ionized molecule.

Electronic supplementary material

The online version of this article (10.1007/s11095-018-2561-8) contains supplementary material, which is available to authorized users.

Study of the Applicability Domain of the QSAR Classification Models by Means of the Rivality and Modelability Indexes

The reliability of a QSAR classification model depends on its capacity to achieve confident predictions of new compounds not considered in the building of the model. The results of this external validation process show the applicability domain (AD) of the QSAR model and, therefore, the robustness of the model to predict the property/activity of new molecules. In this paper we propose the use of the rivality and modelability indexes for the study of the characteristics of the datasets to be correctly modeled by a QSAR algorithm and to predict the reliability of the built model to prognosticate the property/activity of new molecules. The calculation of these indexes has a very low computational cost, not requiring the building of a model, thus being good tools for the analysis of the datasets in the first stages of the building of QSAR classification models. In our study, we have selected two benchmark datasets with similar number of molecules but with very different modelability and we have corroborated the capacity of the predictability of the rivality and modelability indexes regarding the classification models built using Support Vector Machine and Random Forest algorithms with 5-fold cross-validation and leave-one-out techniques. The results have shown the excellent ability of both indexes to predict outliers and the applicability domain of the QSAR classification models. In all cases, these values accurately predicted the statistic parameters of the QSAR models generated by the algorithms.

Data Visualization, Regression, Applicability Domains and Inverse Analysis Based on Generative Topographic Mapping

This paper introduces two generative topographic mapping (GTM) methods that can be used for data visualization, regression analysis, inverse analysis, and the determination of applicability domains (ADs). In GTM-multiple linear regression (GTM-MLR), the prior probability distribution of the descriptors or explanatory variables (X) is calculated with GTM, and the posterior probability distribution of the property/activity or objective variable (y) given X is calculated with MLR; inverse analysis is then performed using the product rule and Bayes' theorem. In GTM-regression (GTMR), X and y are combined and GTM is performed to obtain the joint probability distribution of X and y; this leads to the posterior probability distributions of y given X and of X given y, which are used for regression and inverse analysis, respectively. Simulations using linear and nonlinear datasets and quantitative structure-activity relationship (QSAR) and quantitative structure-property relationship (QSPR) datasets confirm that GTM-MLR and GTMR enable data visualization, regression analysis, and inverse analysis considering appropriate ADs. Python and MATLAB codes for the proposed algorithms are available at https://github.com/hkaneko1985/gtm-generativetopographicmapping.

Computational approaches to chemical hazard assessment

Computational prediction of toxicity has reached new heights as a result of decades of growth in the magnitude and diversity of biological data. Public packages for statistics and machine learning make model creation faster. New theory in machine learning and cheminformatics enables integration of chemical structure, toxicogenomics, simulated and physical data in the prediction of chemical health hazards, and other toxicological information. Our earlier publications have characterized a toxicological dataset of unprecedented scale resulting from the European REACH legislation (Registration Evaluation Authorisation and Restriction of Chemicals). These publications dove into potential use cases for regulatory data and some models for exploiting this data. This article analyzes the options for the identification and categorization of chemicals, moves on to the derivation of descriptive features for chemicals, discusses different kinds of targets modeled in computational toxicology, and ends with a high-level perspective of the algorithms used to create computational toxicology models.

An evaluation of selected (Q)SARs/expert systems for predicting skin sensitisation potential

Predictive testing to characterise substances for their skin sensitisation potential has historically been based on animal models such as the Local Lymph Node Assay (LLNA) and the Guinea Pig Maximisation Test (GPMT). In recent years, EU regulations, have provided a strong incentive to develop non-animal alternatives, such as expert systems software. Here we selected three different types of expert systems: VEGA (statistical), Derek Nexus (knowledge-based) and TIMES-SS (hybrid), and evaluated their performance using two large sets of animal data: one set of 1249 substances from eChemportal and a second set of 515 substances from NICEATM. A model was considered successful at predicting skin sensitisation potential if it had at least the same balanced accuracy as the LLNA and the GPMT had in predicting the other outcomes, which ranged from 79% to 86%. We found that the highest balanced accuracy of any of the expert systems evaluated was 65% when making global predictions. For substances within the domain of TIMES-SS, however, balanced accuracies for the two datasets were found to be 79% and 82%. In those cases where a chemical was within the TIMES-SS domain, the TIMES-SS skin sensitisation hazard prediction had the same confidence as the result from LLNA or GPMT.

PBT assessment under REACH: Screening for low aquatic bioaccumulation with QSAR classifications based on physicochemical properties to replace BCF in vivo testing on fish

Aquatic bioconcentration factors (BCFs) are critical in PBT (persistent, bioaccumulative, toxic) and risk assessment of chemicals. High costs and use of more than 100 fish per standard BCF study (OECD 305) call for alternative methods to replace as much in vivo testing as possible. The BCF waiving scheme is a screening tool combining QSAR classifications based on physicochemical properties related to the distribution (hydrophobicity, ionisation), persistence (biodegradability, hydrolysis), solubility and volatility (Henry's law constant) of substances in water bodies and aquatic biota to predict substances with low aquatic bioaccumulation (nonB, BCF<2000). The BCF waiving scheme was developed with a dataset of reliable BCFs for 998 compounds and externally validated with another 181 substances. It performs with 100% sensitivity (no false negatives), >50% efficacy (waiving potential), and complies with the OECD principles for valid QSARs. The chemical applicability domain of the BCF waiving scheme is given by the structures of the training set, with some compound classes explicitly excluded like organometallics, poly- and perfluorinated compounds, aromatic triphenylphosphates, surfactants. The prediction confidence of the BCF waiving scheme is based on applicability domain compliance, consensus modelling, and the structural similarity with known nonB and B/vB substances. Compounds classified as nonB by the BCF waiving scheme are candidates for waiving of BCF in vivo testing on fish due to low concern with regard to the B criterion. The BCF waiving scheme supports the 3Rs with a possible reduction of >50% of BCF in vivo testing on fish. If the target chemical is outside the applicability domain of the BCF waiving scheme or not classified as nonB, further assessments with in silico, in vitro or in vivo methods are necessary to either confirm or reject bioaccumulative behaviour.

Discussion on Regression Methods Based on Ensemble Learning and Applicability Domains of Linear Submodels

To develop a new ensemble learning method and construct highly predictive regression models in chemoinformatics and chemometrics, applicability domains (ADs) are introduced into the ensemble learning process of prediction. When estimating values of an objective variable using subregression models, only the submodels with ADs that cover a query sample, i.e., the sample is inside the model's AD, are used. By constructing submodels and changing a list of selected explanatory variables, the union of the submodels' ADs, which defines the overall AD, becomes large, and the prediction performance is enhanced for diverse compounds. By analyzing a quantitative structure-activity relationship data set and a quantitative structure-property relationship data set, it is confirmed that the ADs can be enlarged and the estimation performance of regression models is improved compared with traditional methods.

Molecular Descriptors for Structure-Activity Applications: A Hands-On Approach

Molecular descriptors capture diverse parts of the structural information of molecules and they are the support of many contemporary computer-assisted toxicological and chemical applications. After briefly introducing some fundamental concepts of structure-activity applications (e.g., molecular descriptor dimensionality, classical vs. fingerprint description, and activity landscapes), this chapter guides the readers through a step-by-step explanation of molecular descriptors rationale and application. To this end, the chapter illustrates a case study of a recently published application of molecular descriptors for modeling the activity on cytochrome P450.

Applicability Domain: A Step Toward Confident Predictions and Decidability for QSAR Modeling

In the context of human safety assessment through quantitative structure-activity relationship (QSAR) modeling, the concept of applicability domain (AD) has an enormous role to play. The Organization of Economic Co-operation and Development (OECD) for QSAR model validation recommended as principle 3 "A defined domain of applicability" to be present for a predictive QSAR model. The study of AD allows estimating the uncertainty in the prediction for a particular molecule based on how similar it is to the training compounds which are used in the model development. In the current scenario, AD represents an active research topic, and many methods have been designed to estimate the competence of a model and the confidence in its outcome for a given prediction task. Thus, characterization of interpolation space is significant in defining the AD. The diverse set of reported AD methods was constructed through different hypotheses and algorithms. These multiplicities of methodologies mystify the end users and make the comparison of the AD for different models a complex issue to address. We have attempted to summarize in this chapter the important concepts of AD including particulars of the available methods to compute the AD along with their thresholds and criteria for estimating AD through training set interpolation in the descriptor space. The idea about transparent domain and decision domain are also discussed. To help readers determine the AD in their projects, practical examples together with available open source software tools are provided.

CATALOGIC 301C model - validation and improvement

In Europe, REACH legislation encourages the use of alternative in silico methods such as (Q)SAR models. According to the recent progress of Chemical Substances Control Law (CSCL) in Japan, (Q)SAR predictions are also utilized as supporting evidence for the assessment of bioaccumulation potential of chemicals along with read across. Currently, the effective use of read across and QSARs is examined for other hazards, including biodegradability. This paper describes the results of external validation and improvement of CATALOGIC 301C model based on more than 1000 tested new chemical substances of the publication schedule under CSCL. CATALOGIC 301C model meets all REACH requirements to be used for biodegradability assessment. The model formalism built on scientific understanding for the microbial degradation of chemicals has a well-defined and transparent applicability domain. The model predictions are adequate for the evaluation of the ready degradability of chemicals.

Classification of baseline toxicants for QSAR predictions to replace fish acute toxicity studies

Fish acute toxicity studies are required for environmental hazard and risk assessment of chemicals by national and international legislations such as REACH, the regulations of plant protection products and biocidal products, or the GHS (globally harmonised system) for classification and labelling of chemicals. Alternative methods like QSARs (quantitative structure-activity relationships) can replace many ecotoxicity tests. However, complete substitution of in vivo animal tests by in silico methods may not be realistic. For the so-called baseline toxicants, it is possible to predict the fish acute toxicity with sufficient accuracy from log K_ow and, hence, valid QSARs can replace in vivo testing. In contrast, excess toxicants and chemicals not reliably classified as baseline toxicants require further in silico, in vitro or in vivo assessments. Thus, the critical task is to discriminate between baseline and excess toxicants. For fish acute toxicity, we derived a scheme based on structural alerts and physicochemical property thresholds to classify chemicals as either baseline toxicants (=predictable by QSARs) or as potential excess toxicants (=not predictable by baseline QSARs). The step-wise approach identifies baseline toxicants (true negatives) in a precautionary way to avoid false negative predictions. Therefore, a certain fraction of false positives can be tolerated, i.e. baseline toxicants without specific effects that may be tested instead of predicted. Application of the classification scheme to a new heterogeneous dataset for diverse fish species results in 40% baseline toxicants, 24% excess toxicants and 36% compounds not classified. Thus, we can conclude that replacing about half of the fish acute toxicity tests by QSAR predictions is realistic to be achieved in the short-term. The long-term goals are classification criteria also for further groups of toxicants and to replace as many in vivo fish acute toxicity tests as possible with valid QSAR predictions.

Importance of Applicability Domain of QSAR Models

Quantitative Structure-Activity Relationship (QSAR) models have manifold applications in drug discovery, environmental fate modeling, risk assessment, and property prediction of chemicals and pharmaceuticals. One of the principles recommended by the Organization of Economic Co-operation and Development (OECD) for model validation requires defining the Applicability Domain (AD) for QSAR models, which allows one to estimate the uncertainty in the prediction of a compound based on how similar it is to the training compounds, which are used in the model development. The AD is a significant tool to build a reliable QSAR model, which is generally limited in use to query chemicals structurally similar to the training compounds. Thus, characterization of interpolation space is significant in defining the AD. An attempt is made in this chapter to address the important concepts and methodology of the AD as well as criteria for estimating AD through training set interpolation in the descriptor space.

QSAR-Based Studies of Nanomaterials in the Environment

Nanotechnology is a newly emerging field, posing substantial impacts on society, economy, and the environment. In recent years, the development of nanotechnology has led to the design and large-scale production of many new materials and devices with a vast range of applications. However, along with the benefits, the use of nanomaterials raises many questions and generates concerns due to the possible health-risks and environmental impacts. This chapter provides an overview of the Quantitative Structure-Activity Relationships (QSAR) studies performed so far towards predicting nanoparticles' environmental toxicity. Recent progresses on the application of these modeling studies are additionally pointed out. Special emphasis is given to the setup of a QSAR perturbation-based model for the assessment of ecotoxic effects of nanoparticles in diverse conditions. Finally, ongoing challenges that may lead to new and exciting directions for QSAR modeling are discussed.