On the virtues of automated quantitative structure-activity relationship: the new kid on the block

Active and machine learning-enhanced discovery of new FGFR3 inhibitor, Rhapontin, through virtual screening of receptor structures and anti-cancer activity assessment

Introduction: This study bridges traditional remedies and modern pharmacology by exploring the synergy between natural compounds and Ceritinib in treating Non-Small Cell Lung Cancer (NSCLC), aiming to enhance efficacy and reduce toxicities. Methods: Using a combined approach of computational analysis, machine learning, and experimental procedures, we identified and analyzed PD173074, Isoquercitrin, and Rhapontin as potential inhibitors of fibroblast growth factor receptor 3 (FGFR3). Machine learning algorithms guided the initial selection, followed by Quantitative Structure-Activity Relationship (QSAR) modeling and molecular dynamics simulations to evaluate the interaction dynamics and stability of Rhapontin. Physicochemical assessments further verified its drug-like properties and specificity. Results: Our experiments demonstrate that Rhapontin, when combined with Ceritinib, significantly suppresses tumor activity in NSCLC while sparing healthy cells. The molecular simulations and physicochemical evaluations confirm Rhapontin's stability and favorable interaction with FGFR3, highlighting its potential as an effective adjunct in NSCLC therapy. Discussion: The integration of natural compounds with established cancer therapies offers a promising avenue for enhancing treatment outcomes in NSCLC. By combining the ancient wisdom of natural remedies with the precision of modern science, this study contributes to evolving cancer treatment paradigms, potentially mitigating the side effects associated with current therapies.

Discovery of potential epidermal growth factor receptor inhibitors from black pepper for the treatment of lung cancer: an in-silico approach

A tyrosine kinase receptor known as epidermal growth factor receptor (EGFR) is one of the main tumour markers in many cancer types and also plays a crucial role in cell proliferation, differentiation, angiogenesis, and apoptosis, which is a result of the auto-phosphorylations (kinase activity enhancement) that trigger signals involved in different cellular processes. Due to the discovery that non-small cell lung cancer (NSCLC) is a cause of this kinase activity enhancement, so far, several inhibitors have been tested against EGFR, but the side effects of these inhibitors necessitate an urgent measure to come up with an inhibitor that will be more specific to the cancer cells and not affect self-cells. This study was conducted to evaluate the efficacy of 37 compounds derived from Piper nigrum against EGFR using computer-aided drug design. Based on molecular docking, induced-fit docking, calculation of free binding energy, pharmacokinetics, QSAR prediction, and MD simulation. We propose five (5) lead compounds (clarkinol A, isodihydrofutoquinol B, Burchellin, kadsurin B, and lancifolin C) as a novel inhibitor, with clarkinol A demonstrating the highest binding affinity (-7.304 kcal/mol) with EGFR when compared with the standard drug (erlotinib). They also showed significant moderation for parameters investigated for a good pharmacokinetic profile, with a reliable R2 coefficient value predicted using QSAR models. The MD simulation of clarkinol A was found to be stable within the EGFR binding pocket throughout the 75 ns simulation run time. The findings showed that clarkinol A derived from Piper nigrum is worth further investigation and consideration as a possible EGFR inhibitor for the treatment of lung cancer.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00197-1.

Probing the bioconcentration and metabolism disruption of bisphenol A and its analogues in adult female zebrafish from integrated AutoQSAR and metabolomics studies

Plenty of emerging bisphenol A (BPA) substitutes rise to wait for assessment of bioconcentration and metabolism disruption. Computational methods are useful to fill the data gap in chemical risk assessment, such as automated quantitative structure-activity relationship (AutoQSAR). It is not clear how AutoQSAR performs in predicting the bioconcentration factor (BCF) in adult zebrafish. Herein, AutoQSAR was used to predict the logBCFs of BPA, bisphenol AF (BPAF), bisphenol B, bisphenol F and bisphenol S (BPS). For the test set, a linear relationship was shown between the observed and predicted logBCFs with a slope of 0.97. The predicted logBCFs of these five bisphenols were quite close to their experimental data with a slope of 0.94, suggesting better performance than directed message passing neural networks and EPI Suite with a slope of 0.69 and 0.61, respectively. Thus, AutoQSAR is powerful in modeling logBCFs in fish with minimal time and expertise. To link bioconcentration with metabolic effects, female zebrafish were exposed to BPA, BPAF and BPS for metabolomics analysis. BPA caused a significant disturbance in amino acid metabolism, while BPAF and BPS significantly altered another three metabolic pathways, showing chemical-specific responses. BPAF with the highest logBCF elicited the strongest metabolomic responses reflected by the metabolic effect level index, followed by BPA and BPS. Thus, BPAF and BPS elicited higher or similar metabolism disruption compared with BPA in female zebrafish, respectively, reflecting consequences of bioconcentration.

Quantitative structure-activity relationship and machine learning studies of 2-thiazolylhydrazone derivatives with anti-Cryptococcus neoformans activity

Cryptococcus neoformans is a fungus responsible for infections in humans with a significant number of cases in immunosuppressed patients, mainly in underdeveloped countries. In this context, the thiazolylhydrazones are a promising class of compounds with activity against C. neoformans. The understanding of the structure-activity relationship of these derivatives could lead to the design of robust compounds that could be promising drug candidates for fungal infections. Specifically, modern techniques such as 4D-QSAR and machine learning methods were employed in this work to generate two QSAR models (one 2D and one 4D) with high predictive power (r² for the test set equals to 0.934 and 0.831, respectively), and one random forest classification model was reported with Matthews correlation coefficient equals to 1 and 0.62 for internal and external validations, respectively. The physicochemical interpretation of selected models, indicated the importance of aliphatic substituents at the hydrazone moiety to antifungal activity, corroborating experimental data.Communicated by Ramaswamy H. Sarma.

Exploration of the Mechanism of Salvianolic Acid for Injection Against Ischemic Stroke: A Research Based on Computational Prediction and Experimental Validation

Ischemic stroke (IS) is an acute neurological injury that occurs when a vessel supplying blood to the brain is obstructed, which is a leading cause of death and disability. Salvia miltiorrhiza has been used in the treatment of cardiovascular and cerebrovascular diseases for over thousands of years due to its effect activating blood circulation and dissipating blood stasis. However, the herbal preparation is chemically complex and the diversity of potential targets makes difficult to determine its mechanism of action. To gain insight into its mechanism of action, we analyzed “Salvianolic acid for injection” (SAFI), a traditional Chinese herbal medicine with anti-IS effects, using computational systems pharmacology. The potential targets of SAFI, obtained from literature mining and database searches, were compared with IS-associated genes, giving 38 common genes that were related with pathways involved in inflammatory response. This suggests that SAFI might function as an anti-inflammatory agent. Two genes associated with inflammation (PTGS1 and PTGS2), which were inhibited by SAFI, were preliminarily validated in vitro. The results showed that SAFI inhibited PTGS1 and PTGS2 activity in a dose-dependent manner and inhibited the production of prostaglandin E2 induced by lipopolysaccharide in RAW264.7 macrophages and BV-2 microglia. This approach reveals the possible pharmacological mechanism of SAFI acting on IS, and also provides a feasible way to elucidate the mechanism of traditional Chinese medicine (TCM).

Design of Organic Electronic Materials With a Goal-Directed Generative Model Powered by Deep Neural Networks and High-Throughput Molecular Simulations

In recent years, generative machine learning approaches have attracted significant attention as an enabling approach for designing novel molecular materials with minimal design bias and thereby realizing more directed design for a specific materials property space. Further, data-driven approaches have emerged as a new tool to accelerate the development of novel organic electronic materials for organic light-emitting diode (OLED) applications. We demonstrate and validate a goal-directed generative machine learning framework based on a recurrent neural network (RNN) deep reinforcement learning approach for the design of hole transporting OLED materials. These large-scale molecular simulations also demonstrate a rapid, cost-effective method to identify new materials in OLEDs while also enabling expansion into many other verticals such as catalyst design, aerospace, life science, and petrochemicals.

Identification of lead compounds from large natural product library targeting 3C-like protease of SARS-CoV-2 using E-pharmacophore modelling, QSAR and molecular dynamics simulation

COVID-19 is a novel disease caused by SARS-CoV-2 and has made a catastrophic impact on the global economy. As it is, there is no officially FDA approved drug to alleviate the negative impact of SARS-CoV-2 on human health. Numerous drug targets for neutralizing coronavirus infection have been identified, among them is 3-chymotrypsin-like-protease (3CL^pro), a viral protease responsible for the viral replication is chosen for this study. This study aimed at finding novel inhibitors of SARS-CoV-2 3C-like protease from the natural library using computational approaches. A total of 69,000 compounds from natural product library were screened to match a minimum of 3 features from the five sites e-pharmacophore model. Compounds with fitness score of 1.00 and above were consequently filtered by executing molecular docking studies via Glide docking algorithm. Qikprop also predicted the compounds drug-likeness and pharmacokinetic features; besides, the QSAR model built from KPLS analysis with radial as binary fingerprint was used to predict the compounds inhibition properties against SARS-CoV-2 3C-like protease. Fifty ns molecular dynamics (MD) simulation was carried out using GROMACS software to understand the dynamics of binding. Nine (9) lead compounds from the natural products library were discovered; seven among them were found to be more potent than lopinavir based on energies of binding. STOCK1N-98687 with docking score of -9.295 kcal/mol had considerable predicted bioactivity (4.427 µM) against SARS-CoV-2 3C-like protease and satisfactory drug-like features than the experimental drug lopinavir. Post-docking analysis by MM-GBSA confirmed the stability of STOCK1N-98687 bound 3CL^pro crystal structure. MD simulation of STOCKIN-98687 with 3CL^pro at 50 ns showed high stability and low fluctuation of the complex. This study revealed compound STOCK1N-98687 as potential 3CL^pro inhibitor; therefore, a wet experiment is worth exploring to confirm the therapeutic potential of STOCK1N-98687 as an antiviral agent.

Discovery of potent Covid-19 main protease inhibitors using integrated drug-repurposing strategy

The emergence and rapid spreading of novel SARS-CoV-2 across the globe represent an imminent threat to public health. Novel antiviral therapies are urgently needed to overcome this pandemic. Given the significant role of the main protease of Covid-19 for virus replication, we performed a drug-repurposing study using the recently deposited main protease structure, 6LU7. For instance, pharmacophore- and e-pharmacophore-based hypotheses such as AARRH and AARR, respectively, were developed using available small molecule inhibitors and utilized in the screening of the DrugBank repository. Further, a hierarchical docking protocol was implemented with the support of the Glide algorithm. The resultant compounds were then examined for their binding free energy against the main protease of Covid-19 by means of the Prime-MM/GBSA algorithm. Most importantly, the machine learning-based AutoQSAR algorithm was used to predict the antiviral activities of resultant compounds. The hit molecules were also examined for their drug-likeness and toxicity parameters through the QikProp algorithm. Finally, the hit compounds activity against the main protease was validated using molecular dynamics simulation studies. Overall, the present analysis yielded two potential inhibitors (DB02986 and DB08573) that are predicted to bind with the main protease of Covid-19 better than currently used drug molecules such as N3 (cocrystallized native ligand), lopinavir, and ritonavir.

Deciphering the Interactions of Bioactive Compounds in Selected Traditional Medicinal Plants against Alzheimer's Diseases via Pharmacophore Modeling, Auto-QSAR, and Molecular Docking Approaches

Neurodegenerative diseases, for example Alzheimer’s, are perceived as driven by hereditary, cellular, and multifaceted biochemical actions. Numerous plant products, for example flavonoids, are documented in studies for having the ability to pass the blood-brain barrier and moderate the development of such illnesses. Computer-aided drug design (CADD) has achieved importance in the drug discovery world; innovative developments in the aspects of structure identification and characterization, bio-computational science, and molecular biology have added to the preparation of new medications towards these ailments. In this study we evaluated nine flavonoid compounds identified from three medicinal plants, namely T. diversifolia, B. sapida, and I. gabonensis for their inhibitory role on acetylcholinesterase (AChE), butyrylcholinesterase (BChE) and monoamine oxidase (MAO) activity, using pharmacophore modeling, auto-QSAR prediction, and molecular studies, in comparison with standard drugs. The results indicated that the pharmacophore models produced from structures of AChE, BChE and MAO could identify the active compounds, with a recuperation rate of the actives found near 100% in the complete ranked decoy database. Moreso, the robustness of the virtual screening method was accessed by well-established methods including enrichment factor (EF), receiver operating characteristic curve (ROC), Boltzmann-enhanced discrimination of receiver operating characteristic (BEDROC), and area under accumulation curve (AUAC). Most notably, the compounds’ pIC₅₀ values were predicted by a machine learning-based model generated by the AutoQSAR algorithm. The generated model was validated to affirm its predictive model. The best models achieved for AChE, BChE and MAO were models kpls_radial_17 (R² = 0.86 and Q² = 0.73), pls_38 (R² = 0.77 and Q² = 0.72), kpls_desc_44 (R² = 0.81 and Q² = 0.81) and these externally validated models were utilized to predict the bioactivities of the lead compounds. The binding affinity results of the ligands against the three selected targets revealed that luteolin displayed the highest affinity score of −9.60 kcal/mol, closely followed by apigenin and ellagic acid with docking scores of −9.60 and −9.53 kcal/mol, respectively. The least binding affinity was attained by gallic acid (−6.30 kcal/mol). The docking scores of our standards were −10.40 and −7.93 kcal/mol for donepezil and galanthamine, respectively. The toxicity prediction revealed that none of the flavonoids presented toxicity and they all had good absorption parameters for the analyzed targets. Hence, these compounds can be considered as likely leads for drug improvement against the same.

Insight into glycogen synthase kinase-3β inhibitory activity of phyto-constituents from Melissa officinalis: in silico studies

Over activity of Glycogen synthase kinase-3β (GSK-3β), a serine/threonine-protein kinase has been implicated in a number of diseases including stroke, type II diabetes and Alzheimer disease (AD). This study aimed to find novel inhibitors of GSK-3β from phyto-constituents of Melissa officinalis with the aid of computational analysis. Molecular docking, induced-fit docking (IFD), calculation of binding free energy via the MM-GBSA approach and Lipinski's rule of five (RO5) were employed to filter the compounds and determine their druggability. Most importantly, the compounds pIC50 were predicted by machine learning-based model generated by AutoQSAR algorithm. The generated model was validated to affirm its predictive model. The best model obtained was Model kpls_desc_38 (R2 = 0.8467 and Q2 = 0.8069), and this external validated model was utilized to predict the bioactivities of the lead compounds. While a number of characterized compounds from Melissa officinalis showed better docking score, binding free energy alongside adherence to RO5 than co-cystallized ligand, only three compounds (salvianolic acid C, ellagic acid and naringenin) showed more satisfactory pIC50. The results obtained in this study can be useful to design potent inhibitors of GSK-3β.

Exploring structural features of EGFR-HER2 dual inhibitors as anti-cancer agents using G-QSAR approach

Simultaneous inhibition of EGFR and HER2 by dual-targeting inhibitors is an established anti-cancer strategy. Therefore, a recent trend in drug discovery involves understanding the features of such dual inhibitors. In this study, three different G-QSAR models were developed corresponding to individual EGFR, HER2 and the dual-model for both receptors. The dual-model provided site-specific information wherein (i) increasing electronegative character and higher index of saturated carbon at R4 position; (ii) presence of chlorine atom at R2 position; (iii) decreasing alpha modified shape index at R1 and R3 positions; and (iv) less electronegativity at R2 position; were found important for enhancing the dual activity. Also, comparison of dual-model with the EGFR/HER2 individual models revealed that it incorporates the properties of both models and, thus, represents a combination of EGFR/HER2. Further, fragment analysis revealed that R2 and R4 are important for imparting high potency while specificity is decided by R1/R3 fragment. We also checked the predictive ability of the dual-model by determining applicability domain using William's plot. Also, analysis of active molecules showed they show favorable substitutions that agree with the constructed dual-model. Thus, we have been successful in developing a single dual-response QSAR model to get an insight into various structural features influencing EGFR/HER2 activity.