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

Elucidating molecular mechanism and chemical space of chalcones through biological networks and machine learning approaches

We developed a bio-cheminformatics method, exploring disease inhibition mechanisms using machine learning-enhanced quantitative structure-activity relationship (ML-QSAR) models and knowledge-driven neural networks. ML-QSAR models were developed using molecular fingerprint descriptors and the Random Forest algorithm to explore the chemical spaces of Chalcones inhibitors against diverse disease properties, including antifungal, anti-inflammatory, anticancer, antimicrobial, and antiviral effects. We generated and validated robust machine learning-based bioactivity prediction models (https://github.com/RatulChemoinformatics/QSAR) for the top genes. These models underwent ROC and applicability domain analysis, followed by molecular docking studies to elucidate the molecular mechanisms of the molecules. Through comprehensive neural network analysis, crucial genes such as AKT1, HSP90AA1, SRC, and STAT3 were identified. The PubChem fingerprint-based model revealed key descriptors: PubchemFP521 for AKT1, PubchemFP180 for SRC, PubchemFP633 for HSP90AA1, and PubchemFP145 and PubchemFP338 for STAT3, consistently contributing to bioactivity across targets. Notably, chalcone derivatives demonstrated significant bioactivity against target genes, with compound RA1 displaying a predictive pIC50 value of 5.76 against HSP90AA1 and strong binding affinities across other targets. Compounds RA5 to RA7 also exhibited high binding affinity scores comparable to or exceeding existing drugs. These findings emphasize the importance of knowledge-based neural network-based research for developing effective drugs against diverse disease properties. These interactions warrant further in vitro and in vivo investigations to elucidate their potential in rational drug design. The presented models provide valuable insights for inhibitor design and hold promise for drug development. Future research will prioritize investigating these molecules for mycobacterium tuberculosis, enhancing the comprehension of effectiveness in addressing infectious diseases.

Synthesis of thiazolidine-2,4-dione tethered 1,2,3-triazoles as α-amylase inhibitors: In vitro approach coupled with QSAR, molecular docking, molecular dynamics and ADMET studies

A new series of thiazolidine-2,4-dione tethered 1,2,3-triazole derivatives were designed, synthesized and screened for their α-amylase inhibitory potential employing in vitro and in silico approaches. The target compounds were synthesized with the help of Cu (I) catalyzed [3 + 2] cycloaddition of terminal alkyne with numerous azides, followed by unambiguously characterizing the structure by employing various spectroscopic approaches. The synthesized derivatives were assessed for their in vitro α-amylase inhibition and it was found that thiazolidine-2,4-dione derivatives 6e, 6j, 6o, 6u and 6x exhibited comparable inhibition with the standard drug acarbose. The compound 6e with a 7-chloroquinolinyl substituent on the triazole ring exhibited significant inhibition potential with IC50 value of 0.040 μmol mL-1 whereas compound 6c (IC50 = 0.099 μmol mL-1) and 6h (IC50 = 0.098 μmol mL-1) were poor inhibitors. QSAR studies revealed the positively correlating descriptors that aid in the design of novel compounds. Molecular docking was performed to investigate the binding interactions with the active site of the biological receptor and the stability of the complex over a period of 100 ns was examined using molecular dynamics studies. The physiochemical properties and drug-likeliness behavior of the potent derivatives were investigated by carrying out the ADMET studies.

Safer and greener chemicals for the aquatic ecosystem: Chemometric modeling of the prolonged and chronic aquatic toxicity of chemicals on Oryzias latipes

In the modern era, chemicals and their products have been used everywhere like agriculture, healthcare, food, cosmetics, pharmaceuticals, household products, clothing industry, etc. These chemicals find their way to reach the aquatic ecosystem (directly/indirectly) and cause severe chronic and prolonged toxic effects to aquatic species which is also then translated to human beings. Prolonged and chronic toxicity data of many chemicals that are used daily is not available due to high experimentation testing costs, time investment, and the requirement of a large number of animal sacrifices. Thus, in silico approaches (e.g., QSAR (quantitative structure-activity relationship)) are the best alternative for chronic and prolonged toxicity predictions. The present work offers multi-endpoint (five endpoints: chronic_LOEC, prolonged_14D_LC50, prolonged_14D_NOEC, prolonged_21D_LC50, prolonged_21D_NOEC) QSAR models for addressing the prolonged and chronic aquatic toxicity of chemicals toward fish (O. latipes). The statistical results (R2 =0.738-0.869, QLOO2 =0.712-0.831, Q(F1)2 =0.618-0.731) of the developed models show that they were robust, reliable, reproducible, accurate, and predictive. Some of the features that are responsible for prolonged and chronic toxicity of chemicals towards O. latipes are as follows: the presence of substituted benzene, hydrophobicity, unsaturation, electronegativity, the presence of long-chain fragments, the presence of a greater number of atoms at conjugation, and the presence of halogen atoms. On the other hand, hydrophilicity and graph density descriptors retard the aquatic chronic and prolonged toxicity of chemicals toward O. latipes. The PPDB (pesticide properties database) and experimental and investigational classes of drugs from the DrugBank database were also screened using the developed model. Thus, these multi-endpoint models will be helpful for data-gap filling and provide a broad range of applicability. Therefore, this research will aid in the in silico QSAR (quantitative structure-activity relationship) prediction (non-animal testing) of the prolonged and chronic toxicity of untested and new toxic chemicals/drugs/pesticides, design and development of eco-friendly, novel, and safer chemicals, and help to protect the aquatic ecosystem from exposure to toxic and hazardous chemicals.

Machine Learning-Based Approach to Identify Inhibitors of Sterol-14-Alpha Demethylase: A Study on Chagas Disease

Objectives

Chagas Disease, caused by the parasite Trypanosoma cruzi, remains a significant public health concern, particularly in Latin America. The current standard treatment for Chagas Disease, benznidazole, is associated with various side effects, necessitating the search for alternative therapeutic options. In this study, we aimed to identify potential therapeutics for Chagas Disease through a comprehensive computational analysis.

Methods

A library of compounds derived from Cananga odorata was screened using a combination of pharmacophore modeling, structure-based screening, and quantitative structure-activity relationship (QSAR) analysis. The pharmacophore model facilitated the efficient screening of the compound library, while the structure-based screening identified hit compounds with promising inhibitory potential against the target enzyme, sterol-14-alpha demethylase.

Results

The QSAR model predicted the bioactivity of the hit compounds, revealing one compound to exhibit superior activity compared to benznidazole. Evaluation of the physicochemical, pharmacokinetic, toxicity, and medicinal chemistry properties of the hit compounds indicated their drug-like characteristics, oral bioavailability, ease of synthesis, and reduced toxicity profiles.

Conclusion

Overall, our findings present a promising avenue for the discovery of novel therapeutics for Chagas Disease. The identified hit compounds possess favorable drug-like properties and demonstrate potent inhibitory effects against the target enzyme. Further in vitro and in vivo studies are warranted to validate their efficacy and safety profiles.

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.

Unveiling potent inhibitors for schistosomiasis through ligand-based drug design, molecular docking, molecular dynamics simulations and pharmacokinetics predictions

Schistosomiasis is a neglected tropical disease which imposes a considerable and enduring impact on affected regions, leading to persistent morbidity, hindering child development, diminishing productivity, and imposing economic burdens. Due to the emergence of drug resistance and limited management options, there is need to develop additional effective inhibitors for schistosomiasis. In view of this, quantitative structure-activity relationship studies, molecular docking, molecular dynamics simulations, drug-likeness and pharmacokinetics predictions were applied to 39 Schistosoma mansoni Thioredoxin Glutathione Reductase (SmTGR) inhibitors. The chosen QSAR model demonstrated robust statistical parameters, including an R2 of 0.798, R2adj of 0.767, Q2cv of 0.681, LOF of 0.930, R2test of 0.776, and cR2p of 0.746, confirming its reliability. The most active derivative (compound 40) was identified as a lead candidate for the development of new potential non-covalent inhibitors through ligand-based design. Subsequently, 12 novel compounds (40a-40l) were designed with enhanced anti-schistosomiasis activity and binding affinity. Molecular docking studies revealed strong and stable interactions, including hydrogen bonding, between the designed compounds and the target receptor. Molecular dynamics simulations over 100 nanoseconds and MM-PBSA free binding energy (ΔGbind) calculations validated the stability of the two best-designed molecules. Furthermore, drug-likeness and pharmacokinetics prediction analyses affirmed the potential of these designed compounds, suggesting their promise as innovative agents for the treatment of schistosomiasis.

Unveiling critical structural features for effective HDAC8 inhibition: a comprehensive study using quantitative read-across structure-activity relationship (q-RASAR) and pharmacophore modeling

Histone deacetylases constitute a group of enzymes that participate in several biological processes. Notably, inhibiting HDAC8 has become a therapeutic strategy for various diseases. The current inhibitors for HDAC8 lack selectivity and target multiple HDACs. Consequently, there is a growing recognition of the need for selective HDAC8 inhibitors to enhance the effectiveness of therapeutic interventions. In our current study, we have utilized a multi-faceted approach, including Quantitative Structure-Activity Relationship (QSAR) combined with Quantitative Read-Across Structure-Activity Relationship (q-RASAR) modeling, pharmacophore mapping, molecular docking, and molecular dynamics (MD) simulations. The developed q-RASAR model has a high statistical significance and predictive ability (Q2F1:0.778, Q2F2:0.775). The contributions of important descriptors are discussed in detail to gain insight into the crucial structural features in HDAC8 inhibition. The best pharmacophore hypothesis exhibits a high regression coefficient (0.969) and a low root mean square deviation (0.944), highlighting the importance of correctly orienting hydrogen bond acceptor (HBA), ring aromatic (RA), and zinc-binding group (ZBG) features in designing potent HDAC8 inhibitors. To confirm the results of q-RASAR and pharmacophore mapping, molecular docking analysis of the five potent compounds (44, 54, 82, 102, and 118) was performed to gain further insights into these structural features crucial for interaction with the HDAC8 enzyme. Lastly, MD simulation studies of the most active compound (54, mapped correctly with the pharmacophore hypothesis) and the least active compound (34, mapped poorly with the pharmacophore hypothesis) were carried out to validate the observations of the studies above. This study not only refines our understanding of essential structural features for HDAC8 inhibition but also provides a robust framework for the rational design of novel selective HDAC8 inhibitors which may offer insights to medicinal chemists and researchers engaged in the development of HDAC8-targeted therapeutics.

Toxicological evaluation of microbial secondary metabolites in the context of European active substance approval for plant protection products

Risk assessment (RA) of microbial secondary metabolites (SM) is part of the EU approval process for microbial active substances (AS) used in plant protection products (PPP). As the number of potentially produced microbial SM may be high for a certain microbial strain and existing information on the metabolites often are low, data gaps are frequently identified during the RA. Often, RA cannot conclusively clarify the toxicological relevance of the individual substances. This work presents data and RA conclusions on four metabolites, Beauvericin, 2,3-deepoxy-2,3-didehydro-rhizoxin (DDR), Leucinostatin A and Swainsonin in detail as examples for the challenging process of RA. To overcome the problem of incomplete assessment reports, RA of microbial AS for PPP is in need of new approaches. In view of the Next Generation Risk Assessment (NGRA), the combination of literature data, omic-methods, in vitro and in silico methods combined in adverse outcome pathways (AOPs) can be used for an efficient and targeted identification and assessment of metabolites of concern (MoC).

The first report on the assessment of maximum acceptable daily intake (MADI) of pesticides for humans using intelligent consensus predictions

Direct or indirect consumption of pesticides and their related products by humans and other living organisms without safe dosing may pose a health risk. The risk may arise after a short/long time which depends on the nature and amount of chemicals consumed. Therefore, the maximum acceptable daily intake of chemicals must be calculated to prevent these risks. In the present work, regression-based quantitative structure-activity relationship (QSAR) models were developed using 39 pesticides with maximum acceptable daily intake (MADI) for humans as the endpoint. From the statistical results (R2 = 0.674-0.712, QLOO2 = 0.553-0.580, Q(F1)2 = 0.544-0.611, and Q(F2)2 = 0.531-0.599), it can be inferred that the developed models were robust, reliable, reproducible, accurate, and predictive. Intelligent Consensus Prediction (ICP) was employed to improve the external predictivity (Q(F1)2 =0.579-0.657 and Q(F2)2 = 0.563-0.647) of the models. Some of the chemical markers responsible for toxicity enhancement are the presence of unsaturated bonds, lipophilicity, presence of C< (double bond-single bond-single bonded carbon), and the presence of sulphur and phosphate bonds at the topological distances 1 and 6, while the presence of hydrophilic groups and short chain fragments reduces the toxicity. The Pesticide Properties Database (PPDB) (1694 pesticides) was also screened with the developed models. Hence, this research work will be helpful for the toxicity assessment of pesticides before their synthesis, the development of eco-friendly and safer pesticides, and data-gap filling reducing the time, cost, and animal experimentation. Thus, this study might hold promise for future potential MADI assessment of pesticides and provide a meaningful contribution to the field of risk assessment.

Integrating Mechanistic and Toxicokinetic Information in Predictive Models of Cholestasis

Drug development involves the thorough assessment of the candidate's safety and efficacy. In silico toxicology (IST) methods can contribute to the assessment, complementing in vitro and in vivo experimental methods, since they have many advantages in terms of cost and time. Also, they are less demanding concerning the requirements of product and experimental animals. One of these methods, Quantitative Structure-Activity Relationships (QSAR), has been proven successful in predicting simple toxicity end points but has more difficulties in predicting end points involving more complex phenomena. We hypothesize that QSAR models can produce better predictions of these end points by combining multiple QSAR models describing simpler biological phenomena and incorporating pharmacokinetic (PK) information, using quantitative in vitro to in vivo extrapolation (QIVIVE) models. In this study, we applied our methodology to the prediction of cholestasis and compared it with direct QSAR models. Our results show a clear increase in sensitivity. The predictive quality of the models was further assessed to mimic realistic conditions where the query compounds show low similarity with the training series. Again, our methodology shows clear advantages over direct QSAR models in these situations. We conclude that the proposed methodology could improve existing methodologies and could be suitable for being applied to other toxicity end points.

Machine learning driven web-based app platform for the discovery of monoamine oxidase B inhibitors

Monoamine oxidases (MAOs), specifically MAO-A and MAO-B, play important roles in the breakdown of monoamine neurotransmitters. Therefore, MAO inhibitors are crucial for treating various neurodegenerative disorders, including Parkinson's disease (PD), Alzheimer's disease (AD), and amyotrophic lateral sclerosis (ALS). In this study, we developed a novel cheminformatics pipeline by generating three diverse molecular feature-based machine learning-assisted quantitative structural activity relationship (ML-QSAR) models concerning MAO-B inhibition. PubChem fingerprints, substructure fingerprints, and one-dimensional (1D) and two-dimensional (2D) molecular descriptors were implemented to unravel the structural insights responsible for decoding the origin of MAO-B inhibition in 249 non-reductant molecules. Based on a random forest ML algorithm, the final PubChem fingerprint, substructure fingerprint, and 1D and 2D molecular descriptor prediction models demonstrated significant robustness, with correlation coefficients of 0.9863, 0.9796, and 0.9852, respectively. The significant features of each predictive model responsible for MAO-B inhibition were extracted using a comprehensive variance importance plot (VIP) and correlation matrix analysis. The final predictive models were further developed as a web application, MAO-B-pred ( https://mao-b-pred.streamlit.app/ ), to allow users to predict the bioactivity of molecules against MAO-B. Molecular docking and dynamics studies were conducted to gain insight into the atomic-level molecular interactions between the ligand-receptor complexes. These findings were compared with the structural features obtained from the ML-QSAR models, which supported the mechanistic understanding of the binding phenomena. The presented models have the potential to serve as tools for identifying crucial molecular characteristics for the rational design of MAO-B target inhibitors, which may be used to develop effective drugs for neurodegenerative disorders.

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.

Unveiling first report on in silico modeling of aquatic toxicity of organic chemicals to Labeo rohita (Rohu) employing QSAR and q-RASAR

Labeo rohita, a fish species within the Carp family, holds significant dietary and aquacultural importance in South Asian countries. However, the habitats of L. rohita often face exposure to various harmful pesticides and organic compounds originating from industrial and agricultural runoff. It is challenging to individually investigate the effects of each potentially harmful compound. In such cases, in silico techniques like Quantitative Structure-Activity Relationship (QSAR) and quantitative Read-Across Structure-Activity Relationship (q-RASAR) can be employed to construct algorithmic models capable of simultaneously assessing the toxicity of numerous compounds. We utilized the US EPA's ToxValDB database to curate data regarding acute median lethal concentration (LC50) toxicity for L. rohita. The experimental variables included study type (mortality), study duration (ranging from 0.25 h to 4 h), exposure route (static, flowthrough, and renewal), exposure method (drinking water), and types of chemicals (industrial chemicals and pharmaceuticals). Using this dataset, we developed regression-based QSAR and q-RASAR models to predict chemical toxicity to L. rohita based on chemical descriptors. The key descriptors for predicting the toxicity of L. rohita in the regression-based QSAR model include F05[S-Cl], SpMax_EA(ri), s4_relPathLength_2, and SpDiam_AEA(ed). These descriptors can be employed to estimate the toxicity of untested compounds and aid in the development of compounds with lower toxicity based on the presence or absence of these descriptors. Both the QSAR and q-RASAR models serve as valuable tools for understanding the chemicals' structural features responsible for toxicity and for filling gaps in aquatic toxicity data by predicting the toxicity of newly untested compounds in relation to L. rohita. Finally, the developed best model was employed to predict 297 external chemicals, the most toxic substances to L. rohita were identified as cyhalothrin, isobornyl thiocyanatoacetate, and paclobutrzol, while the least toxic ones included ethyl acetate, ethylthiourea, and n-butyric acid.

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.

Deep learning-based design and screening of benzimidazole-pyrazine derivatives as adenosine A2B receptor antagonists

The Adenosine A2B receptor (A2BAR) is considered a novel potential target for the immunotherapy of cancer, and A2BAR antagonists have an inhibitory effect on tumor growth, proliferation, and metastasis. In our previous studies, we identified a class of benzimidazole-pyrazine scaffolds whose derivatives exhibited the antagonistic effect but lacked subtype selectivity towards A2BAR. In this work, we developed a scaffold-based protocol that incorporates a deep generative model and multilayer virtual screening to design benzimidazole-pyrazine derivatives as potential selective A2BAR antagonists. By utilizing a generative model with reported A2BAR antagonists as the training set, we built up a scaffold-focused library of benzimidazole-pyrazine derivatives and processed a virtual screening protocol to discover potential A2BAR antagonists. Finally, five molecules with different Bemis-Murcko scaffolds were identified and exhibited higher binding free energies than the reference molecule 12o. Further computational analysis revealed that the 3-benzyl derivative ABA-1266 presented high selectivity toward A2BAR and showed preferred draggability, providing future potent development of selective A2BAR antagonists.Communicated by Ramaswamy H. Sarma.

A battery of in silico models application for pesticides exerting reproductive health effects: Assessment of performance and prioritization of mechanistic studies

Given the high attention to endocrine disrupting chemicals (EDC), there is an urgent need for the development of rapid and reliable approaches for the screening of large numbers of chemicals with respect to their endocrine disruption potential. This study aimed at the assessment of the correlation between the predicted results of a battery of in silico tools and the reported observed adverse effects from in vivo reproductive toxicity studies. We used VirtualToxLab (VTL) software and the EndocrineDisruptome (ED) online tool to evaluate the binding affinities to nuclear receptors of 17 pesticides, 7 of which were classified as reprotoxic substances under Regulation (EC) No 1272/2008 on the classification, labelling and packaging of substances and mixtures (CLP). Then, we aligned the results of the in silico modelling with data from ToxCast assays and in vivo reproductive toxicity studies. We combined results from different in silico tools in two different ways to improve the characteristics of their predictive performance. Reproductive toxicity can be caused by various mechanisms; however, in this study, we demonstrated that the use of a battery of in silico tools for assessing the binding to nuclear receptors can be useful for identifying hazardous compounds and for prioritizing further studies.

Combination of QSAR Modeling and Hybrid-Based Consensus Scoring to Identify Dual-Targeting Inhibitors of PLK1 and p38γ

Polo-like kinase 1 (PLK1) and p38γ mitogen-activated protein kinase (p38γ) play important roles in cancer pathogenesis by controlling cell cycle progression and are therefore attractive cancer targets. The design of multitarget inhibitors may offer synergistic inhibition of distinct targets and reduce the risk of drug-drug interactions to improve the balance between therapeutic efficacy and safety. We combined deep-learning-based quantitative structure-activity relationship (QSAR) modeling and hybrid-based consensus scoring to screen for inhibitors with potential activity against the targeted proteins. Using this combination strategy, we identified a potent PLK1 inhibitor (compound 4) that inhibited PLK1 activity and liver cancer cell growth in the nanomolar range. Next, we deployed both our QSAR models for PLK1 and p38γ on the Enamine compound library to identify dual-targeting inhibitors against PLK1 and p38γ. Likewise, the identified hits were subsequently subjected to hybrid-based consensus scoring. Using this method, we identified a promising compound (compound 14) that could inhibit both PLK1 and p38γ activities. At nanomolar concentrations, compound 14 inhibited the growth of human hepatocellular carcinoma and hepatoblastoma cells in vitro. This study demonstrates the combined screening strategy to identify novel potential inhibitors for existing targets.

Usage of model combination in computational toxicology

New Approach Methodologies (NAMs) have ushered in a new era in the field of toxicology, aiming to replace animal testing. However, despite these advancements, they are not exempt from the inherent complexities associated with the study's endpoint. In this review, we have identified three major groups of complexities: mechanistic, chemical space, and methodological. The mechanistic complexity arises from interconnected biological processes within a network that are challenging to model in a single step. In the second group, chemical space complexity exhibits significant dissimilarity between compounds in the training and test series. The third group encompasses algorithmic and molecular descriptor limitations and typical class imbalance problems. To address these complexities, this work provides a guide to the usage of a combination of predictive Quantitative Structure-Activity Relationship (QSAR) models, known as metamodels. This combination of low-level models (LLMs) enables a more precise approach to the problem by focusing on different sub-mechanisms or sub-processes. For mechanistic complexity, multiple Molecular Initiating Events (MIEs) or levels of information are combined to form a mechanistic-based metamodel. Regarding the complexity arising from chemical space, two types of approaches were reviewed to construct a fragment-based chemical space metamodel: those with and without structure sharing. Metamodels with structure sharing utilize unsupervised strategies to identify data patterns and build low-level models for each cluster, which are then combined. For situations without structure sharing due to pharmaceutical industry intellectual property, the use of prediction sharing, and federated learning approaches have been reviewed. Lastly, to tackle methodological complexity, various algorithms are combined to overcome their limitations, diverse descriptors are employed to enhance problem definition and balanced dataset combinations are used to address class imbalance issues (methodological-based metamodels). Remarkably, metamodels consistently outperformed classical QSAR models across all cases, highlighting the importance of alternatives to classical QSAR models when faced with such complexities.

A computational investigation of thymidylate synthase inhibitors through a combined approach of 3D-QSAR and pharmacophore modelling

Thymidylate synthase (TS) is a crucial target of cancer drug discovery and is mainly involved in the De novo synthesis of the DNA precursor thymine. In the present study, to generate reliable models and identify a few promising molecules, we combined QSAR modelling with the pharmacophore hypothesis-generating technique. Input molecules were clustered on their similarity, and a cluster of 74 molecules with a pyrimidine moiety was chosen as the set for 3D-QSAR and pharmacophore modelling. Atom-based and field-based 3D-QSAR models were generated and statistically validated with R2 > 0.90 and Q2 > 0.75. The common pharmacophore hypothesis(CPH) generation identified the best six-point model ADHRRR. Using these best models, a library of FDA-approved drugs was screened for activity and filtered via molecular docking, ADME profiling, and molecular dynamics simulations. The top ten promising TS-inhibiting candidates were identified, and their chemical features profitable for TS inhibitors were explored.Communicated by Ramaswamy H. Sarma.

The anti-SARS-CoV-2 activity of novel 9, 10-dihydrophenanthrene derivatives: an insight into molecular docking, ADMET analysis, and molecular dynamics simulation

Originating in Wuhan, the COVID-19 pandemic wave has had a profound impact on the global healthcare system. In this study, we used a 2D QSAR technique, ADMET analysis, molecular docking, and dynamic simulations to sort and evaluate the performance of thirty-nine bioactive analogues of 9,10-dihydrophenanthrene. The primary goal of the study is to use computational approaches to create a greater variety of structural references for the creation of more potent SARS-CoV-2 3Clpro inhibitors. This strategy is to speed up the process of finding active chemicals. Molecular descriptors were calculated using 'PaDEL' and 'ChemDes' software, and then redundant and non-significant descriptors were eliminated by a module in 'QSARINS ver. 2.2.2'. Subsequently, two statistically robust QSAR models were developed by applying multiple linear regression (MLR) methods. The correlation coefficients obtained by the two models are 0.89 and 0.82, respectively. These models were then subjected to internal and external validation tests, Y-randomization, and applicability domain analysis. The best model developed is applied to designate new molecules with good inhibitory activity values against severe acute respiratory syndrome coronavirus 2 (SARS CoV-2). We also examined various pharmacokinetic properties using ADMET analysis. Then, through molecular docking simulations, we used the crystal structure of the main protease of SARS-CoV-2 (3CLpro/Mpro) in a complex with the covalent inhibitor "Narlaprevir" (PDB ID: 7JYC). We also supported our molecular docking predictions with an extended molecular dynamics simulation of a docked ligand-protein complex. We hope that the results obtained in this study can be used as good anti-SARS-CoV-2 inhibitors.

Physicochemical Responsive Integrated Similarity Measure (PRISM) for a Comprehensive Quantitative Perspective of Sample Similarity Dynamically Assessed with NIR Spectra

Determining sample similarity underlies many foundational principles in analytical chemistry. For example, calibration models are unsuitable to predict outliers. Calibration transfer methods assume a moderate degree of sample and measurement dissimilarities between a calibration set and target prediction samples. Classification approaches link target sample similarities to groups of similar class samples. Although similarity is ubiquitous in analytical chemistry and everyday life, quantifying sample similarity is without a straightforward solution, especially when target domain samples are unlabeled and the only known features are measurable, such as spectra (the focus of this paper). The process proposed to assess sample similarity integrates spectral similarity information with contextual considerations among source analyte contents, model, and analyte predictions. This hybrid approach named the physicochemical responsive integrated similarity measure (PRISM) amplifies hidden-but-essential physicochemical properties encoded within respective spectra. PRISM is tested on four near-infrared (NIR) data sets for four diverse application areas to show efficacy. These applications are the assessment of prediction reliability and model updating for model generalizability, outlier detection, and basic matrix matching evaluation. Discussion is provided on adapting PRISM to classification problems. Results indicate that PRISM collects large amounts of similarity information and effectively integrates it to produce a quantitative similarity evaluation between the target sample and a source domain. The approach is also useful for biological samples with additional physiochemical variations. While PRISM is dynamically tested on NIR data, parts of PRISM were previously applied to other data types, and PRISM should be applicable to other measurement systems perturbed by matrix effects.

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

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

Exploring Dimensionality Reduction Techniques for Deep Learning Driven QSAR Models of Mutagenicity

Dimensionality reduction techniques are crucial for enabling deep learning driven quantitative structure-activity relationship (QSAR) models to navigate higher dimensional toxicological spaces, however the use of specific techniques is often arbitrary and poorly explored. Six dimensionality techniques (both linear and non-linear) were hence applied to a higher dimensionality mutagenicity dataset and compared in their ability to power a simple deep learning driven QSAR model, following grid searches for optimal hyperparameter values. It was found that comparatively simpler linear techniques, such as principal component analysis (PCA), were sufficient for enabling optimal QSAR model performances, which indicated that the original dataset was at least approximately linearly separable (in accordance with Cover's theorem). However certain non-linear techniques such as kernel PCA and autoencoders performed at closely comparable levels, while (especially in the case of autoencoders) being more widely applicable to potentially non-linearly separable datasets. Analysis of the chemical space, in terms of XLogP and molecular weight, uncovered that the vast majority of testing data occurred within the defined applicability domain, as well as that certain regions were measurably more problematic and antagonised performances. It was however indicated that certain dimensionality reduction techniques were able to facilitate uniquely beneficial navigations of the chemical space.

Conceptual DFT, machine learning and molecular docking as tools for predicting LD50 toxicity of organothiophosphates

CONTEXT

Several descriptors from conceptual density functional theory (cDFT) and the quantum theory of atoms in molecules (QTAIM) were utilized in Random Forest (RF), LASSO, Ridge, Elastic Net (EN), and Support Vector Machines (SVM) methods to predict the toxicity (LD₅₀) of sixty-two organothiophosphate compounds. The A-RF-G1 and A-RF-G2 models were obtained using the RF method, yielding statistically significant parameters with good performance, as indicated by R² values for the training set (R²_Train) and R² values for the test set (R²_Test), around 0.90.

METHODS

The molecular structure of all organothiophosphates was optimized via the range-separated hybrid functional ωB97XD with the 6-311 +  + G** basis set. Seven hundred and eighty-seven descriptors have been processed using a variety of machine learning algorithms: RF LASSO, Ridge, EN and SVM to generate a predictive model. The properties were obtained with Multiwfn, AIMALL and VMD programs. Docking simulations were performed by using AutoDock 4.2 and LigPlot + programs. All the calculations in this work are carried out in Gaussian 16 program package.

Development of a 2D-QSAR Model for Tissue-to-Plasma Partition Coefficient Value with High Accuracy Using Machine Learning Method, Minimum Required Experimental Values, and Physicochemical Descriptors

BACKGROUND

The demand for physiologically based pharmacokinetic (PBPK) model is increasing currently. New drug application (NDA) of many compounds is submitted with PBPK models for efficient drug development. Tissue-to-plasma partition coefficient (K_p) is a key parameter for the PBPK model to describe differential equations. However, it is difficult to obtain the K_p value experimentally because the measurement of drug concentration in the tissue is much harder than that in plasma.

OBJECTIVE

Instead of experiments, many researchers have sought in silico methods. Today, most of the models for K_p prediction are using in vitro and in vivo parameters as explanatory variables. We thought of physicochemical descriptors that could improve the predictability. Therefore, we aimed to develop the two-dimensional quantitative structure-activity relationship (2D-QSAR) model for K_p using physicochemical descriptors instead of in vivo experimental data as explanatory variables.

METHODS

We compared our model with the conventional models using 20-fold cross-validation according to the published method (Yun et al. J Pharmacokinet Pharmacodyn 41:1-14, 2014). We used random forest algorithm, which is known to be one of the best predictors for the 2D-QSAR model. Finally, we combined minimum in vitro experimental values and physiochemical descriptors. Thus, the prediction method for K_p value using a few in vitro parameters and physicochemical descriptors was developed; this is a multimodal model.

RESULTS

Its accuracy was found to be superior to that of the conventional models. Results of this research suggest that multimodality is useful for the 2D-QSAR model [RMSE and % of two-fold error: 0.66 and 42.2% (Berezohkovsky), 0.52 and 52.2% (Rodgers), 0.65 and 34.6% (Schmitt), 0.44 and 61.1% (published model), 0.41 and 62.1% (traditional model), 0.39 and 64.5% (multimodal model)].

CONCLUSION

We could develop a 2D-QSAR model for K_p value with the highest accuracy using a few in vitro experimental data and physicochemical descriptors.

Artificial intelligence in molecular de novo design: Integration with experiment

In this mini review, we capture the latest progress of applying artificial intelligence (AI) techniques based on deep learning architectures to molecular de novo design with a focus on integration with experimental validation. We will cover the progress and experimental validation of novel generative algorithms, the validation of QSAR models and how AI-based molecular de novo design is starting to become connected with chemistry automation. While progress has been made in the last few years, it is still early days. The experimental validations conducted thus far should be considered proof-of-principle, providing confidence that the field is moving in the right direction.

Modeling study for predicting altered cellular activity induced by nanomaterials based on Dlk1-Dio3 gene expression and structural relationships

Nanomaterials have been widely applied and developed due to its unique physicochemical characteristics, such as their small size. The environmental and biological effects caused by nanomaterials have raised concerns. In particular, some nanometal oxides have obvious biological toxicity and pose a major safety problem. The prediction model established by combining the expression levels of key genes with quantitative structure-activity relationship (QSAR) studies can predict the biotoxicity of nanomaterials by relying on both structural information and gene regulation information. This model can fill the gap of missing mechanisms in QSAR studies. In this study, we exposed A549 cells and BEAS-2B cells to 21 nanometal oxides for 24 h. Cell viability was assessed by measuring absorbance values using the CCK8 assay, and the expression levels of the Dlk1-Dio3 gene cluster were measured. By using the theoretical basis of the nano-QSAR model and the improved principles of the SMILES-based descriptors to combine specific gene expression and structural factors, new models were constructed using Monte Carlo partial least squares (MC-PLS) for the biotoxicity of the nanometal oxides on two different lung cells. The overall quality of the nano-QSAR models constructed by combining specific gene expression and structural parameters for A549 and BEAS-2B cells was better than that of the models constructed based on structural parameters only. The coefficient of determination (R²) of the A549 cell model increased from 0.9044 to 0.9969, and the Root Mean Square Error (RMSE) decreased from 0.1922 to 0.0348. The R² of the BEAS-2B cell model increased from 0.9355 to 0.9705, and the RMSE decreased from 0.1206 to 0.0874. The model validation proved the proposed models have a good prediction, generalization ability and model stability. This study offers a new research perspective for the toxicity assessment of nanometal oxides, contributing to a more systematic safety evaluation of nanomaterials.

A Nano-QSTR model to predict nano-cytotoxicity: an approach using human lung cells data

BACKGROUND

The widespread use of new engineered nanomaterials (ENMs) in industries such as cosmetics, electronics, and diagnostic nanodevices, has been revolutionizing our society. However, emerging studies suggest that ENMs present potentially toxic effects on the human lung. In this regard, we developed a machine learning (ML) nano-quantitative-structure-toxicity relationship (QSTR) model to predict the potential human lung nano-cytotoxicity induced by exposure to ENMs based on metal oxide nanoparticles.

RESULTS

Tree-based learning algorithms (e.g., decision tree (DT), random forest (RF), and extra-trees (ET)) were able to predict ENMs' cytotoxic risk in an efficient, robust, and interpretable way. The best-ranked ET nano-QSTR model showed excellent statistical performance with R² and Q²-based metrics of 0.95, 0.80, and 0.79 for training, internal validation, and external validation subsets, respectively. Several nano-descriptors linked to the core-type and surface coating reactivity properties were identified as the most relevant characteristics to predict human lung nano-cytotoxicity.

CONCLUSIONS

The proposed model suggests that a decrease in the ENMs diameter could significantly increase their potential ability to access lung subcellular compartments (e.g., mitochondria and nuclei), promoting strong nano-cytotoxicity and epithelial barrier dysfunction. Additionally, the presence of polyethylene glycol (PEG) as a surface coating could prevent the potential release of cytotoxic metal ions, promoting lung cytoprotection. Overall, the current work could pave the way for efficient decision-making, prediction, and mitigation of the potential occupational and environmental ENMs risks.

Quantitative structure activity relationship studies of androgen receptor binding affinity of endocrine disruptor chemicals with index of ideality of correlation, their molecular docking, molecular dynamics and ADME studies

Endocrine disrupter chemicals (EDCs) are both natural and man-made chemicals that mimic, block or interfere with human hormonal system. In the present manuscript, QSAR modeling was performed for the androgen disruptors that interfere with biosynthesis, metabolism or action of androgens that causes adverse effects on male reproductive system. A set of 96 EDCs that exhibited affinity towards androgen receptors (Log RBA) in rats were employed for carrying out QSAR studies using Hybrid descriptors (combination of HFG and SMILES) through Monte Carlo Optimization. Using index of ideality of correlation (TF2), five splits were formed and predictability of five models resulting from these splits was assessed by various validation parameters. Models resulted from first split was the top most one with R2validation = 0.7878. Structural attributes responsible for change in endpoint were studied by employing correlation weights of structural attributes. In order to further validate the model, new EDCs were designed using these attributes. In silico molecular modelling studies were performed to assess the detailed interactions with the receptor. The binding energies of all the designed compounds were observed to be better than lead and are in the range of -10.46 to -14.80. Molecular dynamics simulation of 100 ns was performed for ED01 and NED05. The results revealed that the protein-ligand complex bearing NED05 was more stable than lead ED01 exhibiting better interactions with the receptor. Further, in an attempt to assess their metabolism, ADME studies were evaluated using SwissADME. The developed model enables to predict the characteristics of designed compounds in an authentic way.Communicated by Ramaswamy H. Sarma.

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

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

Parsing structural fragments of thiazolidin-4-one based α-amylase inhibitors: A combined approach employing in vitro colorimetric screening and GA-MLR based QSAR modelling supported by molecular docking, molecular dynamics simulation and ADMET studies

α-Amylase (EC.3.2.1.1) is a ubiquitous digestive endoamylase. The abrupt rise in blood glucose levels due to the hydrolysis of carbohydrates by α-amylase at a faster rate is one of the main reasons for type 2 diabetes. The inhibitors prevent the action of digestive enzymes, slowing the digestion of carbs and eventually assisting in the management of postprandial hyperglycemia. In the course of developing α-amylase inhibitors, we have screened 2-aryliminothiazolidin-4-one based analogs for their in vitro α-amylase inhibitory potential and employed various in silico approaches for the detailed exploration of the bioactivity. The DNSA bioassay revealed that compounds 5c, 5e, 5h, 5j, 5m, 5o and 5t were more potent than the reference drug (IC₆₀ value = 22.94 ± 0.24 μg mL^-1). The derivative 5o with -NO₂ group at both the rings was the most potent analog with an IC₆₀ value of 19.67 ± 0.20 μg mL^-1 whereas derivative 5a with unsubstituted aromatic rings showed poor inhibitory potential with an IC₆₀ value of 33.40 ± 0.15 μg mL^-1. The reliable QSAR models were developed using the QSARINS software. The high value of R²_ext = 0.9632 for model IM-9 showed that the built model can be applied to predict the α-amylase inhibitory activity of the untested molecules. A consensus modelling approach was also employed to test the reliability and robustness of the developed QSAR models. Molecular docking and molecular dynamics were employed to validate the bioassay results by studying the conformational changes and interaction mechanisms. A step further, these compounds also exhibited good ADMET characteristics and bioavailability when tested for in silico pharmacokinetics prediction parameters.

Discovery of potential FGFR3 inhibitors via QSAR, pharmacophore modeling, virtual screening and molecular docking studies against bladder cancer

BACKGROUND

Fibroblast growth factor receptor 3 is known as a favorable aim in vast range of cancers, particularly in bladder cancer treatment. Pharmacophore and QSAR modeling approaches are broadly utilized for developing novel compounds for the determination of inhibitory activity versus the biological target. In this study, these methods employed to identify FGFR3 potential inhibitors.

METHODS

To find the potential compounds for bladder cancer targeting, ZINC and NCI databases were screened. Pharmacophore and QSAR modeling of FGFR3 inhibitors were utilized for dataset screening. Then, with regard to several factors such as Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) properties and Lipinski's Rule of Five, the recognized compounds were filtered. In further step, utilizing the flexible docking technique, the obtained compounds interactions with FGFR3 were analyzed.

RESULTS

The best five compounds, namely ZINC09045651, ZINC08433190, ZINC00702764, ZINC00710252 and ZINC00668789 were selected for Molecular Dynamics (MD) studies. Off-targeting of screened compounds was also investigated through CDD search and molecular docking. MD outcomes confirmed docking investigations and revealed that five selected compounds could make steady interactions with the FGFR3 and might have effective inhibitory potencies on FGFR3.

CONCLUSION

These compounds can be considered as candidates for bladder cancer therapy with improved therapeutic properties and less adverse effects.

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

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

Mode of action of p-quinone derivatives with trypanocidal activity studied by experimental and in silico models

Quinones are attractive pharmacological scaffolds for developing new agents for the treatment of different transmissible and non-transmissible human diseases due to their capacity to alter the cell redox homeostasis. The bioactivity and potential mode of action of 19 p-quinone derivatives fused to different aromatic rings (carbo or heterocycles) and harboring distinct substituents were investigated in infective Trypanosoma brucei brucei. All the compounds, except for a furanequinone (EC₅₀₌38 μM), proved to be similarly or even more potent (EC₅₀ = 0.5-5.5 μM) than the clinical drug nifurtimox (EC₅₀ = 5.3 μM). Three furanequinones and one thiazolequinone displayed a higher selectivity than nifurtimox. Two of these selective hits resulted potent inhibitors of T. cruzi proliferation (EC₅₀₌0.8-1.1 μM) but proved inactive against Leishmania infantum amastigotes. Most of the p-quinones induced a rapid and marked intracellular oxidation in T. b. brucei. DFT calculations on the oxidized quinone (Q), semiquinone (Q^•-) and hydroquinone (QH₂) suggest that all quinones have negative ΔG for the formation of Q^•-. Qualitative and quantitative structure-activity relationship analyses in two or three dimensions of different electronic and biophysical descriptors of quinones and their corresponding bioactivities (killing potency and oxidative capacity) were performed. Charge distribution over the quinone ring carbons of Q and Q^.- and the frontier orbitals energies of SUMO (Q^.-) and LUMO (Q) correlate with their oxidative and trypanocidal activity. QSAR analysis also highlighted that both bromine substitution in the p-quinone ring and a bulky phenyl group attached to the furane and thiazole rings (which generates a negative charge due to the π electron system polarized by the nearby heteroatoms) are favorable for activity. By combining experimental and in silico procedures, this study disclosed important information about p-quinones that may help to rationally tune their electronic properties and biological activities.

How the Structure of Per- and Polyfluoroalkyl Substances (PFAS) Influences Their Binding Potency to the Peroxisome Proliferator-Activated and Thyroid Hormone Receptors-An In Silico Screening Study

In this study, we investigated PFAS (per- and polyfluoroalkyl substances) binding potencies to nuclear hormone receptors (NHRs): peroxisome proliferator-activated receptors (PPARs) α, β, and γ and thyroid hormone receptors (TRs) α and β. We have simulated the docking scores of 43 perfluoroalkyl compounds and based on these data developed QSAR (Quantitative Structure-Activity Relationship) models for predicting the binding probability to five receptors. In the next step, we implemented the developed QSAR models for the screening approach of a large group of compounds (4464) from the NORMAN Database. The in silico analyses indicated that the probability of PFAS binding to the receptors depends on the chain length, the number of fluorine atoms, and the number of branches in the molecule. According to the findings, the considered PFAS group bind to the PPARα, β, and γ only with low or moderate probability, while in the case of TR α and β it is similar except that those chemicals with longer chains show a moderately high probability of binding.

Convolutional Neural Network Model Based on 2D Fingerprint for Bioactivity Prediction

Determining and modeling the possible behaviour and actions of molecules requires investigating the basic structural features and physicochemical properties that determine their behaviour during chemical, physical, biological, and environmental processes. Computational approaches such as machine learning methods are alternatives to predicting the physiochemical properties of molecules based on their structures. However, the limited accuracy and high error rates of such predictions restrict their use. In this paper, a novel technique based on a deep learning convolutional neural network (CNN) for the prediction of chemical compounds' bioactivity is proposed and developed. The molecules are represented in the new matrix format Mol2mat, a molecular matrix representation adapted from the well-known 2D-fingerprint descriptors. To evaluate the performance of the proposed methods, a series of experiments were conducted using two standard datasets, namely the MDL Drug Data Report (MDDR) and Sutherland, datasets comprising 10 homogeneous and 14 heterogeneous activity classes. After analysing the eight fingerprints, all the probable combinations were investigated using the five best descriptors. The results showed that a combination of three fingerprints, ECFP4, EPFP4, and ECFC4, along with a CNN activity prediction process, achieved the highest performance of 98% AUC when compared to the state-of-the-art ML algorithms NaiveB, LSVM, and RBFN.

Nano-QSAR modeling for predicting the cytotoxicity of metallic and metal oxide nanoparticles: A review

Given the rapid development of nanotechnology, it is crucial to understand the effects of nanoparticles on living organisms. However, it is laborious to perform toxicological tests on a case-by-case basis. Quantitative structure-activity relationship (QSAR) is an effective computational technique because it saves time, costs, and animal sacrifice. Therefore, this review presents general procedures for the construction and application of nano-QSAR models of metal-based and metal-oxide nanoparticles (MBNPs and MONPs). We also provide an overview of available databases and common algorithms. The molecular descriptors and their roles in the toxicological interpretation of MBNPs and MONPs are systematically reviewed and the future of nano-QSAR is discussed. Finally, we address the growing demand for novel nano-specific descriptors, new computational strategies to address the data shortage, in situ data for regulatory concerns, a better understanding of the physicochemical properties of NPs with bioactivity, and, most importantly, the design of nano-QSAR for real-life environmental predictions rather than laboratory simulations.

Identifications of good and bad structural fragments of hydrazone/2,5-disubstituted-1,3,4-oxadiazole hybrids with correlation intensity index and consensus modelling using Monte Carlo based QSAR studies, their molecular docking and ADME analysis

The application of QSAR along with other in silico tools like molecular docking, and molecular dynamics provide a lot of promise for finding new treatments for life-threatening diseases like Type 2 diabetes mellitus (T2DM). The present study is an attempt to develop Monte Carlo algorithm-based QSAR models using freely available CORAL software. The experimental data on the α-amylase inhibition by a series of benzothiazole-linked hydrazone/2,5-disubstituted-1,3,4-oxadiazole hybrids were selected as endpoint for the model generation. Initially, a total of eight QSAR models were built using correlation intensity index (CII) as a criterion of predictive potential. The model developed from split 6 using CII was the most reliable because of the highest numerical value of the determination coefficient of the validation set (r²_VAL = 0.8739). The important structural fragments responsible for altering the endpoint were also extracted from the best-built model. With the goal of improved prediction quality and lower prediction errors, the validated models were used to build consensus models. Molecular docking was used to know the binding mode and pose of the selected derivatives. Further, to get insight into their metabolism by living beings, ADME studies were investigated using internet freeware, SwissADME.

A Discovery Strategy for Active Compounds of Chinese Medicine Based on the Prediction Model of Compound-Disease Relationship

An accurate characterization of diseases and compounds is the key to predicting the compound-disease relationship (CDR). However, due to the difficulty of a comprehensive description of CDR, the accuracy of traditional drug development models for large-scale CDR prediction is usually unsatisfactory. In order to solve this problem, we propose a new method that integrates the molecular descriptors of compounds and the symptom descriptors of diseases to build a CDR two-dimensional matrix to predict candidate active compounds. The Matlab software draws grayscale images of CDRs, which are used as a benchmark dataset for training convolutional neural network (CNN) models. The trained model is used to predict candidate antitumor active compounds. Among the AlexNet and GoogLeNet models, we selected the GoogLeNet model for the prediction of active compounds in Chinese medicine, and its Acc, Sen, Pre, F-measure, MCC, and AUC are 0.960, 0.956, 0.965, 0.960, 0.920, and 0.964, respectively. In the prediction results of compounds, 1624 candidate CDRs were found in 124 Chinese medicines. Among them, we obtained 31 features of candidate antitumor active compounds. This method provides new insights for the discovery of candidate active compounds in Chinese medicine.

Repurposing FDA approved drugs as possible anti-SARS-CoV-2 medications using ligand-based computational approaches: sum of ranking difference-based model selection

The worldwide burden of coronavirus disease 2019 (COVID-19) is still unremittingly prevailing, with more than 440 million infections and over 5.9 million deaths documented so far since the SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) pandemic. The non-availability of treatment further aggravates the scenario, thereby demanding the exploration of pre-existing FDA-approved drugs for their effectiveness against COVID-19. The current research aims to identify potential anti-SARS-CoV-2 drugs using a computational approach and repurpose them if possible. In the present study, we have collected a set of 44 FDA-approved drugs of different classes from a previously published literature with their potential antiviral activity against COVID-19. We have employed both regression- and classification-based quantitative structure–activity relationship (QSAR) modeling to identify critical chemical features essential for anticoronaviral activity. Multiple models with the consensus algorithm were employed for the regression-based approach to improve the predictions. Additionally, we have employed a machine learning-based read-across approach using Read-Across-v3.1 available from https://sites.google.com/jadavpuruniversity.in/dtc-lab-software/home and linear discriminant analysis for the efficient prediction of potential drug candidate for COVID-19. Finally, the quantitative prediction ability of different modeling approaches was compared using the sum of ranking differences (SRD). Furthermore, we have predicted a true external set of 98 pharmaceuticals using the developed models for their probable anti-COVID activity and their prediction reliability was checked employing the “Prediction Reliability Indicator” tool available from https://dtclab.webs.com/software-tools. Though the present study does not target any protein of viral interaction, the modeling approaches developed can be helpful for identifying or screening potential anti-coronaviral drug candidates.

Monte Carlo based QSGFEAR: prediction of Gibb's free energy of activation at different temperatures using SMILES based descriptors

Monte Carlo optimization based QSGFEAR model development using CII results in the formation of more reliable, robust and predictive models.