Molecular similarity in chemical informatics and predictive toxicity modeling: from quantitative read-across (q-RA) to quantitative read-across structure-activity relationship (q-RASAR) with the application of machine learning.

This article aims to provide a comprehensive critical, yet readable, review of general interest to the chemistry community on molecular similarity as applied to chemical informatics and predictive modeling with a special focus on read-across (RA) and read-across structure-activity relationships (RASAR). Molecular similarity-based computational tools, such as quantitative structure-activity relationships (QSARs) and RA, are routinely used to fill the data gaps for a wide range of properties including toxicity endpoints for regulatory purposes. This review will explore the background of RA starting from how structural information has been used through to how other similarity contexts such as physicochemical, absorption, distribution, metabolism, and elimination (ADME) properties, and biological aspects are being characterized. More recent developments of RA's integration with QSAR have resulted in the emergence of novel models such as ToxRead, generalized read-across (GenRA), and quantitative RASAR (q-RASAR). Conventional QSAR techniques have been excluded from this review except where necessary for context.

Breaking the Barriers: Machine-Learning-Based c-RASAR Approach for Accurate Blood-Brain Barrier Permeability Prediction.

The intricate nature of the blood-brain barrier (BBB) poses a significant challenge in predicting drug permeability, which is crucial for assessing central nervous system (CNS) drug efficacy and safety. This research utilizes an innovative approach, the classification read-across structure-activity relationship (c-RASAR) framework, that leverages machine learning (ML) to enhance the accuracy of BBB permeability predictions. The c-RASAR framework seamlessly integrates principles from both read-across and QSAR methodologies, underscoring the need to consider similarity-related aspects during the development of the c-RASAR model. It is crucial to note that the primary goal of this research is not to introduce yet another model for predicting BBB permeability but rather to showcase the refinement in predicting the BBB permeability of organic compounds through the introduction of a c-RASAR approach. This groundbreaking methodology aims to elevate the accuracy of assessing neuropharmacological implications and streamline the process of drug development. In this study, an ML-based c-RASAR linear discriminant analysis (LDA) model was developed using a dataset of 7807 compounds, encompassing both BBB-permeable and -nonpermeable substances sourced from the B3DB database (freely accessible from https://github.com/theochem/B3DB), for predicting BBB permeability in lead discovery for CNS drugs. The model's predictive capability was then validated using three external sets: one containing 276,518 natural products (NPs) from the LOTUS database (accessible from https://lotus.naturalproducts.net/download) for data gap filling, another comprising 13,002 drug-like/drug compounds from the DrugBank database (available from https://go.drugbank.com/), and a third set of 56 FDA-approved drugs to assess the model's reliability. Further diversifying the predictive arsenal, various other ML-based c-RASAR models were also developed for comparison purposes. The proposed c-RASAR framework emerged as a powerful tool for predicting BBB permeability. This research not only advances the understanding of molecular determinants influencing CNS drug permeability but also provides a versatile computational platform for the rapid assessment of diverse compounds, facilitating informed decision-making in drug development and design.

Green Chemistry in the Synthesis of Pharmaceuticals.

The principles of green chemistry (GC) can be comprehensively implemented in green synthesis of pharmaceuticals by choosing no solvents or green solvents (preferably water), alternative reaction media, and consideration of one-pot synthesis, multicomponent reactions (MCRs), continuous processing, and process intensification approaches for atom economy and final waste reduction. The GC's execution in green synthesis can be performed using a holistic design of the active pharmaceutical ingredient's (API) life cycle, minimizing hazards and pollution, and capitalizing the resource efficiency in the synthesis technique. Thus, the presented review accounts for the comprehensive exploration of GC's principles and metrics, an appropriate implication of those ideas in each step of the reaction schemes, from raw material to an intermediate to the final product's synthesis, and the final execution of the synthesis into scalable industry-based production. For real-life examples, we have discussed the synthesis of a series of established generic pharmaceuticals, starting with the raw materials, and the intermediates of the corresponding pharmaceuticals. Researchers and industries have thoughtfully instigated a green synthesis process to control the atom economy and waste reduction to protect the environment. We have extensively discussed significant reactions relevant for green synthesis, one-pot cascade synthesis, MCRs, continuous processing, and process intensification, which may contribute to the future of green and sustainable synthesis of APIs.

Chronic aquatic toxicity assessment of diverse chemicals on Daphnia magna using QSAR and chemical read-across.

We have modeled here chronic Daphnia toxicity taking pNOEC (negative logarithm of no observed effect concentration in mM) and pEC50 (negative logarithm of half-maximal effective concentration in mM) as endpoints using QSAR and chemical read-across approaches. The QSAR models were developed by strictly obeying the OECD guidelines and were found to be reliable, predictive, accurate, and robust. From the selected features in the developed models, we have found that an increase in lipophilicity and saturation, the presence of electrophilic or electronegative or heavy atoms, the presence of sulphur, amine, and their related functionality, an increase in mean atomic polarizability, and higher number of (thio-) carbamates (aromatic) groups are responsible for chronic toxicity. Therefore, this information might be useful for the development of environmentally friendly and safer chemicals and data-gap filling as well as reducing the use of identified toxic chemicals which have chronic toxic effects on aquatic ecosystems. Approved classes of drugs from DrugBank databases and diverse groups of chemicals from the Chemical and Product Categories (CPDat) database were also assessed through the developed models.

Prediction-Inspired Intelligent Training for the Development of Classification Read-across Structure-Activity Relationship (c-RASAR) Models for Organic Skin Sensitizers: Assessment of Classification Error Rate from Novel Similarity Coefficients.

The advancements in the field of cheminformatics have led to a reduction in animal testing to estimate the activity, property, and toxicity of query chemicals. Read-across structure-activity relationship (RASAR) is an emerging concept that utilizes various similarity functions derived from chemical information to develop highly predictive models. Unlike quantitative structure-activity relationship (QSAR) models, RASAR descriptors of a query compound are computed from its close congeners instead of the compound itself, thus targeting predictions in the model training phase. The objective of the present study is not to propose new QSAR models for skin sensitization but to demonstrate the enhancement in the quality of predictions of the skin-sensitizing potential of organic compounds by developing classification-based RASAR (c-RASAR) models. A diverse, previously curated data set was collected from the literature for which 2D descriptors were computed. The extracted essential features were then used to develop a classification-based linear discriminant analysis (LDA) QSAR model. Furthermore, from the read-across-based predictions, RASAR descriptors were calculated using the basic settings of the hyperparameters for the Laplacian Kernel-based optimum similarity measure. After feature selection, an LDA c-RASAR model was developed, which superseded the prediction quality of the LDA-QSAR model. Various other combinations of RASAR descriptors were also taken to develop additional c-RASAR models, all showing better prediction quality than the LDA QSAR model while using a lower number of descriptors. Various other machine learning c-RASAR models were also developed for comparison purposes. In this work, we have proposed and analyzed three new similarity metrics: gm_class, sm1, and sm2. The first one is an indicator variable used to generate a simple univariate c-RASAR model with good prediction ability, while the remaining two are similarity indices used to analyze possible activity cliffs in the training and test sets and are believed to play an important role in the modelability analysis of data sets.

On Some Novel Similarity-Based Functions Used in the ML-Based q-RASAR Approach for Efficient Quantitative Predictions of Selected Toxicity End Points.

The novel quantitative read-across structure-activity relationship (q-RASAR) approach uses read-across-derived similarity functions in the quantitative structure-activity relationship (QSAR) modeling framework in a unique way for supervised model generation. The aim of this study is to explore how this workflow enhances the external (test set) prediction quality of conventional QSAR models by the incorporation of some novel similarity-based functions as additional descriptors using the same level of chemical information. To establish this, five different toxicity data sets, for which QSAR models were reported previously, have been considered in the q-RASAR modeling exercise, which uses chemical similarity-derived measures. The identical sets of chemical features along with the same compositions of training and test sets as reported previously were used in the present analysis for ease of comparison. The RASAR descriptors were calculated based on a chosen similarity measure with the default setting of relevant hyperparameter(s) and were then clubbed with the original structural and physicochemical descriptors, and the number of selected features was further optimized by employing a grid search technique applied on the respective training sets. These features were then used to develop multiple linear regression (MLR) q-RASAR models that show enhanced predictivity as compared to the QSAR models developed previously. Moreover, various other ML algorithms like support vector machine (SVM), linear SVM, random forest, partial least squares, and ridge regression were also employed using the same feature combinations as used in the MLR models to compare the prediction qualities. The q-RASAR models for five different data sets possess at least one of the RASAR descriptors, RA function, gm, and average similarity, suggesting that these are important determinants of similarities that contribute to the development of predictive q-RASAR models, as also evident from the SHAP analysis of the models.

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.

First report of q-RASAR modeling toward an approach of easy interpretability and efficient transferability.

Quantitative structure-activity relationship (QSAR) and read-across techniques have recently been merged into a new emerging field of read-across structure-activity relationship (RASAR) that uses the chemical similarity concepts of read-across (an unsupervised step) and finally develops a supervised learning model (like QSAR). The RASAR method has so far been used only in case of graded predictions or classification modeling. In this work, we attempt, for the first time, to apply RASAR for quantitative predictions (q-RASAR) using a case study of androgen receptor binding affinity data. We have computed a number of error-based and similarity-based measures such as weighted standard deviation of the predicted values, coefficient of variation of the computed predictions, average similarity level of close training compounds for each query molecule, standard deviation and coefficient of variation of similarity levels, maximum similarity levels to positive and negative close training compounds, a concordance measure indicating similarity to positive, negative or both classes of close training compounds, etc. We have clubbed these additional measures along with the selected chemical descriptors from the previously developed QSAR model and redeveloped new partial least squares models from the training set, and predicted the endpoint using the query data set. Interestingly, these new models outperform the internal and external validation quality of the original QSAR model. In this study, we have also introduced a new similarity-based concordance measure (Banerjee-Roy coefficient) that can significantly contribute to the model quality. A q-RASAR model also has the advantage over read-across predictions in providing easy interpretation and indicating quantitative contributions of important chemical features. The strategy described here should be applicable to other biological/toxicological/property data modeling for enhanced quality of predictions, easy interpretability, and efficient transferability.

Prediction reliability of QSAR models: an overview of various validation tools.

The reliability of any quantitative structure-activity relationship (QSAR) model depends on multiple aspects such as the accuracy of the input dataset, selection of significant descriptors, the appropriate splitting process of the dataset, statistical tools used, and most notably on the measures of validation. Validation, the most crucial step in QSAR model development, confirms the reliability of the developed QSAR models and the acceptability of each step in the model development. The present review deals with various validation tools that involve multiple techniques that improve the model quality and robustness. The double cross-validation tool helps in building improved quality models using different combinations of the same training set in an inner cross-validation loop. This exhaustive method is also integrated for small datasets (< 40 compounds) in another tool, namely the small dataset modeler tool. The main aim of QSAR researchers is to improve prediction quality by lowering the prediction errors for the query compounds. 'Intelligent' selection of multiple models and consensus predictions integrated in the intelligent consensus predictor tool were found to be more externally predictive than individual models. Furthermore, another tool called Prediction Reliability Indicator was explained to understand the quality of predictions for a true external set. This tool uses a composite scoring technique to identify query compounds as 'good' or 'moderate' or 'bad' predictions. We have also discussed a quantitative read-across tool which predicts a chemical response based on the similarity with structural analogues. The discussed tools are freely available from https://dtclab.webs.com/software-tools or http://teqip.jdvu.ac.in/QSAR_Tools/DTCLab/ and https://sites.google.com/jadavpuruniversity.in/dtc-lab-software/home (for read-across).

Combination chemotherapy versus temozolomide for patients with methylated MGMT (m-MGMT) glioblastoma: results of computational biological modeling to predict the magnitude of treatment benefit.

BACKGROUND: A randomized trial in glioblastoma patients with methylated-MGMT (m-MGMT) found an improvement in median survival of 16.7 months for combination therapy with temozolomide (TMZ) and lomustine, however the approach remains controversial and relatively under-utilized. Therefore, we sought to determine whether comprehensive genomic analysis can predict which patients would derive large, intermediate, or negligible benefits from the combination compared to single agent chemotherapy. METHODS: Comprehensive genomic information from 274 newly diagnosed patients with methylated-MGMT glioblastoma (GBM) was downloaded from TCGA. Mutation and copy number changes were input into a computational biologic model to create an avatar of disease behavior and the malignant phenotypes representing hallmark behavior of cancers. In silico responses to TMZ, lomustine, and combination treatment were biosimulated. Efficacy scores representing the effect of treatment for each treatment strategy were generated and compared to each other to ascertain the differential benefit in drug response. RESULTS: Differential benefits for each drug were identified, including strong, modest-intermediate, negligible, and deleterious (harmful) effects for subgroups of patients. Similarly, the benefits of combination therapy ranged from synergy, little or negligible benefit, and deleterious effects compared to single agent approaches. CONCLUSIONS: The benefit of combination chemotherapy is predicted to vary widely in the population. Biosimulation appears to be a useful tool to address the disease heterogeneity, drug response, and the relevance of particular clinical trials observations to individual patients. Biosimulation has potential to spare some patients the experience of over-treatment while identifying patients uniquely situated to benefit from combination treatment. Validation of this new artificial intelligence tool is needed.

Therapeutics for COVID-19: from computation to practices-where we are, where we are heading to.

After the 1918 Spanish Flu pandemic caused by the H1N1 virus, the recent coronavirus disease 2019 (COVID-19) brought us to the time of serious global health catastrophe. Although no proven therapies are identified yet which can offer a definitive treatment of the COVID-19, a series of antiviral, antibacterial, antiparasitic, immunosuppressant drugs have shown clinical benefits based on repurposing theory. However, these studies are made on small number of patients, and, in majority of the cases, have been carried out as nonrandomized trials. As society is running against the time to combat the COVID-19, we present here a comprehensive review dealing with up-to-date information of therapeutics or drug regimens being utilized by physicians to treat COVID-19 patients along with in-depth discussion of mechanism of action of these drugs and their targets. Ongoing vaccine trials, monoclonal antibodies therapy and convalescent plasma treatment are also discussed. Keeping in mind that computational approaches can offer a significant insight to repurposing based drug discovery, an exhaustive discussion of computational modeling studies is performed which can assist target-specific drug discovery.

QSPR modeling of octanol-water partition coefficient and organic carbon normalized sorption coefficient of diverse organic chemicals using Extended Topochemical Atom (ETA) indices.

Octanol-water partition coefficient (logKow) and soil organic carbon content normalized sorption coefficient (logKoc) values are two important physicochemical properties in the context of bioaccumulation and environmental fate of organic compounds and their environmental risk assessment. Simple, interpretable and easy-to-derive extended topochemical atom (ETA) indices obtained from 2D structural representation of compounds were used for quantitative structure-property relationship (QSPR) modeling of these two endpoints. Linear regression based models developed using only ETA indices show encouraging statistical and validation results. Based on the information obtained from developed QSPR models, we may conclude that molecular volume, branching pattern, presence of hydrophobic Cl atoms, cyclicity/fusion, polar environment, electron density, unsaturation content, hydrogen bonding propensity or hydrogen bond donor atoms, local topology, presence of heteroatoms and aromaticity are crucial factors in controlling the logKow and logKoc values of the compounds. The suggested explanatory features for different classes of chemicals or the whole diverse set can help in safer designing of chemicals, which is one of the primary agenda of the "Green Chemistry" program.

Exploring QSAR modeling of toxicity of chemicals on earthworm.

Earthworm provides sustainability towards the agroecosystem which can be degraded day by day by the extensive use of pesticides (e.g., fungicides, insecticides and herbicides). The present study attempts to develop a predictive quantitative structure-activity relationship (QSAR) model for toxicity of pesticides to earthworm in order to give a suitable guidance for designing new analogues with lower toxicity by exploring the important chemical features which are required to develop safer alternatives. The QSAR model was developed by using the negative logarithm of lethal concentration (pLC50) values of pesticides towards earthworm. We have used various 2D descriptors along with extended topochemical atom (ETA) indices as independent variables for the development of the model. The developed partial least squares (PLS) model was subjected to statistical validation tests proving that the model is statistically reliable and robust (R2 = 0.765, Q2 = 0.614, Q2F1 = 0.734, Q2F2 = 0.713). The contributing descriptors in the model suggested that the pesticides may affect the earthworm nucleic acid through various physicochemical interactions including hydrophobicity, hydrogen bonding, electron donor acceptor complex formation, π-π stacking interaction and charge transfer complex formation.

New Workflow for QSAR Model Development from Small Data Sets: Small Dataset Curator and Small Dataset Modeler. Integration of Data Curation, Exhaustive Double Cross-Validation, and a Set of Optimal Model Selection Techniques.

Quantitative structure-activity relationship (QSAR) modeling is a well-known in silico technique with extensive applications in several major fields such as drug design, predictive toxicology, materials science, food science, etc. Handling small-sized datasets due to the lack of experimental data for specialized end points is a crucial task for the QSAR researcher. In the present study, we propose an integrated workflow/scheme capable of dealing with small dataset modeling that integrates dataset curation, "exhaustive" double cross-validation and a set of optimal model selection techniques including consensus predictions. We have developed two software tools, namely, Small Dataset Curator, version 1.0.0, and Small Dataset Modeler, version 1.0.0, to effortlessly execute the proposed workflow. These tools are freely available for download from https://dtclab.webs.com/software-tools . We have performed case studies employing seven diverse datasets to demonstrate the performance of the proposed scheme (including data curation) for small dataset QSAR modeling. The case studies also confirm the usability and stability of the developed software tools.

Consensus QSAR modeling of toxicity of pharmaceuticals to different aquatic organisms: Ranking and prioritization of the DrugBank database compounds.

In the present work, quantitative structure-activity relationship (QSAR) models have been developed for ecotoxicity of pharmaceuticals on four different aquatic species namely Pseudokirchneriella subcapitata, Daphnia magna, Oncorhynchus mykiss and Pimephales promelas using genetic algorithm (GA) for feature selection followed by Partial Least Squares regression technique according to the Organization for Economic Co-operation and Development (OECD) guidelines. Double cross-validation methodology was employed for selecting suitable models. Only 2D descriptors were used for capturing chemical information and model building, whereas validation of the models was performed by considering various stringent internal and external validation metrics. Interestingly, models could be developed even without using any LogP terms in contrary to the usual dependence of toxicity on lipophilicity. However, the current manuscript proposes highly robust and more predictive models employing computed logP descriptors. The applicability domain study was performed in order to set a predefined chemical zone of applicability for the obtained QSAR models, and the test compounds falling outside the domain were not taken for further analysis while making a prioritized list. An additional comparison was made with ECOSAR, an online expert system for toxicity prediction of organic pollutants, in order to prove predictability of the obtained models. The obtained robust consensus models were utilized to predict the toxicity of a large dataset of approximately 9300 drug-like molecules in order to prioritize the existing drug-like substances in accordance to their acute predicted aquatic toxicities following a scaling technique. Finally, prioritized lists of 500 most toxic chemicals obtained by respective consensus models and those predicted from ECOSAR tool have been reported.

Chemometric modeling of aquatic toxicity of contaminants of emerging concern (CECs) in Dugesia japonica and its interspecies correlation with daphnia and fish: QSTR and QSTTR approaches.

The contaminants of emerging concern (CEC) are universally detected in surface water and soil. They can affect the wild life, and their subsequent translocation through the food chain can affect human health, which is an issue of serious concern. Very few amounts of ecotoxicological data are available on the environmental behavior and ecotoxicity of CEC, thus modeling approaches are essential to bridge the existing gap in experimental data. In this present study, we have developed quantitative structure-toxicity relationship (QSTR) models using a data set of 75 compounds for the prediction of aquatic ecotoxicity of CECs on fresh water planarian (Dugesia japonica) by partial least squares (PLS) regression algorithm using simple molecular descriptors selected by genetic algorithm approach. We also explore the correlations between toxicity against D. japonica and those against daphnia (D. magna) and fish (P. promelas), and these were improved on addition of a few molecular descriptors (B08[C-O] and B09[N-O] in case of daphnia and C-006 and H-052 in case of fish) which allowed us to develop predictive interspecies quantitative structure toxicity-toxicity relationship (QSTTR) models, allowing to extrapolate data from one endpoint to another endpoint. The QSTR (Q2LOO ranging from 0.630 to 0.720 and R2pred ranging from 0.723 to 0.798) and QSTTR (Q2LOO = 0.60 and 0.67, R2pred = 0.88 and 0.84) models have desirable statistical qualities and acceptable internal and external validation measures, meeting rigorous criteria of different validation metrics and showing acceptability for regulatory purposes as proposed by Organization for Economic Cooperation and Development (OECD). Consensus predictions were also performed based on multiple models generated in this study by using the "Intelligent Consensus Predictor" (ICP) tool to enhance the prediction quality for external set compounds.

Extrapolating between toxicity endpoints of metal oxide nanoparticles: Predicting toxicity to Escherichia coli and human keratinocyte cell line (HaCaT) with Nano-QTTR.

Synthesis of novel nanoparticles should always be accompanied by a comprehensive assessment of risk to human health and to ecosystem. Application of in silico models is encouraged by regulatory authorities to fill the data gaps related to the properties of nanoparticles affecting the environment and human health. Interspecies toxicity correlations provide a tool for estimation of contaminant's sensitivity with known levels of uncertainty for a diverse pool of species. We propose here first interspecies cytotoxicity correlation models between Escherichia coli (prokaryotic system) and human keratinocyte cell line (HaCaT) (eukaryotic system) to assess the discriminatory features for cytotoxicity of metal oxide nanoparticles. The nano-QTTR models can be employed for extrapolating cytotoxicity to E. coli and human keratinocyte cell line (HaCaT) for metal nanoparticles when the data for the other species are available. Informative illustrations of the contributing mechanisms of toxic action of the metal oxide nanoparticles to the HaCaT cell line as well as to the E. coli are identified from the developed nano quantitative toxicity-toxicity relationship (nano-QTTR) models.

In vitro evaluation and in silico screening of synthetic acetylcholinesterase inhibitors bearing functionalized piperidine pharmacophores.

A series of densely functionalized piperidine (FP) scaffolds was synthesized following a diastereoselective one-pot multicomponent protocol under eco-friendly conditions. The FPs were evaluated in vitro for their acetylcholinesterase (AChE) inhibitory activity, and in silico studies for all the target compounds were carried out using pharmacophore mapping, molecular docking and quantitative structure-activity relationship (QSAR) analysis in order to understand the structural features required for interaction with the AChE enzyme and the key active site residues involved in the intermolecular interactions. Halogenation, nitration or 3,4-methylenedixoxy-substitution at the phenyl ring attached to the 2- and 6-positions of 1,2,5,6-tetrahydropyridine nucleus in compounds 14-17, 19, 20, 24 and 26 greatly enhanced the AChE inhibitory activity. The docking analysis demonstrated that the inhibitors are well-fitted in the active sites. The in silico studies enlighten the future course of studies in modifying the scaffolds for better therapeutic efficacy against the deadly Alzheimer's disease.