Real external predictivity of QSAR models: how to evaluate it? Comparison of different validation criteria and proposal of using the concordance correlation coefficient

Synthesis, biological evaluation, theoretical calculations, QSAR and molecular docking studies of novel arylaminonaphthols as potent antioxidants and BChE inhibitors

A completely green protocol was developed for the synthesis of a series of arylaminonaphthol derivatives in the presence of N-ethylethanolamine (NEEA) as a catalyst under ultrasonic irradiation and solventless conditions. The major assets of this methodology were the use of non-toxic organic medium, available catalyst, mild reaction condition, and good to excellent yield of desired products. All of the synthesized products were screened for their in vitro antioxidant activity using DPPH, ABTS, and Ferric-phenanthroline assays and it was found that most of them are potent antioxidant agents. Also, their butyrylcholinesterase inhibitory activity has been investigated in vitro. All tested compounds exhibited potential inhibitory activity toward BuChE when compared to standard reference drug galantamine, however, compounds 4r, 4u, 4 g and 4x gave higher butyrylcholinesterase inhibitory with IC50 values of 14.78 ± 0.65 µM, 16.18 ± 0.50 µM, 20.00 ± 0.50 µM, and 20.28 ± 0.08 µM respectively. On the other hand, we employed density functional theory (DFT), calculations to analyze molecular geometry and global reactivity descriptors, and MESP analysis to predict electrophilic and nucleophilic attacks. A quantitative structure-activity relationship (QSAR) investigation was conducted on the antioxidant and butyrylcholinesterase properties of 25 arylaminonaphthol derivatives, resulting in robust and satisfactory models. To evaluate their anti-Alzheimer's activity, compounds 4 g, 4q, 4r, 4u, and 4x underwent docking simulations at the active site of the acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), revealing why these compounds displayed superior activity, consistent with the biological findings.

Fluorinated benzoxazinones designed via MIA-QSAR, docking and molecular dynamics as protoporphyrinogen IX oxidase inhibitors

BACKGROUND

Fluorine plays a significant role in agrochemical science because approximately 25% of herbicides licensed worldwide contain this element. In a pool of previously synthesized benzoxazinones, some compounds contained fluorine and demonstrated inhibitory activities against protoporphyrinogen IX oxidase (PPO). Therefore, three data sets of benzoxazinone derivatives with known inhibitory activity against PPO were employed to build a multivariate image analysis applied to a quantitative structure-activity relationships (MIA-QSAR) model to identify improved analogs with at least one fluorine substituent.

RESULTS

The QSAR model was vigorously validated and demonstrated to be highly predictive (r2 = 0.85, q2 = 0.71, and r2 pred = 0.88); thus, the model can provide reliable estimations for the PPO inhibitory activity of unknown derivatives. From these compounds, a couple of N-substituted benzoxazinones that contained the -CH2CHF2 group were found with predicted pKi values larger than 8 (Ki in mol L-1) and higher lipophilicity than the most active data set compounds. In addition, we carried out a systematic investigation of the binding mode of PPO by performing computational docking followed by molecular dynamics simulations. The proposed binding mode was consistent with experimental studies, and several potential key residues were identified.

CONCLUSION

Two new proposed benzoxazinones exhibited better performance than compounds of the data set, and fluorine substituents played pivotal roles in describing the biological activities. © 2024 Society of Chemical Industry.

The rat acute oral toxicity of trifluoromethyl compounds (TFMs): a computational toxicology study combining the 2D-QSTR, read-across and consensus modeling methods

All areas of the modern society are affected by fluorine chemistry. In particular, fluorine plays an important role in medical, pharmaceutical and agrochemical sciences. Amongst various fluoro-organic compounds, trifluoromethyl (CF3) group is valuable in applications such as pharmaceuticals, agrochemicals and industrial chemicals. In the present study, following the strict OECD modelling principles, a quantitative structure-toxicity relationship (QSTR) modelling for the rat acute oral toxicity of trifluoromethyl compounds (TFMs) was established by genetic algorithm-multiple linear regression (GA-MLR) approach. All developed models were evaluated by various state-of-the-art validation metrics and the OECD principles. The best QSTR model included nine easily interpretable 2D molecular descriptors with clear physical and chemical significance. The mechanistic interpretation showed that the atom-type electro-topological state indices, molecular connectivity, ionization potential, lipophilicity and some autocorrelation coefficients are the main factors contributing to the acute oral toxicity of TFMs against rats. To validate that the selected 2D descriptors can effectively characterize the toxicity, we performed the chemical read-across analysis. We also compared the best QSTR model with public OPERA tool to demonstrate the reliability of the predictions. To further improve the prediction range of the QSTR model, we performed the consensus modelling. Finally, the optimum QSTR model was utilized to predict a true external set containing many untested/unknown TFMs for the first time. Overall, the developed model contributes to a more comprehensive safety assessment approach for novel CF3-containing pharmaceuticals or chemicals, reducing unnecessary chemical synthesis whilst saving the development cost of new drugs.

Design, Synthesis, Antitumor, and Antiplasmodial Evaluation of New 7-Chloroquinoline-Benzimidazole Hybrids

Newly synthesized 7-chloro-4-aminoquinoline-benzimidazole hybrids were characterized by NMR and elemental analysis. Compounds were tested for their effects on the growth of the non-tumor cell line MRC-5 (human fetal lung fibroblasts) and carcinoma (HeLa and CaCo-2), leukemia, and lymphoma (Hut78, THP-1, and HL-60) cell lines. The obtained results, expressed as the concentration at which 50% inhibition of cell growth is achieved (IC50 value), show that the tested compounds affect cell growth differently depending on the cell line and the applied dose (IC50 ranged from 0.2 to >100 µM). Also, the antiplasmodial activity of these hybrids was evaluated against two P. falciparum strains (Pf3D7 and PfDd2). The tested compounds showed potent antiplasmodial activity, against both strains, at nanomolar concentrations. Quantitative structure-activity relationship (QSAR) analysis resulted in predictive models for antiplasmodial activity against the 3D7 strain (R2 = 0.886; Rext2 = 0.937; F = 41.589) and Dd2 strain (R2 = 0.859; Rext2 = 0.878; F = 32.525) of P. falciparum. QSAR models identified the structural features of these favorable effects on antiplasmodial activities.

From molecular descriptors to the developmental toxicity prediction of pesticides/veterinary drugs/bio-pesticides against zebrafish embryo: Dual computational toxicological approaches for prioritization

The escalating introduction of pesticides/veterinary drugs into the environment has necessitated a rapid evaluation of their potential risks to ecosystems and human health. The developmental toxicity of pesticides/veterinary drugs was less explored, and much less the large-scale predictions for untested pesticides, veterinary drugs and bio-pesticides. Alternative methods like quantitative structure-activity relationship (QSAR) are promising because their potential to ensure the sustainable and safe use of these chemicals. We collected 133 pesticides and veterinary drugs with half-maximal active concentration (AC50) as the zebrafish embryo developmental toxicity endpoint. The QSAR model development adhered to rigorous OECD principles, ensuring that the model possessed good internal robustness (R2 > 0.6 and QLOO2 > 0.6) and external predictivity (Rtest2 > 0.7, QFn2 >0.7, and CCCtest > 0.85). To further enhance the predictive performance of the model, a quantitative read-across structure-activity relationship (q-RASAR) model was established using the combined set of RASAR and 2D descriptors. Mechanistic interpretation revealed that dipole moment, the presence of C-O fragment at 10 topological distance, molecular size, lipophilicity, and Euclidean distance (ED)-based RA function were main factors influencing toxicity. For the first time, the established QSAR and q-RASAR models were combined to prioritize the developmental toxicity of a vast array of true external compounds (pesticides/veterinary drugs/bio-pesticides) lacking experimental values. The prediction reliability of each query molecule was evaluated by leverage approach and prediction reliability indicator. Overall, the dual computational toxicology models can inform decision-making and guide the design of new pesticides/veterinary drugs with improved safety profiles.

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.

QSPR for the prediction of critical micelle concentration of different classes of surfactants using machine learning algorithms

The determination of the critical micelle concentration (CMC) is a crucial factor when evaluating surfactants, making it an essential tool in studying the properties of surfactants in various industrial fields. In this present research, we assembled a comprehensive set of 593 different classes of surfactants including, anionic, cationic, nonionic, zwitterionic, and Gemini surfactants to establish a link between their molecular structure and the negative logarithmic value of critical micelle concentration (pCMC) utilizing quantitative structure-property relationship (QSPR) methodologies. Statistical analysis revealed that a set of 14 significant Mordred descriptors (SlogP, GATS6d, nAcid, GATS8dv, GATS4dv, PEOE_VSA11, GATS8d, ATS0p, GATS1d, MATS5p, GATS3d, NdssC, GATS6dv and EState_VSA4), along with temperature, served as appropriate inputs. Different machine learning methods, such as multiple linear regression (MLR), random forest regression (RFR), artificial neural network (ANN), and support vector regression (SVM), were employed in this study to build QSPR models. According to the statistical coefficients of QSPR models, SVR with Dragonfly hyperparameter optimization (SVR-DA) was the most accurate in predicting pCMC values, achieving (R2 = 0.9740, Q2 = 0.9739, r‾m2 = 0.9627, and Δrm2 = 0.0244) for the entire dataset.

Toxicity of individual and mixture of organic compounds to P. Phosphoreum and S. Capricornutum using interpretable simple structural parameters

Human and environmental ecosystem beings are exposed to multicomponent compound mixtures but the toxicity nature of compound mixtures is not alike to the individual chemicals. This work introduces four models for the prediction of the negative logarithm of median effective concentration (pEC50) of individual chemicals to marine bacteria Photobacterium Phosphoreum (P. Phosphoreum) and algal test species Selenastrum Capricornutum (S. Capricornutum) as well as their mixtures to P. Phosphoreum, and S. Capricornutum. These models provide the simplest approaches for the forecast of pEC50 of some classes of organic compounds from their interpretable structural parameters. Due to the lack of adequate toxicity data for chemical mixtures, the largest available experimental data of individual chemicals (55 data) and their mixtures (99 data) are used to derive the new correlations. The models of individual chemicals are based on two simple structural parameters but chemical mixture models require further interaction terms. The new model's results are compared with the outputs of the best accessible quantitative structure-activity relationships (QSARs) models. Various statistical parameters are done on the new and comparative complex QSAR models, which confirm the higher reliability and simplicity of the new correlations.

Prediction of acute toxicity for Chlorella vulgaris caused by tire wear particle-derived compounds using quantitative structure-activity relationship models

Tire wear particles (TWPs) enter aquatic ecosystems through various pathways, such as rainwater and urban runoff. Additives in TWPs can harm aquatic organisms in these ecosystems. Therefore, it is essential to investigate their toxicity to aquatic organisms. In our study, we initially recorded the median effective concentrations of 21 TWP-derived compounds on Chlorella vulgaris growth, ranging from 0.04 to 8.60 mg/L. Subsequently, through an extensive review of the literature, we incorporated 112 compounds with specific toxicity endpoints to construct the QSAR model using genetic algorithm and multiple linear regression techniques, followed by the construction of the consensus model and the quantitative read-across structure-activity relationship (q-RASAR) model. Meanwhile, we employed rigorous internal and external validation measures to assess the performance of the model. The results indicated that the developed q-RASAR model exhibited strong adaptation, robustness, and reliable prediction, with q-RASAR indicators of Q2LOO = 0.7673, R2tr = 0.8079, R2test = 0.8610, Q2Fn = 0.8285-0.8614, and CCCtest = 0.9222. Based on an external dataset containing 128 emerging TWP-derived compounds, the model's applicability domain coverage was 90.6 %. The q-RASAR model predicted that the structure of diphenylamine was associated with higher toxicity, possibly liked to the SpMax2_Bhm and LogBCF descriptors. The established model reliably provides prediction and fills a critical data gap. These findings highlight the potential risks posed by emerging TWP-derived compounds to aquatic organisms.

Prediction of mass spectrometry ionization efficiency based on COSMO-RS and machine learning algorithms

Non-targeted analysis of high-resolution mass spectrometry (MS) can identify thousands of compounds, which also gives a huge challenge to their quantification. The aim of this study is to investigate the impact of mass spectrometry ionization efficiency on various compounds in food at different solvent ratios and to develop a predictive model for mass spectrometry ionization efficiency to enable non-targeted quantitative prediction of unknown compounds. This study covered 70 compounds in 14 different mobile phase ratio environments in positive ion mode to analyze the rules of the matrix effect. With the organic phase ratio from low to high, most compounds changed by 1.0 log units in log IE. The addition of formic acid enhanced the signal but also promoted the matrix effect, which often occurred in compounds with strong ionization capacity. It was speculated that the matrix effect was mainly in the form of competitive charge and charged droplet' gasification sites during MS detection. Subsequently, we present a log IE prediction method built using the COSMO-RS software and the artificial neural network (ANN) algorithm to address this difficulty and overcome the shortcomings of previous models, which always ignore the matrix effect. This model was developed following the principles of QSAR modeling recommended by the Organization for Economic Cooperation and Development (OECD). Furthermore, we validated this approach by predicting the log IE of 70 compounds, including those not involved in the log IE model development. The results presented demonstrate that the method we put forward has an excellent prediction accuracy for log IE (R2pred = 0.880), which means that it has the potential to predict the log IE of new compounds without authentic standards.

Toxicity assessment of dioxins and their transformation by-products from inferred degradation pathways

Due to the significant POPs characteristics, dioxins caused concern in public health and environmental protection. Evaluating the toxicity risk of dioxin degradation pathways is critical. OCDD, 1,2,3,4,6,7,8-HpCDD, and 1,2,3,4,6,7,8-HpCDF, which are highly abundant in the environment and have strong biodegradation capabilities, were selected as precursor molecules in this study. Firstly, their transformation pathways were deduced during the metabolism of biometabolism, microbial aerobic, microbial anaerobic, and photodegradation pathways, and density function theory (DFT) was used to calculate the Gibbs free energy to infer the possibility of the occurrence of the transformation pathway. Secondly, the carcinogenic potential of the precursors and their degradation products was evaluated using the TOPKAT modeling method. With the help of the positive indicator (0-1) normalization method and heat map analysis, a significant increase in the toxic effect of some of the transformation products was found, and it was inferred that it was related to the structure of the transformation products. Meanwhile, the strength of the endocrine disrupting effect of dioxin transformation products was quantitatively assessed using molecular docking and subjective assignment methods, and it was found that dioxin transformation products with a higher content of chlorine atoms and molecules similar to those of thyroid hormones exhibited a higher risk of endocrine disruption. Finally, the environmental health risks caused by each degradation pathway were comprehensively assessed with the help of the negative indicator (1-2) standardization method, which provides a theoretical basis for avoiding the toxicity risks caused by dioxin degradation transformation. In addition, the 3D-QSAR model was used to verify the necessity and rationality of this study. This paper provides theoretical support and reference significance for the toxicity assessment of dioxin degradation by-products from inferred degradation pathways.

Application of parallel artificial membrane permeability assay technique and chemometric modeling for blood-brain barrier permeability prediction of protein kinase inhibitors

Aim: This study aims to investigate the passive diffusion of protein kinase inhibitors through the blood-brain barrier (BBB) and to develop a model for their permeability prediction. Materials & methods: We used the parallel artificial membrane permeability assay to obtain logPe values of each of 34 compounds and calculated descriptors for these structures to perform quantitative structure-property relationship modeling, creating different regression models. Results: The logPe values have been calculated for all 34 compounds. Support vector machine regression was considered the most reliable, and CATS2D_09_DA, CATS2D_04_AA, B04[N-S] and F07[C-N] descriptors were identified as the most influential to passive BBB permeability. Conclusion: The quantitative structure-property relationship-support vector machine regression model that has been generated can serve as an efficient method for preliminary screening of BBB permeability of new analogs.

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

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

In silico insights into design of novel VEGFR-2 inhibitors: SMILES-based QSAR modelling, and docking studies on substituted benzo-fused heteronuclear derivatives

Eight QSAR models (M1-M8) were developed from a dataset of 118 benzo-fused heteronuclear derivatives targeting VEGFR-2 by Monte Carlo optimization method of CORALSEA 2023 software. Models were generated with hybrid optimal descriptors using both SMILES and Graphs with zero- and first-order Morgan extended connectivity index from a training set of 103 derivatives. All statistical parameters for model validation were within the prescribed limits, establishing the models to be robust and of excellent quality. Among all models, split-2 of M5 was the best-fit as reflected by   r v a lidation 2 ,   Q v a lidation 2 and MAE . Mechanistic interpretation of this model assisted the identification of structural descriptors as promoters and hinderers for VEGFR-2 inhibition. These descriptors were utilized to design novel VEGFR-2 inhibitors (YS01-YS07) by bringing modifications in compound MS90 in the dataset. Docking of all designed compounds, MS90 and sorafenib with VEGFR-2 binding site revealed favourable binding interactions. Docking score of YS07 was higher than that of MS90 and sorafenib. Molecular dynamics simulation study revealed sustained interactions of YS07 with key amino acids of VEGFR-2 at a run time of 100 ns. This study concludes the development of a best fit QSAR model which can assist the design of new anticancer agents targeting VEGFR-2.

Prediction of critical micelle concentration for per- and polyfluoroalkyl substances

In this study, we focus on the development of Quantitative Structure-Property Relationship (QSPR) models to predict the critical micelle concentration (CMC) for per- and polyfluoroalkyl substances (PFASs). Experimental CMC values for both fluorinated and non-fluorinated compounds were meticulously compiled from existing literature sources. Our approach involved constructing two distinct types of models based on Support Vector Machine (SVM) algorithms applied to the dataset. Type (I) models were trained exclusively on CMC values for fluorinated compounds, while Type (II) models were developed utilizing the entire dataset, incorporating both fluorinated and non-fluorinated compounds. Comparative analyses were conducted against reference data, as well as between the two model types. Encouragingly, both types of models exhibited robust predictive capabilities and demonstrated high reliability. Subsequently, the model having the broadest applicability domain was selected to complement the existing experimental data, thereby enhancing our understanding of PFAS behaviour.

First report on chemometrics-driven multilayered lead prioritization in addressing oxysterol-mediated overexpression of G protein-coupled receptor 183

Contemporary research has convincingly demonstrated that upregulation of G protein-coupled receptor 183 (GPR183), orchestrated by its endogenous agonist, 7α,25-dihydroxyxcholesterol (7α,25-OHC), leads to the development of cancer, diabetes, multiple sclerosis, infectious, and inflammatory diseases. A recent study unveiled the cryo-EM structure of 7α,25-OHC bound GPR183 complex, presenting an untapped opportunity for computational exploration of potential GPR183 inhibitors, which served as our inspiration for the current work. A predictive and validated two-dimensional QSAR model using genetic algorithm (GA) and multiple linear regression (MLR) on experimental GPR183 inhibition data was developed. QSAR study highlighted that structural features like dissimilar electronegative atoms, quaternary carbon atoms, and CH2RX fragment (X: heteroatoms) influence positively, while the existence of oxygen atoms with a topological separation of 3, negatively affects GPR183 inhibitory activity. Post assessment of true external set prediction capability, the MLR model was deployed to screen 12,449 DrugBank compounds, followed by a screening pipeline involving molecular docking, druglikeness, ADMET, protein-ligand stability assessment using deep learning algorithm, molecular dynamics, and molecular mechanics. The current findings strongly evidenced DB05790 as a potential lead for prospective interference of oxysterol-mediated GPR183 overexpression, warranting further in vitro and in vivo validation.

Prioritization of the ecotoxicological hazard of PAHs towards aquatic species spanning three trophic levels using 2D-QSTR, read-across and machine learning-driven modelling approaches

Polycyclic aromatic hydrocarbons (PAHs) represent a common group of environmental pollutants that endanger various aquatic organisms via various pathways. To better prioritize the ecotoxicological hazard of PAHs to aquatic environment, we used 2D descriptors-based quantitative structure-toxicity relationship (QSTR) to assess the toxicity of PAHs toward six aquatic model organisms spanning three trophic levels. According to strict OECD guideline, six easily interpretable, transferable and reproducible 2D-QSTR models were constructed with high robustness and reliability. A mechanistic interpretation unveiled the key structural factors primarily responsible for controlling the aquatic ecotoxicity of PAHs. Furthermore, quantitative read-across and different machine learning approaches were employed to validate and optimize the modelling approach. Importantly, the optimum QSTR models were further applied for predicting the ecotoxicity of hundreds of untested/unknown PAHs gathered from Pesticide Properties Database (PPDB). Especially, we provided a priority list in terms of the toxicity of unknown PAHs to six aquatic species, along with the corresponding mechanistic interpretation. In summary, the models can serve as valuable tools for aquatic risk assessment and prioritization of untested or completely new PAHs chemicals, providing essential guidance for formulating regulatory policies.

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.

Development of a hierarchical support vector regression-based in silico model for the prediction of the cysteine depletion in DPRA

Topical and transdermal treatments have been dramatically growing recently and it is crucial to consider skin sensitization during the drug discovery and development process for these administration routes. Various tests, including animal and non-animal approaches, have been devised to assess the potential for skin sensitization. Furthermore, numerous in silico models have been created, providing swift and cost-effective alternatives to traditional methods such as in vivo, in vitro, and in chemico methods for categorizing compounds. In this study, a quantitative structure-activity relationship (QSAR) model was developed using the innovative hierarchical support vector regression (HSVR) scheme. The aim was to quantitatively predict the potential for skin sensitization by analyzing the percent of cysteine depletion in Direct Peptide Reactivity Assay (DPRA). The results demonstrated accurate, consistent, and robust predictions in the training set, test set, and outlier set. Consequently, this model can be employed to estimate skin sensitization potential of novel or virtual compounds.

QSAR modeling approaches to identify a novel ACE2 inhibitor that selectively bind with the C and N terminals of the ectodomain

Due to the high rates of drug development failure and the massive expenses associated with drug discovery, repurposing existing drugs has become more popular. As a result, we have used QSAR modelling on a large and varied dataset of 657 compounds in an effort to discover both explicit and subtle structural features requisite for ACE2 inhibitory activity, with the goal of identifying novel hit molecules. The QSAR modelling yielded a statistically robust QSAR model with high predictivity (R2tr=0.84, R2ex=0.79), previously undisclosed features, and novel mechanistic interpretations. The developed QSAR model predicted the ACE2 inhibitory activity (PIC50) of 1615 ZINC FDA compounds. This led to the detection of a PIC50 of 8.604 M for the hit molecule (ZINC000027990463). The hit molecule's docking score is -9.67 kcal/mol (RMSD 1.4). The hit molecule revealed 25 interactions with the residue ASP40, which defines the N and C termini of the ectodomain of ACE2. The HIT molecule conducted more than thirty contacts with water molecules and exhibited polar interaction with the ARG522 residue coupled with the second chloride ion, which is 10.4 nm away from the zinc ion. Both molecular docking and QSAR produced comparable findings. Moreover, MD simulation and MMGBSA studies verified docking analysis. The MD simulation showed that the hit molecule-ACE2 receptor complex is stable for 400 ns, suggesting that repurposed hit molecule 3 is a viable ACE2 inhibitor.

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.

Design of new molecules against cervical cancer using DFT, theoretical spectroscopy, 2D/3D-QSAR, molecular docking, pharmacophore and ADMET investigations

Cervical cancer is a major health problem of women. Hormone therapy, via aromatase inhibition, has been proposed as a promising way of blocking estrogen production as well as treating the progression of estrogen-dependent cancer. To overcome the challenging complexities of costly drug design, in-silico strategy, integrating Structure-Based Drug Design (SBDD) and Ligand-Based Drug Design (LBDD), was applied to large representative databases of 39 quinazoline and thioquinazolinone compound derivatives. Quantum chemical and physicochemical descriptors have been investigated using density functional theory (DFT) and MM2 force fields, respectively, to develop 2D-QSAR models, while CoMSIA and CoMFA descriptors were used to build 3D-QSAR models. The robustness and predictive power of the reliable models were verified, via several validation methods, leading to the design of 6 new drug-candidates. Afterwards, 2 ligands were carefully selected using virtual screening methods, taking into account the applicability domain, synthetic accessibility, and Lipinski's criteria. Molecular docking and pharmacophore modelling studies were performed to examine potential interactions with aromatase (PDB ID: 3EQM). Finally, the ADMET properties were investigated in order to select potential drug-candidates against cervical cancer for experimental in vitro and in vivo testing.

The pursuit of accurate predictive models of the bioactivity of small molecules

Property prediction is a key interest in chemistry. For several decades there has been a continued and incremental development of mathematical models to predict properties. As more data is generated and accumulated, there seems to be more areas of opportunity to develop models with increased accuracy. The same is true if one considers the large developments in machine and deep learning models. However, along with the same areas of opportunity and development, issues and challenges remain and, with more data, new challenges emerge such as the quality and quantity and reliability of the data, and model reproducibility. Herein, we discuss the status of the accuracy of predictive models and present the authors' perspective of the direction of the field, emphasizing on good practices. We focus on predictive models of bioactive properties of small molecules relevant for drug discovery, agrochemical, food chemistry, natural product research, and related fields.

In silico soil degradation and ecotoxicity analysis of veterinary pharmaceuticals on terrestrial species: first report

With the aim of persistence property analysis and ecotoxicological impact of veterinary pharmaceuticals on different terrestrial species, different classes of veterinary pharmaceuticals (n = 37) with soil degradation property (DT50) were gathered and subjected to QSAR and q-RASAR model development. The models were developed from 2D descriptors under organization for economic cooperation and development guidelines with the application of multiple linear regressions along with genetic algorithm. All developed QSAR and q-RASAR were statistically significant (Internal = R2adj: 0.721-0.861, Q2LOO: 0.609-0.757, and external = Q2Fn = 0.597-0.933, MAEext = 0.174-0.260). Further, the leverage approach of applicability domain assured the model's reliability. The veterinary pharmaceuticals with no experimental values were classified based on their persistence level. Further, the terrestrial toxicity analysis of persistent veterinary pharmaceuticals was done using toxicity prediction by computer assisted technology and in-house built quantitative structure toxicity relationship models to prioritize the toxic and persistent veterinary pharmaceuticals. This study will be helpful in estimation of persistence and toxicity of existing and upcoming veterinary pharmaceuticals.

Integrated predictive QSAR, Read Across, and q-RASAR analysis for diverse agrochemical phytotoxicity in oat and corn: A consensus-based approach for risk assessment and prioritization

In the modern fast-paced lifestyle, time-efficient and nutritionally rich foods like corn and oat have gained popularity for their amino acids and antioxidant contents. The increasing demand for these cereals necessitates higher production which leads to dependency on agrochemicals, which can pose health risks through residual present in the plant products. To first report the phytotoxicity for corn and oat, our study employs QSAR, quantitative Read-Across and quantitative RASAR (q-RASAR). All developed QSAR and q-RASAR models were equally robust (R2 = 0.680-0.762, Q2Loo = 0.593-0.693, Q2F1 = 0.680-0.860) and find their superiority in either oat or corn model, respectively, based on MAE criteria. AD and PRI had been performed which confirm the reliability and predictability of the models. The mechanistic interpretation reveals that the symmetrical arrangement of electronegative atoms and polar groups directly influences the toxicity of compounds. The final phytotoxicity and prioritization are performed by the consensus approach which results into selection of 15 most toxic compounds for both species.

Toxicity test profile for deep eutectic solvents: A detailed review and future prospects

Deep eutectic solvents (DESs) are preferable in terms of starting materials, storage and synthesis, simplicity, and component material affordability. In several industries ranging from chemical, electrochemical, biological, biotechnology, material science, etc., DES has demonstrated remarkable potential. Despite all these accomplishments, the safety issue with DES must be adequately addressed. Different DES interacts with the cellular membranes differently. It is not possible to classify all DES as easily biodegradable. By expanding the current understanding of the toxicity and biodegradation of DES, interactions between organisms and cellular membranes can be linked. The DES toxicity profile varies according to their concentration, the nature of the individual components, and how they interact with living things. Therefore, the results of this review can serve as a baseline for DES development in the future.

The use of simple structural parameters of organic compounds to assess their PUF-air partition coefficients

A novel approach is introduced for the reliable prediction of PUF-air partition coefficients of organic compounds, which can determine the environmental fate of organic compounds during interactions with air, soil, and water. The biggest accessible measured data of PUF-air partition coefficients for 170 chemicals are used to develop and test the novel model. In comparison to available quantitative structure-property relationship (QSPR) methods for the prediction of PUF-air partition coefficients that need complex descriptors, the here used descriptors are simpler. The assessed various statistical factors of the simple method containing 147 (training) and 23 (test) organic compounds can verify the external and internal cross-validations. Various statistical parameters confirm the high reliability of the novel model as compared with the outputs of complex multiple linear regression (MLR), artificial neural network (ANN) and support vector machine (SVM) methods. The values of R-squared (R2), and root mean square error (RMSE) of the new model are for training/test sets are 0.924/0.894 and 0.374/0.318, respectively. Meanwhile, R2 and RMSE values for three comparative models training/test sets are (i) MLR: 0.848/0.670 (R2) and 0.531/0.573 (RMSE); (ii) ANN: 0.902/0.664 (R2) and 0.425/0.560 (RMSE); (iii) SVM: 0.935/0.794 (R2) and 0.351/0.419 (RMSE). Thus, the new model the simplest approach with higher reliability in comparison to the best available methods.

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.

Aminopyrimidine Derivatives as Multiflavivirus Antiviral Compounds Identified from a Consensus Virtual Screening Approach

Around three billion people are at risk of infection by the dengue virus (DENV) and potentially other flaviviruses. Worldwide outbreaks of DENV, Zika virus (ZIKV), and yellow fever virus (YFV), the lack of antiviral drugs, and limitations on vaccine usage emphasize the need for novel antiviral research. Here, we propose a consensus virtual screening approach to discover potential protease inhibitors (NS3pro) against different flavivirus. We employed an in silico combination of a hologram quantitative structure-activity relationship (HQSAR) model and molecular docking on characterized binding sites followed by molecular dynamics (MD) simulations, which filtered a data set of 7.6 million compounds to 2,775 hits. Lastly, docking and MD simulations selected six final potential NS3pro inhibitors with stable interactions along the simulations. Five compounds had their antiviral activity confirmed against ZIKV, YFV, DENV-2, and DENV-3 (ranging from 4.21 ± 0.14 to 37.51 ± 0.8 μM), displaying aggregator characteristics for enzymatic inhibition against ZIKV NS3pro (ranging from 28 ± 7 to 70 ± 7 μM). Taken together, the compounds identified in this approach may contribute to the design of promising candidates to treat different flavivirus infections.

Application of in-silico drug discovery techniques to discover a novel hit for target-specific inhibition of SARS-CoV-2 Mpro's revealed allosteric binding with MAO-B receptor: A theoretical study to find a cure for post-covid neurological disorder

Several studies have revealed that SARS-CoV-2 damages brain function and produces significant neurological disability. The SARS-CoV-2 coronavirus, which causes COVID-19, may infect the heart, kidneys, and brain. Recent research suggests that monoamine oxidase B (MAO-B) may be involved in metabolomics variations in delirium-prone individuals and severe SARS-CoV-2 infection. In light of this situation, we have employed a variety of computational to develop suitable QSAR model using PyDescriptor and genetic algorithm-multilinear regression (GA-MLR) models (R2 = 0.800-793, Q2LOO = 0.734-0.727, and so on) on the data set of 106 molecules whose anti-SARS-CoV-2 activity was empirically determined. QSAR models generated follow OECD standards and are predictive. QSAR model descriptors were also observed in x-ray-resolved structures. After developing a QSAR model, we did a QSAR-based virtual screening on an in-house database of 200 compounds and found a potential hit molecule. The new hit's docking score (-8.208 kcal/mol) and PIC50 (7.85 M) demonstrated a significant affinity for SARS-CoV-2's main protease. Based on post-covid neurodegenerative episodes in Alzheimer's and Parkinson's-like disorders and MAO-B's role in neurodegeneration, the initially disclosed hit for the SARS-CoV-2 main protease was repurposed against the MAO-B receptor using receptor-based molecular docking, which yielded a docking score of -12.0 kcal/mol. This shows that the compound that inhibits SARS-CoV-2's primary protease may bind allosterically to the MAO-B receptor. We then did molecular dynamic simulations and MMGBSA tests to confirm molecular docking analyses and quantify binding free energy. The drug-receptor complex was stable during the 150-ns MD simulation. The first computational effort to show in-silico inhibition of SARS-CoV-2 Mpro and allosteric interaction of novel inhibitors with MAO-B in post-covid neurodegenerative symptoms and other disorders. The current study seeks a novel compound that inhibits SAR's COV-2 Mpro and perhaps binds MAO-B allosterically. Thus, this study will enable scientists design a new SARS-CoV-2 Mpro that inhibits the MAO-B receptor to treat post-covid neurological illness.

Design and Evaluation of Synthesized Pyrrole Derivatives as Dual COX-1 and COX-2 Inhibitors Using FB-QSAR Approach

This study delves into the intricate dynamics of the inflammatory response, unraveling the pivotal role played by cyclooxygenase (COX) enzymes, particularly COX-1 and COX-2 subtypes. Motivated by the pursuit of advancing scientific knowledge, our contribution to this field is marked by the design and synthesis of novel pyrrole derivatives. Crafted as potential inhibitors of COX-1 and COX-2 enzymes, our goal was to unearth molecules with heightened efficacy in modulating enzyme activity. A meticulous exploration of a synthesis library, housing around 3000 compounds, expedited the identification of potent candidates. Employing advanced docking studies and field-based Quantitative Structure-Activity Relationship (FB-QSAR) analyses enriched our understanding of the complex interactions between synthesized compounds and COX enzymes. Guided by FB-QSAR insights, our synthesis path led to the identification of compounds 4g, 4h, 4l, and 4k as potent COX-2 inhibitors, surpassing COX-1 efficacy. Conversely, compounds 5b and 5e exhibited heightened inhibitory activity against COX-1 relative to COX-2. The utilization of pyrrole derivatives as COX enzyme inhibitors holds promise for groundbreaking advancements in the domain of anti-inflammatory therapeutics, presenting avenues for innovative pharmaceutical exploration.

MASSA Algorithm: an automated rational sampling of training and test subsets for QSAR modeling

QSAR models capable of predicting biological, toxicity, and pharmacokinetic properties were widely used to search lead bioactive molecules in chemical databases. The dataset's preparation to build these models has a strong influence on the quality of the generated models, and sampling requires that the original dataset be divided into training (for model training) and test (for statistical evaluation) sets. This sampling can be done randomly or rationally, but the rational division is superior. In this paper, we present MASSA, a Python tool that can be used to automatically sample datasets by exploring the biological, physicochemical, and structural spaces of molecules using PCA, HCA, and K-modes. The proposed algorithm is very useful when the variables used for QSAR are not available or to construct multiple QSAR models with the same training and test sets, producing models with lower variability and better values for validation metrics. These results were obtained even when the descriptors used in the QSAR/QSPR were different from those used in the separation of training and test sets, indicating that this tool can be used to build models for more than one QSAR/QSPR technique. Finally, this tool also generates useful graphical representations that can provide insights into the data.

Application of QSAR models for acute toxicity of tetrazole compounds administrated orally and intraperitoneally in rat and mouse

Tetrazoles and their derivatives possess various biological activities, such as antibacterial, anti-fungal, and other activities. However, these compounds may induce specific cumulative and toxic effects in living organisms. Therefore, quantitative structure-activity relationship (QSAR) models were constructed to study the acute oral toxicity of tetrazoles in rats and mice. The toxicity data of 111 tetrazole compounds were collected using the ChemIDplus, ChEMBL and ECHA databases as response variables, while the PaDEL-descriptor generated the 2D descriptors as independent variables. The models were developed and validated following the OECD guidelines by the DTC-QSAR tool. Three QSAR models were successfully established for the oral routes of rat and mouse and the intraperitoneal route of mouse, respectively. The scatter plots showed high consistency between the training and test data sets. All the models successfully met the external and internal validation criteria. Most of the descriptors kept in the final models exhibited positive correlations with toxicity, whereas only 6 descriptors exhibited negative associations. Several chemicals were identified as response or structural outliers, based on the standardized residuals and leverage values. In conclusion, the findings of this investigation demonstrate that the proposed QSAR models hold promise in forecasting the acute toxicity of recently developed or synthesized tetrazole compounds, thereby mitigating potential risks to human health and the environment.

Quantitative Structure-Activity Relationship in the Series of 5-Ethyluridine, N2-Guanine, and 6-Oxopurine Derivatives with Pronounced Anti-Herpetic Activity

A quantitative analysis of the relationship between the structure and inhibitory activity against the herpes simplex virus thymidine kinase (HSV-TK) was performed for the series of 5-ethyluridine, N2-guanine, and 6-oxopurines derivatives with pronounced anti-herpetic activity (IC50 = 0.09 ÷ 160,000 μmol/L) using the GUSAR 2019 software. On the basis of the MNA and QNA descriptors and whole-molecule descriptors using the self-consistent regression, 12 statistically significant consensus models for predicting numerical pIC50 values were constructed. These models demonstrated high predictive accuracy for the training and test sets. Molecular fragments of HSV-1 and HSV-2 TK inhibitors that enhance or diminish the anti-herpetic activity are considered. Virtual screening of the ChEMBL database using the developed QSAR models revealed 42 new effective HSV-1 and HSV-2 TK inhibitors. These compounds are promising for further research. The obtained data open up new opportunities for developing novel effective inhibitors of TK.

Application of Machine Learning Methods to Predict the Air Half-Lives of Persistent Organic Pollutants

Persistent organic pollutants (POPs) are ubiquitous and bioaccumulative, posing potential and long-term threats to human health and the ecological environment. Quantitative structure-activity relationship (QSAR) studies play a guiding role in analyzing the toxicity and environmental fate of different organic pollutants. In the current work, five molecular descriptors are utilized to construct QSAR models for predicting the mean and maximum air half-lives of POPs, including specifically the energy of the highest occupied molecular orbital (HOMO_Energy_DMol3), a component of the dipole moment along the z-axis (Dipole_Z), fragment contribution to SAscore (SAscore_Fragments), subgraph counts (SC_3_P), and structural information content (SIC). The QSAR models were achieved through the application of three machine learning methods: partial least squares (PLS), multiple linear regression (MLR), and genetic function approximation (GFA). The determination coefficients (R2) and relative errors (RE) for the mean air half-life of each model are 0.916 and 3.489% (PLS), 0.939 and 5.048% (MLR), 0.938 and 5.131% (GFA), respectively. Similarly, the determination coefficients (R2) and RE for the maximum air half-life of each model are 0.915 and 5.629% (PLS), 0.940 and 10.090% (MLR), 0.939 and 11.172% (GFA), respectively. Furthermore, the mechanisms that elucidate the significant factors impacting the air half-lives of POPs have been explored. The three regression models show good predictive and extrapolation abilities for POPs within the application domain.

Analyzing disinfection by-products yield and mechanisms in UV/Cl2 using response surface methodology and quantitative structure-activity relationship models

The study aimed to investigate the formation of halogenated disinfection byproducts (DBPs) during applying UV/chlorine (UV/Cl2) and unravel the interactive impacts of critical operational parameters and the mechanisms behind DBPs formation. Response surface methodology and quantitative structure-activity relationship models were developed to evaluate the contribution of electrophilic, nucleophilic, and free radical reactions to the formation of DBPs in UV/Cl2. The study found that Cl2 and its interactions dominated the total DBPs and non-Br-DBPs formation, while Br- and the Cl2-Br- interaction played a decisive role in the Br-DBPs formation. The study also observed significant interactions of Br, Cl2, and pH on chloroform, bromodichloromethane, dichloroacetonitrile, 1,1-dichloro-2-propanone, trichloroactic acid, and chlorodibromoacetic acid formations, while no evident interaction on chloral hydrate, dibromochloromethane, trichloroacetone, dibromoacetic acid, and bromodichloroacetic acid formations. The electrophilic substitution of HOBr mainly controlled the formation of trihalomethanes, and the contribution of nucleophilic, electrophilic, and free radical (•OH, Cl•, Cl2•- and ClO•) reactions depended on the molar ratio of Cl2 to Br, and pH-determined hydrolysis rate constants of DBPs and the types of free radicals. Overall, the response surface methodology and quantitative structure-activity relationship models provided a reference for revealing DBPs formation mechanisms in other disinfection processes.

Multivariate QSAR, similarity search and ADMET studies based in a set of methylamine derivatives described as dopamine transporter inhibitors

The dopamine transporter (DAT), responsible for the regulation of dopaminergic neurotransmission, is implicated in the etiology of several neuropsychiatric disorders which, in turn, have contributed to high rates of disability and numerous deaths in recent years, significantly impacting the global health system. Although the research for new drugs for the treatment of neuropsychiatric disorders has evolved in recent years, the availability of DAT-selective drugs that do not generate the same psychostimulant effects observed in drugs of abuse remains scarce. Therefore, we performed a QSAR study based on a dataset of 36 methylamine derivatives described as DAT inhibitors. The model was obtained based only in descriptors derived from 2D structures, and it was validated and generated satisfactory results considering the metrics used for internal and external validation. Subsequently, a virtual screening step also based on 2D similarity was performed, where it was possible to identify a total of 1157 compounds. After a series of reductions of the set using toxicity filters, applicability domain evaluation, and pharmacokinetic properties in silico assessment, seven hit compounds were selected as the most promising to be used, in future studies, as new scaffolds for the development of new DAT inhibitors.

Machine learning-based models for accessing thermal conductivity of liquids at different temperature conditions

Combating global warming-related climate change demands prompt actions to reduce greenhouse gas emissions, particularly carbon dioxide. Biomass-based biofuels represent a promising alternative fossil energy source. To convert biomass into energy, numerous conversion processes are performed at high pressure and temperature conditions, and the design and dimensioning of such processes requires thermophysical property data, particularly thermal conductivity, which are not always available in the literature. In this paper, we proposed the application of Chemoinformatics methodologies to investigate the prediction of thermal conductivity for hydrocarbons and oxygenated compounds. A compilation of experimental data followed by a careful data curation were performed to establish a database. The support vector machine algorithm has been applied to the database leading to models with good predictive abilities. The support vector regression (SVR) model has then been applied to an external set of compounds, i.e. not considered during the training of models. It showed that our SVR model can be used for the prediction of thermal conductivity values for temperatures and/or compounds that are not covered experimentally in the literature.

QSPR models to predict the physical hazards of mixtures: a state of art

Physical hazards of chemical mixtures, associated for example with their fire or explosion risks, are generally characterized using experimental tools. These tests can be expensive, complex, long to perform and even dangerous for operators. Therefore, for several years and especially with the implementation of the REACH regulation, predictive methods like quantitative structure-property relationships have been encouraged as alternatives tests to determine (eco)toxicological but also physical hazards of chemical substances. Initially, these approaches were intended for pure products, by considering a molecular similarity principle. However, additional to those for pure products, QSPR models for mixtures recently appeared and represent an increasing field of research. This study proposes a state of the art of existing QSPR models specifically dedicated to the prediction of the physical hazards of mixtures. Identified models have been analysed on the key elements of model development (experimental data and fields of application, descriptors used, development and validation methods). It draws up an overview of the potential and limitations of current models as well as areas of progress towards enlarged deployment as a complement to experimental characterizations, for example in the search for safer substances (according to safety-by-design concepts).

Bidirectional selection of the functional properties and environmental friendliness of organophosphorus (OP) pesticide derivatives: Design, screening, and mechanism analysis

Organophosphorus pesticides (OPs) are widely used in agricultural production, but the resulting pollution and drug resistance have sparked widespread concern. Therefore, this paper built a model to design OP substitute molecules with high functionality and environmental friendliness, as well as conducted various human health and ecological environment evaluations, synthetic accessibility screening, and easy detection screening. The functionality of the two OP substitute molecules, DIM-100 and DIM-164, increased by 22.79 % and 22.18 %, respectively, and the environmental friendliness increased by 18.07 % and 24.02 %, respectively. The human health risk and ecological, environmental risks were significantly reduced. Both molecules are easy to synthesize, and their detection sensitivity is 9.85 % and 11.24 % higher than that of the target molecule, respectively. Furthermore, significant changes in the distribution of electrons and holes near the C8 and S1 atoms of the OP substitute molecule resulted in easier breakage of the C8-S1 bond, enhancing its photodegradation ability. The charge transfer ability between the atoms of the molecule (as increasing the electron-withdrawing group led to an increase in charge of the P atom) and the volume of the cholinesterase active pocket both affect the functionality of the DIM substitute molecule. That is, the volume of the cholinesterase active pocket of the bee is smaller than that of the brown planthopper and is more affected by the volume of the OP molecule. Furthermore, the mutual verification analysis of the bidirectional selectivity effect of OP substitute molecules between the BayesianRidge model and the 3D-QS(A² + ∀³)R model reveals that the overall charge transfer degree of DIM substitute molecules is the main reason for the increase in the bidirectional selectivity effect.

Ecotoxicological QSAR study of fused/non-fused polycyclic aromatic hydrocarbons (FNFPAHs): Assessment and priority ranking of the acute toxicity to Pimephales promelas by QSAR and consensus modeling methods

Fused/non-fused polycyclic aromatic hydrocarbons (FNFPAHs) have a variety of toxic effects on ecosystems and human body, but the acquisition of their toxicity data is greatly limited by the limited resources available. Here, we followed the EU REACH regulation and used Pimephales promelas as a model organism to investigate the quantitative structure-activity relationship (QSAR) between the FNFPAHs and their toxicity for the aquatic environment for the first time. We developed a single QSAR model (SM1) containing five simple and interpretable 2D molecular descriptors, which met the validation of OECD QSAR-related principles, and analyzed their mechanistic relationships with toxicity in detail. The model had good degree of fitting and robustness, and had better external prediction performance (MAE_test = 0.4219) than ECOSAR model (MAE_test = 0.5614). To further enhance its prediction accuracy, the three qualified single models (SMs) were used for constructing consensus models (CMs), the best one CM2 (MAE_test = 0.3954) had a significantly higher prediction accuracy for test compounds than SM1, and also outperformed the T.E.S.T. consensus model (MAE_test = 0.4233). Subsequently, the toxicity of 252 true external FNFPAHs from Pesticide Properties Database (PPDB) was predicted by SM1, the prediction results showed that 94.84 % compounds were reliably predicted within the model's application domain (AD). We also applied the best CM2 to predict the untested 252 FNFPAHs. Furthermore, we provided a mechanistic analysis and explanation for pesticides ranked as top 10 most toxic FNFPAHs. In summary, all developed QSAR and consensus models can be used as efficient tools for predicting the acute toxicity of unknown FNFPAHs to Pimephales promelas, thus being important for the risk assessment and regulation of FNFPAHs contamination in aquatic environment.

Mechanism analysis and improved molecular modification: Design of high efficiency and environmentally friendly triazole fungicide substitutes

The adverse effects of triazole fungicides (TFs) on the soil and the environmental damage caused by their residues have attracted the attention of the international community. To effectively prevent and control the above problems, this paper designed 72 substitutes of TFs with significantly better molecular functionality (>40%) using Paclobutrazol (PBZ) as the template molecule. Then, the comprehensive scores for environmental effects calculated after normalization by "extreme value method-entropy weight method-weighted average method" was the dependent variable, the structural parameters of TFs molecules was the independent variable (PBZ-214 was the template molecule) to construct the 3D-QSAR model of integrated environmental effects of TFs with high degradability, low bioenrichment, low endocrine disruption effects, and low hepatotoxicity and designed 46 substitutes of TFs with significantly better comprehensive environmental effects (>20%). After confirming the above effects of TFs and assessing human health risk and the universality of biodegradation and endocrine disruption, we screened PBZ-319-175 as the eco-friendly substitute of TF, which had high efficiency (improved functionality) and better environmental effects than those of the target molecule by 51.63% and 36.09%, respectively. Finally, the results of the molecular docking analysis showed that non-bonding interactions (hydrogen bonding, electrostatic, or polar force) predominantly affected the association between PBZ-319-175 and its biodegradable protein, and the hydrophobic effect of the amino acids distributed around PBZ-319-175 played a significant role. Additionally, we determined the microbial degradation path of PBZ-319-175 and found that the steric hindrance of the substituent group after molecular modification promoted its biodegradability. In this study, we enhanced molecular functionality twice and also reduce the major damage of TFs to the environment by performing iterative modifications. This paper provided theoretical support for the development and application of high-performance, eco-friendly substitutes of TFs.

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

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

First report on ecotoxicological QSTR and I-QSTR modeling for the prediction of acute ecotoxicity of diverse organic chemicals against three protozoan species

The recent years have witnessed an upsurge of interest to assess the toxicity of organic chemicals exhibiting harmful impacts on the environment. In this investigation, we have developed regression-based quantitative structure-toxicity relationship (QSTR) models against three protozoan species (Entosiphon sulcantum, Uronema parduczi, and Chilomonas paramecium) using three sets of descriptor combinations such as ETA indices only, non-ETA descriptors only, and both ETA and non-ETA descriptors to examine the key structural features that determine the toxic properties of protozoa. The interspecies models (i-QSTRs) were also generated for efficient data gap-filling of toxicity databases. The statistical results of the validated models in terms of both internal and external validation metrics suggested that the models are statistically reliable and robust. Additionally, using these validated models, we screened the DrugBank database containing 11,300 pharmaceuticals for assessing the ecotoxicological properties. The features appearing in the models suggested that nonpolar characteristics, electronegativity, hydrogen bonding, π-π, and hydrophobic interactions are responsible for chemical toxicity toward protozoan. The validated models may be utilized for the development of eco-friendly drugs & chemicals, data gap-filling of toxicity databases for regulatory purposes and research, as well as to decrease the use of toxic and hazardous chemicals in the environment.

A QSAR Study for Antileishmanial 2-Phenyl-2,3-dihydrobenzofurans †

Leishmaniasis, a parasitic disease that represents a threat to the life of millions of people around the globe, is currently lacking effective treatments. We have previously reported on the antileishmanial activity of a series of synthetic 2-phenyl-2,3-dihydrobenzofurans and some qualitative structure-activity relationships within this set of neolignan analogues. Therefore, in the present study, various quantitative structure-activity relationship (QSAR) models were created to explain and predict the antileishmanial activity of these compounds. Comparing the performance of QSAR models based on molecular descriptors and multiple linear regression, random forest, and support vector regression with models based on 3D molecular structures and their interaction fields (MIFs) with partial least squares regression, it turned out that the latter (i.e., 3D-QSAR models) were clearly superior to the former. MIF analysis for the best-performing and statistically most robust 3D-QSAR model revealed the most important structural features required for antileishmanial activity. Thus, this model can guide decision-making during further development by predicting the activity of potentially new leishmanicidal dihydrobenzofurans before synthesis.

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, g_m, 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.

Study of two combined series of triketones with HPPD inhibitory activity by molecular modelling

Triketones are suitable compounds for 4-hydroxyphenylpyruvate dioxygenase (HPPD) inhibition and are important compounds for eliminating agricultural weeds. We report herein quantitative structure-activity relationship (QSAR) modelling and docking studies for a series of triketone-quinoline hybrids and 2-(aryloxyacetyl)cyclohexane-1,3-diones with the aim of proposing new chemical candidates that exhibit improved performance as herbicides. The QSAR models obtained were reliable and predictive (average r², q², and r²_pred of 0.72, 0.51, and 0.71, respectively). Guided by multivariate image analysis of the PLS regression coefficients and variable importance in projection scores, the substituent effects could be analysed, and a promising derivative with R¹ = H, R² = CN, and R³ = 5,7,8-triCl at the triketone-quinoline scaffold (P18) was proposed. Docking studies demonstrated that π-π stacking interactions and specific interactions between the substituents and amino acid residues in the binding site of the Arabidopsis thaliana HPPD (AtHPPD) enzyme support the desired bioactivity. In addition, compared to a benchmark commercial triketone (mesotrione), the proposed compounds are more lipophilic and less mobile in soil rich in organic matter and are less prone to contaminate groundwater.

A QSAR study to predict the survival motor neuron promoter activity of candidate diaminoquinazoline derivatives for the potential treatment of spinal muscular atrophy

Spinal Muscular Atrophy is a genetic neuromuscular disease that leads to muscle weakness and atrophy and it is characterized by the loss of α-motor neurons in the spinal cord's anterior horn cells. The disease appears due to low levels of the survival motor neuron protein. There are continuing clinical trials for the treatment of Spinal Muscular Atrophy. Quinazoline-based compounds are promising since they were tested on fibroblasts derived from the patients and found to increase the survival motor neuron protein levels. In this study, using multiple linear regression, we generated robust and valid quantitative structure- activity relationship models to predict the survival motor neuron-2 promoter activity of the new candidate compounds using the experimental survival motor neuron-2 promoter activity values of 2,4-diaminoquinazoline derivatives taken from the literature. The novel compounds designed by combining the pyrido[1,2-α]pyrimidin-4-one moeity of the known drug Risdiplam with that of 2,4 - diaminoquinazoline scaffold were predicted to exhibit strong promoter activities.

Development of reliable quantitative structure-toxicity relationship models for toxicity prediction of benzene derivatives using semiempirical descriptors

The Health and environmental hazards of benzene and nitrobenzene (NB) derivatives have remained a topic of interest of researchers. In silico methods for prediction of toxicity of chemicals have proved their worth in accurate forecast of environmental as well as health toxicity and are strongly recommended by regulatory authorities. Two quantitative structure-toxicity relationship (QSTR) models explaining Scenedesmus obliquus toxicity trends among 39 benzene derivatives and Tetrahymena pyriformis toxicity of 103 NB and 392 benzene derivatives are developed using semiempirical quantum chemical parameters. The best constructed QSTR models have good fitting ability (R² = 0.8053, 0.7591, and 0.8283) and robustness (Q²_LOO = 0.7507, 0.7227, and 0.8194; Q²_LMO = 0.7338, 0.7153, and 0.8172). The external predictivity of all the models are quite good (R²_EXT = 0.8256, 0.9349, and 0.8698). Electronegativity, Cosmo volume, total energy, and molecular weight are responsible for the increase and decrease of toxicity of benzene derivatives against S. obliquus while electronegativity, electrophilicity index, the heat of formation, total energy, hydrophobicity, and cosmo volume are responsible for modulation of toxicity of NB and benzene derivatives toward T. pyriformis. These models fulfill the requirements of all the five OECD principles.