On the development and validation of QSAR models

Improved prediction of PFAS partitioning with PPLFERs and QSPRs

Per- and polyfluoroalkyl substances (PFAS) are chemicals of high concern and are undergoing hazard and risk assessment worldwide. Reliable physicochemical property (PCP) data are fundamental to assessments. However, experimental PCP data for PFAS are limited and property prediction tools such as quantitative structure-property relationships (QSPRs) therefore have poor predictive power for PFAS. New experimental data from Endo 2023 are used to improve QSPRs for predicting poly-parameter linear free energy relationship (PPLFER) descriptors for calculating water solubility (SW), vapor pressure (VP) and the octanol-water (KOW), octanol-air (KOA) and air-water (KAW) partition ratios. The new experimental data are only for neutral PFAS, and the QSPRs are only applicable to neutral chemicals. A key PPLFER descriptor for PFAS is the molar volume and this work compares different versions and makes recommendations for obtaining the best PCP predictions. The new models are included in the freely available IFSQSAR package (version 1.1.1), and property predictions are compared to those from the previous IFSQSAR (version 1.1.0) and from QSPRs in the US EPA's EPI Suite (version 4.11) and OPERA (version 2.9) models. The results from the new IFSQSAR models show improvements for predicting PFAS PCPs. The root mean squared error (RMSE) for predicting log KOWversus expected values from quantum chemical calculations was reduced by approximately 1 log unit whereas the RMSE for predicting log KAW and log KOA was reduced by 0.2 log units. IFSQSAR v.1.1.1 has an RMSE one or more log units lower than predictions from OPERA and EPI Suite when compared to expected values of log KOW, log KAW and log KOA for PFAS, except for EPI Suite predictions for log KOW which have a comparable RMSE. Recommendations for future experimental work for PPLFER descriptors for PFAS and future research to improve PCP predictions for PFAS are presented.

The database makes the poison: How the selection of datasets in QSAR models impacts toxicant prediction of higher tier endpoints

As the United States and the European Union continue their steady march towards the acceptance of new approach methodologies (NAMs), we need to ensure that the available tools are fit for purpose. Critics will be well-positioned to caution against NAMs acceptance and adoption if the tools turn out to be inadequate. In this paper, we focus on Quantitative Structure Activity-Relationship (QSAR) models and highlight how the training database affects quality and performance of these models. Our analysis goes to the point of asking, "are the endpoints extracted from the experimental studies in the database trustworthy, or are they false negatives/positives themselves?" We also discuss the impacts of chemistry on QSAR models, including issues with 2-D structure analyses when dealing with isomers, metabolism, and toxicokinetics. We close our analysis with a discussion of challenges associated with translational toxicology, specifically the lack of adverse outcome pathways/adverse outcome pathway networks (AOPs/AOPNs) for many higher tier endpoints. We recognize that it takes a collaborate effort to build better and higher quality QSAR models especially for higher tier toxicological endpoints. Hence, it is critical to bring toxicologists, statisticians, and machine learning specialists together to discuss and solve these challenges to get relevant predictions.

Design of Novel TRPA1 Agonists Based on Structure of Natural Vasodilator Carvacrol-In Vitro and In Silico Studies

Considering the escalating global prevalence and the huge therapeutic demand for the treatment of hypertension, there is a persistent need to identify novel target sites for vasodilator action. This study aimed to investigate the role of TRPA1 channels in carvacrol-induced vasodilation and to design novel compounds based on carvacrol structure with improved activities. In an isolated tissue bath experiment, it was shown that 1 µM of the selective TRPA1 antagonist A967079 significantly (p < 0.001) reduced vasodilation induced by 3 mM of carvacrol. A reliable 3D-QSAR model with good statistical parameters was created (R2 = 0.83; Q2 = 0.59 and Rpred2 = 0.84) using 29 TRPA1 agonists. Obtained results from this model were used for the design of novel TRPA1 activators, and to predict their activity against TRPA1. Predicted pEC50 activities of these molecules range between 4.996 to 5.235 compared to experimental pEC50 of 4.77 for carvacrol. Molecular docking studies showed that designed molecules interact with similar amino acid residues of the TRPA1 channel as carvacrol, with eight compounds showing lower binding energies. In conclusion, carvacrol-induced vasodilation is partly mediated by the activation of TRPA1 channels. Combining different in silico approaches pointed out that the molecule D27 (2-[2-(hydroxymethyl)-4-methylphenyl]acetamide) is the best candidate for further synthesis and experimental evaluation in in vitro conditions.

In-silico mechanistic analysis of adsorption of Iodinated Contrast Media agents on graphene surface

The study aims at assessing the potential of graphene-based adsorbents to reduce environmental impacts of Iodinated Contrast Media Agents (ICMs). We analyze an extensive collection of ICMs. A modeling approach resting on molecular docking and Density Functional Theory simulations is employed to examine the adsorption process at the molecular level. The study also relies on a Quantitative Structure-Activity Relationship (QSAR) modeling framework to correlate molecular properties with the adsorption energy (Ead) of ICMs, thus enabling identification of the key mechanisms underpinning adsorption and of the key factors contributing to it. A collection of distinct QSAR-based models is developed upon relying on Multiple Linear Regression and a standard genetic algorithm method. Having at our disposal multiple models enables us to take into account the uncertainty associated with model formulation. Maximum Likelihood and formal model identification/discrimination criteria (such as Bayesian and/or information theoretic criteria) are then employed to complement the traditional QSAR modeling phase. This has the advantage of (a) providing a rigorous ranking of the alternative models included in the selected set and (b) quantifying the relative degree of likelihood of each of these models through a weight or posterior probability. The resulting workflow of analysis enables one to seamlessly embed DFT and QSAR studies within a theoretical framework of analysis that explicitly takes into account model and parameter uncertainty. Our results suggest that graphene-based surfaces constitute a promising adsorbent for ICMs removal, π-π stacking being the primary mechanism behind ICM adsorption. Furthermore, our findings offer valuable insights into the potential of graphene-based adsorbent materials for effectively removing ICMs from water systems. They contribute to ascertain the significance of various factors (such as, e.g., the distribution of atomic van der Waals volumes, overall molecular complexity, the presence and arrangement of Iodine atoms, and the presence of polar functional groups) on the adsorption process.

Classification Model of Pesticide Toxicity in Americamysis bahia Based on Quantum Chemical Descriptors

A set of quantum chemical descriptors (molecular polarization, heat capacity, entropy, Mulliken net charge of the most positive hydrogen atom, APT charge of the most negative atom and APT charge of the most positive atom with hydrogen summed into heavy atoms) was successfully used to establish the classification models for the toxicity pLC50 of pesticides in Americamysis bahia. The optimal random forest model (Class Model A) yielded predictive accuracy of 100% (training set of 217 pesticides), 95.8% (test set of 72 pesticides) and 99.0% (total set of 289 pesticides), which were very satisfactory, compared with previous classification models reported for the toxicity of compounds in aquatic organisms. Therefore, it is reasonable to apply the quantum chemical descriptors associated with molecular structural information on molecular bulk, chemical reactivity and weak interactions, to develop classification models for the toxicity pLC50 of pesticides in A. bahia.

Probabilistic Reference and 10% Effect Concentrations for Characterizing Inhalation Non-cancer and Developmental/Reproductive Effects for 2,160 Substances

Chemicals assessment and management frameworks rely on regulatory toxicity values, which are based on points of departure (POD) identified following rigorous dose-response assessments. Yet, regulatory PODs and toxicity values for inhalation exposure (i.e., reference concentrations [RfCs]) are available for only ∼200 chemicals. To address this gap, we applied a workflow to determine surrogate inhalation route PODs and corresponding toxicity values, where regulatory assessments are lacking. We curated and selected inhalation in vivo data from the U.S. EPA's ToxValDB and adjusted reported effect values to chronic human equivalent benchmark concentrations (BMCh) following the WHO/IPCS framework. Using ToxValDB chemicals with existing PODs associated with regulatory toxicity values, we found that the 25th %-ile of a chemical's BMCh distribution (POD p 25 BMC h ) could serve as a suitable surrogate for regulatory PODs (Q2 ≥ 0.76, RSE ≤ 0.82 log10 units). We applied this approach to derive POD p 25 BMC h for 2,095 substances with general non-cancer toxicity effects and 638 substances with reproductive/developmental toxicity effects, yielding a total coverage of 2,160 substances. From these POD p 25 BMC h , we derived probabilistic RfCs and human population effect concentrations. With this work, we have expanded the number of chemicals with toxicity values available, thereby enabling a much broader coverage for inhalation risk and impact assessment.

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-α.

Andrographolide: A Diterpenoid from Cymbopogon schoenanthus Identified as a New Hit Compound against Trypanosoma cruzi Using Machine Learning and Experimental Approaches

American Trypanosomiasis, also known as Chagas disease, is caused by the protozoan Trypanosoma cruzi and exhibits limited options for treatment. Natural products offer various structurally complex metabolites with biological activities, including those with anti-T. cruzi potential. The discovery and development of prototypes based on natural products frequently display multiple phases that could be facilitated by machine learning techniques to provide a fast and efficient method for selecting new hit candidates. Using Random Forest and k-Nearest Neighbors, two models were constructed to predict the biological activity of natural products from plants against intracellular amastigotes of T. cruzi. The diterpenoid andrographolide was identified from a virtual screening as a promising hit compound. Hereafter, it was isolated from Cymbopogon schoenanthus and chemically characterized by spectral data analysis. Andrographolide was evaluated against trypomastigote and amastigote forms of T. cruzi, showing IC50 values of 29.4 and 2.9 μM, respectively, while the standard drug benznidazole displayed IC50 values of 17.7 and 5.0 μM, respectively. Additionally, the isolated compound exhibited a reduced cytotoxicity (CC50 = 92.8 μM) against mammalian cells and afforded a selectivity index (SI) of 32, similar to that of benznidazole (SI = 39). From the in silico analyses, we can conclude that andrographolide fulfills many requirements implemented by DNDi to be a hit compound. Therefore, this work successfully obtained machine learning models capable of predicting the activity of compounds against intracellular forms of T. cruzi.

Application of machine learning-based read-across structure-property relationship (RASPR) as a new tool for predictive modelling: Prediction of power conversion efficiency (PCE) for selected classes of organic dyes in dye-sensitized solar cells (DSSCs)

The application of various in-silico-based approaches for the prediction of various properties of materials has been an effective alternative to experimental methods. Recently, the concepts of Quantitative structure-property relationship (QSPR) and read-across (RA) methods were merged to develop a new emerging chemoinformatic tool: read-across structure-property relationship (RASPR). The RASPR method can be applicable to both large and small datasets as it uses various similarity and error-based measures. It has also been observed that RASPR models tend to have an increased external predictivity compared to the corresponding QSPR models. In this study, we have modeled the power conversion efficiency (PCE) of organic dyes used in dye-sensitized solar cells (DSSCs) by using the quantitative RASPR (q-RASPR) method. We have used relatively larger classes of organic dyes-Phenothiazines (n=207), Porphyrins (n=281), and Triphenylamines (n=229) for the modelling purpose. We have divided each of the datasets into training and test sets in 3 different combinations, and with the training sets we have developed three different QSPR models with structural and physicochemical descriptors and validated them with the corresponding test sets. These corresponding modeled descriptors were used to calculate the RASPR descriptors using a Java-based tool RASAR Descriptor Calculator v2.0 (https://sites.google.com/jadavpuruniversity.in/dtc-lab-software/home), and then data fusion was performed by pooling the previously selected structural and physicochemical descriptors with the calculated RASPR descriptors. Further feature selection algorithm was employed to develop the final RASPR PLS models. Here, we also developed different machine learning (ML) models with the descriptors selected in the QSPR PLS and RASPR PLS models, and it was found that models with RASPR descriptors superseded in external predictivity the models with only structural and physicochemical descriptors: RMSEP reduced for phenothiazines from 1.16-1.25 to 1.07-1.18, for porphyrins from 1.60-1.79 to 1.45-1.53, for triphenylamines from 1.27-1.54 to 1.20-1.47.

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.

Correlation intensity index-index of ideality of correlation: A hyphenated target function for furtherance of MAO-B inhibitory activity assessment

Monoamine oxidases are the enzymes involved in the management of brain homeostasis through oxidative deamination of monoamines such as neurotransmitters, tyramine etc. The excessive production of monoamine oxidase-B specifically results in numerous neurodegenerative disorders like Alzheimer's and Parkinson's diseases. Inhibitors of monoamine oxidase-B are applied in the management of these disorders. Here in this article we have developed robust hybrid descriptor based QSAR models related to 123 monoamine oxidase-B inhibitors through CORAL software by means of Monte Carlo optimization method. Three target functions were applied to prepare QSAR models and three splits were made for each target function. The most reliable, robust and better predictive QSAR models were developed with TF3 (correlation intensity index -index of ideality of correlation). Correlation intensity index showed positive effect on QSAR models. The structural features obtained from the QSAR modeling were incorporated in newly designed molecules and exhibited positive effect on their endpoint. Significant binding interactions were represented by these molecules in docking studies. Molecule B5 displayed prominent pIC50 (8.3) and binding affinity (-11.5 kcal mol-1) towards monoamine oxidase-B.

Integration of machine learning in 3D-QSAR CoMSIA models for the identification of lipid antioxidant peptides

The comparative molecular similarity indices analysis (CoMSIA) method is a widely used 3D-quantitative structure-activity relationship (QSAR) approach in the field of medicinal chemistry and drug design. However, relying solely on the Partial Least Square algorithm to build models using numerous CoMSIA indices has, in some cases, led to statistically underperforming models. This issue has also affected 3D-CoMSIA models constructed for the ferric thiocyanate (FTC) dataset from linoleic antioxidant measurements. In this study, a novel modeling routine has been developed incorporating various machine learning (ML) techniques to explore different options for feature selection, model fitting, and tuning algorithms with the ultimate goal of arriving at optimal 3D-CoMSIA models with high predictivity for the FTC activity. Recursive Feature Selection and SelectFromModel techniques were applied for feature selection, resulting in a significant improvement in model fitting and predictivity (R2, RCV2, and R2_test) of 24 estimators. However, these selection methods did not fully address the problem of overfitting and, in some instances, even exacerbated it. On the other hand, hyperparameter tuning for tree-based models resulted in dissimilar levels of model generalization for four tree-based models. GB-RFE coupled with GBR (hyperparameters: learning_rate = 0.01, max_depth = 2, n_estimators = 500, subsample = 0.5) was the only combination that effectively mitigated overfitting and demonstrated superior performance (RCV2 of 0.690, R2_test of 0.759, and R2 of 0.872) compared to the best linear model, PLS (with RCV2 of 0.653, R2_test of 0.575, and R2 of 0.755). Therefore, it was subsequently utilized to screen potential antioxidants among a range of Tryptophyllin L tripeptide fragments, leading to the synthesis and testing of three peptides: F-P-5Htp, F-P-W, and P-5Htp-L. These peptides exhibited promising activity levels, with FTC values of 4.2 ± 0.12, 4.4 ± 0.11, and 1.72 ± 0.15, respectively.

Imagining emotional future events in PTSD: clinical and neurocognitive correlates

Emotional future thinking serves important functions related to goal pursuit and emotion regulation but has been scantly studied in posttraumatic stress disorder (PTSD). The current study sought to characterize emotional future thinking in PTSD and to identify clinical and neurocognitive profiles associated with potential alterations in the level of detail in narratives of imagined future events. Fifty-eight, trauma-exposed, war-zone veterans, who were classified into current PTSD, past PTSD, and no-PTSD groups, were asked to vividly imagine future events in response to positive and negative cue words occurring in the near and distant future. These narratives were scored for internal (i.e., pertaining to the main event) and external (i.e., tangential to the main event) details. Participants also performed neurocognitive tasks of generative ability, working memory, and relational verbal memory. Linear mixed modeling revealed that the current and past PTSD groups generated fewer internal details than the no-PTSD group across positive and negative cue words and across temporal proximity. Partial least squares analysis revealed that symptom severity for all PTSD clusters was inversely associated with production of internal details, albeit with the association relatively weaker for intrusion symptoms. Among the neurocognitive tasks, only relational verbal memory was associated with production of internal details. These findings suggest, as predicted, that functional avoidance may underlie reduced detail generation but also point to potential additional mechanisms to be further investigated. That future event simulation remains overgeneral even when PTSD symptoms abate highlights the importance of addressing alterations in future thinking in this population.

Structure-Activity Relationship Studies of Chalcones and Diarylpentanoids with Antitumor Activity: Potency and Selectivity Optimization

We previously reported that chalcone CM-M345 (1) and diarylpentanoid BP-C4 (2) induced p53-dependent growth inhibitory activity in human cancer cells. Herein, CM-M345 (1) and BP-C4 (2) analogues were designed and synthesized in order to obtain more potent and selective compounds. Compounds 16, 17, 19, 20, and 22-24 caused pronounced in vitro growth inhibitory activity in HCT116 cells (0.09 < GI50 < 3.10 μM). Chemical optimization of CM-M345 (1) led to the identification of compound 36 with increased selectivity for HCT116 cells expressing wild-type p53 compared to its p53-null isogenic derivative and low toxicity to non-tumor HFF-1 cells. The molecular modification of BP-C4 (2) resulted in the discovery of compound 16 with more pronounced antiproliferative activity and being selective for HCT116 cells with p53, as well as 17 with enhanced antiproliferative activity against HCT116 cells and low toxicity to non-tumor cells. Compound 16 behaved as an inhibitor of p53-MDM2 interaction, and compound 17 was shown to induce apoptosis, associated with an increase in cleaved PARP and decreased levels of the anti-apoptotic protein Bcl-2. In silico studies allowed us to predict the druglikeness and ADMET properties for 16 and 17. Docking and molecular dynamics studies predicted that 16 could bind stably to the MDM2 binding pocket.

Thiazolidinedione-triazole conjugates: design, synthesis and probing of the α-amylase inhibitory potential

Aim: The primary objective of this investigation was the synthesis, spectral interpretation and evaluation of the α-amylase inhibition of rationally designed thiazolidinedione-triazole conjugates (7a-7aa). Materials & methods: The designed compounds were synthesized by stirring a mixture of thiazolidine-2,4-dione, propargyl bromide, cinnamaldehyde and azide derivatives in polyethylene glycol-400. The α-amylase inhibitory activity of the synthesized conjugates was examined by integrating in vitro and in silico studies. Results: The investigated derivatives exhibited promising α-amylase inhibitory activity, with IC50 values ranging between 0.028 and 0.088 μmol ml-1. Various computational approaches were employed to get detailed information about the inhibition mechanism. Conclusion: The thiazolidinedione-triazole conjugate 7p, with IC50 = 0.028 μmol ml-1, was identified as the best hit for inhibiting α-amylase.

Garbage in, garbage out: how reliable training data improved a virtual screening approach against SARS-CoV-2 MPro

Introduction: The identification of chemical compounds that interfere with SARS-CoV-2 replication continues to be a priority in several academic and pharmaceutical laboratories. Computational tools and approaches have the power to integrate, process and analyze multiple data in a short time. However, these initiatives may yield unrealistic results if the applied models are not inferred from reliable data and the resulting predictions are not confirmed by experimental evidence. Methods: We undertook a drug discovery campaign against the essential major protease (MPro) from SARS-CoV-2, which relied on an in silico search strategy -performed in a large and diverse chemolibrary- complemented by experimental validation. The computational method comprises a recently reported ligand-based approach developed upon refinement/learning cycles, and structure-based approximations. Search models were applied to both retrospective (in silico) and prospective (experimentally confirmed) screening. Results: The first generation of ligand-based models were fed by data, which to a great extent, had not been published in peer-reviewed articles. The first screening campaign performed with 188 compounds (46 in silico hits and 100 analogues, and 40 unrelated compounds: flavonols and pyrazoles) yielded three hits against MPro (IC₅₀ ≤ 25 μM): two analogues of in silico hits (one glycoside and one benzo-thiazol) and one flavonol. A second generation of ligand-based models was developed based on this negative information and newly published peer-reviewed data for MPro inhibitors. This led to 43 new hit candidates belonging to different chemical families. From 45 compounds (28 in silico hits and 17 related analogues) tested in the second screening campaign, eight inhibited MPro with IC₅₀ = 0.12-20 μM and five of them also impaired the proliferation of SARS-CoV-2 in Vero cells (EC₅₀ 7-45 μM). Discussion: Our study provides an example of a virtuous loop between computational and experimental approaches applied to target-focused drug discovery against a major and global pathogen, reaffirming the well-known "garbage in, garbage out" machine learning principle.

ProfhEX: AI-based platform for small molecules liability profiling

Off-target drug interactions are a major reason for candidate failure in the drug discovery process. Anticipating potential drug's adverse effects in the early stages is necessary to minimize health risks to patients, animal testing, and economical costs. With the constantly increasing size of virtual screening libraries, AI-driven methods can be exploited as first-tier screening tools to provide liability estimation for drug candidates. In this work we present ProfhEX, an AI-driven suite of 46 OECD-compliant machine learning models that can profile small molecules on 7 relevant liability groups: cardiovascular, central nervous system, gastrointestinal, endocrine, renal, pulmonary and immune system toxicities. Experimental affinity data was collected from public and commercial data sources. The entire chemical space comprised 289'202 activity data for a total of 210'116 unique compounds, spanning over 46 targets with dataset sizes ranging from 819 to 18896. Gradient boosting and random forest algorithms were initially employed and ensembled for the selection of a champion model. Models were validated according to the OECD principles, including robust internal (cross validation, bootstrap, y-scrambling) and external validation. Champion models achieved an average Pearson correlation coefficient of 0.84 (SD of 0.05), an R² determination coefficient of 0.68 (SD = 0.1) and a root mean squared error of 0.69 (SD of 0.08). All liability groups showed good hit-detection power with an average enrichment factor at 5% of 13.1 (SD of 4.5) and AUC of 0.92 (SD of 0.05). Benchmarking against already existing tools demonstrated the predictive power of ProfhEX models for large-scale liability profiling. This platform will be further expanded with the inclusion of new targets and through complementary modelling approaches, such as structure and pharmacophore-based models. ProfhEX is freely accessible at the following address: https://profhex.exscalate.eu/ .

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.

Identification of Potential p38γ Inhibitors via In Silico Screening, In Vitro Bioassay and Molecular Dynamics Simulation Studies

Protein kinase p38γ is an attractive target against cancer because it plays a pivotal role in cancer cell proliferation by phosphorylating the retinoblastoma tumour suppressor protein. Therefore, inhibition of p38γ with active small molecules represents an attractive alternative for developing anti-cancer drugs. In this work, we present a rigorous and systematic virtual screening framework to identify potential p38γ inhibitors against cancer. We combined the use of machine learning-based quantitative structure activity relationship modelling with conventional computer-aided drug discovery techniques, namely molecular docking and ligand-based methods, to identify potential p38γ inhibitors. The hit compounds were filtered using negative design techniques and then assessed for their binding stability with p38γ through molecular dynamics simulations. To this end, we identified a promising compound that inhibits p38γ activity at nanomolar concentrations and hepatocellular carcinoma cell growth in vitro in the low micromolar range. This hit compound could serve as a potential scaffold for further development of a potent p38γ inhibitor against cancer.

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.

Development and Validation of Predictive Quantitative Structure-Activity Relationship Models for Estrogenic Activities of Hydroxylated Polychlorinated Biphenyls

Disruption of the endocrine system by hydroxylated polychlorinated biphenyls (OH-PCBs) is hypothesized, among other potential mechanisms, to be mediated via nuclear receptor binding. Due to the high cost and lengthy time required to produce high-quality experimental data, empirical data to support the nuclear receptor binding hypothesis are in short supply. In the present study, two quantitative structure-activity relationship models were developed for predicting the estrogenic activities of OH-PCBs. Findings revealed that model I (for the estrogen receptor α dataset) contained five two-dimensional (2D) descriptors belonging to the classes autocorrelation, Burden modified eigenvalues, chi path, and atom type electrotopological state, whereas model II (for the estrogen receptor β dataset) contained three 2D and three 3D descriptors belonging to the classes autocorrelation, atom type electrotopological state, and Radial Distribution Function descriptors. The internal and external validation metrics reported for models I and II indicate that both models are robust, reliable, and suitable for predicting the estrogenic activities of untested OH-PCB congeners. Environ Toxicol Chem 2023;42:823-834. © 2023 SETAC.

Efficient predictions of cytotoxicity of TiO2-based multi-component nanoparticles using a machine learning-based q-RASAR approach

The availability of experimental nanotoxicity data is in general limited which warrants both the use of in silico methods for data gap filling and exploring novel methods for effective modeling. Read-Across Structure-Activity Relationship (RASAR) is an emerging cheminformatic approach that combines the usefulness of a QSAR model and similarity-based Read-Across predictions. In this work, we have generated simple, interpretable, and transferable quantitative-RASAR (q-RASAR) models which can efficiently predict the cytotoxicity of TiO₂-based multi-component nanoparticles. A data set of 29 TiO₂-based nanoparticles with specific amounts of noble metal precursors was rationally divided into training and test sets, and the Read-Across-based predictions for the test set were generated. The optimized hyperparameters and the similarity approach, which yield the best predictions, were used to calculate the similarity and error-based RASAR descriptors. A data fusion of the RASAR descriptors with the chemical descriptors was done followed by the best subset feature selection. The final set of selected descriptors was used to develop the q-RASAR models, which were validated using the stringent OECD criteria. Finally, a random forest model was also developed with the selected descriptors, which could efficiently predict the cytotoxicity of TiO₂-based multi-component nanoparticles superseding previously reported models in the prediction quality thus showing the merits of the q-RASAR approach. To further evaluate the usefulness of the approach, we have applied the q-RASAR approach also to a second cytotoxicity data set of 34 heterogeneous TiO₂-based nanoparticles which further confirmed the enhancement of external prediction quality of QSAR models after incorporation of RASAR descriptors.

Unsupervised machine learning, QSAR modelling and web tool development for streamlining the lead identification process of antimalarial flavonoids

Identification of lead compounds with the traditional laboratory approach is expensive and time-consuming. Nowadays, in silico techniques have emerged as a promising approach for lead identification. In this study, we aim to develop robust and predictive 2D-QSAR models to identify lead flavonoids by predicting the IC₅₀ against Plasmodium falciparum. We applied machine learning algorithms (Principal component analysis followed by K-means clustering) and Pearson correlation analysis to select 9 molecular descriptors (MDs) for model building. We selected and validated the three best QSAR models after execution of multiple linear regression (MLR) 100 times with different combinations of MDs. The developed models have fulfilled the five principles for QSAR models as specified by the Organization for Economic Co-operation and Development. The outcome of the study is a reliable and sustainable in silico method of IC₅₀ (Mean ± SD) prediction that will positively impact the antimalarial drug development process by reducing the money and time required to identify potential antimalarial lead compounds from the class of flavonoids. We also developed a web tool (JazQSAR, https://etflin.com/news/4) to offer an easily accessible platform for the developed QSAR models.

Insight parameter drug design for human β-tryptase inhibition integrated molecular docking, QSAR, molecular dynamics simulation, and pharmacophore modelling studies of α-keto-[1,2,4]-oxadiazoles

Dengue hemorrhagic fever (DHF) is severe dengue with a hallmark of vascular leakage. β-tryptase has been found to promote vascular leakage in DHF patients, which could be a potential target for DHF treatment. This study aims to develop a theoretical background for designing and selecting human β-tryptase inhibitors through computational studies. Thirty-four α-keto-[1,2,3]-oxadiazoles scaffold-based compounds were used to generate 2D-QSAR models and for molecular docking studies with β-tryptase (PDB Code 4A6L). In addition, molecular dynamics (MD) simulation and molecular mechanics generalised born surface area (MM-GBSA) analysis on the binding of the reported most active compound, compound 11e, towards β-tryptase were performed. Finally, a structure-based pharmacophore model was generated. The selected 2D-QSAR models have statistically proven good models by internal and external validation as well as the y-randomization test. The docking results of compound 11e showed lower CDOCKER energy than the 4A6L co-crystallised ligand and a similar binding pattern as the 4A6L co-crystallised ligand. From molecular dynamics simulation, 4A6L in compound 11e bound state has RMSD below 2 Å throughout the 500 ns simulation, indicating the docked complex is stable. Besides, MM-GBSA analysis suggested the 4A6L-compound 11e docked complex (-66.04 Kcal/mol) is structurally as stable as the 4A6L-native ligand co-crystallized structure (-66.84 Kcal/mol). The best pharmacophore model identified features included hydrogen bond acceptor, ionic interaction, hydrophobic interaction, and aromatic ring, which contribute to the inhibitory potency of a compound. This study supplied insight and knowledge for developing novel chemical compounds with improved inhibition of β-tryptase.Communicated by Ramaswamy H. Sarma.

SIRT2i_Predictor: A Machine Learning-Based Tool to Facilitate the Discovery of Novel SIRT2 Inhibitors

A growing body of preclinical evidence recognized selective sirtuin 2 (SIRT2) inhibitors as novel therapeutics for treatment of age-related diseases. However, none of the SIRT2 inhibitors have reached clinical trials yet. Transformative potential of machine learning (ML) in early stages of drug discovery has been witnessed by widespread adoption of these techniques in recent years. Despite great potential, there is a lack of robust and large-scale ML models for discovery of novel SIRT2 inhibitors. In order to support virtual screening (VS), lead optimization, or facilitate the selection of SIRT2 inhibitors for experimental evaluation, a machine-learning-based tool titled SIRT2i_Predictor was developed. The tool was built on a panel of high-quality ML regression and classification-based models for prediction of inhibitor potency and SIRT1-3 isoform selectivity. State-of-the-art ML algorithms were used to train the models on a large and diverse dataset containing 1797 compounds. Benchmarking against structure-based VS protocol indicated comparable coverage of chemical space with great gain in speed. The tool was applied to screen the in-house database of compounds, corroborating the utility in the prioritization of compounds for costly in vitro screening campaigns. The easy-to-use web-based interface makes SIRT2i_Predictor a convenient tool for the wider community. The SIRT2i_Predictor's source code is made available online.

Pharmacophore Synergism in Diverse Scaffold Clinches in Aurora Kinase B

Aurora kinase B (AKB) is a crucial signaling kinase with an important role in cell division. Therefore, inhibition of AKB is an attractive approach to the treatment of cancer. In the present work, extensive quantitative structure-activity relationships (QSAR) analysis has been performed using a set of 561 structurally diverse aurora kinase B inhibitors. The Organization for Economic Cooperation and Development (OECD) guidelines were used to develop a QSAR model that has high statistical performance (R²_tr = 0.815, Q²_LMO = 0.808, R²_ex = 0.814, CCC_ex = 0.899). The seven-variable-based newly developed QSAR model has an excellent balance of external predictive ability (Predictive QSAR) and mechanistic interpretation (Mechanistic QSAR). The QSAR analysis successfully identifies not only the visible pharmacophoric features but also the hidden features. The analysis indicates that the lipophilic and polar groups-especially the H-bond capable groups-must be present at a specific distance from each other. Moreover, the ring nitrogen and ring carbon atoms play important roles in determining the inhibitory activity for AKB. The analysis effectively captures reported as well as unreported pharmacophoric features. The results of the present analysis are also supported by the reported crystal structures of inhibitors bound to AKB.

White matter microstructure and sleep-wake disturbances in individuals at ultra-high risk of psychosis

Aim

White matter changes in individuals at ultra-high risk for psychosis (UHR) may be involved in the transition to psychosis. Sleep-wake disturbances commonly precede the first psychotic episode and predict development of psychosis. We examined associations between white matter microstructure and sleep-wake disturbances in UHR individuals compared to healthy controls (HC), as well as explored the confounding effect of medication, substance use, and level of psychopathology.

Methods

Sixty-four UHR individuals and 35 HC underwent clinical interviews and diffusion weighted imaging. Group differences on global and callosal mean fractional anisotropy (FA) was tested using general linear modeling. Sleep-wake disturbances were evaluated using the subjective measures disturbed sleep index (DSI) and disturbed awakening index (AWI) from the Karolinska Sleep Questionnaire, supported by objective sleep measures from one-night actigraphy. The primary analyses comprised partial correlation analyses between global FA/callosal FA and sleep-wake measures. Secondary analyses investigated multivariate patterns of covariance between measures of sleep-wake disturbances and FA in 48 white matter regions of interest using partial least square correlations.

Results

Ultra-high risk for psychosis individuals displayed lower global FA (F = 14.56, p < 0.001) and lower callosal FA (F = 11.34, p = 0.001) compared to HC. Subjective sleep-wake disturbances were significantly higher among the UHR individuals (DSI: F = 27.59, p < 0.001, AWI: F = 36.42, p < 0.001). Lower callosal FA was correlated with increased wake after sleep onset (r = -0.34, p = 0.011) and increased sleep fragmentation index (r = -0.31, p = 0.019) in UHR individuals. Multivariate analyses identified a pattern of covariance in regional FA which were associated with DSI and AWI in UHR individuals (p = 0.028), but not in HC. Substance use, sleep medication and antipsychotic medication did not significantly confound these associations. The association with objective sleep-wake measures was sustained when controlling for level of depressive and UHR symptoms, but symptom level confounded the covariation between FA and subjective sleep-wake measures in the multivariate analyses.

Conclusion

Compromised callosal microstructure in UHR individuals was related to objectively observed disruptions in sleep-wake functioning. Lower FA in ventrally located regions was associated with subjectively measured sleep-wake disturbances and was partly explained by psychopathology. These findings call for further investigation of sleep disturbances as a potential treatment target.

Rational Design of Chitin Deacetylase Inhibitors for Sustainable Agricultural Use Based on Molecular Topology

Fungicide resistance is a major concern in modern agriculture; therefore, there is a pressing demand to develop new, greener chemicals. Chitin is a major component of the fungal cell wall and a well-known elicitor of plant immunity. To overcome chitin recognition, fungal pathogens developed different strategies, with chitin deacetylase (CDA) activity being the most conserved. This enzyme is responsible for hydrolyzing the N-acetamido group in N-acetylglucosamine units of chitin to convert it to chitosan, a compound that can no longer be recognized by the plant. In previous works, we observed that treatments with CDA inhibitors, such as carboxylic acids, reduced the symptoms of cucurbit powdery mildew and induced rapid activation of chitin-triggered immunity, indicating that CDA could be an interesting target for fungicide development. In this work, we developed an in silico strategy based on QSAR (quantitative structure-activity relationship) and molecular topology (MT) to discover new, specific, and potent CAD inhibitors. Starting with the chemical structures of few carboxylic acids, with and without disease control activity, three predictive equations based on the MT paradigm were developed to identify a group of potential molecules. Their fungicidal activity was experimentally tested, and their specificity as CDA inhibitors was studied for the three best candidates by molecular docking simulations. To our knowledge, this is the first time that MT has been used for the identification of potential CDA inhibitors to be used against resistant powdery mildew strains. In this sense, we consider of special interest the discovery of molecules capable of stimulating the immune system of plants by triggering a defensive response against fungal species that are highly resistant to fungicides such as powdery mildew.

QSAR models reveal new EPAC-selective allosteric modulators

Exchange proteins directly activated by cAMP (EPAC) are guanine nucleotide exchange factors for the small GTPases, Rap1 and Rap2. They regulate several physiological functions and mitigation of their activity has been suggested as a possible treatment for multiple diseases such as cardiomyopathy, diabetes, chronic pain, and cancer. Several EPAC-specific modulators have been developed, however studies that quantify their structure-activity relationships are still lacking. Here we propose a quantitative structure-activity relationship (QSAR) model for a series of EPAC-specific compounds. The model demonstrated high reproducibility and predictivity and the predictive ability of the model was tested against a series of compounds that were unknown to the model. The compound with the highest predicted affinity was validated experimentally through fluorescence-based competition assays and NMR experiments revealed its mode of binding and mechanism of action as a partial agonist. The proposed QSAR model can, therefore, serve as an effective screening tool to identify promising EPAC-selective drug leads with enhanced potency.

Structure-activity relationship of graphene-related materials: A meta-analysis based on mammalian in vitro toxicity data

To support a safe application of graphene-related materials (GRMs) it is necessary to understand the potential negative impacts they could have on human health, in particular on the lung - one of the most sensitive exposure routes. Machine learning (ML) approaches can help analyse the results of multiple toxicity studies to understand the structure-activity relationship and the effect of experimental conditions, thus supporting predictive nanotoxicology. In this work we collected in vitro cytotoxicity data obtained from studies using lung cells; we then fitted multiple regression models to predict this endpoint based on the material properties and experimental conditions. Moreover, the data set was used to calculate the Benchmark Dose Lower Confidence Interval (BMDL), a dose descriptor widely used in risk assessment. Regression and classification models were applied for the prediction of the BMDL value and BMDL range. The analyses show that both cytotoxicity and the BMDL range can be predicted well (Q² = 0.77 and accuracy = 0.71, respectively). Both physico-chemical characteristics such as the lateral size, number of layers, and functionalization, and experimental conditions such as the assay and media used were important predicting features, confirming the need for thorough characterization and reporting of these parameters.

Repurposing food molecules as a potential BACE1 inhibitor for Alzheimer's disease

Alzheimer's disease (AD) is a severe neurodegenerative disorder of the brain that manifests as dementia, disorientation, difficulty in speech, and progressive cognitive and behavioral impairment. The emerging therapeutic approach to AD management is the inhibition of β-site APP cleaving enzyme-1 (BACE1), known to be one of the two aspartyl proteases that cleave β-amyloid precursor protein (APP). Studies confirmed the association of high BACE1 activity with the proficiency in the formation of β-amyloid-containing neurotic plaques, the characteristics of AD. Only a few FDA-approved BACE1 inhibitors are available in the market, but their adverse off-target effects limit their usage. In this paper, we have used both ligand-based and target-based approaches for drug design. The QSAR study entails creating a multivariate GA-MLR (Genetic Algorithm-Multilinear Regression) model using 552 molecules with acceptable statistical performance (R 2 = 0.82, Q 2 loo = 0.81). According to the QSAR study, the activity has a strong link with various atoms such as aromatic carbons and ring Sulfur, acceptor atoms, sp2-hybridized oxygen, etc. Following that, a database of 26,467 food compounds was primarily used for QSAR-based virtual screening accompanied by the application of the Lipinski rule of five; the elimination of duplicates, salts, and metal derivatives resulted in a truncated dataset of 8,453 molecules. The molecular descriptor was calculated and a well-validated 6-parametric version of the QSAR model was used to predict the bioactivity of the 8,453 food compounds. Following this, the food compounds whose predicted activity (pKi) was observed above 7.0 M were further docked into the BACE1 receptor which gave rise to the Identification of 4-(3,4-Dihydroxyphenyl)-2-hydroxy-1H-phenalen-1-one (PubChem I.D: 4468; Food I.D: FDB017657) as a hit molecule (Binding Affinity = -8.9 kcal/mol, pKi = 7.97 nM, Ki = 10.715 M). Furthermore, molecular dynamics simulation for 150 ns and molecular mechanics generalized born and surface area (MMGBSA) study aided in identifying structural motifs involved in interactions with the BACE1 enzyme. Molecular docking and QSAR yielded complementary and congruent results. The validated analyses can be used to improve a drug/lead candidate's inhibitory efficacy against the BACE1. Thus, our approach is expected to widen the field of study of repurposing nutraceuticals into neuroprotective as well as anti-cancer and anti-viral therapeutic interventions.

Influence of chemical structure of organic micropollutants on the degradability with ozonation

The increasing environmental problems due to various organic micropollutants in water cause the search of suitable additional water treatment methods. Gaining experimental data for the large amount and variety of pollutants would consume a lot of time as well as economic and ecologic resources. An alternative approach is predictive quantitative structure-property relationship (QSPR) modeling, which establishes a correlation between the structural properties of a molecules with a biological, physical, or chemical property. Therefore, in this study, QSPR modeling has been conducted using extensive validation techniques and statistical test to investigate the structural influence on the degradability of organic micropollutants with ozonation. In contrast to most of the other studies, the underlying dataset - rate constants for 92 organic molecules - were obtained under standardized conditions with defined experimental parameters. QSPR modeling was executed using a combination of the software PaDEL for descriptor calculation and QSARINS for the modeling process respecting all five OECD-requirements for applicable QSAR/QSPR-models. The final model was selected using a multi-criteria decision-making tool to evaluate the model quality based on all calculated statistical quality parameters. The model included 10 selected descriptors and fingerprints and showed good regression abilities, predictive power, and stability (R² = 0.8221, CCC_tr = 0.9024, Q²_loo = 0.7436, R²_ext = 0.8420, Q²_F1 = 0.8104). The applicability domain of the QSPR model was defined and an interpretation of selected model descriptors has been connected to previous experimental studies. A significant influence of the interpretable descriptors was put into experimental context and compared with previous studies and models. For example, the molar refractivity as a measure of size and polarizability of a molecule and the occurrence of important substructures such as a formamide group seem to decrease the removal rate constant. The contribution of lone electrons entering into resonance as well as the occurrence of fused rings were identified as influences for the increase of the degradability of micropollutants by ozonation.

Predicting protection capacities of pyrimidine-based corrosion inhibitors for mild steel/HCl interface using linear and nonlinear QSPR models

Pyrimidine compounds have proven to be effective and efficient additives capable of protecting mild steel in acidic media. This class of organic compounds often functions as adsorption-type inhibitors of corrosion by forming a protective layer on the metallic substrate. The present study reports a computational study of forty pyrimidine compounds that have been investigated as sustainable inhibitors of mild steel corrosion in molar HCl solution. Quantitative structure property relationship was conducted using linear (multiple linear regression) and nonlinear (artificial neural network) models. Standardization method was employed in variable selection yielding five top chemical descriptors utilized for model development along with the inhibitor concentration. Multiple linear regression model yielded a fair predictive model. Artificial neural network model developed using k-fold cross-validation method provided a comprehensive insight into the corrosion protection mechanism of studied pyrimidine-based corrosion inhibitors. Using a multilayer perceptron with Levenberg-Marquardt algorithm, the study obtained the optimal model having a MSE of 8.479, RMSE of 2.912, MAD of 1.791, and MAPE of 2.648. The optimal neural network model was further utilized to forecast the protection capacities of nine non-synthesized pyrimidine derivatives. The predicted inhibition efficiencies ranged from 89 to 98%, revealing the significance of the considered chemical descriptors, the predictive capacity of the developed model, and the potency of the theoretical inhibitors.

QSAR, Molecular Docking, MD Simulation and MMGBSA Calculations Approaches to Recognize Concealed Pharmacophoric Features Requisite for the Optimization of ALK Tyrosine Kinase Inhibitors as Anticancer Leads

ALK tyrosine kinase ALK TK is an important target in the development of anticancer drugs. In the present work, we have performed a QSAR analysis on a dataset of 224 molecules in order to quickly predict anticancer activity on query compounds. Double cross validation assigns an upward plunge to the genetic algorithm−multi linear regression (GA-MLR) based on robust univariate and multivariate QSAR models with high statistical performance reflected in various parameters like, fitting parameters; R2 = 0.69−0.87, F = 403.46−292.11, etc., internal validation parameters; Q2LOO = 0.69−0.86, Q2LMO = 0.69−0.86, CCCcv = 0.82−0.93, etc., or external validation parameters Q2F1 = 0.64−0.82, Q2F2 = 0.63−0.82, Q2F3 = 0.65−0.81, R2ext = 0.65−0.83 including RMSEtr < RMSEcv. The present QSAR evaluation successfully identified certain distinct structural features responsible for ALK TK inhibitory potency, such as planar Nitrogen within four bonds from the Nitrogen atom, Fluorine atom within five bonds beside the non-ring Oxygen atom, lipophilic atoms within two bonds from the ring Carbon atoms. Molecular docking, MD simulation, and MMGBSA computation results are in consensus with and complementary to the QSAR evaluations. As a result, the current study assists medicinal chemists in prioritizing compounds for experimental detection of anticancer activity, as well as their optimization towards more potent ALK tyrosine kinase inhibitor.

Predicting the Surface Tension of Deep Eutectic Solvents: A Step Forward in the Use of Greener Solvents

Deep eutectic solvents (DES) are an important class of green solvents that have been developed as an alternative to toxic solvents. However, the large-scale industrial application of DESs requires fine-tuning their physicochemical properties. Among others, surface tension is one of such properties that have to be considered while designing novel DESs. In this work, we present the results of a detailed evaluation of Quantitative Structure-Property Relationships (QSPR) modeling efforts designed to predict the surface tension of DESs, following the Organization for Economic Co-operation and Development (OECD) guidelines. The data set used comprises a large number of structurally diverse binary DESs and the models were built systematically through rigorous validation methods, including ‘mixtures-out’- and ‘compounds-out’-based data splitting. The most predictive individual QSPR model found is shown to be statistically robust, besides providing valuable information about the structural and physicochemical features responsible for the surface tension of DESs. Furthermore, the intelligent consensus prediction strategy applied to multiple predictive models led to consensus models with similar statistical robustness to the individual QSPR model. The benefits of the present work stand out also from its reproducibility since it relies on fully specified computational procedures and on publicly available tools. Finally, our results not only guide the future design and screening of novel DESs with a desirable surface tension but also lays out strategies for efficiently setting up silico-based models for binary mixtures.

Prediction of Combined Sorbent and Catalyst Materials for SE-SMR, Using QSPR and Multitask Learning

The process of sorption enhanced steam methane reforming (SE-SMR) is an emerging technology for the production of low carbon hydrogen. The development of a suitable catalytic material, as well as a CO₂ adsorbent with high capture capacity, has slowed the upscaling of this process to date. In this study, to aid the development of a combined sorbent catalyst material (CSCM) for SE-SMR, a novel approach involving quantitative structure–property relationship analysis (QSPR) has been proposed. Through data-mining, two databases have been developed for the prediction of the last cycle capacity (g_CO2/g_sorbent) and methane conversion (%). Multitask learning (MTL) was applied for the prediction of CSCM properties. Patterns in the data of this study have also yielded further insights; colored scatter plots were able to show certain patterns in the input data, as well as suggestions on how to develop an optimal material. With the results from the actual vs predicted plots collated, raw materials and synthesis conditions were proposed that could lead to the development of a CSCM that has good performance with respect to both the last cycle capacity and the methane conversion.

Perceiving the Concealed and Unreported Pharmacophoric Features of the 5-Hydroxytryptamine Receptor Using Balanced QSAR Analysis

The 5-hydroxytryptamine receptor 6 (5-HT6) has gained attention as a target for developing therapeutics for Alzheimer’s disease, schizophrenia, cognitive dysfunctions, anxiety, and depression, to list a few. In the present analysis, a larger and diverse dataset of 1278 molecules covering a broad chemical and activity space was used to identify visual and concealed structural features associated with binding affinity for 5-HT6. For this, quantitative structure–activity relationships (QSAR) and molecular docking analyses were executed. This led to the development of a statistically robust QSAR model with a balance of excellent predictivity (R²_tr = 0.78, R²_ex = 0.77), the identification of unreported aspects of known features, and also novel mechanistic interpretations. Molecular docking and QSAR provided similar as well as complementary results. The present analysis indicates that the partial charges on ring carbons present within four bonds from a sulfur atom, the occurrence of sp3-hybridized carbon atoms bonded with donor atoms, and a conditional occurrence of lipophilic atoms/groups from nitrogen atoms, which are prominent but unreported pharmacophores that should be considered while optimizing a molecule for 5-HT6. Thus, the present analysis led to identification of some novel unreported structural features that govern the binding affinity of a molecule. The results could be beneficial in optimizing the molecules for 5-HT6.

Ligand-based and structure-based studies to develop predictive models for SARS-CoV-2 main protease inhibitors through the 3d-qsar.com portal

The main protease (Mpro) of SARS-Cov-2 is the essential enzyme for maturation of functional proteins implicated in viral replication and transcription. The peculiarity of its specific cleavage site joint with its high degree of conservation among all coronaviruses promote it as an attractive target to develop broad-spectrum inhibitors, with high selectivity and tolerable safety profile. Herein is reported a combination of three-dimensional quantitative structure-activity relationships (3-D QSAR) and comparative molecular binding energy (COMBINE) analysis to build robust and predictive ligand-based and structure-based statistical models, respectively. Models were trained on experimental binding poses of co-crystallized Mpro-inhibitors and validated on available literature data. By means of deep optimization both models' goodness and robustness reached final statistical values of r2/q2 values of 0.97/0.79 and 0.93/0.79 for the 3-D QSAR and COMBINE approaches respectively, and an overall predictiveness values of 0.68 and 0.57 for the SDEPPRED and AAEP metrics after application to a test set of 60 compounds covered by the training set applicability domain. Despite the different nature (ligand-based and structure-based) of the employed methods, their outcome fully converged. Furthermore, joint ligand- and structure-based structure-activity relationships were found in good agreement with nirmatrelvir chemical features properties, a novel oral Mpro-inhibitor that has recently received U.S. FDA emergency use authorization (EUA) for the oral treatment of mild-to-moderate COVID-19 infected patients. The obtained results will guide future rational design and/or virtual screening campaigns with the aim of discovering new potential anti-coronavirus lead candidates, minimizing both time and financial resources. Moreover, as most of calculation were performed through the well-established web portal 3d-qsar.com the results confirm the portal as a useful tool for drug design.

Integration of Ligand-Based and Structure-Based Methods for the Design of Small-Molecule TLR7 Antagonists

Toll-like receptor 7 (TLR7) is activated in response to the binding of single-stranded RNA. Its over-activation has been implicated in several autoimmune disorders, and thus, it is an established therapeutic target in such circumstances. TLR7 small-molecule antagonists are not yet available for therapeutic use. We conducted a ligand-based drug design of new TLR7 antagonists through a concerted effort encompassing 2D-QSAR, 3D-QSAR, and pharmacophore modelling of 54 reported TLR7 antagonists. The developed 2D-QSAR model depicted an excellent correlation coefficient (R²_training: 0.86 and R²_test: 0.78) between the experimental and estimated activities. The ligand-based drug design approach utilizing the 3D-QSAR model (R²_training: 0.95 and R²_test: 0.84) demonstrated a significant contribution of electrostatic potential and steric fields towards the TLR7 antagonism. This consolidated approach, along with a pharmacophore model with high correlation (R_training: 0.94 and R_test: 0.92), was used to design quinazoline-core-based hTLR7 antagonists. Subsequently, the newly designed molecules were subjected to molecular docking onto the previously proposed binding model and a molecular dynamics study for a better understanding of their binding pattern. The toxicity profiles and drug-likeness characteristics of the designed compounds were evaluated with in silico ADMET predictions. This ligand-based study contributes towards a better understanding of lead optimization and the future development of potent TLR7 antagonists.

Mechanistic Analysis of Chemically Diverse Bromodomain-4 Inhibitors Using Balanced QSAR Analysis and Supported by X-ray Resolved Crystal Structures

Bromodomain-4 (BRD-4) is a key enzyme in post-translational modifications, transcriptional activation, and many other cellular processes. Its inhibitors find their therapeutic usage in cancer, acute heart failure, and inflammation to name a few. In the present study, a dataset of 980 molecules with a significant diversity of structural scaffolds and composition was selected to develop a balanced QSAR model possessing high predictive capability and mechanistic interpretation. The model was built as per the OECD (Organisation for Economic Co-operation and Development) guidelines and fulfills the endorsed threshold values for different validation parameters (R²_tr = 0.76, Q²_LMO = 0.76, and R²_ex = 0.76). The present QSAR analysis identified that anti-BRD-4 activity is associated with structural characters such as the presence of saturated carbocyclic rings, the occurrence of carbon atoms near the center of mass of a molecule, and a specific combination of planer or aromatic nitrogen with ring carbon, donor, and acceptor atoms. The outcomes of the present analysis are also supported by X-ray-resolved crystal structures of compounds with BRD-4. Thus, the QSAR model effectively captured salient as well as unreported hidden pharmacophoric features. Therefore, the present study successfully identified valuable novel pharmacophoric features, which could be beneficial for the future optimization of lead/hit compounds for anti-BRD-4 activity.

Prediction Models for Brain Distribution of Drugs Based on Biomimetic Chromatographic Data

The development of high-throughput approaches for the valid estimation of brain disposition is of great importance in the early drug screening of drug candidates. However, the complexity of brain tissue, which is protected by a unique vasculature formation called the blood−brain barrier (BBB), complicates the development of robust in silico models. In addition, most computational approaches focus only on brain permeability data without considering the crucial factors of plasma and tissue binding. In the present study, we combined experimental data obtained by HPLC using three biomimetic columns, i.e., immobilized artificial membranes, human serum albumin, and α1-acid glycoprotein, with molecular descriptors to model brain disposition of drugs. Kp,uu,brain, as the ratio between the unbound drug concentration in the brain interstitial fluid to the corresponding plasma concentration, brain permeability, the unbound fraction in the brain, and the brain unbound volume of distribution, was collected from literature. Given the complexity of the investigated biological processes, the extracted models displayed high statistical quality (R2 > 0.6), while in the case of the brain fraction unbound, the models showed excellent performance (R2 > 0.9). All models were thoroughly validated, and their applicability domain was estimated. Our approach highlighted the importance of phospholipid, as well as tissue and protein, binding in balance with BBB permeability in brain disposition and suggests biomimetic chromatography as a rapid and simple technique to construct models with experimental evidence for the early evaluation of CNS drug candidates.

Quantitative Retention (Structure)–Activity Relationships in Predicting the Pharmaceutical and Toxic Properties of Potential Pesticides

The micellar liquid chromatography technique and quantitative retention (structure)–activity relationships method were used to predict properties of carbamic and phenoxyacetic acids derivatives, newly synthesized in our laboratory and considered as potential pesticides. Important properties of the test substances characterizing their potential significance as pesticides as well as threats to humans were considered: the volume of distribution, the unbonded fractions, the blood–brain distribution, the rate of skin and cell permeation, the dermal absorption, the binding to human serum albumin, partitioning between water and plants’ cuticles, and the lethal dose. Pharmacokinetic and toxicity parameters were predicted as functions of the solutes’ lipophilicities and the number of hydrogen bond donors, the number of hydrogen bond acceptors, and the number of rotatable bonds. The equations that were derived were evaluated statistically and cross-validated. Important features of the molecular structure influencing the properties of the tested substances were indicated. The QSAR models that were developed had high predictive ability and high reliability in modeling the properties of the molecules that were tested. The investigations highlighted the applicability of combined chromatographic technique and QS(R)ARs in modeling the important properties of potential pesticides and reducing unethical animal testing.

A novel model for prediction of stability constants of the thiosemicarbazone ligands with different types of toxic heavy metal ions using structural parameters and multivariate linear regression method

A novel model is presented for reliable estimation of the stability constants of the thiosemicarbazone ligands with different types of toxic heavy metal ions (log β₁₁) in an aqueous solution, which has wide usage in environmental safety and ecotoxicology applications. The biggest reported data of log β₁₁ for 120 metalthiosemicarbazone complexes are used for deriving and testing the novel model. In contrast to available methods where they need the two-dimensional (2D) and three-dimensional (3D) complex molecular descriptors as well as expert users and computer codes, the novel correlation uses four additive and two non-additive structural parameters of thiosemicarbazone ligands. The calculated results of the novel correlation are compared with the outputs of the genetic algorithm with multivariate linear regression method (GA-MLR) as one of the best existing methods, which requires seven complex descriptors. The estimated results for 78 of training as well as 42 of two different test sets were established by external and internal validations. The values of statistical parameters comprising average deviation, average absolute deviation, average absolute relative deviation, absolute maximum deviation, and the coefficient of determination for 73 data of training set of New model/GA-MLR are 0.04/ - 0.25, 1.06/1.31, 14.4/18.7, 3.18/7.92, and 0.830/0.652, respectively. Thus, the predicted results of the new model are worthy as compared to the complex GA-MLR model. Moreover, assessments of various statistical parameters confirm that the new model provides great reliability, goodness-of-fit, accuracy, and precision.

QSAR analysis on a large and diverse set of potent phosphoinositide 3-kinase gamma (PI3Kγ) inhibitors using MLR and ANN methods

Phosphorylation of PI3Kγ as a member of lipid kinases-enzymes, plays a crucial role in regulating immune cells through the generation of intracellular signals. Deregulation of this pathway is involved in several tumors. In this research, diverse sets of potent and selective isoform-specific PI3Kγ inhibitors whose drug-likeness was confirmed based on Lipinski’s rule of five were used in the modeling process. Genetic algorithm (GA)-based multivariate analysis was employed on the half-maximal inhibitory concentration (IC₅₀) of them. In this way, multiple linear regression (MLR) and artificial neural network (ANN) algorithm, were used to QSAR models construction on 245 compounds with a wide range of pIC₅₀ (5.23–9.32). The stability and robustness of the models have been evaluated by external and internal validation methods (R² 0.623–0.642, RMSE 0.464–0.473, F 40.114, Q²_LOO 0.600, and R²_y-random 0.011). External verification using a wide variety of structures out of the training and test sets show that ANN is superior to MLR. The descriptors entered into the model are in good agreement with the X-ray structures of target-ligand complexes; so the model is interpretable. Finally, Williams plot-based analysis was applied to simultaneously compare the inhibitory activity and structural similarity of training, test and validation sets.

Identification of Potent Aldose Reductase Inhibitors as Antidiabetic (Anti-hyperglycemic) agents using QSAR Based Virtual Screening, Molecular Docking, MD Simulation and MMGBSA Approaches

The aldose reductase (AR) enzyme is an important target enzyme in the development of therapeutics against hyperglycaemia induced health complications such as retinopathy, etc. In the present study, a quantitative structure activity relationship (QSAR) evaluation of a dataset of 226 reported AR inhibitor (ARi) molecules is performed using a genetic algorithm – multi linear regression (GA-MLR) technique. Multi-criteria decision making (MCDM) analysis furnished two five variables based QSAR models with acceptably high performance reflected in various statistical parameters such as, R² = 0.79–0.80, Q²_LOO = 0.78–0.79, Q²_LMO = 0.78–0.79. The QSAR model analysis revealed some of the molecular features that play crucial role in deciding inhibitory potency of the molecule against AR such as; hydrophobic Nitrogen within 2 Å of the center of mass of the molecule, non-ring Carbon separated by three and four bonds from hydrogen bond donor atoms, number of sp2 hybridized Oxygen separated by four bonds from sp2 hybridized Carbon atoms, etc. 14 in silico generated hits, using a compound 18 (a most potent ARi from present dataset with pIC₅₀ = 8.04 M) as a template, on QSAR based virtual screening (QSAR-VS) furnished a scaffold 5 with better ARi activity (pIC₅₀ = 8.05 M) than template compound 18. Furthermore, molecular docking of compound 18 (Docking Score = –7.91 kcal/mol) and scaffold 5 (Docking Score = –8.08 kcal/mol) against AR, divulged that they both occupy the specific pocket(s) in AR receptor binding sites through hydrogen bonding and hydrophobic interactions. Molecular dynamic simulation (MDS) and MMGBSA studies right back the docking results by revealing the fact that binding site residues interact with scaffold 5 and compound 18 to produce a stable complex similar to co-crystallized ligand's conformation. The QSAR analysis, molecular docking, and MDS results are all in agreement and complementary. QSAR-VS successfully identified a more potent novel ARi and can be used in the development of therapeutic agents to treat diabetes.