A Review on Applications of Computational Methods in Drug Screening and Design

Targeting bacterial growth in biofilm conditions: rational design of novel inhibitors to mitigate clinical and food contamination using QSAR

Antimicrobial resistance (AMR) is a pressing global issue exacerbated by the abuse of antibiotics and the formation of bacterial biofilms, which cause up to 80% of human bacterial infections. This study presents a computational strategy to address AMR by developing three novel quantitative structure-activity relationship (QSAR) models based on molecular topology to identify potential anti-biofilm and antibacterial agents. The models aim to determine the chemo-topological pattern of Gram (+) antibacterial, Gram (-) antibacterial, and biofilm formation inhibition activity. The models were applied to the virtual screening of a commercial chemical database, resulting in the selection of 58 compounds. Subsequent in vitro assays showed that three of these compounds exhibited the most promising antibacterial activity, with potential applications in enhancing food and medical device safety.

Addressing docking pose selection with structure-based deep learning: Recent advances, challenges and opportunities

Molecular docking is a widely used technique in drug discovery to predict the binding mode of a given ligand to its target. However, the identification of the near-native binding pose in docking experiments still represents a challenging task as the scoring functions currently employed by docking programs are parametrized to predict the binding affinity, and, therefore, they often fail to correctly identify the ligand native binding conformation. Selecting the correct binding mode is crucial to obtaining meaningful results and to conveniently optimizing new hit compounds. Deep learning (DL) algorithms have been an area of a growing interest in this sense for their capability to extract the relevant information directly from the protein-ligand structure. Our review aims to present the recent advances regarding the development of DL-based pose selection approaches, discussing limitations and possible future directions. Moreover, a comparison between the performances of some classical scoring functions and DL-based methods concerning their ability to select the correct binding mode is reported. In this regard, two novel DL-based pose selectors developed by us are presented.

Comprehensive applications of the artificial intelligence technology in new drug research and development

Purpose

Target-based strategy is a prevalent means of drug research and development (R&D), since targets provide effector molecules of drug action and offer the foundation of pharmacological investigation. Recently, the artificial intelligence (AI) technology has been utilized in various stages of drug R&D, where AI-assisted experimental methods show higher efficiency than sole experimental ones. It is a critical need to give a comprehensive review of AI applications in drug R &D for biopharmaceutical field.

Methods

Relevant literatures about AI-assisted drug R&D were collected from the public databases (Including Google Scholar, Web of Science, PubMed, IEEE Xplore Digital Library, Springer, and ScienceDirect) through a keyword searching strategy with the following terms [("Artificial Intelligence" OR "Knowledge Graph" OR "Machine Learning") AND ("Drug Target Identification" OR "New Drug Development")].

Results

In this review, we first introduced common strategies and novel trends of drug R&D, followed by characteristic description of AI algorithms widely used in drug R&D. Subsequently, we depicted detailed applications of AI algorithms in target identification, lead compound identification and optimization, drug repurposing, and drug analytical platform construction. Finally, we discussed the challenges and prospects of AI-assisted methods for drug discovery.

Conclusion

Collectively, this review provides comprehensive overview of AI applications in drug R&D and presents future perspectives for biopharmaceutical field, which may promote the development of drug industry.

Discovery of novel and highly potential inhibitors of glycogen synthase kinase 3-beta (GSK-3β) through structure-based pharmacophore modeling, virtual computational screening, docking and in silico ADMET analysis

The protein Glycogen Synthase Kinase 3-Beta (GSK-3β), is a promising therapeutic target for treating various diseases such as neurodegenerative disorders, diabetes, inflammation and cancer. This study aims to investigate the potential of compounds targeting inflammation or carbohydrate metabolism to selectively inhibit GSK3β by binding to its ATP site. To achieve this goal, we filtered a database of 49367 molecules involved in carbohydrate metabolism or targeting inflammation using various computational analyses, including pharmacophore modeling, molecular docking, dynamic simulation, prime MM-GBSA calculation, and in silico ADME studies. We generated a pharmacophore model (hypo S: AADDHRR) using two different crystallographic complexes of GSK3β and evaluated the model's performance in identifying hits using various parameters, including EF, GH, ROC, AUC and BEDROC. Subsequently, we performed various dockings (HTVS, SP, XP and IFD) for the retrieved hits and found that, 5 out of the top 10 ranked compounds had the scaffold of pyrazolidine 3,5-dione, which has never been reported to inhibit kinases. We also conducted ADMET studies to and concluded that compound N6 exhibited the best pharmacokinetic profile passing the blood-brain barrier, possessing high lipophilicity and a high coefficient of skin permeability in the intestines, along with good bioavailability and low toxicity risk assessment. Dynamic simulation were also performed indicating that compounds N6 derived from pyrazolidine 3,5-dione demonstrated better binding potential for GSK3β during the simulation period. Therefore, we propose that compounds derived from pyrazolidine-3,5-dione, which modulate the activity of lysosomal alpha-glucosidase could serve as a novel scaffold for the selective inhibition of GSK-3β.Communicated by Ramaswamy H. Sarma.

In silico analysis of potential inhibitors for breast cancer targeting 17beta-hydroxysteroid dehydrogenase type 1 (17beta-HSD1) catalyses

Breast cancer (BC) is still one of the major issues in world health, especially for women, which necessitates innovative therapeutic strategies. In this study, we investigated the efficacy of retinoic acid derivatives as inhibitors of 17beta-hydroxysteroid dehydrogenase type 1 (17beta-HSD1), which plays a crucial role in the biosynthesis and metabolism of oestrogen and thereby influences the progression of BC and, the main objective of this investigation is to identify the possible drug candidate against BC through computational drug design approach including PASS prediction, molecular docking, ADMET profiling, molecular dynamics simulations (MD) and density functional theory (DFT) calculations. The result has reported that total eight derivatives with high binding affinity and promising pharmacokinetic properties among 115 derivatives. In particular, ligands 04 and 07 exhibited a higher binding affinity with values of -9.9 kcal/mol and -9.1 kcal/mol, respectively, than the standard drug epirubicin hydrochloride, which had a binding affinity of -8.2 kcal/mol. The stability of the ligand-protein complexes was further confirmed by MD simulations over a 100-ns trajectory, which included assessments of hydrogen bonds, root mean square deviation (RMSD), root mean square Fluctuation (RMSF), dynamic cross-correlation matric (DCCM) and principal component analysis. The study emphasizes the need for experimental validation to confirm the therapeutic utility of these compounds. This study enhances the computational search for new BC drugs and establishes a solid foundation for subsequent experimental and clinical research.

Synthesis and characterization of gold(I) thiolate derivatives and bimetallic complexes for HIV inhibition

Introduction: The human immunodeficiency virus (HIV) remains a significant global health concern, with a reported high infection rate of 38.4 million cases globally; an estimated 2 million new infections and approximately 700,000 HIV/AIDS-related deaths were reported in 2021. Despite the advent of anti-retroviral therapy (ART), HIV/AIDS persists as a chronic disease. To combat this, several studies focus on developing inhibitors targeting various stages of the HIV infection cycle, including HIV-1 protease. This study aims to synthesize and characterize novel glyco diphenylphosphino metal complexes with potential HIV inhibitory properties. Method: A series of new gold(I) thiolate derivatives and three bimetallic complexes, incorporating amino phosphines and thiocarbohydrate as auxiliary ligands, were synthesized using procedures described by Jiang, et al. (2009) and Coetzee et al. (2007). Structural elucidation and purity assessment of the synthesized compounds (1-11) were conducted using micro-analysis, NMR, and infrared spectrometry. Results and Discussion: Using molecular modeling techniques, three of the metal complexes were identified as potential HIV protease inhibitors, exhibiting strong binding affinity interactions with binding pocket residues. These inhibitors demonstrated an ability to inhibit the flexibility of the flap regions of the HIV protease, similar to the known HIV protease inhibitor, darunavir. This study sheds light on the promising avenues for the development of novel therapeutic agents against HIV/AIDS.

A Perspective on Interdicting in Protein Misfolding for Therapeutic Drug Design: Modulating the Formation of Nonlocal Contacts in α-Synuclein as a Strategy against Parkinson's Disease

In recent work we proposed that interdiction in the earliest contact-formation events along the folding pathway of key viral proteins could provide a novel avenue for therapeutic drug design. In this Perspective we explore the potential applicability of the protein folding interdiction strategy in the realm of neurodegenerative diseases with a specific focus on synucleinopathies. In order to fulfill this goal we review the interdiction proposal and its practical challenges, and we present new results concerning design strategies for possible peptide drugs that could be useful in preventing α-synuclein aggregation.

Acceleration of chondrogenic differentiation utilizing biphasic core-shell alginate sulfate electrospun nanofibrous scaffold

Engineering new biomedical materials with tailored physicochemical, mechanical and biological virtues in order to differentiate stem cells into chondrocytes for cartilage regeneration has garnered much scientific interest. In this study, core/shell nanofibrous scaffold based on poly(ɛ-caprolactone) (PCL) as a core material and alginate sulfate (AlgS)-poly(vinyl alcohol) (PVA) blend as shell materials (AlgS-PVA/PCL) was fabricated by emulsion electrospinning. In this vein, the influence of AlgS to PVA ratio (30:70, 50:50), organic to aqueous phase ratio (1:2, 1:3 and 1:5) and acid concentration (0, 10, 20, 30, 40 and 50 %) on nanofibers morphology were investigated. SEM images depicted that AlgS to PVA ratio of 30:70 and 50:50, organic to aqueous phase ratio of 1:3 and 1:5 and acid concentration of 30 % led to uniform, bead-free fibrous mats. AlgS-PVA/PCL scaffolds with AlgS to PVA ratio of 30:70 and organic to aqueous phase ratio of 1:3, showed admirable mechanical features, high porosity (>90 %) with desirable swelling ratio in wet condition. In vitro assays indicated that the AlgS-PVA/PCL scaffold surface had desirable interaction with stem cells and promotes cells attachment, proliferation and differentiation. Thus, we envision that this salient structure could be an intriguing construction as a cartilage tissue-engineered scaffold.

Design of vilazodone-donepezil chimeric derivatives as acetylcholinesterase inhibitors by QSAR, molecular docking and molecular dynamics simulations

Alzheimer's disease (AD) is a disease that affects the cognitive abilities of older adults, and it is one of the biggest global medical challenges of the 21st century. Acetylcholinesterase (AChE) can increase acetylcholine concentrations and improve cognitive function in patients, and is a potential target to develop small molecule inhibitors for the treatment of Alzheimer's disease (AD). In this study, 29 vilazodone-donepezil chimeric derivatives are systematically studied using 3D-QSAR modeling, and a robust and reliable Topomer CoMFA model was obtained with: q2 = 0.720, r2 = 0.991, F = 287.234, N = 6, and SEE = 0.098. Based on the established model and combined with the ZINC20 database, 33 new compounds with ideal inhibitory activity are successfully designed. Molecular docking and ADMET property prediction also show that these newly designed compounds have a good binding ability to the target protein and can meet the medicinal conditions. Subsequently, four new compounds with good comprehensive ability are selected for molecular dynamics simulation, and the simulation results confirm that the newly designed compounds have a certain degree of reliability and stability. This study provides guidance for vilazodone-donepezil chimeric derivatives as a potential AChE inhibitor and has certain theoretical value.

POxload: Machine Learning Estimates Drug Loadings of Polymeric Micelles

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

Pharmaceutical Potential of Remedial Plants and Helminths for Treating Inflammatory Bowel Disease

Research is increasingly revealing that inflammation significantly contributes to various diseases, particularly inflammatory bowel disease (IBD). IBD is a major medical challenge due to its chronic nature, affecting at least one in a thousand individuals in many Western countries, with rising incidence in developing nations. Historically, indigenous people have used natural products to treat ailments, including IBD. Ethnobotanically guided studies have shown that plant-derived extracts and compounds effectively modulate immune responses and reduce inflammation. Similarly, helminths and their products offer unique mechanisms to modulate host immunity and alleviate inflammatory responses. This review explored the pharmaceutical potential of Aboriginal remedial plants and helminths for treating IBD, emphasizing recent advances in discovering anti-inflammatory small-molecule drug leads. The literature from Scopus, MEDLINE Ovid, PubMed, Google Scholar, and Web of Science was retrieved using keywords such as natural product, small molecule, cytokines, remedial plants, and helminths. This review identified 55 important Aboriginal medicinal plants and 9 helminth species that have been studied for their anti-inflammatory properties using animal models and in vitro cell assays. For example, curcumin, berberine, and triptolide, which have been isolated from plants; and the excretory-secretory products and their protein, which have been collected from helminths, have demonstrated anti-inflammatory activity with lower toxicity and fewer side effects. High-throughput screening, molecular docking, artificial intelligence, and machine learning have been engaged in compound identification, while clustered regularly interspaced short palindromic repeats (CRISPR) gene editing and RNA sequencing have been employed to understand molecular interactions and regulations. While there is potential for pharmaceutical application of Aboriginal medicinal plants and gastrointestinal parasites in treating IBD, there is an urgent need to qualify these plant and helminth therapies through reproducible clinical and mechanistic studies.

One-Step Synthesis, Crystallography, and Acute Toxicity of Two Boron-Carbohydrate Adducts That Induce Sedation in Mice

Boronic acids form diester bonds with cis-hydroxyl groups in carbohydrates. The formation of these adducts could impair the physical and chemical properties of precursors, even their biological activity. Two carbohydrate derivatives from d-fructose and d-arabinose and phenylboronic acid were synthesized in a straightforward one-step procedure and chemically characterized via spectroscopy and X-ray diffraction crystallography. Additionally, an acute toxicity test was performed to determine their lethal dose 50 (LD50) values by using Lorke's method. Analytical chemistry assays confirmed the formation of adducts by the generation of diester bonds with the β-d-pyranose of carbohydrates, including signals corresponding to the formation of new bonds, such as the stretching of B-O bonds. NMR spectra yielded information about the stereoselectivity in the synthesis reaction: Just one signal was found in the range for the anomeric carbon in the 13C NMR spectra of both adducts. The acute toxicity tests showed that the LD50 value for both compounds was 1265 mg/kg, while the effective dose 50 (ED50) for sedation was 531 mg/kg. However, differences were found in the onset and lapse of sedation. For example, the arabinose derivative induced sedation for more than 48 h at 600 mg/kg, while the fructose derivative induced sedation for less than 6 h at the same dose without the death of the mice. Thus, we report for the first time two boron-containing carbohydrate derivatives inducing sedation after intraperitoneal administration. They are bioactive and highly safe agents. Further biological evaluation is desirable to explore their medical applications.

Drug-target interaction predictions with multi-view similarity network fusion strategy and deep interactive attention mechanism

MOTIVATION

Accurately identifying the drug-target interactions (DTIs) is one of the crucial steps in the drug discovery and drug repositioning process. Currently, many computational-based models have already been proposed for DTI prediction and achieved some significant improvement. However, these approaches pay little attention to fuse the multi-view similarity networks related to drugs and targets in an appropriate way. Besides, how to fully incorporate the known interaction relationships to accurately represent drugs and targets is not well investigated. Therefore, there is still a need to improve the accuracy of DTI prediction models.

RESULTS

In this study, we propose a novel approach that employs Multi-view similarity network fusion strategy and deep Interactive attention mechanism to predict Drug-Target Interactions (MIDTI). First, MIDTI constructs multi-view similarity networks of drugs and targets with their diverse information and integrates these similarity networks effectively in an unsupervised manner. Then, MIDTI obtains the embeddings of drugs and targets from multi-type networks simultaneously. After that, MIDTI adopts the deep interactive attention mechanism to further learn their discriminative embeddings comprehensively with the known DTI relationships. Finally, we feed the learned representations of drugs and targets to the multilayer perceptron model and predict the underlying interactions. Extensive results indicate that MIDTI significantly outperforms other baseline methods on the DTI prediction task. The results of the ablation experiments also confirm the effectiveness of the attention mechanism in the multi-view similarity network fusion strategy and the deep interactive attention mechanism.

AVAILABILITY AND IMPLEMENTATION

https://github.com/XuLew/MIDTI.

DTKGIN: Predicting drug-target interactions based on knowledge graph and intent graph

Knowledge graph intent graph attention mechanism Predicting drug-target interactions (DTIs) plays a crucial role in drug discovery and drug development. Considering the high cost and risk of biological experiments, developing computational approaches to explore the interactions between drugs and targets can effectively reduce the time and cost of drug development. Recently, many methods have made significant progress in predicting DTIs. However, existing approaches still suffer from the high sparsity of DTI datasets and the cold start problem. In this paper, we develop a new model to predict drug-target interactions via a knowledge graph and intent graph named DTKGIN. Our method can effectively capture biological environment information for targets and drugs by mining their associated relations in the knowledge graph and considering drug-target interactions at a fine-grained level in the intent graph. DTKGIN learns the representation of drugs and targets from the knowledge graph and the intent graph. Then the probabilities of interactions between drugs and targets are obtained through the inner product of the representation of drugs and targets. Experimental results show that our proposed method outperforms other state-of-the-art methods in 10-fold cross-validation, especially in cold-start experimental settings. Furthermore, the case studies demonstrate the effectiveness of DTKGIN in predicting potential drug-target interactions. The code is available on GitHub: https://github.com/Royluoyi123/DTKGIN.

PharmaCore: The Automatic Generation of 3D Structure-Based Pharmacophore Models from Protein/Ligand Complexes

In this work, we present PharmaCore: a new, completely automatic workflow aimed at generating three-dimensional (3D) structure-based pharmacophore models toward any target of interest. The proposed approach relies on using cocrystallized ligands to create the input files for generating the pharmacophore hypotheses, integrating not only the three-dimensional structural information on the ligand but also data concerning the binding mode of these molecules put in the protein cavity. We developed a Python library that, starting from the specific UniProt ID of the protein under investigation as the only element that requires user intervention, subsequently collects and aligns the corresponding structures bearing a known ligand in a fully automated fashion, bringing them all into the same coordinate system. The protocol includes a final phase in which the aligned small molecules are used to produce the pharmacophore hypotheses directly onto the protein structure using a specific software, e.g., Phase (Schrödinger LLC). To validate the entire procedure and highlight the possible applications in the field of drug discovery and repositioning, we first generated pharmacophores for soluble epoxide hydrolase (sEH) and compared with already-published ones. Then, we reproduced the binding profile of a reported selective binder of ATAD2 bromodomain (AM879), testing it against a panel of 1741 pharmacophores related to 16 epigenetic proteins and automatically generated with PharmaCore, finally disclosing putative unprecedented off-targets. The computational predictions were successfully validated with AlphaScreen assays, highlighting the applicability of the proposed workflow in drug discovery and repositioning. Finally, the process was also validated on tankyrase 2 and SARS-CoV-2 MPro, confirming the robustness of PharmaCore.

Screening approach by a combination of computational and in vitro experiments: identification of fluvastatin sodium as a potential SIRT2 inhibitor

BACKGROUND

Sirtuins (SIRTs) are NAD+-dependent deacetylases that play various roles in numerous pathophysiological processes, holding promise as therapeutic targets worthy of further investigation. Among them, the SIRT2 subtype is closely associated with tumorigenesis and malignancies. Dysregulation of SIRT2 activation can regulate the expression levels of related genes in cancer cells, leading to tumor occurrence and metastasis.

METHODS

In this study, we used computer simulations to screen for novel SIRT2 inhibitors from the FDA database, based on which 10 compounds with high docking scores and good interactions were selected for in vitro anti-pancreatic cancer metastasis testing and enzyme binding inhibition experiments. The results showed that fluvastatin sodium may possess inhibitory activity against SIRT2. Subsequently, fluvastatin sodium was subjected to molecular docking experiments with various SIRT isoforms, and the combined results from Western blotting experiments indicated its potential as a SIRT2 inhibitor. Next, molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations were performed, revealing the binding mode of fluvastatin sodium at the SIRT2 active site, further validating the stability and interaction of the ligand-protein complex under physiological conditions.

RESULTS

Overall, this study provides a systematic virtual screening workflow for the discovery of SIRT2 activity inhibitors, identifies the potential inhibitory effect of fluvastatin sodium as a lead compound on SIRT2, and opens up a new direction for developing highly active and selectively targeted SIRT2 inhibitors.

Subtractive genomics study for the identification of therapeutic targets against Cronobacter sakazakii: A threat to infants

Cronobacter sakazakii is an opportunistic pathogen that has been associated with severe infection in neonates such as necrotizing enterocolitis (NEC), neonatal meningitis, and bacteremia. This pathogen can survive in a relatively dry environment, especially in powdered infant formula (PIF). Unfortunately, conventional drugs that were once effective against C. sakazakii are gradually losing their efficacy due to rising antibiotic resistance. In this study, a subtractive genomic approach was followed in order to identify potential therapeutic targets in the pathogen. The whole proteome of the pathogen was filtered through a step-by-step process, which involved removing paralogous proteins, human homologs, sequences that are less essential for survival, proteins with shared metabolic pathways, and proteins that are located in cells other than the cytoplasmic membrane. As a result, nine novel drug targets were identified. Further, the analysis also unveiled that the FDA-approved drug Terbinafine can be repurposed against the Glutathione/l-cysteine transport system ATP-binding/permease protein CydC of C. sakazakii. Moreover, molecular docking and dynamics studies of Terbinafine and CydC suggested that this drug can be used to treat C. sakazakii infection in neonates. However, for clinical purposes further in vitro and in vivo studies are necessary.

An Overview of SARS-CoV-2 Potential Targets, Inhibitors, and Computational Insights to Enrich the Promising Treatment Strategies

The rapid spread of the SARS-CoV-2 virus has emphasized the urgent need for effective therapies to combat COVID-19. Investigating the potential targets, inhibitors, and in silico approaches pertinent to COVID-19 are of utmost need to develop novel therapeutic agents and reprofiling of existing FDA-approved drugs. This article reviews the viral enzymes and their counter receptors involved in the entry of SARS-CoV-2 into host cells, replication of genomic RNA, and controlling the host cell physiology. In addition, the study provides an overview of the computational techniques such as docking simulations, molecular dynamics, QSAR modeling, and homology modeling that have been used to find the FDA-approved drugs and other inhibitors against SARS-CoV-2. Furthermore, a comprehensive overview of virus-based and host-based druggable targets from a structural point of view, together with the reported therapeutic compounds against SARS-CoV-2 have also been presented. The current study offers future perspectives for research in the field of network pharmacology investigating the large unexplored molecular libraries. Overall, the present in-depth review aims to expedite the process of identifying and repurposing drugs for researchers involved in the field of COVID-19 drug discovery.

The Histone Deacetylase Family: Structural Features and Application of Combined Computational Methods

Histone deacetylases (HDACs) are crucial in gene transcription, removing acetyl groups from histones. They also influence the deacetylation of non-histone proteins, contributing to the regulation of various biological processes. Thus, HDACs play pivotal roles in various diseases, including cancer, neurodegenerative disorders, and inflammatory conditions, highlighting their potential as therapeutic targets. This paper reviews the structure and function of the four classes of human HDACs. While four HDAC inhibitors are currently available for treating hematological malignancies, numerous others are undergoing clinical trials. However, their non-selective toxicity necessitates ongoing research into safer and more efficient class-selective or isoform-selective inhibitors. Computational methods have aided the discovery of HDAC inhibitors with the desired potency and/or selectivity. These methods include ligand-based approaches, such as scaffold hopping, pharmacophore modeling, three-dimensional quantitative structure-activity relationships, and structure-based virtual screening (molecular docking). Moreover, recent developments in the field of molecular dynamics simulations, combined with Poisson-Boltzmann/molecular mechanics generalized Born surface area techniques, have improved the prediction of ligand binding affinity. In this review, we delve into the ways in which these methods have contributed to designing and identifying HDAC inhibitors.

TransGEM: a molecule generation model based on Transformer with gene expression data

MOTIVATION

It is difficult to generate new molecules with desirable bioactivity through ligand-based de novo drug design, and receptor-based de novo drug design is constrained by disease target information availability. The combination of artificial intelligence and phenotype-based de novo drug design can generate new bioactive molecules, independent from disease target information. Gene expression profiles can be used to characterize biological phenotypes. The Transformer model can be utilized to capture the associations between gene expression profiles and molecular structures due to its remarkable ability in processing contextual information.

RESULTS

We propose TransGEM (Transformer-based model from gene expression to molecules), which is a phenotype-based de novo drug design model. A specialized gene expression encoder is used to embed gene expression difference values between diseased cell lines and their corresponding normal tissue cells into TransGEM model. The results demonstrate that the TransGEM model can generate molecules with desirable evaluation metrics and property distributions. Case studies illustrate that TransGEM model can generate structurally novel molecules with good binding affinity to disease target proteins. The majority of genes with high attention scores obtained from TransGEM model are associated with the onset of the disease, indicating the potential of these genes as disease targets. Therefore, this study provides a new paradigm for de novo drug design, and it will promote phenotype-based drug discovery.

AVAILABILITY AND IMPLEMENTATION

The code is available at https://github.com/hzauzqy/TransGEM.

G-K BertDTA: A graph representation learning and semantic embedding-based framework for drug-target affinity prediction

Developing new drugs is costly, time-consuming, and risky. Drug-target affinity (DTA), indicating the binding capability between drugs and target proteins, is a crucial indicator for drug development. Accurately predicting interaction strength between new drug-target pairs by analyzing previous experiments aids in screening potential drug molecules, repurposing them, and developing safe and effective medicines. Existing computational models for DTA prediction rely on strings or single-graph neural networks, lacking consideration of protein structure and molecular semantic information, leading to limited accuracy. Our experiments demonstrate that string-based methods may overlook protein conformations, causing a high root mean square error (RMSE) of 3.584 in affinity due to a lack of spatial context. Single graph networks also underperform on topology features, with a 6% lower confidence interval (CI) for activity classification. Absent semantic information also limits generalization across diverse compounds, resulting in 18% increment in RMSE and 5% in misclassifications within quantifications study, restricting potential drug discovery. To address these limitations, we propose G-K BertDTA, a novel framework for accurate DTA prediction incorporating protein features, molecular semantic features, and molecular structural information. In this proposed model, we represent drugs as graphs, with a GIN employed to learn the molecular topological information. For the extraction of protein structural features, we utilize a DenseNet architecture. A knowledge-based BERT semantic model is incorporated to obtain rich pre-trained semantic embeddings, thereby enhancing the feature information. We extensively evaluated our proposed approach on the publicly available benchmark datasets (i.e., KIBA and Davis), and experimental results demonstrate the promising performance of our method, which consistently outperforms previous state-of-the-art approaches. Code is available at https://github.com/AmbitYuki/G-K-BertDTA.

A Multi-Drug Concentration Gradient Mixing Chip: A Novel Platform for High-Throughput Drug Combination Screening

Combinatorial drug therapy has emerged as a critically important strategy in medical research and patient treatment and involves the use of multiple drugs in concert to achieve a synergistic effect. This approach can enhance therapeutic efficacy while simultaneously mitigating adverse side effects. However, the process of identifying optimal drug combinations, including their compositions and dosages, is often a complex, costly, and time-intensive endeavor. To surmount these hurdles, we propose a novel microfluidic device capable of simultaneously generating multiple drug concentration gradients across an interlinked array of culture chambers. This innovative setup allows for the real-time monitoring of live cell responses. With minimal effort, researchers can now explore the concentration-dependent effects of single-agent and combination drug therapies. Taking neural stem cells (NSCs) as a case study, we examined the impacts of various growth factors-epithelial growth factor (EGF), platelet-derived growth factor (PDGF), and fibroblast growth factor (FGF)-on the differentiation of NSCs. Our findings indicate that an overdose of any single growth factor leads to an upsurge in the proportion of differentiated NSCs. Interestingly, the regulatory effects of these growth factors can be modulated by the introduction of additional growth factors, whether singly or in combination. Notably, a reduced concentration of these additional factors resulted in a decreased number of differentiated NSCs. Our results affirm that the successful application of this microfluidic device for the generation of multi-drug concentration gradients has substantial potential to revolutionize drug combination screening. This advancement promises to streamline the process and accelerate the discovery of effective therapeutic drug combinations.

New Theobromine Apoptotic Analogue with Anticancer Potential Targeting the EGFR Protein: Computational and In Vitro Studies

AIM

The aim of this study was to design and examine a novel epidermal growth factor receptor (EGFR) inhibitor with apoptotic properties by utilizing the essential structural characteristics of existing EGFR inhibitors as a foundation.

METHOD

The study began with the natural alkaloid theobromine and developed a new semisynthetic derivative (T-1-PMPA). Computational ADMET assessments were conducted first to evaluate its anticipated safety and general drug-likeness. Deep density functional theory (DFT) computations were initially performed to validate the three-dimensional (3D) structure and reactivity of T-1-PMPA. Molecular docking against the EGFR proteins was conducted to investigate T-1-PMPA's binding affinity and inhibitory potential. Additional molecular dynamics (MD) simulations over 200 ns along with MM-GPSA, PLIP, and principal component analysis of trajectories (PCAT) experiments were employed to verify the binding and inhibitory properties of T-1-PMPA. Afterward, T-1-PMPA was semisynthesized to validate the proposed design and in silico findings through several in vitro examinations.

RESULTS

DFT studies indicated T-1-PMPA's reactivity using electrostatic potential, global reactive indices, and total density of states. Molecular docking, MD simulations, MM-GPSA, PLIP, and ED suggested the binding and inhibitory properties of T-1-PMPA against the EGFR protein. The in silico ADMET predicted T-1-PMPA's safety and general drug-likeness. In vitro experiments demonstrated that T-1-PMPA effectively inhibited EGFRWT and EGFR790m, with IC50 values of 86 and 561 nM, respectively, compared to Erlotinib (31 and 456 nM). T-1-PMPA also showed significant suppression of the proliferation of HepG2 and MCF7 malignant cell lines, with IC50 values of 3.51 and 4.13 μM, respectively. The selectivity indices against the two cancer cell lines indicated the overall safety of T-1-PMPA. Flow cytometry confirmed the apoptotic effects of T-1-PMPA by increasing the total percentage of apoptosis to 42% compared to 31, and 3% in Erlotinib-treated and control cells, respectively. The qRT-PCR analysis further supported the apoptotic effects by revealing significant increases in the levels of Casp3 and Casp9. Additionally, T-1-PMPA controlled the levels of TNFα and IL2 by 74 and 50%, comparing Erlotinib's values (84 and 74%), respectively.

CONCLUSION

In conclusion, our study's findings suggest the potential of T-1-PMPA as a promising apoptotic anticancer lead compound targeting the EGFR.

Computational biology approaches for drug repurposing

The drug discovery and development (DDD) process greatly relies on the data available in various forms to generate hypotheses for novel drug design. The complex and heterogeneous nature of biological data makes it difficult to utilize or gather meaningful information as such. Computational biology techniques have provided us with opportunities to better understand biological systems through refining and organizing large amounts of data into actionable and systematic purviews. The drug repurposing approach has been utilized to overcome the expansive time periods and costs associated with traditional drug development. It deals with discovering new uses of already approved drugs that have an established safety and efficacy profile, thereby, requiring them to go through fewer development phases. Thus, drug repurposing through computational biology provides a systematic approach to drug development and overcomes the constraints of traditional processes. The current chapter covers the basics, approaches and tools of computational biology that can be employed to effectively develop repurposing profile of already approved drug molecules.

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

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

Discovery of novel BRD4-BD2 inhibitors via in silico approaches: QSAR techniques, molecular docking, and molecular dynamics simulations

Bromodomain-containing protein 4(BRD4) plays an important role in the occurrence and development of various malignant tumors, which has attracted the attention of scientific research institutions and pharmaceutical companies. The structural modification of most currently available BRD4 inhibitors is relatively simple, but the drug effectiveness is limited. Research has found that the inhibition of BD1 may promote the differentiation of oligodendrocyte progenitor cell; however, the inhibition of BD2 will not cause this outcome. Therefore, newly potential drugs which target BRD4-BD2 need further research. Herein, we initially built QSAR models out of 49 compounds using HQSAR, CoMFA, CoMSIA, and Topomer CoMFA technology. All of the models have shown suitable reliabilities (q2 = 0.778, 0.533, 0.640, 0.702, respectively) and predictive abilities (r2pred = 0.716, 0.6289, 0.6153, 0.7968, respectively) for BRD4-BD2 inhibitors. On the basis of QSAR results and the search of the R-group in the topomer search module, we designed 20 new compounds with high activity that showed appropriate docking score and suitable ADMET. Docking studies and MD simulation were carried out to reveal the amino acid residues (Asn351, Cys347, Tyr350, Pro293, and Asp299) at the active site of BRD4-BD2. Free energy calculations and free energy landscapes verified the stable binding results and indicated stable conformations of the complexes. These theoretical studies provide guidance and theoretical basis for designing and developing novel BRD4-BD2 inhibitors.

Automatic generation of functional peptides with desired bioactivity and membrane permeability using Bayesian optimization

Peptides are potentially useful modalities of drugs; however, cell membrane permeability is an obstacle in peptide drug discovery. The identification of bioactive peptides for a therapeutic target is also challenging because of the huge amino acid sequence patterns of peptides. In this study, we propose a novel computational method, PEptide generation system using Neural network Trained on Amino acid sequence data and Gaussian process-based optimizatiON (PENTAGON), to automatically generate new peptides with desired bioactivity and cell membrane permeability. In the algorithm, we mapped peptide amino acid sequences onto the latent space constructed using a variational autoencoder and searched for peptides with desired bioactivity and cell membrane permeability using Bayesian optimization. We used our proposed method to generate peptides with cell membrane permeability and bioactivity for each of the nine therapeutic targets, such as the estrogen receptor (ER). Our proposed method outperformed a previously developed peptide generator in terms of similarity to known active peptide sequences and the length of generated peptide sequences.

Implementation of IFPTML Computational Models in Drug Discovery Against Flaviviridae Family

The Flaviviridae family consists of single-stranded positive-sense RNA viruses, which contains the genera Flavivirus, Hepacivirus, Pegivirus, and Pestivirus. Currently, there is an outbreak of viral diseases caused by this family affecting millions of people worldwide, leading to significant morbidity and mortality rates. Advances in computational chemistry have greatly facilitated the discovery of novel drugs and treatments for diseases associated with this family. Chemoinformatic techniques, such as the perturbation theory machine learning method, have played a crucial role in developing new approaches based on ML models that can effectively aid drug discovery. The IFPTML models have shown its capability to handle, classify, and process large data sets with high specificity. The results obtained from different models indicates that this methodology is proficient in processing the data, resulting in a reduction of the false positive rate by 4.25%, along with an accuracy of 83% and reliability of 92%. These values suggest that the model can serve as a computational tool in assisting drug discovery efforts and the development of new treatments against Flaviviridae family diseases.

Brain Organoids: A Game-Changer for Drug Testing

Neurological disorders are the second cause of death and the leading cause of disability worldwide. Unfortunately, no cure exists for these disorders, but the actual therapies are only able to ameliorate people's quality of life. Thus, there is an urgent need to test potential therapeutic approaches. Brain organoids are a possible valuable tool in the study of the brain, due to their ability to reproduce different brain regions and maturation stages; they can be used also as a tool for disease modelling and target identification of neurological disorders. Recently, brain organoids have been used in drug-screening processes, even if there are several limitations to overcome. This review focuses on the description of brain organoid development and drug-screening processes, discussing the advantages, challenges, and limitations of the use of organoids in modeling neurological diseases. We also highlighted the potential of testing novel therapeutic approaches. Finally, we examine the challenges and future directions to improve the drug-screening process.

Identification of bioactive compounds of Zanthoxylum armatum as potential inhibitor of pyruvate kinase M2 (PKM2): Computational and virtual screening approaches

PKM2 (Pyruvate kinase M2) is the isoform of pyruvate kinase which is known to catalyse the last step of glycolysis that is responsible for energy production. This specific isoform is known to be highly expressed in certain cancerous conditions. Considering the role of this protein in various cancer conditions, we used PKM2 as a target protein to identify the potential compounds against this target. In this study, we have examined 96 compounds of Zanthoxylum armatum using an array of computational and in silico tools. The compounds were assessed for toxicity then their anticancer potential was predicted. The virtual screening was done with molecular docking followed by a detailed examination using molecular dynamics simulation. The majority of the compounds showed a higher probability of being antineoplastic. Based on toxicity, predicted anticancer potential, binding affinity, and binding site, three compounds (nevadensin, asarinin, and kaempferol) were selected as hit compounds. The binding energy of these compounds with PKM2 ranged from -7.7 to -8.3 kcal/mol and all hit compounds interact at the active site of the protein. The selected hit compounds formed a stable complex with PKM2 when simulated under physiological conditions. The dynamic analysis showed that these compounds remained attached to the active site till the completion of molecular simulation. MM-PBSA analysis showed that nevadensin exhibited a higher affinity towards PKM2 compared to asarinin and kaempferol. These compounds need to be assessed properties in vivo and in vitro to validate their efficacy.

Genetic Algorithm-Based Receptor Ligand: A Genetic Algorithm-Guided Generative Model to Boost the Novelty and Drug-Likeness of Molecules in a Sampling Chemical Space

Deep learning-based de novo molecular design has recently gained significant attention. While numerous DL-based generative models have been successfully developed for designing novel compounds, the majority of the generated molecules lack sufficiently novel scaffolds or high drug-like profiles. The aforementioned issues may not be fully captured by commonly used metrics for the assessment of molecular generative models, such as novelty, diversity, and quantitative estimation of the drug-likeness score. To address these limitations, we proposed a genetic algorithm-guided generative model called GARel (genetic algorithm-based receptor-ligand interaction generator), a novel framework for training a DL-based generative model to produce drug-like molecules with novel scaffolds. To efficiently train the GARel model, we utilized dense net to update the parameters based on molecules with novel scaffolds and drug-like features. To demonstrate the capability of the GARel model, we used it to design inhibitors for three targets: AA2AR, EGFR, and SARS-Cov2. The results indicate that GARel-generated molecules feature more diverse and novel scaffolds and possess more desirable physicochemical properties and favorable docking scores. Compared with other generative models, GARel makes significant progress in balancing novelty and drug-likeness, providing a promising direction for the further development of DL-based de novo design methodology with potential impacts on drug discovery.

Unveiling the molecular interaction of hepatitis B virus inhibitor, entecavir with human serum albumin through computational, microscopic and spectroscopic approaches

Molecular docking, molecular dynamics (MD) simulation, atomic force microscopy (AFM) and multi-spectroscopic techniques were selected to unveil the molecular association between the hepatitis B virus (HBV) inhibitor, entecavir (ETR), and the major blood plasma transporter, human serum albumin (HSA). The entire docking and simulation analyses recognized ETR binding to subdomain IIA (Site I) of HSA through hydrogen bonds, hydrophobic and van der Waals forces while maintaining the complex's stability throughout the 100 ns. A gradual lessening in the Stern-Volmer quenching constant (K sv ) with rising temperatures registered ETR-induced quenching of HBV fluorescence as static quenching, thus advising complexation between ETR and HSA. The further advocation of this conclusion was seen from a larger value of the biomolecular quenching rate constant ((kq ) > 1010 M-1s-1), changes in the spectra (UV-Vis absorption) of HSA following ETR inclusion and ETR-induced swelling of HSA in the AFM results. The ETR appeared to bind to HSA with moderate affinity (K a = 1.87 - 1.19 × 10 4   M-1) at 290, 300 and 310 K. Significant alterations in the protein's secondary and tertiary structures, including changes in the protein's Tyr/Trp microenvironment, were also detected by circular dichroism and three-dimensional fluorescence spectra when the protein was bound to ETR. The findings of the drug displacement study backed the docking results of Site I as ETR's preferred binding site in HSA.Communicated by Ramaswamy H. Sarma.

In silico drug repurposing by combining machine learning classification model and molecular dynamics to identify a potential OGT inhibitor

O-linked N-acetylglucosamine (O-GlcNAc) is a unique intracellular post-translational glycosylation at the hydroxyl group of serine or threonine residues in nuclear, cytoplasmic and mitochondrial proteins. The enzyme O-GlcNAc transferase (OGT) is responsible for adding GlcNAc, and anomalies in this process can lead to the development of diseases associated with metabolic imbalance, such as diabetes and cancer. Repurposing approved drugs can be an attractive tool to discover new targets reducing time and costs in the drug design. This work focuses on drug repurposing to OGT targets by virtual screening of FDA-approved drugs through consensus machine learning (ML) models from an imbalanced dataset. We developed a classification model using docking scores and ligand descriptors. The SMOTE approach to resampling the dataset showed excellent statistical values in five of the seven ML algorithms to create models from the training set, with sensitivity, specificity and accuracy over 90% and Matthew's correlation coefficient greater than 0.8. The pose analysis obtained by molecular docking showed only H-bond interaction with the OGT C-Cat domain. The molecular dynamics simulation showed the lack of H-bond interactions with the C- and N-catalytic domains allowed the drug to exit the binding site. Our results showed that the non-steroidal anti-inflammatory celecoxib could be a potentially OGT inhibitor.

GTAMP-DTA: Graph transformer combined with attention mechanism for drug-target binding affinity prediction

Drug target affinity prediction (DTA) is critical to the success of drug development. While numerous machine learning methods have been developed for this task, there remains a necessity to further enhance the accuracy and reliability of predictions. Considerable bias in drug target binding prediction may result due to missing structural information or missing information. In addition, current methods focus only on simulating individual non-covalent interactions between drugs and proteins, thereby neglecting the intricate interplay among different drugs and their interactions with proteins. GTAMP-DTA combines special Attention mechanisms, assigning each atom or amino acid an attention vector. Interactions between drug forms and protein forms were considered to capture information about their interactions. And fusion transformer was used to learn protein characterization from raw amino acid sequences, which were then merged with molecular map features extracted from SMILES. A self-supervised pre-trained embedding that uses pre-trained transformers to encode drug and protein attributes is introduced in order to address the lack of labeled data. Experimental results demonstrate that our model outperforms state-of-the-art methods on both the Davis and KIBA datasets. Additionally, the model's performance undergoes evaluation using three distinct pooling layers (max-pooling, mean-pooling, sum-pooling) along with variations of the attention mechanism. GTAMP-DTA shows significant performance improvements compared to other methods.

3D-QSAR pharmacophore modeling, virtual screening, molecular docking, MD simulations, in vitro and in vivo studies to identify potential anti-hyperplasia drugs

Psoriasis is a common immune-mediated skin condition characterized by aberrant keratinocytes and cell proliferation. The purpose of this study was to explore the FDA-approved drugs by 3D-QSAR pharmacophore model and evaluate their efficiency by in-silico, in vitro, and in vivo psoriasis animal model. A 3D-QSAR pharmacophore model was developed by utilizing HypoGen algorithm using the structural features of 48 diaryl derivatives with diverse molecular patterns. The model was validated by a test set of 27 compounds, by cost analysis method, and Fischer's randomization test. The correlation coefficient of the best model (Hypo2) was 0.9601 for the training set while it was 0.805 for the test set. The selected model was taken as a 3D query for the virtual screening of over 3000 FDA-approved drugs. Compounds mapped with the pharmacophore model were further screened through molecular docking. The hits that showed the best docking results were screened through in silico skin toxicity approach. Top five hits were selected for the MD simulation studies. Based on MD simulations results, the best two hit molecules, that is, ebastine (Ebs) and mebeverine (Mbv) were selected for in vitro and in vivo antioxidant studies performed in mice. TNF-α and COX pro-inflammatory mediators, biochemical assays, histopathological analyses, and immunohistochemistry observations confirmed the anti-inflammatory response of the selected drugs. Based on these findings, it appeared that Ebs can effectively treat psoriasis-like skin lesions and down-regulate inflammatory responses which was consistent with docking predictions and could potentially be employed for further research on inflammation-related skin illnesses such as psoriasis.

De Novo Drug Design Using Transformer-Based Machine Translation and Reinforcement Learning of an Adaptive Monte Carlo Tree Search

The discovery of novel therapeutic compounds through de novo drug design represents a critical challenge in the field of pharmaceutical research. Traditional drug discovery approaches are often resource intensive and time consuming, leading researchers to explore innovative methods that harness the power of deep learning and reinforcement learning techniques. Here, we introduce a novel drug design approach called drugAI that leverages the Encoder-Decoder Transformer architecture in tandem with Reinforcement Learning via a Monte Carlo Tree Search (RL-MCTS) to expedite the process of drug discovery while ensuring the production of valid small molecules with drug-like characteristics and strong binding affinities towards their targets. We successfully integrated the Encoder-Decoder Transformer architecture, which generates molecular structures (drugs) from scratch with the RL-MCTS, serving as a reinforcement learning framework. The RL-MCTS combines the exploitation and exploration capabilities of a Monte Carlo Tree Search with the machine translation of a transformer-based Encoder-Decoder model. This dynamic approach allows the model to iteratively refine its drug candidate generation process, ensuring that the generated molecules adhere to essential physicochemical and biological constraints and effectively bind to their targets. The results from drugAI showcase the effectiveness of the proposed approach across various benchmark datasets, demonstrating a significant improvement in both the validity and drug-likeness of the generated compounds, compared to two existing benchmark methods. Moreover, drugAI ensures that the generated molecules exhibit strong binding affinities to their respective targets. In summary, this research highlights the real-world applications of drugAI in drug discovery pipelines, potentially accelerating the identification of promising drug candidates for a wide range of diseases.

Identifying novel inhibitors targeting Exportin-1 for the potential treatment of COVID-19

The nuclear export protein 1 (XPO1) mediates the nucleocytoplasmic transport of proteins and ribonucleic acids (RNAs) and plays a prominent role in maintaining cellular homeostasis. XPO1 has emerged as a promising therapeutic approach to interfere with the lifecycle of many viruses. In our earlier study, we proved the inhibition of XPO1 as a therapeutic strategy for managing SARS-COV-2 and its variants. In this study, we have utilized pharmacophore-assisted computational methods to identify prominent XPO1 inhibitors. After several layers of screening, a few molecules were shortlisted for further experimental validation on the in vitro SARS-CoV-2 cell infection model. It was observed that these compounds reduced spike positivity, suggesting inhibition of SARS-COV-2 infection. The outcome of this study could be considered further for developing novel antiviral therapeutic strategies against SARS-CoV-2.

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

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

Mutational analysis of SARS-CoV-2 ORF6-KPNA2 binding interface and identification of potent small molecule inhibitors to recuse the host immune system

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) surfaced on 31 December, 2019, and was identified as the causative agent of the global COVID-19 pandemic, leading to a pneumonia-like disease. One of its accessory proteins, ORF6, has been found to play a critical role in immune evasion by interacting with KPNA2 to antagonize IFN signaling and production pathways, resulting in the inhibition of IRF3 and STAT1 nuclear translocation. Since various mutations have been observed in ORF6, therefore, a comparative binding, biophysical, and structural analysis was used to reveal how these mutations affect the virus's ability to evade the human immune system. Among the identified mutations, the V9F, V24A, W27L, and I33T, were found to have a highly destabilizing effect on the protein structure of ORF6. Additionally, the molecular docking analysis of wildtype and mutant ORF6 and KPNA2 revealed the docking score of - 53.72 kcal/mol for wildtype while, -267.90 kcal/mol, -258.41kcal/mol, -254.51 kcal/mol and -268.79 kcal/mol for V9F, V24A, W27L, and I33T respectively. As compared to the wildtype the V9F showed a stronger binding affinity with KPNA2 which is further verified by the binding free energy (-42.28 kcal/mol) calculation. Furthermore, to halt the binding interface of the ORF6-KPNA2 complex, we used a computational molecular search of potential natural products. A multi-step virtual screening of the African natural database identified the top 5 compounds with best docking scores of -6.40 kcal/mol, -6.10 kcal/mol, -6.09 kcal/mol, -6.06 kcal/mol, and -6.03 kcal/mol for tophit1-5 respectively. Subsequent all-atoms simulations of these top hits revealed consistent dynamics, indicating their stability and their potential to interact effectively with the interface residues. In conclusion, our study represents the first attempt to establish a foundation for understanding the heightened infectivity of new SARS-CoV-2 variants and provides a strong impetus for the development of novel drugs against them.

Exploring the natural products chemical space to abrogate the F3L-dsRNA interface of monkeypox virus to enhance the immune responses using molecular screening and free energy calculations

Amid the ongoing monkeypox outbreak, there is an urgent need for the rapid development of effective therapeutic interventions capable of countering the immune evasion mechanisms employed by the monkeypox virus (MPXV). The evasion strategy involves the binding of the F3L protein to dsRNA, resulting in diminished interferon (IFN) production. Consequently, our current research focuses on utilizing virtual drug screening techniques to target the RNA binding domain of the F3L protein. Out of the 954 compounds within the South African natural compound database, only four demonstrated notable docking scores: -6.55, -6.47, -6.37, and -6.35 kcal/mol. The dissociation constant (KD) analysis revealed a stronger binding affinity of the top hits 1-4 (-5.34, -5.32, -5.29, and -5.36 kcal/mol) with the F3L in the MPXV. All-atom simulations of the top-ranked hits 1 to 4 consistently exhibited stable dynamics, suggesting their potential to interact effectively with interface residues. This was further substantiated through analyses of parameters such as radius of gyration (Rg), Root Mean Square Fluctuation, and hydrogen bonding. Cumulative assessments of binding free energy confirmed the top-performing candidates among all the compounds, with values of -35.90, -52.74, -28.17, and -32.11 kcal/mol for top hits 1-4, respectively. These results indicate that compounds top hit 1-4 could hold significant promise for advancing innovative drug therapies, suggesting their suitability for both in vivo and in vitro experiments.

Investigation and drug design for novel molecules from natural products as inhibitors for controlling multiple myeloma disease using in-silico tools

Multiple myeloma (MM) is a disease that causes plasma cell growth in the bone marrow and immune globulin buildup in blood and urine. Despite recent advances in MM therapy, many still die due to its high mortality rate. A study using computational simulations analyzed 100 natural ingredients from the SANC database to determine if they inhibited the IgH domain, a known cause of multiple myeloma. Natural component Diospyrin inhibited the IgH enzyme with the best binding energy of -10.3 kcal/mol and three carbon-hydrogen bonds, followed by Parviflorone F complex with a binding energy of -10.1 kcal/mol and two conventional-hydrogen bonds. As a result, the Molecular Dynamic simulation was used to test the stability of the two complexes. During the simulation, the Diospyrin molecule dissociated from the protein at roughly 67.5 ns, whereas the Parviflorone F molecule stayed attached to the protein throughout. The latter was the subject of the investigation. The analysis of the production run data revealed that the Parviflorone F molecule exhibits a variety of conformations within the binding pocket while keeping a relatively constant distance from the protein's center of mass. The analysis of the production run data revealed that the Parviflorone F molecule exhibited a variety of conformations within the binding pocket while keeping a relatively constant distance from the protein's center of mass. The root mean square deviation (RMSD) plots for both the protein and complex showed a stable and steady average value of 4.4 Å for the first 82 nanoseconds of manufacture. As a result, the average value increased to 8.3 Å. Furthermore, the components of the binding free energy, as computed by MM-GBSA, revealed that the mean binding energy of the Parviflorone F molecule was -23.88 kcal/mol. Finally, after analyzing all of the examination data, Parviflorone F was identified as a powerful inhibitor of the IgH domain and hence of the MM disease, which requires further in-vivo conformation.Communicated by Ramaswamy H. Sarma.

Identification of a Histone Deacetylase 8 Inhibitor through Drug Screenings Based on Machine Learning

Histone deacetylase 8 (HDAC8) is a zinc-dependent HDAC that catalyzes the deacetylation of nonhistone proteins. It is involved in cancer development and HDAC8 inhibitors are promising candidates as anticancer agents. However, most reported HDAC8 inhibitors contain a hydroxamic acid moiety, which often causes mutagenicity. Therefore, we used machine learning for drug screening and attempted to identify non-hydroxamic acids as HDAC8 inhibitors. In this study, we established a prediction model based on the random forest (RF) algorithm for screening HDAC8 inhibitors because it exhibited the best predictive accuracy in the training dataset, including data generated by the synthetic minority over-sampling technique (SMOTE). Using the trained RF-SMOTE model, we screened the Osaka University library for compounds and selected 50 virtual hits. However, the 50 hits in the first screening did not show HDAC8-inhibitory activity. In the second screening, using the RF-SMOTE model, which was established by retraining the dataset including 50 inactive compounds, we identified non-hydroxamic acid 12 as an HDAC8 inhibitor with an IC50 of 842 nM. Interestingly, its IC50 values for HDAC1 and HDAC3-inhibitory activity were 38 and 12 µM, respectively, showing that compound 12 has high HDAC8 selectivity. Using machine learning, we expanded the chemical space for HDAC8 inhibitors and identified non-hydroxamic acid 12 as a novel HDAC8 selective inhibitor.

Structure-based Virtual Screening from Natural Products as Inhibitors of SARS-CoV-2 Spike Protein and ACE2 Receptor Binding and their Biological Evaluation In vitro

BACKGROUND

In the last years, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused more than 760 million infections and 6.9 million deaths. Currently, remains a public health problem with limited pharmacological treatments. Among the virus drug targets, the SARS-CoV-2 spike protein attracts the development of new anti-SARS-CoV-2 agents.

OBJECTIVE

The aim of this work was to identify new compounds derived from natural products (BIOFACQUIM and Selleckchem databases) as potential inhibitors of the spike receptor binding domain (RBD)-ACE2 binding complex.

METHODS

Molecular docking, molecular dynamics simulations, and ADME-Tox analysis were performed to screen and select the potential inhibitors. ELISA-based enzyme assay was done to confirm our predictive model.

RESULTS

Twenty compounds were identified as potential binders of RBD of the spike protein. In vitro assay showed compound B-8 caused 48% inhibition at 50 μM, and their binding pattern exhibited interactions via hydrogen bonds with the key amino acid residues present on the RBD.

CONCLUSION

Compound B-8 can be used as a scaffold to develop new and more efficient antiviral drugs.

Elucidation of binding mechanism of rhodanine derivative P4OC on bovine serum albumin

Rhodanine is an important scaffold in medicinal chemistry and it act as potent anticancer agent and other pharmacological effects. In pharmacokinetics and pharmacodynamics studies of the drug, the drug binding properties on serum protein is crucial for producing better drug. This study was designed to explore the binding interactions between the Rhodanine derivative (P4OC) on Bovine Serum Albumin (BSA). The interactions between P4OC and BSA were investigated using biophysical approach and molecular docking. The quenching mechanism and binding constants of P4OC on BSA were determined by biophysical approach through fluorescence spectroscopic experiments. Circular dichroism (CD) spectroscopy was used to study the secondary structural changes of BSA upon P4OC binding. The fluorescence experiments of P4OC binding on BSA show good drug binding with static quenching constants using stern Volmer plot and found the quenching constant value KP4OC = 1.12762 × 1013 M-1 with corresponding binding free energy (ΔG) -2.303 kcal/mol. The molecular displacement fluorescence emission on BSA-P4OC complex by site specific markers shows that P4OC binds at I A sub-domain of BSA further confirmed peak shift by synchronous fluorescence of P4OC on BSA with tyrosine, tryptophan and phenylalanine amino acids. Increasing concentration of P4OC on BSA found secondary structural changes, the percentage of α-helix was decreased as well increase percentage of β-sheet and random coil. The binding of P4OC to BSA was computationally studied by molecular docking methods. Thus, results obtained are in excellent agreement with experimental and theoretical results with respect to the binding mechanism and binding constant of P4OC on BSA. We concluded that, the rhodanine derivative P4OC possesses good drug binding properties on BSA. Further P4OC may be evaluated its potential pharmacological activities on clinical trial.Communicated by Ramaswamy H. Sarma.

Recent Advances and Techniques for Identifying Novel Antibacterial Targets

BACKGROUND

With the emergence of drug-resistant bacteria, the development of new antibiotics is urgently required. Target-based drug discovery is the most frequently employed approach for the drug development process. However, traditional drug target identification techniques are costly and time-consuming. As research continues, innovative approaches for antibacterial target identification have been developed which enabled us to discover drug targets more easily and quickly.

METHODS

In this review, methods for finding drug targets from omics databases have been discussed in detail including principles, procedures, advantages, and potential limitations. The role of phage-driven and bacterial cytological profiling approaches is also discussed. Moreover, current article demonstrates the advancements being made in the establishment of computational tools, machine learning algorithms, and databases for antibacterial target identification.

RESULTS

Bacterial drug targets successfully identified by employing these aforementioned techniques are described as well.

CONCLUSION

The goal of this review is to attract the interest of synthetic chemists, biologists, and computational researchers to discuss and improve these methods for easier and quicker development of new drugs.

Machine learning in Alzheimer's disease drug discovery and target identification

Alzheimer's disease (AD) stands as a formidable neurodegenerative ailment that poses a substantial threat to the elderly population, with no known curative or disease-slowing drugs in existence. Among the vital and time-consuming stages in the drug discovery process, disease modeling and target identification hold particular significance. Disease modeling allows for a deeper comprehension of disease progression mechanisms and potential therapeutic avenues. On the other hand, target identification serves as the foundational step in drug development, exerting a profound influence on all subsequent phases and ultimately determining the success rate of drug development endeavors. Machine learning (ML) techniques have ushered in transformative breakthroughs in the realm of target discovery. Leveraging the strengths of large dataset analysis, multifaceted data processing, and the exploration of intricate biological mechanisms, ML has become instrumental in the quest for effective AD treatments. In this comprehensive review, we offer an account of how ML methodologies are being deployed in the pursuit of drug discovery for AD. Furthermore, we provide an overview of the utilization of ML in uncovering potential intervention strategies and prospective therapeutic targets for AD. Finally, we discuss the principal challenges and limitations currently faced by these approaches. We also explore the avenues for future research that hold promise in addressing these challenges.

Enhanced Sampling in Molecular Dynamics Simulations: How Many MD Snapshots can be Needed to Reproduce the Biological Behavior?

Since its early days in the 19th century, medicinal chemistry has concentrated its efforts on the treatment of diseases, using tools from areas such as chemistry, pharmacology, and molecular biology. The understanding of biological mechanisms and signaling pathways is crucial information for the development of potential agents for the treatment of diseases mainly because they are such complex processes. Given the limitations that the experimental approach presents, computational chemistry is a valuable alternative for the study of these systems and their behavior. Thus, classical molecular dynamics, based on Newton's laws, is considered a technique of great accuracy, when appropriated force fields are used, and provides satisfactory contributions to the scientific community. However, as many configurations are generated in a large MD simulation, methods such as Statistical Inefficiency and Optimal Wavelet Signal Compression Algorithm are great tools that can reduce the number of subsequent QM calculations. Accordingly, this review aims to briefly discuss the importance and relevance of medicinal chemistry allied to computational chemistry as well as to present a case study where, through a molecular dynamics simulation of AMPK protein (50 ns) and explicit solvent (TIP3P model), a minimum number of snapshots necessary to describe the oscillation profile of the protein behavior was proposed. For this purpose, the RMSD calculation, together with the sophisticated OWSCA method was used to propose the minimum number of snapshots.

The Discovery of Novel Agents against Staphylococcus aureus by Targeting Sortase A: A Combination of Virtual Screening and Experimental Validation

Staphylococcus aureus (S. aureus), commonly known as "superbugs", is a highly pathogenic bacterium that poses a serious threat to human health. There is an urgent need to replace traditional antibiotics with novel drugs to combat S. aureus. Sortase A (SrtA) is a crucial transpeptidase involved in the adhesion process of S. aureus. The reduction in virulence and prevention of S. aureus infections have made it a significant target for antimicrobial drugs. In this study, we combined virtual screening with experimental validation to identify potential drug candidates from a drug library. Three hits, referred to as Naldemedine, Telmisartan, and Azilsartan, were identified based on docking binding energy and the ratio of occupied functional sites of SrtA. The stability analysis manifests that Naldemedine and Telmisartan have a higher binding affinity to the hydrophobic pockets. Specifically, Telmisartan forms stable hydrogen bonds with SrtA, resulting in the highest binding energy. Our experiments prove that the efficiency of adhesion and invasion by S. aureus can be decreased without significantly affecting bacterial growth. Our work identifies Telmisartan as the most promising candidate for inhibiting SrtA, which can help combat S. aureus infection.

Annonaceous acetogenins as promising DNA methylation inhibitors to prevent and treat leukemogenesis - an in silico approach

Leukemia is a haematological malignancy affecting blood and bone marrow, ranking 10th among the other common cancers. DNA methylation is an epigenetic dysregulation that plays a critical role in leukemogenesis. DNA methyltransferases (DNMTs) such as DNMT1, DNMT3A and DNMT3B are the key enzymes catalysing DNA methylation. Inhibition of DNMT1 with secondary metabolites from medicinal plants helps reverse DNA methylation. The present study focuses on inhibiting DNMT1 protein (PDB ID: 3PTA) with annonaceous acetogenins through in-silico studies. The docking and molecular dynamic (MD) simulation study was carried out using Schrödinger Maestro and Desmond, respectively. These compounds' drug likeliness, ADMET properties and bioactivity scores were analysed. About 76 different acetogenins were chosen for this study, among which 17 showed the highest binding energy in the range of -8.312 to -10.266 kcal/mol. The compounds with the highest negative binding energy were found to be annohexocin (-10.266 kcal/mol), isoannonacinone (-10.209 kcal/mol) and annonacin (-9.839 kcal/mol). MD simulation results reveal that annonacin remains stable throughout the simulation time of 100 ns and also binds to the catalytic domain of DNMT1 protein. From the above results, it can be concluded that annonacin has the potential to inhibit the DNA methylation process and prevent leukemogenesis.Communicated by Ramaswamy H. Sarma.