Online chemical modeling environment (OCHEM): web platform for data storage, model development and publishing of chemical information

Investigation of in silico studies for cytochrome P450 isoforms specificity

Cytochrome P450 (CYP450) enzymes comprise a highly diverse superfamily of heme-thiolate proteins that responsible for catalyzing over 90 % of enzymatic reactions associated with xenobiotic metabolism in humans. Accurately predicting whether chemicals are substrates or inhibitors of different CYP450 isoforms can aid in pre-selecting hit compounds for the drug discovery process, chemical toxicology studies, and patients treatment planning. In this work, we investigated in silico studies on CYP450s specificity over past twenty years, categorizing these studies into structure-based and ligand-based approaches. Subsequently, we utilized 100 of the most frequently prescribed drugs to test eleven machine learning-based prediction models which were published between 2015 and 2024. We analyzed various aspects of the evaluated models, such as their datasets, algorithms, and performance. This will give readers with a comprehensive overview of these prediction models and help them choose the most suitable one to do prediction. We also provide our insights for future research trend in both structure-based and ligand-based approaches in this field.

Molecular designing of potential environmentally friendly PFAS based on deep learning and generative models

Perfluoroalkyl and polyfluoroalkyl substances (PFAS) are widely used across a spectrum of industrial and consumer goods. Nonetheless, their persistent nature and tendency to accumulate in biological systems pose substantial environmental and health threats. Consequently, striking a balance between maximizing product efficiency and minimizing environmental and health risks by tailoring the molecular structure of PFAS has become a pivotal challenge in the fields of environmental chemistry and sustainable development. To address this issue, a computational workflow was proposed for designing an environmentally friendly PFAS by incorporating deep learning (DL) and molecular generative models. The hybrid DL architecture MolHGT+ based on heterogeneous graph neural network with transformer-like attention was applied to predict the surface tension, bioaccumulation, and hepatotoxicity of the molecules. Through virtual screening of the PFAS master database using MolHGT+, the findings indicate that incorporating the siloxane group and betaine fragment can effectively decrease both the bioaccumulation and hepatotoxicity of PFAS while preserving low surface tension. In addition, molecular generative models were employed to create a structurally diverse pool of novel PFASs with the aforementioned hit molecules serving as the initial template structures. Overall, our study presents a promising AI-driven method for advancing the development of environmentally friendly PFAS.

Predicting fluorescence to singlet oxygen generation quantum yield ratio for BODIPY dyes using QSPR and machine learning

Functional dyes that are capable of both bright fluorescence and efficient singlet oxygen generation are crucial for theranostic techniques, which integrate fluorescence imaging and photodynamic therapy (PDT). The development of new functional dyes for theranostics is often costly and time-consuming due to laborious synthesis and post-synthetic screening of large libraries of compounds. In this work, we describe machine learning methods suitable for simultaneous prediction of fluorescence and photosensitizing ability of heavy-atom-free boron dipyrromethene (BODIPY) compounds. We analysed the ratio between fluorescence quantum yield (ΦFl) and singlet oxygen quantum yield (ΦΔ) for over 70 BODIPY structures in polar (acetonitrile) and non-polar (toluene) solvents, which mimic hydrophilic and hydrophobic cell environments, respectively. QSPR models were developed based on more than 5000 calculated molecular descriptors, including quantum chemical and topological descriptors. We applied multiple linear regression (MLR), support vector regression (SVR), and random forest regression (RFR) methods for model building and optimization. The resulting models demonstrated robust statistical parameters (R2 = 0.73-0.91) for both polar and non-polar media. The relative contributions of the descriptors to the models were assessed, identifying Eig03_EA(dm), F01[C-N], and TDB06p as the most influential. These results demonstrate that QSPR machine learning methods are effective in predicting key photochemical parameters of BODIPY photosensitizers, thereby potentially streamlining the development of theranostic agents.

Development of a machine learning-based target-specific scoring function for structure-based binding affinity prediction for human dihydroorotate dehydrogenase inhibitors

Human dihydroorotate dehydrogenase (hDHODH) is a flavin mononucleotide-dependent enzyme that can limit de novo pyrimidine synthesis, making it a therapeutic target for diseases such as autoimmune disorders and cancer. In this study, using the docking structures of complexes generated by AutoDock Vina, we integrate interaction features and ligand features, and employ support vector regression to develop a target-specific scoring function for hDHODH (TSSF-hDHODH). The Pearson correlation coefficient values of TSSF-hDHODH in the cross-validation and external validation are 0.86 and 0.74, respectively, both of which are far superior to those of classic scoring function AutoDock Vina and random forest (RF) based generic scoring function RF-Score. TSSF-hDHODH is further used for the virtual screening of potential inhibitors in the FDA-Approved & Pharmacopeia Drug Library. In conjunction with the results from molecular dynamics simulations, crizotinib is identified as a candidate for subsequent structural optimization. This study can be useful for the discovery of hDHODH inhibitors and the development of scoring functions for additional targets.

Computationally accelerated identification of P-glycoprotein inhibitors

Overexpression of the polyspecific efflux transporter, P-glycoprotein (P-gp, MDR1, ABCB1 ), is a major mechanism by which cancer cells acquire multidrug resistance (MDR), the resistance to diverse chemotherapeutic drugs. Inhibiting drug transport by P-gp can resensitize cancer cells to chemotherapy, but there are no P-gp inhibitors available to patients. Clinically unsuccessful P-gp inhibitors tend to bind at the pump's transmembrane drug binding domains and are often P-gp transport substrates, resulting in lowered intracellular concentration of the drug and altered pharmacokinetics. In prior work, we used computationally accelerated drug discovery to identify novel P-gp inhibitors that target the pump's cytoplasmic nucleotide binding domains. Our first-draft study provided conclusive evidence that the nucleotide binding domains of P-gp are viable targets for drug discovery. Here we develop an enhanced, computationally accelerated drug discovery pipeline that expands upon our prior work by iteratively screening compounds against multiple conformations of P-gp with molecular docking. Targeted molecular dynamics simulations with our homology model of human P-gp were used to generate docking receptors in conformations mimicking a putative drug transport cycle. We offset the increased computational complexity using custom Tanimoto chemical datasets, which maximize the chemical diversity of ligands screened by docking. Using our expanded, virtual-assisted pipeline, we identified nine novel P-gp inhibitors that reverse MDR in two types of P-gp overexpressing human cancer cell lines, reflecting a 13.4% hit rate. Of these inhibitors, all were non-toxic to non-cancerous human cells, and six were not likely to be transport substrates of P-gp. Our novel P-gp inhibitors are chemically diverse and are good candidates for lead optimization. Our results demonstrate that the nucleotide binding domains of P-gp are an underappreciated target in the effort to reverse P-gp-mediated multidrug resistance in cancer.

Prediction of Cytochrome P450 Substrates Using the Explainable Multitask Deep Learning Models

Cytochromes P450 (P450s or CYPs) are the most important phase I metabolic enzymes in the human body and are responsible for metabolizing ∼75% of the clinically used drugs. P450-mediated metabolism is also closely associated with the formation of toxic metabolites and drug-drug interactions. Therefore, it is of high importance to predict if a compound is the substrate of a given P450 in the early stage of drug development. In this study, we built the multitask learning models to simultaneously predict the substrates of five major drug-metabolizing P450 enzymes, namely, CYP3A4, 2C9, 2C19, 2D6, and 1A2, based on the collected substrate data sets. Compared to the single-task model and conventional machine learning models, the multitask fingerprints and graph neural networks model achieved superior performance with the average AUC values of 90.8% on the test set. Notably, the multitask model demonstrated its good performance on the small amount of substrate data sets such as CYP1A2, 2C9, and 2C19. In addition, the Shapley additive explanation and the attention mechanism were used to reveal specific substructures associated with P450 substrates, which were further confirmed and complemented by the substructure mining tool and the literature.

In Vitro Screening for ToxCast Chemicals Binding to Thyroxine-Binding Globulin

Thyroid hormone (TH) carrier proteins play an important role in distributing TH to target tissue as well as maintaining the balance of free versus bound TH in the blood. Interference with the TH carrier proteins has been identified as a potential mechanism of thyroid system disruption. To address the lack of data regarding chemicals binding to these carrier proteins and displacing TH, a fluorescence-based in vitro screening assay was utilized to screen over 1,400 chemicals from the U.S. EPA's ToxCast phase1_v2, phase 2, and e1k libraries for competitive binding to one of the carrier proteins, thyroxine-binding globulin. Initial screening at a single high concentration of 100 μM identified 714 chemicals that decreased signal of the bound fluorescent ligand by 20% or higher. Of these, 297 produced 50% or greater reduction in fluorescence and were further tested in concentration-response (0.004 to 150 μM) to determine relative potency. Ten chemicals were found to have EC50 values <1 μM, 63 < 10 μM, and 141 chemicals between 10 and 100 μM. Utilization of this assay contributes to expanding the number of in vitro assays available for identifying chemicals with the potential to disrupt TH homeostasis. These results support ranking and prioritization of chemicals to be tested in vivo to aid in the development of a framework for predicting in vivo effects from in vitro high-throughput data.

Screening the ToxCast Chemical Libraries for Binding to Transthyretin

Transthyretin (TTR) is one of the serum binding proteins responsible for transport of thyroid hormones (TH) to target tissue and for maintaining the balance of available TH. Chemical binding to TTR and subsequent displacement of TH has been identified as an end point in screening chemicals for potential disruption of the thyroid system. To address the lack of data regarding chemicals binding to TTR, we optimized an in vitro assay utilizing the fluorescent probe 8-anilino-1-napthalenesulfonic acid (ANSA) and the human protein TTR to screen over 1500 chemicals from the U.S. EPA's ToxCast ph1_v2, ph2, and e1k libraries utilizing a tiered approach. Testing of a single high concentration (target 100 μM) resulted in 888 chemicals with 20% or greater activity based on displacement of ANSA from TTR. Of these, 282 chemicals had activity of 85% or greater and were further tested in 12-point concentration-response with target concentrations ranging from 0.015 to 100 μM. An EC50 was obtained for 276 of these 301 chemicals. To date, this is the largest set of chemicals screened for binding to TTR. Utilization of this assay is a significant contribution toward expanding the suite of in vitro assays used to identify chemicals with the potential to disrupt thyroid hormone homeostasis.

Toxicity of Bromo-DragonFLY as a New Psychoactive Substance: Application of In Silico Methods for the Prediction of Key Toxicological Parameters Important to Clinical and Forensic Toxicology

Bromo-DragonFLY is a synthetic new psychoactive substance (NPS) that has gained attention due to its powerful and long-lasting hallucinogenic effects, legal status, and widespread availability. This study aimed to use various in silico toxicology methods to predict key toxicological parameters for Bromo-DragonFLY, including acute toxicity (LD50), genotoxicity, cardiotoxicity, health effects, and the potential for endocrine disruption. The results indicate significant acute toxicity with noticeable variations across different species, a low likelihood of genotoxic potential suggesting potential DNA damage, and a notable risk of cardiotoxicity associated with inhibition of the hERG channel. Evaluation of endocrine disruption suggests a low probability of Bromo-DragonFLY interacting with the estrogen receptor α (ER-α), indicating minimal estrogenic activity. These insights from in silico investigations are important for advancing our understanding of this NPS in forensic and clinical toxicology. These initial toxicological examinations establish a foundation for future research efforts and contribute to developing risk assessment and management strategies for using and misusing NPS.

In Silico Design and In Vitro Assessment of Bicyclic Trifluoromethylated Pyrroles as New Antibacterial and Antifungal Agents

QSAR studies on the number of compounds tested as S. aureus inhibitors were performed using an interactive Online Chemical Database and Modeling Environment (OCHEM) web platform. The predictive ability of the developed consensus QSAR model was q2=0.79±0.02. The consensus prediction for the external evaluation set afforded high predictive power (q2=0.82±0.03). The models were applied to screen a virtual chemical library with anti-S. aureus activity. Six promising new bicyclic trifluoromethylated pyrroles were identified, synthesized and evaluated in vitro against S. aureus, E. coli, and A. baumannii for their antibacterial activity and against C. albicans, C. krusei and C. glabrata for their antifungal activity. The synthesized compounds were characterized by 1H, 19F, and 13C NMR and elemental analysis. The antimicrobial activity assessment indicated that trifluoromethylated pyrroles 9 and 11 demonstrated the greatest antibacterial and antifungal effects against all the tested pathogens, especially against multidrug-resistant strains. The acute toxicity of the compounds to Daphnia magna ranged from 1.21 to 33.39 mg/L (moderately and slightly toxic). Based on the docking results, it can be suggested that the antibacterial and antifungal effects of the compounds can be explained by the inhibition of bacterial wall component synthesis.

Classification models for predicting the bioactivity of pan-TRK inhibitors and SAR analysis

Tropomyosin receptor kinases (TRKs) are important broad-spectrum anticancer targets. The oncogenic rearrangement of the NTRK gene disrupts the extracellular structural domain and epitopes for therapeutic antibodies, making small-molecule inhibitors essential for treating NTRK fusion-driven tumors. In this work, several algorithms were used to construct descriptor-based and nondescriptor-based models, and the models were evaluated by outer 10-fold cross-validation. To find a model with good generalization ability, the dataset was partitioned by random and cluster-splitting methods to construct in- and cross-domain models, respectively. Among the 48 models built, the model with the combination of the deep neural network (DNN) algorithm and extended connectivity fingerprints 4 (ECFP4) descriptors achieved excellent performance in both dataset divisions. The results indicate that the DNN algorithm has a strong generalization prediction ability, and the richness of features plays a vital role in predicting unknown spatial molecules. Additionally, we combined the clustering results and decision tree models of fingerprint descriptors to perform structure-activity relationship analysis. It was found that nitrogen-containing aromatic heterocyclic and benzo heterocyclic structures play a crucial role in enhancing the activity of TRK inhibitors.

QSARtuna: An Automated QSAR Modeling Platform for Molecular Property Prediction in Drug Design

Machine-learning (ML) and deep-learning (DL) approaches to predict the molecular properties of small molecules are increasingly deployed within the design-make-test-analyze (DMTA) drug design cycle to predict molecular properties of interest. Despite this uptake, there are only a few automated packages to aid their development and deployment that also support uncertainty estimation, model explainability, and other key aspects of model usage. This represents a key unmet need within the field, and the large number of molecular representations and algorithms (and associated parameters) means it is nontrivial to robustly optimize, evaluate, reproduce, and deploy models. Here, we present QSARtuna, a molecule property prediction modeling pipeline, written in Python and utilizing the Optuna, Scikit-learn, RDKit, and ChemProp packages, which enables the efficient and automated comparison between molecular representations and machine learning models. The platform was developed by considering the increasingly important aspect of model uncertainty quantification and explainability by design. We provide details for our framework and provide illustrative examples to demonstrate the capability of the software when applied to simple molecular property, reaction/reactivity prediction, and DNA encoded library enrichment classification. We hope that the release of QSARtuna will further spur innovation in automatic ML modeling and provide a platform for education of best practices in molecular property modeling. The code for the QSARtuna framework is made freely available via GitHub.

ADMETlab 3.0: an updated comprehensive online ADMET prediction platform enhanced with broader coverage, improved performance, API functionality and decision support

ADMETlab 3.0 is the second updated version of the web server that provides a comprehensive and efficient platform for evaluating ADMET-related parameters as well as physicochemical properties and medicinal chemistry characteristics involved in the drug discovery process. This new release addresses the limitations of the previous version and offers broader coverage, improved performance, API functionality, and decision support. For supporting data and endpoints, this version includes 119 features, an increase of 31 compared to the previous version. The updated number of entries is 1.5 times larger than the previous version with over 400 000 entries. ADMETlab 3.0 incorporates a multi-task DMPNN architecture coupled with molecular descriptors, a method that not only guaranteed calculation speed for each endpoint simultaneously, but also achieved a superior performance in terms of accuracy and robustness. In addition, an API has been introduced to meet the growing demand for programmatic access to large amounts of data in ADMETlab 3.0. Moreover, this version includes uncertainty estimates in the prediction results, aiding in the confident selection of candidate compounds for further studies and experiments. ADMETlab 3.0 is publicly for access without the need for registration at: https://admetlab3.scbdd.com.

Graph theoretical descriptors differentiate d-Mannose isomers in the principal component proposed feature space: A computational approach

The intricate nature of carbohydrates, particularly monosaccharides, stems from the existence of several chiral centers within their tertiary structures. Predicting and characterizing the molecular geometries and electrostatic landscapes of these substances is difficult due to their complex electrical properties. Moreover, these structures can display a substantial degree of conformational flexibility due to the presence of many rotatable bonds. Moreover, identifying and distinguishing between D and L enantiomers of monosaccharides presents a significant analytical obstacle since there is a need for empirically measurable properties that can distinguish them. This work uses Principal Component Analysis (PCA) to explore the chemical information included in 3D descriptors in order to comprehend the conformational space of d-Mannose stereoisomers. The isomers may be discriminated by utilizing 3D matrix-based indices, geometrical descriptors, and RDF descriptors. The isomers can be distinguished by descriptors, such as the Harary-like index from the reciprocal squared geometrical matrix (H_RG), Harary-like index from Coulomb matrix (H_Coulomb), Wiener-like index from Coulomb matrix (Wi_Coulomb), Wiener-like index from geometrical matrix (Wi_G), Graph energy from Coulomb matrix (SpAbs_Coulomb), Spectral absolute deviation from Coulomb matrix (SpAD_Coulomb), and Spectral positive sum from Coulomb matrix (SpPos_Coulomb). Among these descriptors, the first two, H_RG and H_Coulomb, perform the best in differentiation among the 3D-Matrix-Based Descriptors (3D-MBD) class. The results obtained from this study provide a significant chemical insight into the structural characteristics of the compounds inside the graph theoretical framework. These findings are likely to serve as the basis for developing new methods for analytical experiments.

Physicochemical modelling of the retention mechanism of temperature-responsive polymeric columns for HPLC through machine learning algorithms

Temperature-responsive liquid chromatography (TRLC) offers a promising alternative to reversed-phase liquid chromatography (RPLC) for environmentally friendly analytical techniques by utilizing pure water as a mobile phase, eliminating the need for harmful organic solvents. TRLC columns, packed with temperature-responsive polymers coupled to silica particles, exhibit a unique retention mechanism influenced by temperature-induced polymer hydration. An investigation of the physicochemical parameters driving separation at high and low temperatures is crucial for better column manufacturing and selectivity control. Assessment of predictability using a dataset of 139 molecules analyzed at different temperatures elucidated the molecular descriptors (MDs) relevant to retention mechanisms. Linear regression, support vector regression (SVR), and tree-based ensemble models were evaluated, with no standout performer. The precision, accuracy, and robustness of models were validated through metrics, such as r and mean absolute error (MAE), and statistical analysis. At 45 ∘ C , logP predominantly influenced retention, akin to reversed-phase columns, while at5 ∘ C , complex interactions with lipophilic and negative MDs, along with specific functional groups, dictated retention. These findings provide deeper insights into TRLC mechanisms, facilitating method development and maximizing column potential.

New Pyrimidinone Bearing Aminomethylenes and Schiff Bases as Potent Antioxidant, Antibacterial, SARS-CoV-2, and COVID-19 Main Protease MPro Inhibitors: Design, Synthesis, Bioactivities, and Computational Studies

New 2-thioxopyrimidinone derivatives (A1-A10) were synthesized in 87-96% yields via a simple three-component condensation reaction. These compounds were screened extensively through in vitro assays for antioxidant and antibacterial investigations. The DPPH assays resulted in the excellent potency of A6-A10 as antioxidants with IC50 values of 0.83 ± 0.125, 0.90 ± 0.77, 0.36 ± 0.063, 1.4 ± 0.07, and 1.18 ± 0.06 mg/mL, which were much better than 1.79 ± 0.045 mg/mL for the reference ascorbic acid. These compounds exhibited better antibacterial potency against Klebsiella with IC50 values of 2 ± 7, 1.32 ± 8.9, 1.19 ± 11, 1.1 ± 12, and 1.16 ± 11 mg/mL for A6-A10. High-throughput screenings (HTS) of these motifs were carried out including investigation of drug-like behaviors, physiochemical property evaluation, and structure-related studies involving DFT and metabolic transformation trends. The radical scavenging ability of the synthesized motifs was validated through molecular docking studies through ligand-protein binding against human inducible nitric oxide synthase (HINOS) PDB ID: 4NOS, and the results were promising. Furthermore, the antiviral capability of the compounds was examined by in silico studies using two viral proteins PDB ID: 6Y84 and PDB ID: 6LU7. Binding poses of ligands were discussed, and amino acids in the protein binding pockets were investigated, where the tested compounds showed much better binding affinities than the standard inhibitors, proving to be suitable leads for antiviral drug discovery. The stabilities of the molecular docked complexes in real systems were validated by molecular dynamics simulations.

OLB-AC: toward optimizing ligand bioactivities through deep graph learning and activity cliffs

MOTIVATION

Deep graph learning (DGL) has been widely employed in the realm of ligand-based virtual screening. Within this field, a key hurdle is the existence of activity cliffs (ACs), where minor chemical alterations can lead to significant changes in bioactivity. In response, several DGL models have been developed to enhance ligand bioactivity prediction in the presence of ACs. Yet, there remains a largely unexplored opportunity within ACs for optimizing ligand bioactivity, making it an area ripe for further investigation.

RESULTS

We present a novel approach to simultaneously predict and optimize ligand bioactivities through DGL and ACs (OLB-AC). OLB-AC possesses the capability to optimize ligand molecules located near ACs, providing a direct reference for optimizing ligand bioactivities with the matching of original ligands. To accomplish this, a novel attentive graph reconstruction neural network and ligand optimization scheme are proposed. Attentive graph reconstruction neural network reconstructs original ligands and optimizes them through adversarial representations derived from their bioactivity prediction process. Experimental results on nine drug targets reveal that out of the 667 molecules generated through OLB-AC optimization on datasets comprising 974 low-activity, noninhibitor, or highly toxic ligands, 49 are recognized as known highly active, inhibitor, or nontoxic ligands beyond the datasets' scope. The 27 out of 49 matched molecular pairs generated by OLB-AC reveal novel transformations not present in their training sets. The adversarial representations employed for ligand optimization originate from the gradients of bioactivity predictions. Therefore, we also assess OLB-AC's prediction accuracy across 33 different bioactivity datasets. Results show that OLB-AC achieves the best Pearson correlation coefficient (r2) on 27/33 datasets, with an average improvement of 7.2%-22.9% against the state-of-the-art bioactivity prediction methods.

AVAILABILITY AND IMPLEMENTATION

The code and dataset developed in this work are available at github.com/Yueming-Yin/OLB-AC.

Toxicity of the New Psychoactive Substance (NPS) Clephedrone (4-Chloromethcathinone, 4-CMC): Prediction of Toxicity Using In Silico Methods for Clinical and Forensic Purposes

This study reports the first application of in silico methods to assess the toxicity of 4-chloromethcathinone (4-CMC), a novel psychoactive substance (NPS). Employing advanced toxicology in silico tools, it was possible to predict crucial aspects of the toxicological profile of 4-CMC, including acute toxicity (LD50), genotoxicity, cardiotoxicity, and its potential for endocrine disruption. The obtained results indicate significant acute toxicity with species-specific variability, moderate genotoxic potential suggesting the risk of DNA damage, and a notable cardiotoxicity risk associated with hERG channel inhibition. Endocrine disruption assessment revealed a low probability of 4-CMC interacting with estrogen receptor alpha (ER-α), suggesting minimal estrogenic activity. These insights, derived from in silico studies, are critical in advancing the understanding of 4-CMC properties in forensic and clinical toxicology. These initial toxicological findings provide a foundation for future research and aid in the formulation of risk assessment and management strategies in the context of the use and abuse of NPSs.

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

Enhancing Multi-species Liver Microsomal Stability Prediction through Artificial Intelligence

Liver microsomal stability, a crucial aspect of metabolic stability, significantly impacts practical drug discovery. However, current models for predicting liver microsomal stability are based on limited molecular information from a single species. To address this limitation, we constructed the largest public database of compounds from three common species: human, rat, and mouse. Subsequently, we developed a series of classification models using both traditional descriptor-based and classic graph-based machine learning (ML) algorithms. Remarkably, the best-performing models for the three species achieved Matthews correlation coefficients (MCCs) of 0.616, 0.603, and 0.574, respectively, on the test set. Furthermore, through the construction of consensus models based on these individual models, we have demonstrated their superior predictive performance in comparison with the existing models of the same type. To explore the similarities and differences in the properties of liver microsomal stability among multispecies molecules, we conducted preliminary interpretative explorations using the Shapley additive explanations (SHAP) and atom heatmap approaches for the models and misclassified molecules. Additionally, we further investigated representative structural modifications and substructures that decrease the liver microsomal stability in different species using the matched molecule pair analysis (MMPA) method and substructure extraction techniques. The established prediction models, along with insightful interpretation information regarding liver microsomal stability, will significantly contribute to enhancing the efficiency of exploring practical drugs for development.

Advancing Adverse Drug Reaction Prediction with Deep Chemical Language Model for Drug Safety Evaluation

The accurate prediction of adverse drug reactions (ADRs) is essential for comprehensive drug safety evaluation. Pre-trained deep chemical language models have emerged as powerful tools capable of automatically learning molecular structural features from large-scale datasets, showing promising capabilities for the downstream prediction of molecular properties. However, the performance of pre-trained chemical language models in predicting ADRs, especially idiosyncratic ADRs induced by marketed drugs, remains largely unexplored. In this study, we propose MoLFormer-XL, a pre-trained model for encoding molecular features from canonical SMILES, in conjunction with a CNN-based model to predict drug-induced QT interval prolongation (DIQT), drug-induced teratogenicity (DIT), and drug-induced rhabdomyolysis (DIR). Our results demonstrate that the proposed model outperforms conventional models applied in previous studies for predicting DIQT, DIT, and DIR. Notably, an analysis of the learned linear attention maps highlights amines, alcohol, ethers, and aromatic halogen compounds as strongly associated with the three types of ADRs. These findings hold promise for enhancing drug discovery pipelines and reducing the drug attrition rate due to safety concerns.

Expanding Predictive Capacities in Toxicology: Insights from Hackathon-Enhanced Data and Model Aggregation

In the realm of predictive toxicology for small molecules, the applicability domain of QSAR models is often limited by the coverage of the chemical space in the training set. Consequently, classical models fail to provide reliable predictions for wide classes of molecules. However, the emergence of innovative data collection methods such as intensive hackathons have promise to quickly expand the available chemical space for model construction. Combined with algorithmic refinement methods, these tools can address the challenges of toxicity prediction, enhancing both the robustness and applicability of the corresponding models. This study aimed to investigate the roles of gradient boosting and strategic data aggregation in enhancing the predictivity ability of models for the toxicity of small organic molecules. We focused on evaluating the impact of incorporating fragment features and expanding the chemical space, facilitated by a comprehensive dataset procured in an open hackathon. We used gradient boosting techniques, accounting for critical features such as the structural fragments or functional groups often associated with manifestations of toxicity.

Excito-repellent and Pesticide-Likeness Properties of Essential Oils on Carpophilus dimidiatus (Fabricius) (Nitidulidae) and Oryzaephilus mercator (L.) (Silvanidae)

Carpophilus dimidiatus (Fabricius) (Nitidulidae) and Oryzaephilus mercator (L.) (Silvanidae) are insect pests that cause severe damage in important walnut growing regions in the northwest of Argentina. The current management approaches for these pests involve the use of unsafe phosphorus pesticides whose overuse have led to farmworker poisoning, pest resistance issues, and environmental contamination. Plant extracts, particularly essential oils, are an alternative source of insect control agents. Excito-repellent essential oils can be used to develop ecofriendly tools for managing the pest population without affecting quality and visual appearance of the stored walnuts. Laboratory studies were conducted to assess the excito-repellent effects of C. dimidiatus and O. mercator of 12 essential oils derived from aromatic plants used as food additives and traditional medicine in Argentina: Aloysia citrodora (AC), Aloysia gratissima (AG), Aloysia gratissima var. Gratissima (AGG), Blepharocalyx salicipholius (BS), Hyptis mutabilis (HM), Lippia junelliana (LJ), Lippia turbinata (LT), Mentha x piperita (MP), Minthostachys mollis (MM), Minthostachys verticillata (MV), Origanum vulgare(OV), and Rosmarinus officinalis (RO). The most bioactive EOs (ERijk ≥ 70%) were Aloysia gratissima var. Gratissima (AGG), Minthostachys verticillata, and Lippia junelliana. Their bioactivity profile and chemical space, characterized from GC-MS measures, Generalized Estimating Equations, and Hierarchical Cluster Analysis, revealed that they are mixtures of very functionalized molecules with physicochemical properties similar to those of insecticides with low residual property that enter the insect body through the respiratory system by inhalation. The AGG, MV, and LJ oils are promising as protective agents of walnut products. In our laboratory, studies of their formulations for use in integrated pest management programs are still ongoing.

In Silico Design and Experimental Validation of Novel Oxazole Derivatives Against Varicella zoster virus

Varicella zoster virus (VZV) infection causes severe disease such as chickenpox, shingles, and postherpetic neuralgia, often leading to disability. Reactivation of latent VZV is associated with a decrease in specific cellular immunity in the elderly and in patients with immunodeficiency. However, due to the limited efficacy of existing therapy and the emergence of antiviral resistance, it has become necessary to develop new and effective antiviral drugs for the treatment of diseases caused by VZV, particularly in the setting of opportunistic infections. The goal of this work is to identify potent oxazole derivatives as anti-VZV agents by machine learning, followed by their synthesis and experimental validation. Predictive QSAR models were developed using the Online Chemical Modeling Environment (OCHEM). Data on compounds exhibiting antiviral activity were collected from the ChEMBL and uploaded in the OCHEM database. The predictive ability of the models was tested by cross-validation, giving coefficient of determination q2 = 0.87-0.9. The validation of the models using an external test set proves that the models can be used to predict the antiviral activity of newly designed and known compounds with reasonable accuracy within the applicability domain (q2 = 0.83-0.84). The models were applied to screen a virtual chemical library with expected activity of compounds against VZV. The 7 most promising oxazole derivatives were identified, synthesized, and tested. Two of them showed activity against the VZV Ellen strain upon primary in vitro antiviral screening. The synthesized compounds may represent an interesting starting point for further development of the oxazole derivatives against VZV. The developed models are available online at OCHEM http://ochem.eu/article/145978 and can be used to virtually screen for potential compounds with anti-VZV activity.

Will we ever be able to accurately predict solubility?

Accurate prediction of thermodynamic solubility by machine learning remains a challenge. Recent models often display good performances, but their reliability may be deceiving when used prospectively. This study investigates the origins of these discrepancies, following three directions: a historical perspective, an analysis of the aqueous solubility dataverse and data quality. We investigated over 20 years of published solubility datasets and models, highlighting overlooked datasets and the overlaps between popular sets. We benchmarked recently published models on a novel curated solubility dataset and report poor performances. We also propose a workflow to cure aqueous solubility data aiming at producing useful models for bench chemist. Our results demonstrate that some state-of-the-art models are not ready for public usage because they lack a well-defined applicability domain and overlook historical data sources. We report the impact of factors influencing the utility of the models: interlaboratory standard deviation, ionic state of the solute and data sources. The herein obtained models, and quality-assessed datasets are publicly available.

Alcohol-free synthesis, biological assessment, in vivo toxicological evaluation, and in silico analysis of novel silane quaternary ammonium compounds differing in structure and chain length as promising disinfectants

Quaternary ammonium compounds (QACs) are commonly used as disinfectants for industrial, medical, and residential applications. However, adverse health outcomes have been reported. Therefore, biocompatible disinfectants must be developed to reduce these adverse effects. In this context, QACs with various alkyl chain lengths (C12-C18) were synthesized by reacting QACs with the counterion silane. The antimicrobial activities of the novel compounds against four strains of microorganisms were assessed. Several in vivo assays were conducted on Drosophila melanogaster to determine the toxicological outcomes of Si-QACs, followed by computational analyses (molecular docking, simulation, and prediction of skin sensitization). The in vivo results were combined using a cheminformatics approach to understand the descriptors responsible for the safety of Si-QAC. Si-QAC-2 was active against all tested bacteria, with minimal inhibitory concentrations ranging from 13.65 to 436.74 ppm. Drosophila exposed to Si-QAC-2 have moderate-to-low toxicological outcomes. The molecular weight, hydrophobicity/lipophilicity, and electron diffraction properties were identified as crucial descriptors for ensuring the safety of the Si-QACs. Furthermore, Si-QAC-2 exhibited good stability and notable antiviral potential with no signs of skin sensitization. Overall, Si-QAC-2 (C14) has the potential to be a novel disinfectant.

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.

The openOCHEM consensus model is the best-performing open-source predictive model in the First EUOS/SLAS joint compound solubility challenge

The EUOS/SLAS challenge aimed to facilitate the development of reliable algorithms to predict the aqueous solubility of small molecules using experimental data from 100 K compounds. In total, hundred teams took part in the challenge to predict low, medium and highly soluble compounds as measured by the nephelometry assay. This article describes the winning model, which was developed using the publicly available Online CHEmical database and Modeling environment (OCHEM) available on the website https://ochem.eu/article/27. We describe in detail the assumptions and steps used to select methods, descriptors and strategy which contributed to the winning solution. In particular we show that consensus based on 28 models calculated using descriptor-based and representation learning methods allowed us to obtain the best score, which was higher than those based on individual approaches or consensus models developed using each individual approach. A combination of diverse models allowed us to decrease both bias and variance of individual models and to calculate the highest score. The model based on Transformer CNN contributed the best individual score thus highlighting the power of Natural Language Processing (NLP) methods. The inclusion of information about aleatoric uncertainty would be important to better understand and use the challenge data by the contestants.

Shaping the Future of Obesity Treatment: In Silico Multi-Modeling of IP6K1 Inhibitors for Obesity and Metabolic Dysfunction

Recent research has uncovered a promising approach to addressing the growing global health concern of obesity and related disorders. The inhibition of inositol hexakisphosphate kinase 1 (IP6K1) has emerged as a potential therapeutic strategy. This study employs multiple ligand-based in silico modeling techniques to investigate the structural requirements for benzisoxazole derivatives as IP6K1 inhibitors. Firstly, we developed linear 2D Quantitative Structure-Activity Relationship (2D-QSAR) models to ensure both their mechanistic interpretability and predictive accuracy. Then, ligand-based pharmacophore modeling was performed to identify the essential features responsible for the compounds' high activity. To gain insights into the 3D requirements for enhanced potency against the IP6K1 enzyme, we employed multiple alignment techniques to set up 3D-QSAR models. Given the absence of an available X-ray crystal structure for IP6K1, a reliable homology model for the enzyme was developed and structurally validated in order to perform structure-based analyses on the selected dataset compounds. Finally, molecular dynamic simulations, using the docked poses of these compounds, provided further insights. Our findings consistently supported the mechanistic interpretations derived from both ligand-based and structure-based analyses. This study offers valuable guidance on the design of novel IP6K1 inhibitors. Importantly, our work exclusively relies on non-commercial software packages, ensuring accessibility for reproducing the reported models.

Screening of the Antagonistic Activity of Potential Bisphenol A Alternatives toward the Androgen Receptor Using Machine Learning and Molecular Dynamics Simulation

Over the past few decades, extensive research has indicated that exposure to bisphenol A (BPA) increases the health risks in humans. Toxicological studies have demonstrated that BPA can bind to the androgen receptor (AR), resulting in endocrine-disrupting effects. In recent investigations, many alternatives to BPA have been detected in various environmental media as major pollutants. However, related experimental evaluations of BPA alternatives have not been systematically implemented for the assessment of chemical safety and the effects of structural characteristics on the antagonistic activity of the AR. To promote the green development of BPA alternatives, high-throughput toxicological screening is fundamental for prioritizing chemical tests. Therefore, we proposed a hybrid deep learning architecture that combines molecular descriptors and molecular graphs to predict AR antagonistic activity. Compared to previous models, this hybrid architecture can extract substantial chemical information from various molecular representations to improve the model's generalization ability for BPA alternatives. Our predictions suggest that lignin-derivable bisguaiacols, as alternatives to BPA, are likely to be nonantagonist for AR compared to bisphenol analogues. Additionally, molecular dynamics (MD) simulations identified the dihydrotestosterone-bound pocket, rather than the surface, as the major binding site of bisphenol analogues. The conformational changes of key helix H12 from an agonistic to an antagonistic conformation can be evaluated qualitatively by accelerated MD simulations to explain the underlying mechanism. Overall, our computational study is helpful for toxicological screening of BPA alternatives and the design of environmentally friendly BPA alternatives.

Predicting Retention Time in Unified-Hydrophilic-Interaction/Anion-Exchange Liquid Chromatography High-Resolution Tandem Mass Spectrometry (Unified-HILIC/AEX/HRMS/MS) for Comprehensive Structural Annotation of Polar Metabolome

The accuracy of the structural annotation of unidentified peaks obtained in metabolomic analysis using liquid chromatography/tandem mass spectrometry (LC/MS/MS) can be enhanced using retention time (RT) information as well as precursor and product ions. Unified-hydrophilic-interaction/anion-exchange liquid chromatography high-resolution tandem mass spectrometry (unified-HILIC/AEX/HRMS/MS) has been recently developed as an innovative method ideal for nontargeted polar metabolomics. However, the RT prediction for unified-HILIC/AEX has not been developed because of the complex separation mechanism characterized by the continuous transition of the separation modes from HILIC to AEX. In this study, we propose an RT prediction model of unified-HILIC/AEX/HRMS/MS, which enables the comprehensive structural annotation of polar metabolites. With training data for 203 polar metabolites, we ranked the feature importance using a random forest among 12,420 molecular descriptors (MDs) and constructed an RT prediction model with 26 selected MDs. The accuracy of the RT model was evaluated using test data for 51 polar metabolites, and 86.3% of the ΔRTs (difference between measured and predicted RTs) were within ±1.50 min, with a mean absolute error of 0.80 min, indicating high RT prediction accuracy. Nontargeted metabolomic data from the NIST SRM 1950-Metabolites in frozen human plasma were analyzed using the developed RT model and in silico MS/MS prediction, resulting in a successful structural estimation of 216 polar metabolites, in addition to the 62 identified based on standards. The proposed model can help accelerate the structural annotation of unknown hydrophilic metabolites, which is a key issue in metabolomic research.

Artificial Intelligence for Drug Discovery: Are We There Yet?

Drug discovery is adapting to novel technologies such as data science, informatics, and artificial intelligence (AI) to accelerate effective treatment development while reducing costs and animal experiments. AI is transforming drug discovery, as indicated by increasing interest from investors, industrial and academic scientists, and legislators. Successful drug discovery requires optimizing properties related to pharmacodynamics, pharmacokinetics, and clinical outcomes. This review discusses the use of AI in the three pillars of drug discovery: diseases, targets, and therapeutic modalities, with a focus on small-molecule drugs. AI technologies, such as generative chemistry, machine learning, and multiproperty optimization, have enabled several compounds to enter clinical trials. The scientific community must carefully vet known information to address the reproducibility crisis. The full potential of AI in drug discovery can only be realized with sufficient ground truth and appropriate human intervention at later pipeline stages.

Anticancer activity features of imidazole-based ionic liquids and lysosomotropic detergents: in silico and in vitro studies

Long-chain imidazole-based ionic liquids (compounds 2, 4, 9) and lysosomotropic detergents (compounds 7, 3, 8) with potent anticancer activity were synthesized. Their inhibitory activities against neuroblastoma and leukaemia cell lines were predicted by the new in silico QSAR models. The cytotoxic activities of the synthesized imidazole derivatives were investigated on the SK-N-DZ (human neuroblastoma) and K-562 (human chronic myeloid leukaemia) cell lines. Compounds 2 and 7 showed the highest in vitro cytotoxic effect on both cancer cell lines. The docking procedure of compounds 2 and 7 into the NAD+ coenzyme binding site of deacetylase Sirtuin-1 (SIRT-1) showed the formation of protein-ligand complexes with calculated binding energies of - 8.0 and - 8.1 kcal/mol, respectively. The interaction of SIRT1 with compounds 2, 7 and 9 and the interaction of Bromodomain-containing protein 4 (BRD4) with compounds 7 and 9 were also demonstrated by thermal shift assay. Compounds 2, 4, 7 and 9 inhibited SIRT1 deacetylase activity in the SIRT-Glo assay. Compounds 7 and 9 showed a moderate inhibitory activity against Aurora kinase A. In addition, compounds 3, 4, 8 and 9 inhibited the Janus kinase 2 activity. The results obtained showed that long-chain imidazole derivatives exhibited cytotoxic activities on K562 leukaemia and SK-N-DZ neuroblastoma cell lines. Furthermore, these compounds inhibited a panel of molecular targets involved in leukaemia and neuroblastoma tumorigenesis. All these results suggest that both long-chain imidazole-based ionic liquids and lysosomotropic detergents may be an effective alternative for the treatment of neuroblastoma and chronic myeloid leukemia and merit further investigation.

Combining Experimental and Computational Methods to Produce Conjugates of Anticholinesterase and Antioxidant Pharmacophores with Linker Chemistries Affecting Biological Activities Related to Treatment of Alzheimer's Disease

Effective therapeutics for Alzheimer's disease (AD) are in great demand worldwide. In our previous work, we responded to this need by synthesizing novel drug candidates consisting of 4-amino-2,3-polymethylenequinolines conjugated with butylated hydroxytoluene via fixed-length alkylimine or alkylamine linkers (spacers) and studying their bioactivities pertaining to AD treatment. Here, we report significant extensions of these studies, including the use of variable-length spacers and more detailed biological characterizations. Conjugates were potent inhibitors of acetylcholinesterase (AChE, the most active was 17d IC50 15.1 ± 0.2 nM) and butyrylcholinesterase (BChE, the most active was 18d: IC50 5.96 ± 0.58 nM), with weak inhibition of off-target carboxylesterase. Conjugates with alkylamine spacers were more effective cholinesterase inhibitors than alkylimine analogs. Optimal inhibition for AChE was exhibited by cyclohexaquinoline and for BChE by cycloheptaquinoline. Increasing spacer length elevated the potency against both cholinesterases. Structure-activity relationships agreed with docking results. Mixed-type reversible AChE inhibition, dual docking to catalytic and peripheral anionic sites, and propidium iodide displacement suggested the potential of hybrids to block AChE-induced β-amyloid (Aβ) aggregation. Hybrids also exhibited the inhibition of Aβ self-aggregation in the thioflavin test; those with a hexaquinoline ring and C8 spacer were the most active. Conjugates demonstrated high antioxidant activity in ABTS and FRAP assays as well as the inhibition of luminol chemiluminescence and lipid peroxidation in mouse brain homogenates. Quantum-chemical calculations explained antioxidant results. Computed ADMET profiles indicated favorable blood-brain barrier permeability, suggesting the CNS activity potential. Thus, the conjugates could be considered promising multifunctional agents for the potential treatment of AD.

Wide-scope targeted analysis of bioactive lipids in human plasma by LC/MS/MS

Quantitative information on blood metabolites can be used in developing advanced medical strategies such as early detection and prevention of disease. Monitoring bioactive lipids such as steroids, bile acids, and PUFA metabolites could be a valuable indicator of health status. However, a method for simultaneously measuring these bioactive lipids has not yet been developed. Here, we report a LC/MS/MS method that can simultaneously measure 144 bioactive lipids, including steroids, bile acids, and PUFA metabolites, from human plasma, and a sample preparation method for these targets. Protein removal by methanol precipitation and purification of bioactive lipids by solid-phase extraction improved the recovery of the targeted compounds in human plasma samples, demonstrating the importance of sample preparation methods for a wide range of bioactive lipid analyses. Using the developed method, we studied the plasma from healthy human volunteers and confirmed the presence of bioactive lipid molecules associated with sex differences and circadian rhythms. The developed method of bioactive lipid analysis can be applied to health monitoring and disease biomarker discovery in precision medicine.

Finding structural requirements of structurally diverse α-glucosidase and α-amylase inhibitors through validated and predictive 2D-QSAR and 3D-QSAR analyses

Diabetes mellitus (DM) is a chronic metabolic disorder characterized by hyperglycemic state. The α-glucosidase and α-amylase are considered two major targets for the management of Type 2 DM due to their ability of metabolizing carbohydrates into simpler sugars. In the current study, cheminformatics analyses were performed to develop validated and predictive models with a dataset of 187 α-glucosidase and α-amylase dual inhibitors. Separate linear, interpretable and statistically robust 2D-QSAR models were constructed with datasets containing the activities of α-glucosidase and α-amylase inhibitors with an aim to explain the crucial structural and physicochemical attributes responsible for higher activity towards these targets. Consequently, some descriptors of the models pointed out the importance of specific structural moieties responsible for the higher activities for these targets and on the other hand, properties such as ionization potential and mass of the compounds as well as number of hydrogen bond donors in molecules were found to be crucial in determining the binding potentials of the dataset compounds. Statistically significant 3D-QSAR models were developed with both α-glucosidase and α-amylase inhibition datapoints to estimate the importance of 3D electrostatic and steric fields for improved potentials towards these two targets. Molecular docking performed with selected compounds with homology model of α-glucosidase and X-ray crystal structure of α-amylase largely supported the interpretations obtained from the cheminformatic analyses. The current investigation should serve as important guidelines for the design of future α-glucosidase and α-amylase inhibitors. Besides, the current investigation is entirely performed by using non-commercial open-access tools to ensure easy accessibility and reproducibility of the investigation which may help researchers throughout the world to work more on drug design and discovery.

Mining Chemogenomic Spaces for Prediction of Drug-Target Interactions

The pipeline of drug discovery consists of a number of processes; drug-target interaction determination is one of the salient steps among them. Computational prediction of drug-target interactions can facilitate in reducing the search space of experimental wet lab-based verifications steps, thus considerably reducing time and other resources dedicated to the drug discovery pipeline. While machine learning-based methods are more widespread for drug-target interaction prediction, network-centric methods are also evolving. In this chapter, we focus on the process of the drug-target interaction prediction from the perspective of using machine learning algorithms and the various stages involved for developing an accurate predictor.

GUIDEMOL: A Python graphical user interface for molecular descriptors based on RDKit

GUIDEMOL is a Python computer program based on the RDKit software to process molecular structures and calculate molecular descriptors with a graphical user interface using the tkinter package. It can calculate descriptors already implemented in RDKit as well as grid representations of 3D molecular structures using the electrostatic potential or voxels. The GUIDEMOL app provides easy access to RDKit tools for chemoinformatics users with no programming skills and can be adapted to calculate other descriptors or to trigger other procedures. A command line interface (CLI) is also provided for the calculation of grid representations. The source code is available at https://github.com/jairesdesousa/guidemol.

Autonomous, multiproperty-driven molecular discovery: From predictions to measurements and back

A closed-loop, autonomous molecular discovery platform driven by integrated machine learning tools was developed to accelerate the design of molecules with desired properties. We demonstrated two case studies on dye-like molecules, targeting absorption wavelength, lipophilicity, and photooxidative stability. In the first study, the platform experimentally realized 294 unreported molecules across three automatic iterations of molecular design-make-test-analyze cycles while exploring the structure-function space of four rarely reported scaffolds. In each iteration, the property prediction models that guided exploration learned the structure-property space of diverse scaffold derivatives, which were realized with multistep syntheses and a variety of reactions. The second study exploited property models trained on the explored chemical space and previously reported molecules to discover nine top-performing molecules within a lightly explored structure-property space.

Bioactive Molecules from the Innate Immunity of Ascidians and Innovative Methods of Drug Discovery: A Computational Approach Based on Artificial Intelligence

The study of bioactive molecules of marine origin has created an important bridge between biological knowledge and its applications in biotechnology and biomedicine. Current studies in different research fields, such as biomedicine, aim to discover marine molecules characterized by biological activities that can be used to produce potential drugs for human use. In recent decades, increasing attention has been paid to a particular group of marine invertebrates, the Ascidians, as they are a source of bioactive products. We describe omics data and computational methods relevant to identifying the mechanisms and processes of innate immunity underlying the biosynthesis of bioactive molecules, focusing on innovative computational approaches based on Artificial Intelligence. Since there is increasing attention on finding new solutions for a sustainable supply of bioactive compounds, we propose that a possible improvement in the biodiscovery pipeline might also come from the study and utilization of marine invertebrates' innate immunity.

Machine Learning Enables Accurate Prediction of Quinone Formation during Drug Metabolism

Metabolism helps in the elimination of drugs from the human body by making them more hydrophilic. Sometimes, drugs can be bioactivated to highly reactive metabolites or intermediates during metabolism. These reactive metabolites are often responsible for the toxicities associated with the drugs. Identification of reactive metabolites of drug candidates can be very helpful in the initial stages of drug discovery. Quinones are soft electrophiles that are generated as reactive intermediates during metabolism. Quinones make up more than 40% of the reactive metabolites. In this work, a reliable data set of 510 molecules was used to develop machine learning and deep learning-based predictive models to predict the formation of quinone-type metabolites. For representing molecules, two-dimensional (2D) descriptors, PubChem fingerprints, electro-topological state (E-state) fingerprints, and metabolic reactivity-based descriptors were used. Developed models were compared to the existing Xenosite web server using the untouched test set of 102 molecules. The best model achieved an accuracy of 86.27%, while the Xenosite server could achieve an accuracy of only 52.94% on the test set. Descriptor analysis revealed that the presence of greater numbers of polar moieties in a molecule can prevent the formation of quinone-type metabolites. In addition, the presence of a nitrogen atom in an aromatic ring and the presence of metabolophores V51, V52, and V53 (SMARTCyp descriptors) decrease the probability of quinone formation. Finally, a tool based on the best machine learning models was developed, which is accessible at http://14.139.57.41/quinonepred/.

Predicting absolute aqueous solubility by applying a machine learning model for an artificially liquid-state as proxy for the solid-state

In this study, we use machine learning algorithms with QM-derived COSMO-RS descriptors, along with Morgan fingerprints, to predict the absolute solubility of drug-like compounds. The QM-derived descriptors account for the molecular properties of the solute, i.e., the solute-solute interactions in an artificial-liquid-state (super-cooled liquid), and the solute-solvent interactions in solution. We employ two main approaches to predict solubility: (i) a hypothetical pathway that involves melting the solute at room temperature T = T¯ ([Formula: see text]) and mixing the artificially liquid solute into the solvent ([Formula: see text]). In this approach [Formula: see text] is predicted using machine learning models, and the [Formula: see text] is obtained from COSMO-RS calculations; (ii) direct solubility prediction using machine learning algorithms. The models were trained on a large number of Bayer in-house compounds for which water solubility data is available at physiological pH of 6.5 and ambient temperature. We also evaluated our models using external datasets from a solubility challenge. Our models present great improvements compared to the absolute solubility prediction with the QSAR model for the artificial liquid state as implemented in the COSMOtherm software, for both in-house and external datasets. We are furthermore able to demonstrate the superiority of QM-derived descriptors compared to cheminformatics descriptors. We finally present low-cost alternative models using fragment-based COSMOquick calculations with only marginal reduction in the quality of predicted solubility.

In Silico, in Vitro and in Vivo Study of Substituted Imidazolidinone Sulfonamides as Antibacterial Agents

New substituted imidazolidinone sulfonamides have been developed using a rational drug design strategy. Predictive QSAR models for the search of new antibacterials were created using the OCHEM platform. Regression models were applied to verify a virtual chemical library of new imidazolidinone derivatives designed to have antibacterial activity. A number of substituted imidazolidinone sulfonamides as effective antibacterial agents were identified by QSAR prediction, synthesized and characterized by spectral and elemental, and tested in vitro. Six studied compounds have shown the highest in vitro antibacterial activity against Gram-negative E. coli and Gram-positive S. aureus multidrug-resistant strains. The in vivo acute toxicity of these imidazolidinone sulfonamides based on the LC50 value ranged from 16.01 to 44.35 mg/L (slightly toxic compounds class). The results of molecular docking suggest that the antibacterial mechanism of the compounds can be associated with the inhibition of post-translational modification processes of bacterial peptides and proteins.

The Trends and Future Prospective of In Silico Models from the Viewpoint of ADME Evaluation in Drug Discovery

Drug discovery and development are aimed at identifying new chemical molecular entities (NCEs) with desirable pharmacokinetic profiles for high therapeutic efficacy. The plasma concentrations of NCEs are a biomarker of their efficacy and are governed by pharmacokinetic processes such as absorption, distribution, metabolism, and excretion (ADME). Poor ADME properties of NCEs are a major cause of attrition in drug development. ADME screening is used to identify and optimize lead compounds in the drug discovery process. Computational models predicting ADME properties have been developed with evolving model-building technologies from a simplified relationship between ADME endpoints and physicochemical properties to machine learning, including support vector machines, random forests, and convolution neural networks. Recently, in the field of in silico ADME research, there has been a shift toward evaluating the in vivo parameters or plasma concentrations of NCEs instead of using predictive results to guide chemical structure design. Another research hotspot is the establishment of a computational prediction platform to strengthen academic drug discovery. Bioinformatics projects have produced a series of in silico ADME models using free software and open-access databases. In this review, we introduce prediction models for various ADME parameters and discuss the currently available academic drug discovery platforms.

5-Nitroisoxazoles in SNAr Reactions: A Novel Chemo- and Regioselective Approach to Isoxazole-Based Bivalent Ligands of AMPA Receptors

An efficient regioselective approach to novel functionalized bis(isoxazoles) with a variety of aromatic and aliphatic linkers was elaborated, based on the heterocyclization reaction of electrophilic alkenes under the treatment with tetranitromethane-triethylamine complex affording 3-EWG-5-nitroisoxazoles. The subsequent SNAr reactions of 5-nitroisoxazoles with various O,O-, N,N- and S,S-bis(nucleophiles) provide a wide range of bis(isoxazole) derivatives in good isolated yields. Employing an elaborated method, a series of novel bis(3-EWG-isoxazoles) as the promising allosteric modulators of AMPA receptors were designed and synthesized. The effect of the compounds on the kainate-induced currents was studied in the patch clamp experiments, revealing modulator properties for several of them. The best positive modulator potency was found for dimethyl 5,5'-(ethane-1,2-diylbis(sulfanediyl))bis(isoxazole-3-carboxylate), which potentiated the kainate-induced currents in a wide concentration range (10-12-10-6 M) with maximum potentiation of 77% at 10-10 M. The results were rationalized using molecular docking and molecular dynamics simulations of modulator complexes with the dimeric ligand-binding domain of the GluA2 AMPA receptor. The predicted physicochemical, ADMET, and PAINS properties confirmed that the AMPA receptor modulators based on the bis(isoxazole) scaffold may serve as potential lead compounds for the development of neuroprotective drugs.

Recent Studies of Artificial Intelligence on In Silico Drug Absorption

Absorption is an important area of research in pharmacochemistry and drug development, because the drug has to be absorbed before any drug effects can occur. Furthermore, the ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) profile of drugs can be directly and considerably altered by modulating factors affecting absorption. Many drugs in development fail because of poor absorption. The research and continuous efforts of researchers in recent years have brought many successes and promises in drug absorption property prediction, especially in silico, which helps to reduce the time and cost significantly for screening undesirable drug candidates. In this report, we explicitly provide an overview of recent in silico studies on predicting absorption properties, especially from 2019 to the present, using artificial intelligence. Additionally, we have collected and investigated public databases that support absorption prediction research. On those grounds, we also proposed the challenges and development directions of absorption prediction in the future. We hope this review can provide researchers with valuable guidelines on absorption prediction to facilitate the development of newer approaches in drug discovery.

Modeling and insights into the structural characteristics of endocrine-disrupting chemicals

Endocrine-disrupting chemicals (EDCs) can cause serious harm to human health and the environment; therefore, it is important to rapidly and correctly identify EDCs. Different computational models have been proposed for the prediction of EDCs over the past few decades, but the reported models are not always easily available, and few studies have investigated the structural characteristics of EDCs. In the present study, we have developed a series of artificial intelligence models targeting EDC receptors: the androgen receptor (AR); estrogen receptor (ER); and pregnane X receptor (PXR). The consensus models achieved good predictive results for validation sets with balanced accuracy values of 87.37%, 90.13%, and 79.21% for AR, ER, and PXR binding assays, respectively. Analysis of the physical-chemical properties suggested that several chemical properties were significantly (p < 0.05) different between EDCs and non-EDCs. We also identified structural alerts that can indicate an EDC, which were integrated into the web server SApredictor. These models and structural characteristics can provide useful tools and information in the discrimination and mechanistic understanding of EDCs in drug discovery and environmental risk assessment.

Modeling and Insights into the Structural Characteristics of Chemical Mitochondrial Toxicity

Mitochondria are the energy metabolism center of cells and are involved in a number of other processes, such as cell differentiation and apoptosis, signal transduction, and regulation of cell cycle and cell proliferation. It is of great significance to evaluate the mitochondrial toxicity of drugs and other chemicals. In the present study, we aimed to propose easily available artificial intelligence (AI) models for the prediction of chemical mitochondrial toxicity and investigate the structural characteristics with the analysis of molecular properties and structural alerts. The consensus model achieved good predictive results with high total accuracy at 87.21% for validation sets. The models can be accessed freely via https://ochem.eu/article/158582. Besides, several commonly used chemical properties were significantly different between chemicals with and without mitochondrial toxicity. We also detected the structural alerts (SAs) responsible for mitochondrial toxicity and integrated them into the web-server SApredictor (www.sapredictor.cn). The study may provide useful tools for in silico estimation of mitochondrial toxicity and be helpful to understand the mechanisms of mitochondrial toxicity.

LogD7.4 prediction enhanced by transferring knowledge from chromatographic retention time, microscopic pKa and logP

Lipophilicity is a fundamental physical property that significantly affects various aspects of drug behavior, including solubility, permeability, metabolism, distribution, protein binding, and toxicity. Accurate prediction of lipophilicity, measured by the logD7.4 value (the distribution coefficient between n-octanol and buffer at physiological pH 7.4), is crucial for successful drug discovery and design. However, the limited availability of data for logD modeling poses a significant challenge to achieving satisfactory generalization capability. To address this challenge, we have developed a novel logD7.4 prediction model called RTlogD, which leverages knowledge from multiple sources. RTlogD combines pre-training on a chromatographic retention time (RT) dataset since the RT is influenced by lipophilicity. Additionally, microscopic pKa values are incorporated as atomic features, providing valuable insights into ionizable sites and ionization capacity. Furthermore, logP is integrated as an auxiliary task within a multitask learning framework. We conducted ablation studies and presented a detailed analysis, showcasing the effectiveness and interpretability of RT, pKa, and logP in the RTlogD model. Notably, our RTlogD model demonstrated superior performance compared to commonly used algorithms and prediction tools. These results underscore the potential of the RTlogD model to improve the accuracy and generalization of logD prediction in drug discovery and design. In summary, the RTlogD model addresses the challenge of limited data availability in logD modeling by leveraging knowledge from RT, microscopic pKa, and logP. Incorporating these factors enhances the predictive capabilities of our model, and it holds promise for real-world applications in drug discovery and design scenarios.

Quinoline Hydrazone Derivatives as New Antibacterials against Multidrug Resistant Strains

To develop novel antimicrobial agents a series of 2(4)-hydrazone derivatives of quinoline were designed, synthesized and tested. QSAR models of the antibacterial activity of quinoline derivatives were developed by the OCHEM web platform using different machine learning methods. A virtual set of quinoline derivatives was verified with a previously published classification model of anti-E. coli activity and screened using the regression model of anti-S. aureus activity. Selected and synthesized 2(4)-hydrazone derivatives of quinoline exhibited antibacterial activity against the standard and antibiotic-resistant S. aureus and E. coli strains in the range from 15 to 30 mm by the diameter of growth inhibition zones. Molecular docking showed the complex formation of the studied compounds into the catalytic domain of dihydrofolate reductase with an estimated binding affinity from -8.4 to -9.4 kcal/mol.