QSAR modeling: where have you been? Where are you going to?

Explainable artificial intelligence in the design of selective carbonic anhydrase I-II inhibitors via molecular fingerprinting

Inhibiting the enzymes carbonic anhydrase I (CA I) and carbonic anhydrase II (CA II) presents a potential avenue for addressing nervous system ailments such as glaucoma and Alzheimer's disease. Our study explored harnessing explainable artificial intelligence (XAI) to unveil the molecular traits inherent in CA I and CA II inhibitors. The PubChem molecular fingerprints of these inhibitors, sourced from the ChEMBL database, were subjected to detailed XAI analysis. The study encompassed training 10 regression models using IC50 values, and their efficacy was gauged using metrics including R2, RMSE, and time taken. The Decision Tree Regressor algorithm emerged as the optimal performer (R2: 0.93, RMSE: 0.43, time-taken: 0.07). Furthermore, the PFI method unveiled key molecular features for CA I inhibitors, notably PubChemFP432 (C(O)N) and PubChemFP6978 (C(O)O). The SHAP analysis highlighted the significance of attributes like PubChemFP539 (C(O)NCC), PubChemFP601 (C(O)OCC), and PubChemFP432 (C(O)N) in CA I inhibitiotable n. Likewise, features for CA II inhibitors encompassed PubChemFP528(C(O)OCCN), PubChemFP791 (C(O)OCCC), PubChemFP696 (C(O)OCCCC), PubChemFP335 (C(O)NCCN), PubChemFP580 (C(O)NCCCN), and PubChemFP180 (C(O)NCCC), identified through SHAP analysis. The sulfonamide group (S), aromatic ring (A), and hydrogen bonding group (H) exert a substantial impact on CA I and CA II enzyme activities and IC50 values through the XAI approach. These insights into the CA I and CA II inhibitors are poised to guide future drug discovery efforts, serving as a beacon for innovative therapeutic interventions.

Topological regression as an interpretable and efficient tool for quantitative structure-activity relationship modeling

Quantitative structure-activity relationship (QSAR) modeling is a powerful tool for drug discovery, yet the lack of interpretability of commonly used QSAR models hinders their application in molecular design. We propose a similarity-based regression framework, topological regression (TR), that offers a statistically grounded, computationally fast, and interpretable technique to predict drug responses. We compare the predictive performance of TR on 530 ChEMBL human target activity datasets against the predictive performance of deep-learning-based QSAR models. Our results suggest that our sparse TR model can achieve equal, if not better, performance than the deep learning-based QSAR models and provide better intuitive interpretation by extracting an approximate isometry between the chemical space of the drugs and their activity space.

Chemical hazard assessment toward safer electrolytes for lithium-ion batteries

Commercialization of rechargeable lithium-ion (Li-ion) batteries has revolutionized the design of portable electronic devices and is facilitating the current transition to electric vehicles. The technological specifications of Li-ion batteries continue to evolve through the introduction of various high-risk liquid electrolyte chemicals, yet critical evaluation of the physical, environmental, and human health hazards of these substances is lacking. Using the GreenScreen for Safer Chemicals approach, we conducted a chemical hazard assessment (CHA) of 103 electrolyte chemicals categorized into seven chemical groups: salts, carbonates, esters, ethers, sulfoxides-sulfites-sulfones, overcharge protection additives, and flame-retardant additives. To minimize data gaps, we focused on six toxicity and hazard data sources, including three empirical and three nonempirical predictive data sources. Furthermore, we investigated the structural similarities among selected electrolyte chemicals using the ChemMine tool and the simplified molecular input line entry system inputs from PubChem to evaluate whether chemicals with similar structures exhibit similar toxicity. The results demonstrate that salts, overcharge protection additives, and flame-retardant additives contain the most toxic components in the electrolyte solutions. Furthermore, carbonates, esters, and ethers account for most flammability hazards in Li-ion batteries. This study supports the complementary use of quantitative structure-activity relationship models to minimize data gaps and inconsistencies in CHA. Integr Environ Assess Manag 2024;00:1-14. © 2024 The Author(s). Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Impact of 2-hydroxypropyl-β-cyclodextrin inclusion complex formation on dopamine receptor-ligand interaction - A case study

The octanol-water distribution coefficient (logP), used as a measure of lipophilicity, plays a major role in the drug design and discovery processes. While average logP values remain unchanged in approved oral drugs since 1983, current medicinal chemistry trends towards increasingly lipophilic compounds that require adapted analytical workflows and drug delivery systems. Solubility enhancers like cyclodextrins (CDs), especially 2-hydroxypropyl-β-CD (2-HP-β-CD), have been studied in vitro and in vivo investigating their ADMET (adsorption, distribution, metabolism, excretion and toxicity)-related properties. However, data is scarce regarding the applicability of CD inclusion complexes (ICs) in vitro compared to pure compounds. In this study, dopamine receptor (DR) ligands were used as a case study, utilizing a combined in silico/in vitro workflow. Media-dependent solubility and IC stoichiometry were investigated using HPLC. NMR was used to observe IC formation-caused chemical shift deviations while in silico approaches utilizing basin hopping global minimization were used to propose putative IC binding modes. A cell-based in vitro homogeneously time-resolved fluorescence (HTRF) assay was used to quantify ligand binding affinity at the DR subtype 2 (D2R). While all ligands showed increased solubility using 2-HP-β-CD, they differed regarding IC stoichiometry and receptor binding affinity. This case study shows that IC-formation was ligand-dependent and sometimes altering in vitro binding. Therefore, IC complex formation can't be recommended as a general means of improving compound solubility for in vitro studies as they may alter ligand binding.

Development of a novel representation of drug 3D structures and enhancement of the TSR-based method for probing drug and target interactions

Understanding the mechanisms underlying interactions between drugs and target proteins is critical for drug discovery. In our earlier studies, we introduced the Triangular Spatial Relationship (TSR)-based algorithm, which enables the representation of a protein's 3D structure as a vector of integers (TSR keys). These TSR keys correspond to substructures of the 3D structure of a protein and are computed based on the triangles constructed by all possible triples of Cα atoms within the protein. In this study, we report on a new TSR-based algorithm for probing drug and target interactions. Specifically, we have extended the previous algorithm in three novel directions: TSR keys for representing the 3D structure of a drug or a ligand, cross TSR keys between drugs and their targets and intra-residual TSR keys for phosphorylated amino acids. The outcomes illustrate the key contributions as follows: (i) The TSR-based method, which uses the TSR keys as features, is unique in its capability to interpret hierarchical relationships of drugs as well as drug - target complexes using common and specific TSR keys. (ii) The method can distinguish not only the binding sites from the rest of the protein structures, but also the binding sites of primary targets from those of off-targets. (iii) The method has the potential to correlate the 3D structures of drugs with their functions. (iv) Representation of 3D structures by TSR keys has its unique advantage in terms of ease of making searching for similar substructures across structure datasets easier. In summary, this study presents a novel computational methodology, with significant advantages, for providing insights into the mechanism underlying drug and target interactions.

In silico prediction of ocular toxicity of compounds using explainable machine learning and deep learning approaches

The accurate identification of chemicals with ocular toxicity is of paramount importance in health hazard assessment. In contemporary chemical toxicology, there is a growing emphasis on refining, reducing, and replacing animal testing in safety evaluations. Therefore, the development of robust computational tools is crucial for regulatory applications. The performance of predictive models is heavily reliant on the quality and quantity of data. In this investigation, we amalgamated the most extensive dataset (4901 compounds) sourced from governmental GHS-compliant databases and literature to develop binary classification models of chemical ocular toxicity. We employed 12 molecular representations in conjunction with six machine learning algorithms and two deep learning algorithms to create a series of binary classification models. The findings indicated that the deep learning method GCN outperformed the machine learning models in cross-validation, achieving an impressive AUC of 0.915. However, the top-performing machine learning model (RF-Descriptor) demonstrated excellent performance with an AUC of 0.869 on the test set and was therefore selected as the best model. To enhance model interpretability, we conducted the SHAP method and attention weights analysis. The two approaches offered visual depictions of the relevance of key descriptors and substructures in predicting ocular toxicity of chemicals. Thus, we successfully struck a delicate balance between data quality and model interpretability, rendering our model valuable for predicting and comprehending potential ocular-toxic compounds in the early stages of drug discovery.

Solubilization of drugs using beta-cyclodextrin: Experimental data and modeling

Many drug candidates fail to complete the entire drug development process because of poor physicochemical properties. Solubility is an important physicochemical property which plays a vital role in various stages of drug discovery and development. Several methods have been proposed to enhance the solubility of drugs, and complex formation with cyclodextrins is among them. Beta-cyclodextrin (βCD) is a common excipient for solubilization of drugs. The aim of this study is to develop the mechanistic QSPR models to predict the solubility enhancement of a drug in the presence of βCD. In this study, the solubility enhancement of some drugs in the presence of 10mM βCD at 25°C was experimentally determined or collected from the literature. Two different models to predict the solubilization by βCD were developed by binary logistic regression using structural properties of drugs with more than 80% accuracy. Polar surface area and excess molar refraction are the main parameters for estimating solubilization by βCD. Moreover, other descriptors related to hydrophobicity and the capability of hydrogen bonding formation of molecules could improve the accuracy of the established models.

Multitarget anti-parasitic activities of isoquinoline alkaloids isolated from Hippeastrum aulicum (Amaryllidaceae)

BACKGROUND

Chagas disease and leishmaniasis affect a significant portion of the Latin American population and still lack efficient treatments. In this context, natural products emerge as promising compounds for developing more effective therapies, aiming to mitigate side effects and drug resistance. Notably, species from the Amaryllidaceae family emerge as potential reservoirs of antiparasitic agents due to the presence of diverse biologically active alkaloids.

PURPOSE

To assess the anti-Trypanosoma cruzi and anti-Leishmania infantum activity of five isolated alkaloids from Hippeastrum aulicum Herb. (Amaryllidaceae) against different life stages of the parasites using in silico and in vitro assays. Furthermore, molecular docking was employed to evaluate the interaction of the most active alkaloids.

METHODS

Five natural isoquinoline alkaloids isolated in suitable quantities for in vitro testing underwent preliminary in silico analysis to predict their potential efficacy against Trypanosoma cruzi (amastigote and trypomastigote forms) and Leishmania infantum (amastigote and promastigote forms). The in vitro antiparasitic activity and mammalian cytotoxicity were investigated with a subsequent comparison of both analysis (in silico and in vitro) findings. Additionally, this study employed the molecular docking technique, utilizing cruzain (T. cruzi) and sterol 14α-demethylase (CYP51, L. infantum) as crucial biological targets for parasite survival, specifically focusing on compounds that exhibited promising activities against both parasites.

RESULTS

Through computational techniques, it was identified that the alkaloids haemanthamine (1) and lycorine (8) were the most active against T. cruzi (amastigote and trypomastigote) and L. infantum (amastigote and promastigote), while also revealing unprecedented activity of alkaloid 7‑methoxy-O-methyllycorenine (6). The in vitro analysis confirmed the in silico tests, in which compound 1 presented the best activities against the promastigote and amastigote forms of L. infantum with half-maximal inhibitory concentration (IC50) 0.6 µM and 1.78 µM, respectively. Compound 8 exhibited significant activity against the amastigote form of T. cruzi (IC50 7.70 µM), and compound 6 demonstrated activity against the trypomastigote forms of T. cruzi and amastigote of L. infantum, with IC50 values of 89.55 and 86.12 µM, respectively. Molecular docking analyses indicated that alkaloids 1 and 8 exhibited superior interaction energies compared to the inhibitors.

CONCLUSION

The hitherto unreported potential of compound 6 against T. cruzi trypomastigotes and L. infantum amastigotes is now brought to the forefront. Furthermore, the acquired dataset signifies that the isolated alkaloids 1 and 8 from H. aulicum might serve as prototypes for subsequent structural refinements aimed at the exploration of novel leads against both T. cruzi and L. infantum parasites.

Comprehensive evaluation framework for intervention on health effects of ambient temperature

Despite the existence of many interventions to mitigate or adapt to the health effects of climate change, their effectiveness remains unclear. Here, we introduce the Comprehensive Evaluation Framework for Intervention on Health Effects of Ambient Temperature to evaluate study designs and effects of intervention studies. The framework comprises three types of interventions: proactive, indirect, and direct, and four categories of indicators: classification, methods, scope, and effects. We trialed the framework by an evaluation of existing intervention studies. The evaluation revealed that each intervention has its own applicable characteristics in terms of effectiveness, feasibility, and generalizability scores. We expanded the framework's potential by offering a list of intervention recommendations in different scenarios. Future applications are then explored to establish models of the relationship between study designs and intervention effects, facilitating effective interventions to address the health effects of ambient temperature under climate change.

POxload: Machine Learning Estimates Drug Loadings of Polymeric Micelles

Block copolymers, composed of poly(2-oxazoline)s and poly(2-oxazine)s, can serve as drug delivery systems; they form micelles that carry poorly water-soluble drugs. Many recent studies have investigated the effects of structural changes of the polymer and the hydrophobic cargo on drug loading. In this work, we combine these data to establish an extended formulation database. Different molecular properties and fingerprints are tested for their applicability to serve as formulation-specific mixture descriptors. A variety of classification and regression models are built for different descriptor subsets and thresholds of loading efficiency and loading capacity, with the best models achieving overall good statistics for both cross- and external validation (balanced accuracies of 0.8). Subsequently, important features are dissected for interpretation, and the DrugBank is screened for potential therapeutic use cases where these polymers could be used to develop novel formulations of hydrophobic drugs. The most promising models are provided as an open-source software tool for other researchers to test the applicability of these delivery systems for potential new drug candidates.

Py-CoMFA, docking, and molecular dynamics simulations of Leishmania (L.) amazonensis arginase inhibitors

Leishmaniasis is a disease caused by a protozoan of the genus Leishmania, affecting millions of people, mainly in tropical countries, due to poor social conditions and low economic development. First-line chemotherapeutic agents involve highly toxic pentavalent antimonials, while treatment failure is mainly due to the emergence of drug-resistant strains. Leishmania arginase (ARG) enzyme is vital in pathogenicity and contributes to a higher infection rate, thus representing a potential drug target. This study helps in designing ARG inhibitors for the treatment of leishmaniasis. Py-CoMFA (3D-QSAR) models were constructed using 34 inhibitors from different chemical classes against ARG from L. (L.) amazonensis (LaARG). The 3D-QSAR predictions showed an excellent correlation between experimental and calculated pIC50 values. The molecular docking study identified the favorable hydrophobicity contribution of phenyl and cyclohexyl groups as substituents in the enzyme allosteric site. Molecular dynamics simulations of selected protein-ligand complexes were conducted to understand derivatives' interaction modes and affinity in both active and allosteric sites. Two cinnamide compounds, 7g and 7k, were identified, with similar structures to the reference 4h allosteric site inhibitor. These compounds can guide the development of more effective arginase inhibitors as potential antileishmanial drugs.

Interpreting Neural Network Models for Toxicity Prediction by Extracting Learned Chemical Features

Neural network models have become a popular machine-learning technique for the toxicity prediction of chemicals. However, due to their complex structure, it is difficult to understand predictions made by these models which limits confidence. Current techniques to tackle this problem such as SHAP or integrated gradients provide insights by attributing importance to the input features of individual compounds. While these methods have produced promising results in some cases, they do not shed light on how representations of compounds are transformed in hidden layers, which constitute how neural networks learn. We present a novel technique to interpret neural networks which identifies chemical substructures in training data found to be responsible for the activation of hidden neurons. For individual test compounds, the importance of hidden neurons is determined, and the associated substructures are leveraged to explain the model prediction. Using structural alerts for mutagenicity from the Derek Nexus expert system as ground truth, we demonstrate the validity of the approach and show that model explanations are competitive with and complementary to explanations obtained from an established feature attribution method.

FEP Augmentation as a Means to Solve Data Paucity Problems for Machine Learning in Chemical Biology

In the realm of medicinal chemistry, the primary objective is to swiftly optimize a multitude of chemical properties of a set of compounds to yield a clinical candidate poised for clinical trials. In recent years, two computational techniques, machine learning (ML) and physics-based methods, have evolved substantially and are now frequently incorporated into the medicinal chemist's toolbox to enhance the efficiency of both hit optimization and candidate design. Both computational methods come with their own set of limitations, and they are often used independently of each other. ML's capability to screen extensive compound libraries expediently is tempered by its reliance on quality data, which can be scarce especially during early-stage optimization. Contrarily, physics-based approaches like free energy perturbation (FEP) are frequently constrained by low throughput and high cost by comparison; however, physics-based methods are capable of making highly accurate binding affinity predictions. In this study, we harnessed the strength of FEP to overcome data paucity in ML by generating virtual activity data sets which then inform the training of algorithms. Here, we show that ML algorithms trained with an FEP-augmented data set could achieve comparable predictive accuracy to data sets trained on experimental data from biological assays. Throughout the paper, we emphasize key mechanistic considerations that must be taken into account when aiming to augment data sets and lay the groundwork for successful implementation. Ultimately, the study advocates for the synergy of physics-based methods and ML to expedite the lead optimization process. We believe that the physics-based augmentation of ML will significantly benefit drug discovery, as these techniques continue to evolve.

AutoDock-SS: AutoDock for Multiconformational Ligand-Based Virtual Screening

Ligand-based virtual screening (LBVS) can be pivotal for identifying potential drug leads, especially when the target protein's structure is unknown. However, current LBVS methods are limited in their ability to consider the ligand conformational flexibility. This study presents AutoDock-SS (Similarity Searching), which adapts protein-ligand docking for use in LBVS. AutoDock-SS integrates novel ligand-based grid maps and AutoDock-GPU into a novel three-dimensional LBVS workflow. Unlike other approaches based on pregenerated conformer libraries, AutoDock-SS's built-in conformational search optimizes conformations dynamically based on the reference ligand, thus providing a more accurate representation of relevant ligand conformations. AutoDock-SS supports two modes: single and multiple ligand queries, allowing for the seamless consideration of multiple reference ligands. When tested on the Directory of Useful Decoys─Enhanced (DUD-E) data set, AutoDock-SS surpassed alternative 3D LBVS methods, achieving a mean AUROC of 0.775 and an EF1% of 25.72 in single-reference mode. The multireference mode, evaluated on the augmented DUD-E+ data set, demonstrated superior accuracy with a mean AUROC of 0.843 and an EF1% of 34.59. This enhanced performance underscores AutoDock-SS's ability to treat compounds as conformationally flexible while considering the ligand's shape, pharmacophore, and electrostatic potential, expanding the potential of LBVS methods.

Prediction of the First-Pass Metabolism of a Drug After Oral Intake Based on Structural Parameters and Physicochemical Properties

BACKGROUND AND OBJECTIVE

The oral first-pass metabolism is a crucial factor that plays a key role in a drug's pharmacokinetic profile. Prediction of the oral first-pass metabolism based on chemical structural parameters can be useful in the drug-design process. Developing an orally administered drug with an acceptable pharmacokinetic profile is necessary to reduce the cost and time associated with evaluating the extent of the first-pass metabolism of a candidate compound in preclinical studies. The aim of this study is to estimate the first-pass metabolism of an orally administered drug.

METHODS

A set of compounds with reported first-pass metabolism data were collected. Moreover, human intestinal absorption percentage and oral bioavailability data were extracted from the literature to propose a classification system that split the drugs up based on their first-pass metabolism extents. Various structural parameters were calculated for each compound. The relations of the structural and physicochemical values of each compound to the class the compound belongs to were obtained using logistic regression.

RESULTS

Initial analysis showed that compounds with logD7.4 > 1 or a rugosity factor of > 1.5 are more likely to have high first-pass metabolism. Four different models that can predict the oral first-pass metabolism with acceptable error were introduced. The overall accuracies of the models were in the range of 72% (for models with simple descriptors) to 78% (for models with complex descriptors). Although the models with simple descriptors have lower accuracies compared to complex models, they are more interpretable and easier for researchers to utilize.

CONCLUSION

A novel classification of drugs based on the extent of the oral first-pass metabolism was introduced, and mechanistic models were developed to assign candidate compounds to the appropriate proposed classes.

A Chemical Structure and Machine Learning Approach to Assess the Potential Bioactivity of Endogenous Metabolites and Their Association with Early Childhood Systemic Inflammation

Metabolomics has gained much attention due to its potential to reveal molecular disease mechanisms and present viable biomarkers. This work uses a panel of untargeted serum metabolomes from 602 children from the COPSAC2010 mother-child cohort. The annotated part of the metabolome consists of 517 chemical compounds curated using automated procedures. We created a filtering method for the quantified metabolites using predicted quantitative structure-bioactivity relationships for the Tox21 database on nuclear receptors and stress response in cell lines. The metabolites measured in the children's serums are predicted to affect specific targeted models, known for their significance in inflammation, immune function, and health outcomes. The targets from Tox21 have been used as targets with quantitative structure-activity relationships (QSARs). They were trained for ~7000 structures, saved as models, and then applied to the annotated metabolites to predict their potential bioactivities. The models were selected based on strict accuracy criteria surpassing random effects. After application, 52 metabolites showed potential bioactivity based on structural similarity with known active compounds from the Tox21 set. The filtered compounds were subsequently used and weighted by their bioactive potential to show an association with early childhood hs-CRP levels at six months in a linear model supporting a physiological adverse effect on systemic low-grade inflammation.

QuanDB: a quantum chemical property database towards enhancing 3D molecular representation learning

Previous studies have shown that the three-dimensional (3D) geometric and electronic structure of molecules play a crucial role in determining their key properties and intermolecular interactions. Therefore, it is necessary to establish a quantum chemical (QC) property database containing the most stable 3D geometric conformations and electronic structures of molecules. In this study, a high-quality QC property database, called QuanDB, was developed, which included structurally diverse molecular entities and featured a user-friendly interface. Currently, QuanDB contains 154,610 compounds sourced from public databases and scientific literature, with 10,125 scaffolds. The elemental composition comprises nine elements: H, C, O, N, P, S, F, Cl, and Br. For each molecule, QuanDB provides 53 global and 5 local QC properties and the most stable 3D conformation. These properties are divided into three categories: geometric structure, electronic structure, and thermodynamics. Geometric structure optimization and single point energy calculation at the theoretical level of B3LYP-D3(BJ)/6-311G(d)/SMD/water and B3LYP-D3(BJ)/def2-TZVP/SMD/water, respectively, were applied to ensure highly accurate calculations of QC properties, with the computational cost exceeding 107 core-hours. QuanDB provides high-value geometric and electronic structure information for use in molecular representation models, which are critical for machine-learning-based molecular design, thereby contributing to a comprehensive description of the chemical compound space. As a new high-quality dataset for QC properties, QuanDB is expected to become a benchmark tool for the training and optimization of machine learning models, thus further advancing the development of novel drugs and materials. QuanDB is freely available, without registration, at https://quandb.cmdrg.com/ .

Concurrent Optimizations of Efficacy and Blood-Brain Barrier Permeability in New Macrocyclic LRRK2 Inhibitors for Potential Parkinson's Disease Therapeutics

The elevated activity of leucine-rich repeat kinase 2 (LRRK2) is implicated in the pathogenesis of Parkinson's disease (PD). The quest for effective LRRK2 inhibitors has been impeded by the formidable challenge of crossing the blood-brain barrier (BBB). We leveraged structure-based de novo design and developed robust three-dimensional quantitative structure-activity relationship (3D-QSAR) models to predict BBB permeability, enhancing the likelihood of the inhibitor's brain accessibility. Our strategy involved the synthesis of macrocyclic molecules by linking the two terminal nitrogen atoms of HG-10-102-01 with an alkyl chain ranging from 2 to 4 units, laying the groundwork for innovative LRRK2 inhibitor designs. Through meticulous computational and synthetic optimization of both biochemical efficacy and BBB permeability, 9 out of 14 synthesized candidates demonstrated potent low-nanomolar inhibition and significant BBB penetration. Further assessments of in vitro and in vivo effectiveness, coupled with pharmacological profiling, highlighted 8 as the promising new lead compound for PD therapeutics.

LGGA-MPP: Local Geometry-Guided Graph Attention for Molecular Property Prediction

Molecular property prediction is a fundamental task of drug discovery. With the rapid development of deep learning, computational approaches for predicting molecular properties are experiencing increasing popularity. However, these existing methods often ignore the 3D information on molecules, which is critical in molecular representation learning. In the past few years, several self-supervised learning (SSL) approaches have been proposed to exploit the geometric information by using pre-training on 3D molecular graphs and fine-tuning on 2D molecular graphs. Most of these approaches are based on the global geometry of molecules, and there is still a challenge in capturing the local structure and local interpretability. To this end, we propose local geometry-guided graph attention (LGGA), which integrates local geometry into the attention mechanism and message-passing of graph neural networks (GNNs). LGGA introduces a novel method to model molecules, enhancing the model's ability to capture intricate local structural details. Experiments on various data sets demonstrate that the integration of local geometry has a significant impact on the improved results, and our model outperforms the state-of-the-art methods for molecular property prediction, establishing its potential as a promising tool in drug discovery and related fields.

Predicting Elimination of Small-Molecule Drug Half-Life in Pharmacokinetics Using Ensemble and Consensus Machine Learning Methods

Half-life is a significant pharmacokinetic parameter included in the excretion phase of absorption, distribution, metabolism, and excretion. It is one of the key factors for the successful marketing of drug candidates. Therefore, predicting half-life is of great significance in drug design. In this study, we employed eXtreme Gradient Boosting (XGboost), randomForest (RF), gradient boosting machine (GBM), and supporting vector machine (SVM) to build quantitative structure-activity relationship (QSAR) models on 3512 compounds and evaluated model performance by using root-mean-square error (RMSE), R2, and mean absolute error (MAE) metrics and interpreted features by SHapley Additive exPlanation (SHAP). Furthermore, we developed consensus models through integrating four individual models and validated their performance using a Y-randomization test and applicability domain analysis. Finally, matched molecular pair analysis was used to extract the transformation rules. Our results revealed that XGboost outperformed other individual models (RMSE = 0.176, R2 = 0.845, MAE = 0.141). The consensus model integrating all four models continued to enhance prediction performance (RMSE = 0.172, R2 = 0.856, MAE = 0.138). We evaluated the reliability, robustness, and generalization ability via Y-randomization test and applicability domain analysis. Meanwhile, we utilized SHAP to interpret features and employed matched molecular pair analysis to extract chemical transformation rules that provide suggestions for optimizing drug structure. In conclusion, we believe that the consensus model developed in this study serve as a reliable tool to evaluate half-life in drug discovery, and the chemical transformation rules concluded in this study could provide valuable suggestions in drug discovery.

Exploring the potential of 2-arylbenzimidazole scaffolds as novel α-amylase inhibitors: QSAR, molecular docking, simulation and pharmacokinetic studies

Previous studies have shown that 2-arylbenzimidazole derivatives have a strong anti-diabetic effect. To further explore this potential, we develop new analogues of the compound using ligand-based drug design and tested their inhibitory and binding properties through QSAR analyses, molecular docking, dynamic simulations and pharmacokinetic studies. By using quantitative structure activity relationship and ligand-based modification, a highly precise predictive model and design of potent compounds was developed from the derivatives of 2-arylbenzimidazoles. Molecular docking and simulation studies were then conducted to identify the optimal binding poses and pharmacokinetic profiles of the newly generated therapeutic drugs. DFT was employed to optimize the chemical structures of 2-arylbenzimidazole derivatives using B3LYP/6-31G* as the basis set. The model with the highest R2trng set, R2adj, Q2cv, and R2test sets (0.926, 0.912, 0.903, and 0.709 respectively) was chosen to predict the inhibitory activities of the derivatives. Five analogues designed using ligand-based strategy had higher activity than the hit molecule. Additionally, the designed molecules had more favorable MolDock scores than the hit molecule and acarbose and simulation studies confirm on their stability and binding affinities towards the protein. The ADME and druglikeness properties of the analogues indicated that they are safe to consume orally and have a high potential for total clearance. The results of this study showed that the suggested analogues could act as α-amylase inhibitors, which could be used as a basis for the creation of new drugs to treat type 2 diabetes mellitus.

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

Gene Target Prediction of Environmental Chemicals Using Coupled Matrix-Matrix Completion

Human exposure to toxic chemicals presents a huge health burden. Key to understanding chemical toxicity is knowledge of the molecular target(s) of the chemicals. Because a comprehensive safety assessment for all chemicals is infeasible due to limited resources, a robust computational method for discovering targets of environmental exposures is a promising direction for public health research. In this study, we implemented a novel matrix completion algorithm named coupled matrix-matrix completion (CMMC) for predicting direct and indirect exposome-target interactions, which exploits the vast amount of accumulated data regarding chemical exposures and their molecular targets. Our approach achieved an AUC of 0.89 on a benchmark data set generated using data from the Comparative Toxicogenomics Database. Our case studies with bisphenol A and its analogues, PFAS, dioxins, PCBs, and VOCs show that CMMC can be used to accurately predict molecular targets of novel chemicals without any prior bioactivity knowledge. Our results demonstrate the feasibility and promise of computationally predicting environmental chemical-target interactions to efficiently prioritize chemicals in hazard identification and risk assessment.

Recent advances in the area of plant-based anti-cancer drug discovery using computational approaches

Phytocompounds are a well-established source of drug discovery due to their unique chemical and functional diversities. In the area of cancer therapeutics, several phytocompounds have been used till date to design and develop new drugs. One of the desired interests of pharmaceutical companies and researchers globally is that new anti-cancer leads are discovered, for which phytocompounds can be considered a valuable source. Simultaneously, in recent years, the growth of computational approaches like virtual screening (VS), molecular dynamics (MD), pharmacophore modelling, Quantitative structure-activity relationship (QSAR), Absorption Distribution Metabolism Excretion and Toxicity (ADMET), network biology, and machine learning (ML) has gained importance due to their efficiency, reduced time-consuming nature, and cost-effectiveness. Therefore, the present review amalgamates the information on plant-based molecules identified for cancer lead discovery from in silico approaches. The mandate of this review is to discuss studies published in the last 5-6 years that aim to identify the phytomolecules as leads against cancer with the help of traditional computational approaches as well as newer techniques like network pharmacology and ML. This review also lists the databases and webservers available in the public domain for phytocompounds related information that can be harnessed for drug discovery. It is expected that the present review would be useful to pharmacologists, medicinal chemists, molecular biologists, and other researchers involved in the development of natural products (NPs) into clinically effective lead molecules.

Transformers enable accurate prediction of acute and chronic chemical toxicity in aquatic organisms

Environmental hazard assessments are reliant on toxicity data that cover multiple organism groups. Generating experimental toxicity data is, however, resource-intensive and time-consuming. Computational methods are fast and cost-efficient alternatives, but the low accuracy and narrow applicability domains have made their adaptation slow. Here, we present a AI-based model for predicting chemical toxicity. The model uses transformers to capture toxicity-specific features directly from the chemical structures and deep neural networks to predict effect concentrations. The model showed high predictive performance for all tested organism groups-algae, aquatic invertebrates and fish-and has, in comparison to commonly used QSAR methods, a larger applicability domain and a considerably lower error. When the model was trained on data with multiple effect concentrations (EC50/EC10), the performance was further improved. We conclude that deep learning and transformers have the potential to markedly advance computational prediction of chemical toxicity.

Relationships between the Photodegradation Reaction Rate and Structural Properties of Polymer Systems

The development of reusable polymeric materials inspires an attempt to combine renewable biomass with upcycling to form a biorenewable closed system. It has been reported that 2,5-furandicarboxylic acid (FDCA) can be recovered for recycling when incorporated as monomers into photodegradable polymeric systems. Here, we develop a procedure to better understand the photodegradation reactions combining density functional theory (DFT) based time-dependent excited-state molecular dynamics (TDESMD) studies with machine learning-based quantitative structure-activity relationships (QSAR) methodology. This procedure allows for the unveiling of hidden structural features between active orbitals that affect the rate of photodegradation and is coined InfoTDESMD. Findings show that electrotopological features are influential factors affecting the rate of photodegradation in differing environments. Additionally, statistical validations and knowledge-based analysis of descriptors are conducted to further understand the structural features' influence on the rate of photodegradation of polymeric materials.

AI-enhanced chemical paradigm: From molecular graphs to accurate prediction and mechanism

The development of accurate and interpretable models for predicting reaction constants of organic compounds with hydroxyl radicals is vital for advancing quantitative structure-activity relationships (QSAR) in pollutant degradation. Methods like molecular descriptors, molecular fingerprinting, and group contribution methods have limitations, as traditional machine learning struggles to capture all intramolecular information simultaneously. To address this, we established an integrated graph neural network (GNN) with approximately 12 million learnable parameters. GNN represents atoms as nodes and chemical bonds as edges, thus transforming molecules into a graph structures, effectively capturing microscopic properties while depicting atom connectivity in non-Euclidean space. Our datasets comprise 1401 pollutants to develop an integrated GNN model with Bayesian optimization, the model achieves root mean square errors of 0.165, 0.172, and 0.189 on the training, validation, and test datasets, respectively. Furthermore, we assess molecular structure similarity using molecular fingerprint to enhance the model's applicability. Afterwards, we propose a gradient weight mapping method for model explainability, uncovering the key functional groups in chemical reactions in artificial intelligence perspective, which would boost chemistry through artificial intelligence extreme arithmetic power.

Monte Carlo Optimization Method Based QSAR Modeling of Placental Barrier Permeability

PURPOSE

In order to ensure that drug administration is safe during pregnancy, it is crucial to have the possibility to predict the placental permeability of drugs in humans. The experimental method which is most widely used for the said purpose is in vitro human placental perfusion, though the approach is highly expensive and time consuming. Quantitative structure-activity relationship (QSAR) modeling represents a powerful tool for the assessment of the drug placental transfer, and can be successfully employed to be an alternative in in vitro experiments.

METHODS

The conformation-independent QSAR models covered in the present study were developed through the use of the SMILES notation descriptors and local molecular graph invariants. What is more, the Monte Carlo optimization method, was used in the test sets and the training sets as the model developer with three independent molecular splits.

RESULTS

A range of different statistical parameters was used to validate the developed QSAR model, including the standard error of estimation, mean absolute error, root-mean-square error (RMSE), correlation coefficient, cross-validated correlation coefficient, Fisher ratio, MAE-based metrics and the correlation ideality index. Once the mentioned statistical methods were employed, an excellent predictive potential and robustness of the developed QSAR model was demonstrated. In addition, the molecular fragments, which are derived from the SMILES notation descriptors accounting for the decrease or increase in the investigated activity, were revealed.

CONCLUSION

The presented QSAR modeling can be an invaluable tool for the high-throughput screening of the placental permeability of drugs.

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.

Beyond Amyloid: A Machine Learning-Driven Approach Reveals Properties of Potent GSK-3β Inhibitors Targeting Neurofibrillary Tangles

Current treatments for Alzheimer's disease (AD) focus on slowing memory and cognitive decline, but none offer curative outcomes. This study aims to explore and curate the common properties of active, drug-like molecules that modulate glycogen synthase kinase 3β (GSK-3β), a well-documented kinase with increased activity in tau hyperphosphorylation and neurofibrillary tangles-hallmarks of AD pathology. Leveraging quantitative structure-activity relationship (QSAR) data from the PubChem and ChEMBL databases, we employed seven machine learning models: logistic regression (LogR), k-nearest neighbors (KNN), random forest (RF), support vector machine (SVM), extreme gradient boosting (XGB), neural networks (NNs), and ensemble majority voting. Our goal was to correctly predict active and inactive compounds that inhibit GSK-3β activity and identify their key properties. Among the six individual models, the NN demonstrated the highest performance with a 79% AUC-ROC on unbalanced external validation data, while the SVM model was superior in accurately classifying the compounds. The SVM and RF models surpassed NN in terms of Kappa values, and the ensemble majority voting model demonstrated slightly better accuracy to the NN on the external validation data. Feature importance analysis revealed that hydrogen bonds, phenol groups, and specific electronic characteristics are important features of molecular descriptors that positively correlate with active GSK-3β inhibition. Conversely, structural features like imidazole rings, sulfides, and methoxy groups showed a negative correlation. Our study highlights the significance of structural, electronic, and physicochemical descriptors in screening active candidates against GSK-3β. These predictive features could prove useful in therapeutic strategies to understand the important properties of GSK-3β candidate inhibitors that may potentially benefit non-amyloid-based AD treatments targeting neurofibrillary tangles.

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.

A Primer on 2D Descriptors in Selectivity Modeling for Asymmetric Catalysis

Machine learning has permeated all fields of research, including chemistry, and is now an integral part of the design of novel compounds with desired properties. In the field of asymmetric catalysis, the preference still lies with models based on a physical understanding of the catalysis phenomenon and the electronic and steric properties of catalysts. However, such models require quantum chemical calculations and are thus limited by their computational cost. Here, we highlight the recent advances in modeling catalyst selectivity by using the 2D structures of catalysts and substrates. While these have a less explicit mechanistic connection to the modeled property, 2D descriptors, such as topological indices, molecular fingerprints, and fragments, offer the tremendous advantages of low cost and high speed of calculations. This makes them optimal for the in-silico screening of large amounts of data. We provide an overview of common quantitative structure-property relationship workflow, model building and validation techniques, applications of these methodologies in asymmetric catalysis design, and an outlook on improving the understanding of 2D-based models.

Harmonization Risks and Rewards: Nano-QSAR for Agricultural Nanomaterials

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

Rethinking the applicability domain analysis in QSAR models

Notwithstanding the wide adoption of the OECD principles (or best practices) for QSAR modeling, disparities between in silico predictions and experimental results are frequent, suggesting that model predictions are often too optimistic. Of these OECD principles, the applicability domain (AD) estimation has been recognized in several reports in the literature to be one of the most challenging, implying that the actual reliability measures of model predictions are often unreliable. Applying tree-based error analysis workflows on 5 QSAR models reported in the literature and available in the QsarDB repository, i.e., androgen receptor bioactivity (agonists, antagonists, and binders, respectively) and membrane permeability (highest membrane permeability and the intrinsic permeability), we demonstrate that predictions erroneously tagged as reliable (AD prediction errors) overwhelmingly correspond to instances in subspaces (cohorts) with the highest prediction error rates, highlighting the inhomogeneity of the AD space. In this sense, we call for more stringent AD analysis guidelines which require the incorporation of model error analysis schemes, to provide critical insight on the reliability of underlying AD algorithms. Additionally, any selected AD method should be rigorously validated to demonstrate its suitability for the model space over which it is applied. These steps will ultimately contribute to more accurate estimations of the reliability of model predictions. Finally, error analysis may also be useful in "rational" model refinement in that data expansion efforts and model retraining are focused on cohorts with the highest error rates.

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.

Creating New Antifoulants Using the Tools and Tactics of Medicinal Chemistry

The unwanted accumulation of marine micro- and macroorganisms such as algae and barnacles on submerged man-made structures and vessel hulls is a major challenge for any marine operation. Known as biofouling, this problem leads to reduced hydrodynamic efficiency, significantly increased fuel usage, microbially induced corrosion, and, if not managed appropriately, eventual loss of both performance and structural integrity. Ship hull biofouling in the international maritime transport network conservatively accounts for 0.6% of global carbon emissions, highlighting the global scale and the importance of this problem. Improved antifouling strategies to limit surface colonization are paramount for essential activities such as shipping, aquaculture, desalination, and the marine renewable energy sector, representing both a multibillion dollar cost and a substantial practical challenge. From an ecological perspective, biofouling is a primary contributor to the global spread of invasive marine species, which has extensive implications for the marine environment.Historically, heavy metal-based toxic biocides have been used to control biofouling. However, their unwanted collateral ecological damage on nontarget species and bioaccumulation has led to recent global bans. With expanding human activities within aquaculture and offshore energy, it is both urgent and apparent that environmentally friendly surface protection remains key for maintaining the function of both moving and stationary marine structures. Biofouling communities are typically a highly complex network of both micro- and macroorganisms, representing a broad section of life from bacteria to macrophytes and animals. Given this diversity, it is unrealistic to expect that a single antifouling "silver bullet" will prevent colonization with the exception of generally toxic biocides. For that reason, modern and future antifouling solutions are anticipated to rely on novel coating technologies and "combination therapies" where mixtures of narrow-spectrum bioactive components are used to provide coverage across fouling species. In contrast to the existing cohort of outdated, toxic antifouling strategies, such as copper- and tributyltin-releasing paints, modern drug discovery techniques are increasingly being employed for the rational design of effective yet safe alternatives. The challenge for a medicinal chemistry approach is to effectively account for the large taxonomic diversity among fouling organisms combined with a lack of well-defined conserved molecular targets within most taxa.The current Account summarizes our work employing the tools of modern medicinal chemistry to discover, modify, and develop optimized and scalable antifouling solutions based on naturally occurring antifouling and repelling compounds from both marine and terrestrial sources. Inspiration for rational design comes from targeted studies on allelopathic natural products, natural repelling peptides, and secondary metabolites from sessile marine organisms with clean exteriors, which has yielded several efficient and promising antifouling leads.

Integrating QSAR modelling and deep learning in drug discovery: the emergence of deep QSAR

Quantitative structure-activity relationship (QSAR) modelling, an approach that was introduced 60 years ago, is widely used in computer-aided drug design. In recent years, progress in artificial intelligence techniques, such as deep learning, the rapid growth of databases of molecules for virtual screening and dramatic improvements in computational power have supported the emergence of a new field of QSAR applications that we term 'deep QSAR'. Marking a decade from the pioneering applications of deep QSAR to tasks involved in small-molecule drug discovery, we herein describe key advances in the field, including deep generative and reinforcement learning approaches in molecular design, deep learning models for synthetic planning and the application of deep QSAR models in structure-based virtual screening. We also reflect on the emergence of quantum computing, which promises to further accelerate deep QSAR applications and the need for open-source and democratized resources to support computer-aided drug design.

The role and potential of computer-aided drug discovery strategies in the discovery of novel antimicrobials

Antimicrobial resistance (AMR) has become more of a concern in recent decades, particularly in infections associated with global public health threats. The development of new antibiotics is crucial to ensuring infection control and eradicating AMR. Although drug discovery and development are essential processes in the transformation of a drug candidate from the laboratory to the bedside, they are often very complicated, expensive, and time-consuming. The pharmaceutical sector is continuously innovating strategies to reduce research costs and accelerate the development of new drug candidates. Computer-aided drug discovery (CADD) has emerged as a powerful and promising technology that renews the hope of researchers for the faster identification, design, and development of cheaper, less resource-intensive, and more efficient drug candidates. In this review, we discuss an overview of AMR, the potential, and limitations of CADD in AMR drug discovery, and case studies of the successful application of this technique in the rapid identification of various drug candidates. This review will aid in achieving a better understanding of available CADD techniques in the discovery of novel drug candidates against resistant pathogens and other infectious agents.

Clustering of atoms relative to vector space in the Z-matrix coordinate system and 'graphical fingerprint' analysis of 3D pharmacophore structure

The behavior of a molecule within its environment is governed by chemical fields present in 3D space. However, beyond local descriptors in 3D, the conformations a molecule assumes, and the resulting clusters also play a role in influencing structure-activity models. This study focuses on the clustering of atoms according to the vector space of four atoms aligned in the Z-Matrix Reference system for molecular similarity. Using 3D-QSAR analysis, it was aimed to determine the pharmacophore groups as interaction points in the binding region of the β2-adrenoceptor target of fenoterol stereoisomers. Different types of local reactive descriptors of ligands have been used to elucidate points of interaction with the target. Activity values for ligand-receptor interaction energy were determined using the Levenberg-Marquardt algorithm. Using the Molecular Comparative Electron Topology method, the 3D pharmacophore model (3D-PhaM) was obtained after aligning and superimposing the molecules and was further validated by the molecular docking method. Best guesses were calculated with a non-output validation (LOO-CV) method. Finally, the data were calculated using the 'graphic fingerprint' technique. Based on the eLKlopman (Electrostatic LUMO Klopman) descriptor, the Q2 value of this derivative set was calculated as 0.981 and the R2ext value is calculated as 0.998.

Reliable method for predicting the binding affinity of RNA-small molecule interactions using machine learning

Ribonucleic acids (RNAs) play important roles in cellular regulation. Consequently, dysregulation of both coding and non-coding RNAs has been implicated in several disease conditions in the human body. In this regard, a growing interest has been observed to probe into the potential of RNAs to act as drug targets in disease conditions. To accelerate this search for disease-associated novel RNA targets and their small molecular inhibitors, machine learning models for binding affinity prediction were developed specific to six RNA subtypes namely, aptamers, miRNAs, repeats, ribosomal RNAs, riboswitches and viral RNAs. We found that differences in RNA sequence composition, flexibility and polar nature of RNA-binding ligands are important for predicting the binding affinity. Our method showed an average Pearson correlation (r) of 0.83 and a mean absolute error of 0.66 upon evaluation using the jack-knife test, indicating their reliability despite the low amount of data available for several RNA subtypes. Further, the models were validated with external blind test datasets, which outperform other existing quantitative structure-activity relationship (QSAR) models. We have developed a web server to host the models, RNA-Small molecule binding Affinity Predictor, which is freely available at: https://web.iitm.ac.in/bioinfo2/RSAPred/.

Developing a Novel Read-Across Concept for Ecotoxicological Risk Assessment of Phosphate Chemicals: A Case Study

This study introduces a novel concept approach for a read-across assessment, considering species sensitivity differences among phosphate chemicals within structurally similar compound groups. Twenty-five organic chemicals, with a log Kow of 5 or less, were categorized into three functional groups based on acetylcholinesterase (AChE) inhibition as a specific mode of action (MOA). The short-term aquatic toxicity data (LC50) for fish, crustaceans, and insects were collected from the U.S. EPA Ecotoxicology (ECOTOX) Knowledgebase. A geometric mean calculation method was applied for multiple toxic endpoints. Performance metrics for the new read-across concept, including correlation coefficient, bias, precision, and accuracy, were calculated. Overall, a slightly higher overestimation (49.2%) than underestimation (48.4%) in toxicity predictions was observed in two case studies. In Case study I, a strong positive correlation (r = 0.93) between the predicted and known toxicity values of target chemicals was observed, while in Case study II, with limited information on species and their ecotoxicity, showed a moderate correlation (r = 0.75). Overall, the bias and precision for Case study I were 0.32 ± 0.01, while Case study II showed 0.65 ± 0.06; however, the relative bias (%) increased from 37.65% (Case study I) to 91.94% (Case study II). Bland-Altman plots highlight the mean differences of 1.33 (Case study I) and 1.24 (Case study II), respectively. The new read-across concept, focusing on AChE inhibition and structural similarity, demonstrated good reliability, applicability, and accuracy with minimal bias. Future studies are needed to evaluate various types of chemical substances, diverse modes of action, functional groups, toxic endpoints, and test species to ensure overall comprehensiveness and robustness in toxicity predictions.

Artificial Intelligence and Machine Learning Methods to Evaluate Cardiotoxicity following the Adverse Outcome Pathway Frameworks

Cardiovascular disease is a leading global cause of mortality. The potential cardiotoxic effects of chemicals from different classes, such as environmental contaminants, pesticides, and drugs can significantly contribute to effects on health. The same chemical can induce cardiotoxicity in different ways, following various Adverse Outcome Pathways (AOPs). In addition, the potential synergistic effects between chemicals further complicate the issue. In silico methods have become essential for tackling the problem from different perspectives, reducing the need for traditional in vivo testing, and saving valuable resources in terms of time and money. Artificial intelligence (AI) and machine learning (ML) are among today's advanced approaches for evaluating chemical hazards. They can serve, for instance, as a first-tier component of Integrated Approaches to Testing and Assessment (IATA). This study employed ML and AI to assess interactions between chemicals and specific biological targets within the AOP networks for cardiotoxicity, starting with molecular initiating events (MIEs) and progressing through key events (KEs). We explored methods to encode chemical information in a suitable way for ML and AI. We started with commonly used approaches in Quantitative Structure-Activity Relationship (QSAR) methods, such as molecular descriptors and different types of fingerprint. We then increased the complexity of encoders, incorporating graph-based methods, auto-encoders, and character embeddings employed in neural language processing. We also developed a multimodal neural network architecture, capable of considering the complementary nature of different chemical representations simultaneously. The potential of this approach, compared to more conventional architectures designed to handle a single encoder, becomes apparent when the amount of data increases.

Anticancer Drug Discovery Based on Natural Products: From Computational Approaches to Clinical Studies

Globally, malignancies cause one out of six mortalities, which is a serious health problem. Cancer therapy has always been challenging, apart from major advances in immunotherapies, stem cell transplantation, targeted therapies, hormonal therapies, precision medicine, and palliative care, and traditional therapies such as surgery, radiation therapy, and chemotherapy. Natural products are integral to the development of innovative anticancer drugs in cancer research, offering the scientific community the possibility of exploring novel natural compounds against cancers. The role of natural products like Vincristine and Vinblastine has been thoroughly implicated in the management of leukemia and Hodgkin's disease. The computational method is the initial key approach in drug discovery, among various approaches. This review investigates the synergy between natural products and computational techniques, and highlights their significance in the drug discovery process. The transition from computational to experimental validation has been highlighted through in vitro and in vivo studies, with examples such as betulinic acid and withaferin A. The path toward therapeutic applications have been demonstrated through clinical studies of compounds such as silvestrol and artemisinin, from preclinical investigations to clinical trials. This article also addresses the challenges and limitations in the development of natural products as potential anti-cancer drugs. Moreover, the integration of deep learning and artificial intelligence with traditional computational drug discovery methods may be useful for enhancing the anticancer potential of natural products.

Unraveling the Push-Pull Effect in Acenes, Polyenes and Polyynes

Substituent effects (SEs) are fundamental for predicting molecular reactivity, while polyene, polyyne and acene derivatives are precursors to compounds with diverse applications. Computations were performed for Y-R-X systems, where reaction sites Y=NO2 and O- , substituents X=NO2 , CN, Cl, H, OH, NH2 , and spacers R=polyene, polyyne (n=1-5, 10 repeating units) and acene (up to tetracene). The cSAR (charge of the substituent active region) approach allowed to present, for the first time, quantitative relations describing the spacer's electron-donating and withdrawing properties as a function of n and the spacer type. The electronic properties of the X substituents depend on the type of spacer, its length and the Y group, which is an example of the reverse SE. To describe how the SE between Y and X weakens with n, two approaches were compared: cSAR and SESE (SE stabilization energy). The EDDB (electron density of delocalized bonds) characterize changes in electron delocalization in spacers due to the SE. A new approach - EDDB differential maps - allow to extract the effect of X substitution on the electron delocalization. The charges at spacer's C atoms correlate with cSAR; changes in the slopes confirm the charge transfer by resonance.

Bioactivity predictions and virtual screening using machine learning predictive model

Recently, there has been significant attention on machine learning algorithms for predictive modeling. Prediction models for enzyme inhibitors are limited, and it is essential to account for chemical biases while developing them. The lack of repeatability in available models and chemical bias issues constrain drug discovery and development. A new prediction model for enzyme inhibitors has been developed, and the model efficacy was checked using Dipeptidyl peptidase 4 (DPP-4) inhibitors. A Python script was prepared and can be provided for personal use upon request. Among various machine learning algorithms, it was found that Random Forest offers the best accuracy. Two models were compared, one with diverse training and test data and the other with a random split. It was concluded that machine learning predictive models based on the Murcko scaffold can address chemical bias concerns. In-silico screening of the Drug Bank database identified two molecules against DPP-4, which are previously proven hit molecules. The approach was further validated through molecular docking studies and molecular dynamics simulations, demonstrating the credibility and relevance of the developed model for future investigations and potential translation into clinical applications.Communicated by Ramaswamy H. Sarma.

Viscosity of Ionic Liquids: Theories and Models

Ionic liquids (ILs) offer a wide range of promising applications due to their unique and designable properties compared to conventional solvents. Further development and application of ILs require correlating/predicting their pressure-viscosity-temperature behavior. In this review, we firstly introduce methods for calculation of thermodynamic inputs of viscosity models. Next, we introduce theories, theoretical and semi-empirical models coupling various theories with EoSs or activity coefficient models, and empirical and phenomenological models for viscosity of pure ILs and IL-related mixtures. Our modelling description is followed immediately by model application and performance. Then, we propose simple predictive equations for viscosity of IL-related mixtures and systematically compare performances of the above-mentioned theories and models. In concluding remarks, we recommend robust predictive models for viscosity at atmospheric pressure as well as proper and consistent theories and models for P-η-T behavior. The work that still remains to be done to obtain the desired theories and models for viscosity of ILs and IL-related mixtures is also presented. The present review is structured from pure ILs to IL-related mixtures and aims to summarize and quantitatively discuss the recent advances in theoretical and empirical modelling of viscosity of ILs and IL-related mixtures.

Computational Nanotoxicology Models for Environmental Risk Assessment of Engineered Nanomaterials

Although engineered nanomaterials (ENMs) have tremendous potential to generate technological benefits in numerous sectors, uncertainty on the risks of ENMs for human health and the environment may impede the advancement of novel materials. Traditionally, the risks of ENMs can be evaluated by experimental methods such as environmental field monitoring and animal-based toxicity testing. However, it is time-consuming, expensive, and impractical to evaluate the risk of the increasingly large number of ENMs with the experimental methods. On the contrary, with the advancement of artificial intelligence and machine learning, in silico methods have recently received more attention in the risk assessment of ENMs. This review discusses the key progress of computational nanotoxicology models for assessing the risks of ENMs, including material flow analysis models, multimedia environmental models, physiologically based toxicokinetics models, quantitative nanostructure-activity relationships, and meta-analysis. Several challenges are identified and a perspective is provided regarding how the challenges can be addressed.

Curated mode-of-action data and effect concentrations for chemicals relevant for the aquatic environment

Chemicals in the aquatic environment can be harmful to organisms and ecosystems. Knowledge on effect concentrations as well as on mechanisms and modes of interaction with biological molecules and signaling pathways is necessary to perform chemical risk assessment and identify toxic compounds. To this end, we developed criteria and a pipeline for harvesting and summarizing effect concentrations from the US ECOTOX database for the three aquatic species groups algae, crustaceans, and fish and researched the modes of action of more than 3,300 environmentally relevant chemicals in literature and databases. We provide a curated dataset ready to be used for risk assessment based on monitoring data and the first comprehensive collection and categorization of modes of action of environmental chemicals. Authorities, regulators, and scientists can use this data for the grouping of chemicals, the establishment of meaningful assessment groups, and the development of in vitro and in silico approaches for chemical testing and assessment.

Machine Learning-Assisted Computational Screening of Adhesive Molecules Derived from Dihydroxyphenyl Alanine

Marine mussels adhere to virtually any surface via 3,4-dihydroxyphenyl-L-alanines (L-DOPA), an amino acid largely contained in their foot proteins. The biofriendly, water-repellent, and strong adhesion of L-DOPA are unparalleled by any synthetic adhesive. Inspired by this, we computationally designed diverse derivatives of DOPA and studied their potential as adhesives or coating materials. We used first-principles calculations to investigate the adsorption of the DOPA derivatives on graphite. The presence of an electron-withdrawing group, such as nitrogen dioxide, strengthens the adsorption by increasing the π-π interaction between DOPA and graphite. To quantify the distribution of electron charge and to gain insights into the charge distribution at interfaces, we performed Bader charge analysis and examined charge density difference plots. We developed a quantitative structure-property relationship (QSPR) model using an artificial neural network (ANN) to predict the adsorption energy. Using the three-dimensional and quantum mechanical electrostatic potential of a molecule as a descriptor, the present quantum NN model shows promising performance as a predictive QSPR model.