Drug Design by Pharmacophore and Virtual Screening Approach

Structure-guided identification and characterization of potent inhibitors targeting PhoP and MtrA to combat mycobacteria

Mycobacteria are causative agents of tuberculosis (TB), which is a global health concern. Drug-resistant TB strains are rapidly emerging, thereby necessitating the urgent development of new drugs. Two-component signal transduction systems (TCSs) are signaling pathways involved in the regulation of various bacterial behaviors and responses to environmental stimuli. Applying specific inhibitors of TCSs can disrupt bacterial signaling, growth, and virulence, and can help combat drug-resistant TB. We conducted a comprehensive pharmacophore-based inhibitor screening and biochemical and biophysical examinations to identify, characterize, and validate potential inhibitors targeting the response regulators PhoP and MtrA of mycobacteria. The constructed pharmacophore model Phar-PR-n4 identified effective inhibitors of formation of the PhoP-DNA complex: ST132 (IC50 = 29 ± 1.6 µM) and ST166 (IC50 = 18 ± 1.3 µM). ST166 (KD = 18.4 ± 4.3 μM) and ST132 (KD = 14.5 ± 0.1 μM) strongly targeted PhoP in a slow-on, slow-off manner. The inhibitory potency and binding affinity of ST166 and ST132 for MtrAC were comparable to those of PhoP. Structural analyses and molecular dynamics simulations revealed that ST166 and ST132 mainly interact with the α8-helix and C-terminal β-hairpin of PhoP, with functionally essential residue hotspots for structure-based inhibitor optimization. Moreover, ST166 has in vitro antibacterial activity against Macrobacterium marinum. Thus, ST166, with its characteristic 1,2,5,6-tetrathiocane and terminal sulphonic groups, has excellent potential as a candidate for the development of novel antimicrobial agents to combat pathogenic mycobacteria.

Ligand-based pharmacophore modelling, structure optimisation, and biological evaluation for the identification of 2-heteroarylthio-N-arylacetamides as novel HSP90 C-terminal inhibitors

Targeting Heat shock protein 90 (HSP90) C-terminus is an important strategy to develop HSP90 inhibitors without inducing heat shock response. The development of C-terminal inhibitors, however, is hampered by a lack of understanding regarding the interaction between the HSP90 C-terminus and the present inhibitors. We collected seven classical and structurally diverse HSP90 C-terminal inhibitors and constructed a ligand-based pharmacophore model. The subsequent virtual screening and structural optimisation led to the identification of 2-heteroarylthio-N-arylacetamides as novel HSP90 C-terminal inhibitors. 9 and 27 exhibited strong antitumour activity in vitro by inhibiting proliferation and inducing apoptosis in multiple cancer cell lines. These compounds disrupted the interaction between HSP90 C-terminus and peptidylprolyl isomerase D, exerting a stronger inhibitory effect than novobiocin. 27 significantly induced the degradation of HSP90 clients without triggering heat shock response. In an in vivo study using 4T1 mice breast cancer models, 9 showed a potent antitumour effect without obvious toxicity.

Virtual screening and molecular growth guide the design of inhibitors for the influenza virus drug-resistant mutant M2-V27A/S31N

The influenza A virus matrix protein 2 (AM2) protein is a proton-gated, proton-selective ion channel essential for influenza replication that has been identified as an antiviral target. The drug-resistance of the M2-V27A/S31N strain, which has been growing more prevalent in recent years and has the potential to spread globally, prevents current amantadine inhibitors from having the desired impact. In this study, we compiled the most common influenza A virus strains from 2001-2020 from the U.S. National Center for Biotechnology Information database and hypothesized that M2-V27A/S31N would become a common strain. The lead compound ZINC299830590 was screened for M2-V27A/S31N in the ZINC15 database using a pharmacophore model and molecular descriptors. This lead compound was then optimized by molecular growth, which allowed us to identify important amino acid residues and create interactions with them to produce compound 4. Molecular dynamics simulation showed that the complex of compound 4 and M2-V27A/S31N had certain degrees of stability and flexibility. The binding free energy of compound 4 was calculated using the MM/PB(GB)SA method and totaled -106.525 kcal/mol. Finally, physicochemical and pharmacokinetic profiles were predicted using the Absorption, Distribution, Metabolism, Excretion, and Toxicity model, which indicated the good bioavailability of compound 4. These results provide the basis for further in vivo and in vitro studies to demonstrate that compound 4 is a promising drug candidate against M2-V27A/S31N.Communicated by Ramaswamy H. Sarma.

Pharmacophore-based virtual screening, molecular docking and molecular dynamics studies for the discovery of novel neuraminidase inhibitors

The in silico evaluation of 27 p-aminosalicylic acid derivatives, also referred to as neuraminidase inhibitors was the focus of the current study. To search and predict new potential neuraminidase inhibitors, this study was based on the ligand-based pharmacophore modeling, 3D QSAR, molecular docking, ADMET and MD simulation studies. The data was generated from recently reported inhibitors and divided into two groups, one of these group has 17 compounds for training and the second group has 10 compounds for testing purpose. The generated pharmacophore has known as ADDPR_4 was found statistically significant 3D-QSAR model owing the high trust scores (R2 = 0.974, Q2 = 0.905, RMSE = 0.23). Morever external validation was also employed to evaluate the prediction capacity of the built pharmacophore model (R2pred = 0.905). In addition, in silico ADMET, analyses were employed to evaluate the obtained hits for drug likeness properties. The stability of formed complexes was further evaluated using molecular dynamics. Top two hits showed stable complexes with Neuraminidase based on calculated total binding energy by MM-PBSA.Communicated by Ramaswamy H. Sarma.

Multi-target drugs for Alzheimer's disease

Alzheimer's disease (AD), a leading cause of dementia, increasingly challenges our healthcare systems and society. Traditional therapies aimed at single targets have fallen short owing to the complex, multifactorial nature of AD that necessitates simultaneous targeting of various disease mechanisms for clinical success. Therefore, targeting multiple pathologies at the same time could provide a synergistic therapeutic effect. The identification of new disease targets beyond the classical hallmarks of AD offers a fertile ground for the design of new multi-target drugs (MTDs), and building on existing compounds have the potential to yield in successful disease modifying therapies. This review discusses the evolving landscape of MTDs, focusing on their potential as AD therapeutics. Analysis of past and current trials of compounds with multi-target activity underscores the capacity of MTDs to offer synergistic therapeutic effects, and the flourishing genetic understanding of AD will inform and inspire the development of MTD-based AD therapies.

In silico assessment of biocompatibility and toxicity: molecular docking and dynamics simulation of PMMA-based dental materials for interim prosthetic restorations

AIM

This study aimed to comprehensively assess the biocompatibility and toxicity profiles of poly(methyl methacrylate) (PMMA) and its monomeric unit, methyl methacrylate (MMA), crucial components in dental materials for interim prosthetic restorations.

METHODOLOGY

Molecular docking was employed to predict the binding affinities, energetics, and steric features of MMA and PMMA with selected receptors involved in bone metabolism and tissue development, including RANKL, Fibronectin, BMP9, NOTCH2, and other related receptors. The HADDOCK standalone version was utilized for docking calculations, employing a Lamarckian genetic algorithm to explore the conformational space of ligand-receptor interactions. Furthermore, molecular dynamics (MD) simulations over 100 nanoseconds were conducted using the GROMACS package to evaluate dynamic actions and structural stability. The LigandScout was utilized for pharmacophore modeling, which employs a shape-based screening approach to identify potential ligand binding sites on protein targets.

RESULTS

The molecular docking studies elucidated promising interactions between PMMA and MMA with key biomolecular targets relevant to dental applications. MD simulation results provided strong evidence supporting the structural stability of PMMA complexes over time. Pharmacophore modeling highlighted the significance of carbonyl and hydroxyl groups as pharmacophoric features, indicating compounds with favorable biocompatibility profiles.

CONCLUSION

This study underscores the potential of PMMA in dental applications, emphasizing its structural stability, molecular interactions, and safety considerations. These findings lay a foundation for future advancements in dental biomaterials, guiding the design and optimization of materials for enhanced biocompatibility. Future directions include experimental validation of computational findings and the development of PMMA-based dental materials with improved biocompatibility and clinical performance.

Unveiling G-protein coupled receptor kinase-5 inhibitors for chronic degenerative diseases: Multilayered prioritization employing explainable machine learning-driven multi-class QSAR, ligand-based pharmacophore and free energy-inspired molecular simulation

GRK5 holds a pivotal role in cellular signaling pathways, with its overexpression in cardiomyocytes, neuronal cells, and tumor cells strongly associated with various chronic degenerative diseases, which highlights the urgent need for potential inhibitors. In this study, multiclass classification-based QSAR models were developed using diverse machine learning algorithms. These models were built from curated compounds with experimentally derived GRK5 inhibitory activity. Additionally, a pharmacophore model was constructed using active compounds from the dataset. Among the models, the SVM-based approach proved most effective and was initially used to screen DrugBank compounds within the applicability domain. Compounds showing significant GRK5 inhibitory potential underwent evaluation for key pharmacophoric features. Prospective compounds were subjected to molecular docking to assess binding affinity towards GRK5's key active site amino acid residues. Stability at the binding site was analyzed through 200 ns molecular dynamics simulations. MM-GBSA analysis quantified individual free energy components contributing to the total binding energy with respect to binding site residues. Metadynamics analysis, including PCA, FEL, and PDF, provided crucial insights into conformational changes of both apo and holo forms of GRK5 at defined energy states. The study identifies DB02844 (S-Adenosyl-1,8-Diamino-3-Thiooctane) and DB13155 (Esculin) as promising GRK5 inhibitors, warranting further in vitro and in vivo validation studies.

Unraveling potential neuroprotective mechanisms of herbal medicine for Alzheimer's diseases through comprehensive molecular docking analyses

Alzheimer's disease (AD) continues to be a worldwide health concern, demanding innovative therapeutic approaches. This study investigates the neuroprotective potential of herbal compounds by scrutinizing their interactions with Beta-Secretase-1 (BACE1). Through comprehensive molecular docking analyses, three compounds, Masticadienonic acid (ΔG: -9.6 kcal/mol), Hederagenin (ΔG: -9.3 kcal/mol), and Anthocyanins (ΔG: -8.1 kcal/mol), emerge as promising BACE1 ligands, displaying low binding energies and strong affinities. ADME parameter predictions, drug-likeness assessments, and toxicity analyses reveal favorable pharmacokinetic profiles for these compounds. Notably, Masticadienonic Acid exhibits optimal drug-likeness (-3.3736) and negligible toxicity concerns. Hederagenin (drug-likeness: -5.3272) and Anthocyanins (drug-likeness: -6.2041) also demonstrate promising safety profiles. Furthermore, pharmacophore modeling elucidates the compounds' unique interaction landscapes within BACE1's active site. Masticadienonic acid showcases seven hydrophobic interactions and a hydrogen bond acceptor interaction with Thr232. Hederagenin exhibits a specific hydrogen bond acceptor interaction with Trp76, emphasizing its selective binding. Anthocyanins reveal a multifaceted engagement, combining hydrophobic contacts and hydrogen bond interactions with key residues. In conclusion, Masticadienonic acid, Hederagenin, and Anthocyanins stand out as promising candidates for further experimental validation, presenting a synergistic balance of efficacy and safety in combating AD through BACE1 inhibition.

Synthesis, Biological and Molecular Docking Studies of Thiazole-Thiadiazole derivatives as potential Anti-Tuberculosis Agents

Tuberculosis remains a global health threat, with increasing infection rates and mortality despite existing anti-TB drugs. The present work focuses on the research findings regarding the development and evaluation of thiadiazole-linked thiazole derivatives as potential anti-tuberculosis agents. We present the synthesis data and confirm the compound structures using spectroscopic techniques. The current study reports twelve thiazole-thiadiazole compounds (5 a-5 l) for their anti-tuberculosis and related bioactivities. This paper emphasizes compounds 5 g, 5 i, and 5 l, which exhibited promising MIC values, leading to further in silico and interaction analysis. Pharmacophore mapping data included in the present analysis identified tubercular ThyX as potential drug targets. The compounds were evaluated for anti-tubercular activity using standard methods, revealing significant MIC values, particularly compound 5 l, with the best MIC value of 7.1285 μg/ml. Compounds 5 g and 5 i also demonstrated moderate to good MIC values against M. tuberculosis (H37Ra). Structural inspection of the docked poses revealed interactions such as hydrogen bonds, halogen bonds, and interactions containing Pi electron cloud, shedding light on conserved interactions with residues like Arg 95, Cys 43, His 69, and Arg 87 from the tubercular ThyX enzyme.

Unveiling potential repurposed drug candidates for Plasmodium falciparum through in silico evaluation: A synergy of structure-based approaches, structure prediction, and molecular dynamics simulations

The rise of drug resistance in Plasmodium falciparum, rendering current treatments ineffective, has hindered efforts to eliminate malaria. To address this issue, the study employed a combination of Systems Biology approach and a structure-based pharmacophore method to identify a target against P. falciparum. Through text mining, 448 genes were extracted, and it was discovered that plasmepsins, found in the Plasmodium genus, play a crucial role in the parasite's survival. The metabolic pathways of these proteins were determined using the PlasmoDB genomic database and recreated using CellDesigner 4.4.2. To identify a potent target, Plasmepsin V (PF13_0133) was selected and examined for protein-protein interactions (PPIs) using the STRING Database. Topological analysis and global-based methods identified PF13_0133 as having the highest centrality. Moreover, the static protein knockout PPIs demonstrated the essentiality of PF13_0133 in the modeled network. Due to the unavailability of the protein's crystal structure, it was modeled and subjected to a molecular dynamics simulation study. The structure-based pharmacophore modeling utilized the modeled PF13_0133 (PfPMV), generating 10 pharmacophore hypotheses with a library of active and inactive compounds against PfPMV. Through virtual screening, two potential candidates, hesperidin and rutin, were identified as potential drugs which may be repurposed as potential anti-malarial agents.

Insights into the computer-aided drug design and discovery based on anthraquinone scaffold for cancer treatment: A systematic review

BACKGROUND

In the search for better anticancer drugs, computer-aided drug design (CADD) techniques play an indispensable role in facilitating the lengthy and costly drug discovery process especially when natural products are involved. Anthraquinone is one of the most widely-recognized natural products with anticancer properties. This review aimed to systematically assess and synthesize evidence on the utilization of CADD techniques centered on the anthraquinone scaffold for cancer treatment.

METHODS

The conduct and reporting of this review were done in accordance to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA) 2020 guideline. The protocol was registered in the "International prospective register of systematic reviews" database (PROSPERO: CRD42023432904) and also published recently. The search strategy was designed based on the combination of concept 1 "CADD or virtual screening", concept 2 "anthraquinone" and concept 3 "cancer". The search was executed in PubMed, Scopus, Web of Science and MedRxiv on 30 June 2023.

RESULTS

Databases searching retrieved a total of 317 records. After deduplication and applying the eligibility criteria, the final review ended up with 32 articles in which 3 articles were found by citation searching. The CADD methods used in the studies were either structure-based alone (69%) or combined with ligand-based methods via parallel (9%) or sequential (22%) approaches. Molecular docking was performed in all studies, with Glide and AutoDock being the most popular commercial and public software used respectively. Protein data bank was used in most studies to retrieve the crystal structure of the targets of interest while the main ligand databases were PubChem and Zinc. The utilization of in-silico techniques has enabled a deeper dive into the structural, biological and pharmacological properties of anthraquinone derivatives, revealing their remarkable anticancer properties in an all-rounded fashion.

CONCLUSION

By harnessing the power of computational tools and leveraging the natural diversity of anthraquinone compounds, researchers can expedite the development of better drugs to address the unmet medical needs in cancer treatment by improving the treatment outcome for cancer patients.

Anti-Obesity Therapeutic Targets Studied In Silico and In Vivo: A Systematic Review

In the age of information technology and the additional computational search tools and software available, this systematic review aimed to identify potential therapeutic targets for obesity, evaluated in silico and subsequently validated in vivo. The systematic review was initially guided by the research question "What therapeutic targets have been used in in silico analysis for the treatment of obesity?" and structured based on the acronym PECo (P, problem; E, exposure; Co, context). The systematic review protocol was formulated and registered in PROSPERO (CRD42022353808) in accordance with the Preferred Reporting Items Checklist for Systematic Review and Meta-Analysis Protocols (PRISMA-P), and the PRISMA was followed for the systematic review. The studies were selected according to the eligibility criteria, aligned with PECo, in the following databases: PubMed, ScienceDirect, Scopus, Web of Science, BVS, and EMBASE. The search strategy yielded 1142 articles, from which, based on the evaluation criteria, 12 were included in the systematic review. Only seven these articles allowed the identification of both in silico and in vivo reassessed therapeutic targets. Among these targets, five were exclusively experimental, one was exclusively theoretical, and one of the targets presented an experimental portion and a portion obtained by modeling. The predominant methodology used was molecular docking and the most studied target was Human Pancreatic Lipase (HPL) (n = 4). The lack of methodological details resulted in more than 50% of the papers being categorized with an "unclear risk of bias" across eight out of the eleven evaluated criteria. From the current systematic review, it seems evident that integrating in silico methodologies into studies of potential drug targets for the exploration of new therapeutic agents provides an important tool, given the ongoing challenges in controlling obesity.

A Computational Approach for the Discovery of Novel DNA Methyltransferase Inhibitors

Nowadays, the explosion of knowledge in the field of epigenetics has revealed new pathways toward the treatment of multifactorial diseases, rendering the key players of the epigenetic machinery the focus of today's pharmaceutical landscape. Among epigenetic enzymes, DNA methyltransferases (DNMTs) are first studied as inhibition targets for cancer treatment. The increasing clinical interest in DNMTs has led to advanced experimental and computational strategies in the search for novel DNMT inhibitors. Considering the importance of epigenetic targets as a novel and promising pharmaceutical trend, the present study attempted to discover novel inhibitors of natural origin against DNMTs using a combination of structure and ligand-based computational approaches. Particularly, a pharmacophore-based virtual screening was performed, followed by molecular docking and molecular dynamics simulations in order to establish an accurate and robust selection methodology. Our screening protocol prioritized five natural-derived compounds, derivatives of coumarins, flavones, chalcones, benzoic acids, and phenazine, bearing completely diverse chemical scaffolds from FDA-approved "Epi-drugs". Their total DNMT inhibitory activity was evaluated, revealing promising results for the derived hits with an inhibitory activity ranging within 30-45% at 100 µM of the tested compounds.

Integration of the Butina algorithm and ensemble learning strategies for the advancement of a pharmacophore ligand-based model: an in silico investigation of apelin agonists

Introduction: 3D pharmacophore models describe the ligand's chemical interactions in their bioactive conformation. They offer a simple but sophisticated approach to decipher the chemically encoded ligand information, making them a valuable tool in drug design. Methods: Our research summarized the key studies for applying 3D pharmacophore models in virtual screening for 6,944 compounds of APJ receptor agonists. Recent advances in clustering algorithms and ensemble methods have enabled classical pharmacophore modeling to evolve into more flexible and knowledge-driven techniques. Butina clustering categorizes molecules based on their structural similarity (indicated by the Tanimoto coefficient) to create a structurally diverse training dataset. The learning method combines various individual pharmacophore models into a set of pharmacophore models for pharmacophore space optimization in virtual screening. Results: This approach was evaluated on Apelin datasets and afforded good screening performance, as proven by Receiver Operating Characteristic (AUC score of 0.994 ± 0.007), enrichment factor of (EF1% of 50.07 ± 0.211), Güner-Henry score of 0.956 ± 0.015, and F-measure of 0.911 ± 0.031. Discussion: Although one of the high-scoring models achieved statistically superior results in each dataset (AUC of 0.82; an EF1% of 19.466; GH of 0.131 and F1-score of 0.071), the ensemble learning method including voting and stacking method balanced the shortcomings of each model and passed with close performance measures.

A combination of virtual screening, molecular dynamics simulation, MM/PBSA, ADMET, and DFT calculations to identify a potential DPP4 inhibitor

DPP4 inhibitors can control glucose homeostasis by increasing the level of GLP-1 incretins hormone due to dipeptidase mimicking. Despite the potent effects of DPP4 inhibitors, these compounds cause unwanted toxicity attributable to their effect on other enzymes. As a result, it seems essential to find novel and DPP4 selective compounds. In this study, we introduce a potent and selective DPP4 inhibitor via structure-based virtual screening, molecular docking, molecular dynamics simulation, MM/PBSA calculations, DFT analysis, and ADMET profile. The screened compounds based on similarity with FDA-approved DPP4 inhibitors were docked towards the DPP4 enzyme. The compound with the highest docking score, ZINC000003015356, was selected. For further considerations, molecular docking studies were performed on selected ligands and FDA-approved drugs for DPP8 and DPP9 enzymes. Molecular dynamics simulation was run during 200 ns and the analysis of RMSD, RMSF, Rg, PCA, and hydrogen bonding were performed. The MD outputs showed stability of the ligand-protein complex compared to available drugs in the market. The total free binding energy obtained for the proposed DPP4 inhibitor was more negative than its co-crystal ligand (N7F). ZINC000003015356 confirmed the role of the five Lipinski rule and also, have low toxicity parameter according to properties. Finally, DFT calculations indicated that this compound is sufficiently soft.

A robust computational quest: Discovering potential hits to improve the treatment of pyrazinamide-resistant Mycobacterium tuberculosis

The rise of pyrazinamide (PZA)-resistant strains of Mycobacterium tuberculosis (MTB) poses a major challenge to conventional tuberculosis (TB) treatments. PZA, a cornerstone of TB therapy, must be activated by the mycobacterial enzyme pyrazinamidase (PZase) to convert its active form, pyrazinoic acid, which targets the ribosomal protein S1. Resistance, often associated with mutations in the RpsA protein, complicates treatment and highlights a critical gap in the understanding of structural dynamics and mechanisms of resistance, particularly in the context of the G97D mutation. This study utilizes a novel integration of computational techniques, including multiscale biomolecular and molecular dynamics simulations, physicochemical and medicinal chemistry predictions, quantum computations and virtual screening from the ZINC and Chembridge databases, to elucidate the resistance mechanism and identify lead compounds that have the potential to improve treatment outcomes for PZA-resistant MTB, namely ZINC15913786, ZINC20735155, Chem10269711, Chem10279789 and Chem10295790. These computational methods offer a cost-effective, rapid alternative to traditional drug trials by bypassing the need for organic subjects while providing highly accurate insight into the binding sites and efficacy of new drug candidates. The need for rapid and appropriate drug development emphasizes the need for robust computational analysis to justify further validation through in vitro and in vivo experiments.

Deciphering the Lexicon of Protein Targets: A Review on Multifaceted Drug Discovery in the Era of Artificial Intelligence

Understanding protein sequence and structure is essential for understanding protein-protein interactions (PPIs), which are essential for many biological processes and diseases. Targeting protein binding hot spots, which regulate signaling and growth, with rational drug design is promising. Rational drug design uses structural data and computational tools to study protein binding sites and protein interfaces to design inhibitors that can change these interactions, thereby potentially leading to therapeutic approaches. Artificial intelligence (AI), such as machine learning (ML) and deep learning (DL), has advanced drug discovery and design by providing computational resources and methods. Quantum chemistry is essential for drug reactivity, toxicology, drug screening, and quantitative structure-activity relationship (QSAR) properties. This review discusses the methodologies and challenges of identifying and characterizing hot spots and binding sites. It also explores the strategies and applications of artificial-intelligence-based rational drug design technologies that target proteins and protein-protein interaction (PPI) binding hot spots. It provides valuable insights for drug design with therapeutic implications. We have also demonstrated the pathological conditions of heat shock protein 27 (HSP27) and matrix metallopoproteinases (MMP2 and MMP9) and designed inhibitors of these proteins using the drug discovery paradigm in a case study on the discovery of drug molecules for cancer treatment. Additionally, the implications of benzothiazole derivatives for anticancer drug design and discovery are deliberated.

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

Sunset Yellow protects against oxidative damage and exhibits chemoprevention in chemically induced skin cancer model

Skin cancer and other skin-related inflammatory pathologies are rising due to heightened exposure to environmental pollutants and carcinogens. In this context, natural products and repurposed compounds hold promise as novel therapeutic and preventive agents. Strengthening the skin's antioxidant defense mechanisms is pivotal in neutralizing reactive oxygen species (ROS) and mitigating oxidative stress. Sunset Yellow (SY) exhibits immunomodulatory characteristics, evidenced by its capacity to partially inhibit the secretion of proinflammatory cytokines, regulate immune cell populations, and modulate the activation of lymphocytes. This study aimed to investigate the antioxidant and anti-genotoxic properties of SY using in-silico, in vitro, and physiochemical test systems, and to further explore its potential role in 7,12-dimethylbenz(a) anthracene (DMBA)/ 12-o-tetradecanoylphorbol-13-acetate (TPA)-induced two-stage skin carcinogenesis. In vitro experiments showed that pre-treatment of SY significantly enhanced the cell viability of HaCaT cells when exposed to tertiary-Butyl Hydrogen Peroxide (tBHP). This increase was accompanied by reduced ROS levels, restoration of mitochondrial membrane potential, and notable reduction in DNA damage in (SY + tBHP) treated cells. Mechanistic investigations using DPPH chemical antioxidant activity test and potentiometric titrations confirmed SY's antioxidant properties, with a standard reduction potential (E o ) of 0.211 V. Remarkably, evaluating the effect of topical application of SY in DMBA/TPA-induced two-step skin carcinogenesis model revealed dose-dependent decreases in tumor latency, incidence, yield, and burden over 21-weeks. Furthermore, computational analysis and experimental validations identified GSK3β, KEAP1 and EGFR as putative molecular targets of SY. Collectively, our findings reveal that SY enhances cellular antioxidant defenses, exhibits anti-genotoxic effects, and functions as a promising chemopreventive agent.

Multi-method computational evaluation of the inhibitors against leucine-rich repeat kinase 2 G2019S mutant for Parkinson's disease

Leucine-rich repeat kinase 2 G2019S mutant (LRRK2 G2019S) is a potential target for Parkinson's disease therapy. In this work, the computational evaluation of the LRRK2 G2019S inhibitors was conducted via a combined approach which contains a preliminary screening of a large database of compounds via similarity and pharmacophore, a secondary selection via structure-based affinity prediction and molecular docking, and a rescoring treatment for the final selection. MD simulations and MM/GBSA calculations were performed to check the agreement between different prediction methods for these inhibitors. 331 experimental ligands were collected, and 170 were used to build the structure-activity relationship. Eight representative ligand structural models were employed in similarity searching and pharmacophore screening over 14 million compounds. The process for selecting proper molecular descriptors provides a successful sample which can be used as a general strategy in QSAR modelling. The rescoring used in this work presents an alternative useful treatment for ranking and selection.

Computational resources and chemoinformatics for translational health research

The integration of computational resources and chemoinformatics has revolutionized translational health research. It has offered a powerful set of tools for accelerating drug discovery. This chapter overviews the computational resources and chemoinformatics methods used in translational health research. The resources and methods can be used to analyze large datasets, identify potential drug candidates, predict drug-target interactions, and optimize treatment regimens. These resources have the potential to transform the drug discovery process and foster personalized medicine research. We discuss insights into their various applications in translational health and emphasize the need for addressing challenges, promoting collaboration, and advancing the field to fully realize the potential of these tools in transforming healthcare.

Identification of novel compounds and repurposing of FDA drugs for 1-deoxy-D-xylulose 5-phosphate reductoisomerase enzyme of Plasmodium falciparum to combat malaria resistance

The rise of Plasmodium falciparum resistance to Artemisinin-based combination therapies (ACTs) is a significant concern in the fight against malaria. This situation calls for the search for novel anti-malarial candidates. 1-deoxy-D-xylulose 5-phosphate reductoisomerase (IspC) is a potential target involved in various cellular processes in P. falciparum (Pf). We screened ∼0.69 billion novel compounds from the ZINC20 library and repurposed ∼1400 FDA drugs using computational drug discovery methods against PfIspC. Following our computational pipeline, we found five novel ZINC20 compounds (Z-2, Z-3, Z-10, Z-13, and Z-14) and three FDA drugs (Aliskiren, Ceftolozane, and Ombitasvir) that showed striking docking energy (ranging from -8.405 to -10.834 kcal/mol), and strong interactions with key binding site residues (Ser269, Ser270, Ser306, Asn311, Lys312, and Met360) of PfIspC. The novel anti-malarial compounds also exhibited favorable pharmacokinetics and physicochemical properties. Furthermore, through molecular dynamics simulation, we observed the stable dynamics of PfIspC-inhibitor complexes and the influence of inhibitor binding on the protein's conformational arrangements. Notably, the binding free energy estimation confirmed high binding affinity (varied from -11.68 to -33.16 kcal/mol) of these compounds for PfIspC. Our findings could contribute to the ongoing efforts in combating malaria and invite experimental-lab researchers for validation.

Inhibitor design for TMPRSS2: insights from computational analysis of its backbone hydrogen bonds using a simple descriptor

Transmembrane protease serine 2 (TMPRSS2) is an important drug target due to its role in the infection mechanism of coronaviruses including SARS-CoV-2. Current understanding regarding the molecular mechanisms of known inhibitors and insights required for inhibitor design are limited. This study investigates the effect of inhibitor binding on the intramolecular backbone hydrogen bonds (BHBs) of TMPRSS2 using the concept of hydrogen bond wrapping, which is the phenomenon of stabilization of a hydrogen bond in a solvent environment as a result of being surrounded by non-polar groups. A molecular descriptor which quantifies the extent of wrapping around BHBs is introduced for this. First, virtual screening for TMPRSS2 inhibitors is performed by molecular docking using the program DOCK 6 with a Generalized Born surface area (GBSA) scoring function. The docking results are then analyzed using this descriptor and its relationship to the solvent-accessible surface area term ΔGsa of the GBSA score is demonstrated with machine learning regression and principal component analysis. The effect of binding of the inhibitors camostat, nafamostat, and 4-guanidinobenzoic acid (GBA) on the wrapping of important BHBs in TMPRSS2 is also studied using molecular dynamics. For BHBs with a large increase in wrapping groups due to these inhibitors, the radial distribution function of water revealed that certain residues involved in these BHBs, like Gln438, Asp440, and Ser441, undergo preferential desolvation. The findings offer valuable insights into the mechanisms of these inhibitors and may prove useful in the design of new inhibitors.

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

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

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.

Characterization of prevalent tyrosine kinase inhibitors and their challenges in glioblastoma treatment

Glioblastoma multiforme (GBM) is a highly aggressive malignant primary tumor in the central nervous system. Despite extensive efforts in radiotherapy, chemotherapy, and neurosurgery, there remains an inadequate level of improvement in treatment outcomes. The development of large-scale genomic and proteomic analysis suggests that GBMs are characterized by transcriptional heterogeneity, which is responsible for therapy resistance. Hence, knowledge about the genetic and epigenetic heterogeneity of GBM is crucial for developing effective treatments for this aggressive form of brain cancer. Tyrosine kinases (TKs) can act as signal transducers, regulate important cellular processes like differentiation, proliferation, apoptosis and metabolism. Therefore, TK inhibitors (TKIs) have been developed to specifically target these kinases. TKIs are categorized into allosteric and non-allosteric inhibitors. Irreversible inhibitors form covalent bonds, which can lead to longer-lasting effects. However, this can also increase the risk of off-target effects and toxicity. The development of TKIs as therapeutics through computer-aided drug design (CADD) and bioinformatic techniques enhance the potential to improve patients' survival rates. Therefore, the continued exploration of TKIs as drug targets is expected to lead to even more effective and specific therapeutics in the future.

Targeting ion channels with ultra-large library screening for hit discovery

Ion channels play a crucial role in a variety of physiological and pathological processes, making them attractive targets for drug development in diseases such as diabetes, epilepsy, hypertension, cancer, and chronic pain. Despite the importance of ion channels in drug discovery, the vastness of chemical space and the complexity of ion channels pose significant challenges for identifying drug candidates. The use of in silico methods in drug discovery has dramatically reduced the time and cost of drug development and has the potential to revolutionize the field of medicine. Recent advances in computer hardware and software have enabled the screening of ultra-large compound libraries. Integration of different methods at various scales and dimensions is becoming an inevitable trend in drug development. In this review, we provide an overview of current state-of-the-art computational chemistry methodologies for ultra-large compound library screening and their application to ion channel drug discovery research. We discuss the advantages and limitations of various in silico techniques, including virtual screening, molecular mechanics/dynamics simulations, and machine learning-based approaches. We also highlight several successful applications of computational chemistry methodologies in ion channel drug discovery and provide insights into future directions and challenges in this field.

Advances in Drug Discovery and Design using Computer-aided Molecular Modeling

Computer-aided molecular modeling is a rapidly emerging technology that is being used to accelerate the discovery and design of new drug therapies. It involves the use of computer algorithms and 3D structures of molecules to predict interactions between molecules and their behavior in the body. This has drastically improved the speed and accuracy of drug discovery and design. Additionally, computer-aided molecular modeling has the potential to reduce costs, increase the quality of data, and identify promising targets for drug development. Through the use of sophisticated methods, such as virtual screening, molecular docking, pharmacophore modeling, and quantitative structure-activity relationships, scientists can achieve higher levels of efficacy and safety for new drugs. Moreover, it can be used to understand the activity of known drugs and simplify the process of formulating, optimizing, and predicting the pharmacokinetics of new and existing drugs. In conclusion, computer-aided molecular modeling is an effective tool to rapidly progress drug discovery and design by predicting the interactions between molecules and anticipating the behavior of new drugs in the body.

Design, synthesis, and pharmacological evaluation of indazole carboxamides of N-substituted pyrrole derivatives as soybean lipoxygenase inhibitors

In this paper, we attempted to develop a novel class of compounds against lipoxygenase, a key enzyme in the biosynthesis of leukotrienes implicated in a series of inflammatory diseases. Given the absence of appropriate human 5-lipoxygenase crystallographic data, solved soybean lipoxygenase-1 and -3 structures were used as a template to generate an accurate pharmacophore model which was further used for virtual screening purposes. Eight compounds (1-8) have been derived from the in-house library consisting of N-substituted pyrroles conjugated with 5- or 6-indazole moieties through a carboxamide linker. This study led to the discovery of hit molecule 8 bearing a naphthyl group with the IC50 value of 22 μM according to soybean lipoxygenase in vitro assay. Isosteric replacement of naphthyl ring with quinoline moieties and reduction of carbonyl carboxamide group resulted in compounds 9-12 and 13, respectively. Compound 12 demonstrated the most promising enzyme inhibition. In addition, compounds 8 and 12 were found to reduce the carrageenan-induced paw edema in vivo by 52.6 and 49.8%, respectively. In view of the encouraging outcomes concerning their notable in vitro and in vivo anti-inflammatory activities, compounds 8 and 12 could be further optimized for the discovery of novel 5-lipoxygenase inhibitors in future. A structure-based 3D pharmacophore model was used in the virtual screening of in-house library to discover novel potential 5-lipoxygenase inhibitors.

Molecular Docking Study of Phytosterols in Lygodium microphyllum Towards SIRT1 and AMPK, the in vitro Brine Shrimp Toxicity Test, and the Phenols and Sterols Levels in the Extract

Background

Lygodium microphyllum is a fern plant with various pharmacological activities, and phytosterols were reported contained in the n-hexane and ethyl acetate extract of this plant. Phytosterols are known to inhibit steatosis, oxidative stress, and inflammation. Sirtuin 1 (SIRT1) and adenosine monophosphate-activated protein kinase (AMPK) are the key proteins that control lipogenesis. However, information about L. microphyllum on SIRT1 and AMPK is still lacking.

Purpose

This study aims to investigate the binding mode of phytosterols in L. microphyllum extract towards AMPK and SIRT1, and the toxicity of the extract against brine shrimp (Artemia salina) larvae, and to determine the phenols and sterols levels in the extract.

Methods

The molecular docking was performed towards SIRT1 and AMPK using AutoDock v4.2.6, the toxicity of the extract was assayed against brine shrimp (Artemia salina) larvae, and the phytosterols were analyzed by employing a thin layer chromatography densitometry, and the total phenols were by spectrophotometry.

Results

The molecular docking study revealed that β-sitosterol and stigmasterol could occupy the active allosteric-binding site of SIRT1 and AMPK by binding to important residues similar to the protein's activators. The cold extraction of the plant yields 15.86% w/w. Phytochemical screening revealed the presence of phenols, steroids, flavonoids, alkaloids, and saponins. The total phenols are equivalent to 126 mg gallic acid (GAE)/g dry extract, the total sterols are 954.04 µg/g, and the β-sitosterol level is 283.55 µg/g. The LC50 value of the extract towards A. salina larvae is 203.704 ppm.

Conclusion

Lygodium microphyllum extract may have the potential to be further explored for its pharmacology activities, particularly in the discovery of plant-based anti-dyslipidemic drug candidates. However, further studies are needed to confirm their roles in alleviating lipid disorders.

Computational 3D Modeling-Based Identification of Inhibitors Targeting Cysteine Covalent Bond Catalysts for JAK3 and CYP3A4 Enzymes in the Treatment of Rheumatoid Arthritis

This work aimed to find new inhibitors of the CYP3A4 and JAK3 enzymes, which are significant players in autoimmune diseases such as rheumatoid arthritis. Advanced computer-aided drug design techniques, such as pharmacophore and 3D-QSAR modeling, were used. Two strong 3D-QSAR models were created, and their predictive power was validated by the strong correlation (R2 values > 80%) between the predicted and experimental activity. With an ROC value of 0.9, a pharmacophore model grounded in the DHRRR hypothesis likewise demonstrated strong predictive ability. Eight possible inhibitors were found, and six new inhibitors were designed in silico using these computational models. The pharmacokinetic and safety characteristics of these candidates were thoroughly assessed. The possible interactions between the inhibitors and the target enzymes were made clear via molecular docking. Furthermore, MM/GBSA computations and molecular dynamics simulations offered insightful information about the stability of the binding between inhibitors and CYP3A4 or JAK3. Through the integration of various computational approaches, this study successfully identified potential inhibitor candidates for additional investigation and efficiently screened compounds. The findings contribute to our knowledge of enzyme-inhibitor interactions and may help us create more effective treatments for autoimmune conditions like rheumatoid arthritis.

In Silico Analysis of SARS-CoV-2 Non-Structural Proteins Reveals an Interaction with the Host's Heat Shock Proteins That May Contribute to Viral Replications and Development

The non-structural protein 2 (NSP2) is an RNA-binding protein involved in coronavirus genome replication, and it often decreases human immune response to promote viral invasion and development. It is believed that the NSP2 associates itself with polyamines and heat shock proteins inside the host cell to proceed with viral development. This study aimed to investigate how the SARS-CoV-2 virus' key non-structural proteins (NSP2) utilize polyamines and heat shock proteins using a molecular docking approach and molecular dynamics (MD). ClusPro and HADDOCK servers were used for the docking and Discovery Studio, chimera, and PyMOL were used for analysis. Docking of the heat shock proteins 40 (HSP40), 70 (HSP70), and 90 (HSP90) with SARS-CoV-2 NSP2 resulted in 32, 28, and 19 interactions, respectively. Molecular dynamics revealed Arg458, Asn508, Met297, Arg301, and Trp417 as active residues, and pharmacophore modeling indicated ZINC395648, ZINC01150525, and ZINC85324008 from the zinc database as possible inhibitors for this NSP2.

Quantitative Structure-Activity Relationship in the Series of 5-Ethyluridine, N2-Guanine, and 6-Oxopurine Derivatives with Pronounced Anti-Herpetic Activity

A quantitative analysis of the relationship between the structure and inhibitory activity against the herpes simplex virus thymidine kinase (HSV-TK) was performed for the series of 5-ethyluridine, N2-guanine, and 6-oxopurines derivatives with pronounced anti-herpetic activity (IC50 = 0.09 ÷ 160,000 μmol/L) using the GUSAR 2019 software. On the basis of the MNA and QNA descriptors and whole-molecule descriptors using the self-consistent regression, 12 statistically significant consensus models for predicting numerical pIC50 values were constructed. These models demonstrated high predictive accuracy for the training and test sets. Molecular fragments of HSV-1 and HSV-2 TK inhibitors that enhance or diminish the anti-herpetic activity are considered. Virtual screening of the ChEMBL database using the developed QSAR models revealed 42 new effective HSV-1 and HSV-2 TK inhibitors. These compounds are promising for further research. The obtained data open up new opportunities for developing novel effective inhibitors of TK.

Synthesis, molecular docking, and in vivo antidiabetic evaluation of new benzylidene-2,4-thiazolidinediones as partial PPAR-γ agonists

Peroxisome proliferator-activated receptor-γ (PPAR-γ) partial agonists or antagonists, also termed as selective PPAR-γ modulators, are more beneficial than full agonists because they can avoid the adverse effects associated with PPAR-γ full agonists, such as weight gain and congestive heart disorders, while retaining the antidiabetic efficiency. In this study, we designed and synthesized new benzylidene-thiazolidine-2,4-diones while keeping the acidic thiazolidinedione (TZD) ring at the center, which is in contrast with the typical pharmacophore of PPAR-γ agonists. Five compounds (5a-e) were designed and synthesized in moderate to good yields and were characterized using spectral techniques. The in vivo antidiabetic efficacy of the synthesized compounds was assessed on streptozotocin-induced diabetic mice using standard protocols, and their effect on weight gain was also studied. Molecular docking and molecular dynamics (MD) simulation studies were performed to investigate the binding interactions of the title compounds with the PPAR-γ receptor and to establish their binding mechanism. Antidiabetic activity results revealed that compounds 5d and 5e possess promising antidiabetic activity comparable with the standard drug rosiglitazone. No compound showed considerable effect on the body weight of animals after 21 days of administration, and the findings showed statistical difference (p < 0.05 to p < 0.0001) among the diabetic control and standard drug rosiglitazone groups. In molecular docking study, compounds 5c and 5d exhibited higher binding energies (- 10.1 and - 10.0 kcal/mol, respectively) than the native ligand, non-thiazolidinedione PPAR-γ partial agonist (nTZDpa) (- 9.8 kcal/mol). MD simulation further authenticated the stability of compound 5c-PPAR-γ complex over the 150 ns duration. The RMSD, RMSF, rGyr, SASA, and binding interactions of compound 5c-PPAR-γ complex were comparable to those of native ligand nTZDpa-PPAR-γ complex, suggesting that the title compounds have the potential to be developed as partial PPAR-γ agonists.

Repurposing Drugs for Inhibition against ALDH2 via a 2D/3D Ligand-Based Similarity Search and Molecular Simulation

Aldehyde dehydrogenase-2 (ALDH2) is a crucial enzyme participating in intracellular aldehyde metabolism and is acknowledged as a potential therapeutic target for the treatment of alcohol use disorder and other addictive behaviors. Using previously reported ALDH2 inhibitors of Daidzin, CVT-10216, and CHEMBL114083 as reference molecules, here we perform a ligand-based virtual screening of world-approved drugs via 2D/3D similarity search methods, followed by the assessments of molecular docking, toxicity prediction, molecular simulation, and the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis. The 2D molecular fingerprinting of ECFP4 and FCFP4 and 3D molecule-shape-based USRCAT methods show good performances in selecting compounds with a strong binding behavior with ALDH2. Three compounds of Zeaxanthin (q = 0), Troglitazone (q = 0), and Sequinavir (q = +1 e) are singled out as potential inhibitors; Zeaxanthin can only be hit via USRCAT. These drugs displayed a stronger binding strength compared to the reported potent inhibitor CVT-10216. Sarizotan (q = +1 e) and Netarsudil (q = 0/+1 e) displayed a strong binding strength with ALDH2 as well, whereas they displayed a shallow penetration into the substrate-binding tunnel of ALDH2 and could not fully occupy it. This likely left a space for substrate binding, and thus they were not ideal inhibitors. The MM-PBSA results indicate that the selected negatively charged compounds from the similarity search and Vina scoring are thermodynamically unfavorable, mainly due to electrostatic repulsion with the receptor (q = -6 e for ALDH2). The electrostatic attraction with positively charged compounds, however, yielded very strong binding results with ALDH2. These findings reveal a deficiency in the modeling of electrostatic interactions (in particular, between charged moieties) in the virtual screening via the 2D/3D similarity search and molecular docking with the Vina scoring system.

Targeting allosteric binding site in methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) to identify natural product inhibitors via structure-based computational approach

Cancer has been viewed as one of the deadliest diseases worldwide. Among various types of cancer, breast cancer is the most common type of cancer in women. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is a promising druggable target and is overexpressed in cancerous cells, like, breast cancer. We conducted structure-based modeling on the allosteric site of the enzyme. Targeting the allosteric site avoids the problem of drug resistance. Pharmacophore modeling, molecular docking, HYDE assessment, drug-likeness, ADMET predictions, simulations, and free-energy calculations were performed. The RMSD, RMSF, RoG, SASA, and Hydrogen-bonding studies showed that seven candidates displayed stable behaviour. As per the literature, average superimposed simulated structures revealed a similar protein conformational change in the αE'-βf' loop, causing its displacement away from the allosteric site. The MM-PBSA showed tight binding of six compounds with the allosteric pocket. The effect of inhibitors interacting in the allosteric site causes a decrease in the binding energy of J49 (active-site inhibitor), suggesting the effect of allosteric binding. The PCA and FEL analysis revealed the significance of the docked compounds in the stable behaviour of the complexes. The outcome can contribute to the development of potential natural products with drug-like properties that can inhibit the MTHFD2 enzyme.

Identification of novel dual acting ligands targeting the adenosine A2A and serotonin 5-HT1A receptors

GPCRs are a family of transmembrane receptors that are profoundly linked to various neurological disorders, among which is Parkinson's disease (PD). PD is the second most ubiquitous neurological disorder after Alzheimer's disease, characterized by the depletion of dopamine in the central nervous system due to the impairment of dopaminergic neurons, leading to involuntary movements or dyskinesia. The current standard of care for PD is Levodopa, a dopamine precursor, yet the chronic use of this agent can exacerbate motor symptoms. Recent studies have investigated the effects of combining A2AR antagonist and 5-HT1A agonist on dyskinesia and motor complications in animal models of PD. It has been proved that the drug combination has significantly improved involuntary movements while maintaining motor activity, highlighting as a result new lines of therapy for PD treatments, through the regulation of both receptors. Using a combination of ligand-based pharmacophore modelling, virtual screening, and molecular dynamics simulation, this study intends on identifying potential dual-target compounds from IBScreen. Results showed that the selected models displayed good enrichment metrics with a near perfect receiver operator characteristic (ROC) and Area under the accumulation curve (AUAC) values, signifying that the models are both specific and sensitive. Molecular docking and ADMET analysis revealed that STOCK2N-00171 could be potentially active against A2AR and 5-HT1A. Post-MD analysis confirmed that the ligand exhibits a stable behavior throughout the simulation while maintaining crucial interactions. These results imply that STOCK2N-00171 can serve as a blueprint for the design of novel and effective dual-acting ligands targeting A2AR and 5-HT1A.Communicated by Ramaswamy H. Sarma.

Discovery and pharmacophoric characterization of chemokine network inhibitors using phage-display, saturation mutagenesis and computational modelling

CC and CXC-chemokines are the primary drivers of chemotaxis in inflammation, but chemokine network redundancy thwarts pharmacological intervention. Tick evasins promiscuously bind CC and CXC-chemokines, overcoming redundancy. Here we show that short peptides that promiscuously bind both chemokine classes can be identified from evasins by phage-display screening performed with multiple chemokines in parallel. We identify two conserved motifs within these peptides and show using saturation-mutagenesis phage-display and chemotaxis studies of an exemplar peptide that an anionic patch in the first motif and hydrophobic, aromatic and cysteine residues in the second are functionally necessary. AlphaFold2-Multimer modelling suggests that the peptide occludes distinct receptor-binding regions in CC and in CXC-chemokines, with the first and second motifs contributing ionic and hydrophobic interactions respectively. Our results indicate that peptides with broad-spectrum anti-chemokine activity and therapeutic potential may be identified from evasins, and the pharmacophore characterised by phage display, saturation mutagenesis and computational modelling.

An Updated Review on Developing Small Molecule Kinase Inhibitors Using Computer-Aided Drug Design Approaches

Small molecule kinase inhibitors (SMKIs) are of heightened interest in the field of drug research and development. There are 79 (as of July 2023) small molecule kinase inhibitors that have been approved by the FDA and hundreds of kinase inhibitor candidates in clinical trials that have shed light on the treatment of some major diseases. As an important strategy in drug design, computer-aided drug design (CADD) plays an indispensable role in the discovery of SMKIs. CADD methods such as docking, molecular dynamic, quantum mechanics/molecular mechanics, pharmacophore, virtual screening, and quantitative structure-activity relationship have been applied to the design and optimization of small molecule kinase inhibitors. In this review, we provide an overview of recent advances in CADD and SMKIs and the application of CADD in the discovery of SMKIs.

Identification of Promising Sulfonamide Chalcones as Inhibitors of SARS-CoV-2 3CLpro through Structure-Based Virtual Screening and Experimental Approaches

3CLpro is a viable target for developing antiviral therapies against the coronavirus. With the urgent need to find new possible inhibitors, a structure-based virtual screening approach was developed. This study recognized 75 pharmacologically bioactive compounds from our in-house library of 1052 natural product-based compounds that satisfied drug-likeness criteria and exhibited good bioavailability and membrane permeability. Among these compounds, three promising sulfonamide chalcones were identified by combined theoretical and experimental approaches, with SWC423 being the most suitable representative compound due to its competitive inhibition and low cytotoxicity in Vero E6 cells (EC50 = 0.89 ± 0.32 μM; CC50 = 25.54 ± 1.38 μM; SI = 28.70). The binding and stability of SWC423 in the 3CLpro active site were investigated through all-atom molecular dynamics simulation and fragment molecular orbital calculation, indicating its potential as a 3CLpro inhibitor for further SARS-CoV-2 therapeutic research. These findings suggested that inhibiting 3CLpro with a sulfonamide chalcone such as SWC423 may pave the effective way for developing COVID-19 treatments.

Using macromolecular electron densities to improve the enrichment of active compounds in virtual screening

The quest for effective virtual screening algorithms is hindered by the scarcity of training data, calling for innovative approaches. This study presents the use of experimental electron density (ED) data for improving active compound enrichment in virtual screening, supported by ED's ability to reflect the time-averaged behavior of ligands and solvents in the binding pocket. Experimental ED-based grid matching score (ExptGMS) was developed to score compounds by measuring the degree of matching between their binding conformations and a series of multi-resolution experimental ED grids. The efficiency of ExptGMS was validated using both in silico tests with the Directory of Useful Decoys-Enhanced dataset and wet-lab tests on Covid-19 3CLpro-inhibitors. ExptGMS improved the active compound enrichment in top-ranked molecules by approximately 20%. Furthermore, ExptGMS identified four active inhibitors of 3CLpro, with the most effective showing an IC50 value of 1.9 µM. We also developed an online database containing experimental ED grids for over 17,000 proteins to facilitate the use of ExptGMS for academic users.

Computational exploration of natural compounds targeting Staphylococcus aureus: inhibiting AgrA promoter binding for antimicrobial intervention

Staphylococcus aureus is a highly virulent nosocomial pathogen that poses a significant threat to individuals exposed to healthcare settings. Due to its sophisticated machinery for producing virulence factors, S. aureus can cause severe and potentially fatal infections in humans. This study focuses on the response regulator AgrA, which plays a crucial role in regulating the production of virulence factors in S. aureus. The objective is to identify natural compounds that can inhibit the binding of AgrA to its promoter site, thus inhibiting the expression of virulence genes. To achieve this, a pharmacophore model was generated using known drugs and applied to screen the ZINC natural product database. The resulting compounds were subjected to molecular docking-based virtual screening against the C-terminal DNA binding domain of AgrA. Three compounds, namely ZINC000077269178, ZINC000051012304, and ZINC000004266026, were shortlisted based on their strong affinity for key residues involved in DNA binding and transcription initiation. Subsequently, the unbound and ligand-bound complexes were subjected to a 200 ns molecular dynamics simulation to assess their conformational stability. Various analyses, including RMSD, RMSF, Rg, SASA, Principal Component Analysis, and Gibbs free energy landscape, were conducted on the simulation trajectory. The RMSD profile indicated similar fluctuations in both bound and unbound structures, while the Rg profile demonstrated the compactness of the protein without any unfolding during the simulation. Furthermore, Principal component analysis revealed that ligand binding reduced the overall atomic motion of the protein whereas free energy landscape suggested the energy variations obtained in complexes.Communicated by Ramaswamy H. Sarma.

Computer-Aided Drug Design of Novel Derivatives of 2-Amino-7,9-dihydro-8H-purin-8-one as Potent Pan-Janus JAK3 Inhibitors

Rheumatoid arthritis (RA) remains one of the most prevalent autoimmune diseases worldwide. Janus kinase 3 (JAK3) is an essential enzyme for treating autoimmune diseases, including RA. Molecular modeling techniques play a crucial role in the search for new drugs by reducing time delays. In this study, the 3D-QSAR approach is employed to predict new JAK3 inhibitors. Two robust models, both field-based with R2 = 0.93, R = 0.96, and Q2 = 87, and atom-based with R2 = 0.94, R = 0.97, and Q2 = 86, yielded good results by identifying groups that may readily direct their interaction. A reliable pharmacophore model, DHRRR1, was provided in this work to enable the clear characterization of chemical features, leading to the design of 13 inhibitors with their pIC50 values. The DHRRR1 model yielded a validation result with a ROC value of 0.87. Five promising inhibitors were selected for further study based on an ADMET analysis of their pharmacokinetic properties and covalent docking (CovDock). Compared to the FDA-approved drug tofacitinib, the pharmaceutical features, binding affinity and stability of the inhibitors were analyzed through CovDock, 300 ns molecular dynamics simulations, free energy binding calculations and ADMET predictions. The results show that the inhibitors have strong binding affinity, stability and favorable pharmaceutical properties. The newly predicted molecules, as JAK3 inhibitors for the treatment of RA, are promising candidates for use as drugs.

Identification of inhibitors for Agr quorum sensing system of Staphylococcus aureus by machine learning, pharmacophore modeling, and molecular dynamics approaches

CONTEXT

Staphylococcus aureus is a highly pathogenic organism that is the most common cause of postoperative complications as well as severe infections like bacteremia and infective endocarditis. By mediating the formation of biofilms and the expression of virulent genes, the quorum sensing (QS) mechanism is a major contributor to the development of these diseases. By hindering its QS network, an innovative approach to avoiding this bacterial infection is taken. Targeting the AgrA of the Agr system serves as beneficial in holding the top position in the QS system cascade.

METHODS

Using known AgrA inhibitors, the machine learning algorithms (artificial neural network, naïve Bayes, random forest, and support vector machine) and pharmacophore model were developed. The potential lead compounds were screened against the Zinc and COCONUT databases using the best pharmacophore hypothesis. The hits were then subjected second screening process using the best machine learning model. The predicted active compounds were then reranked based on the docking score. The stability of AgrA-lead compounds was studied using molecular dynamics approaches, and an ADME profile was also carried out. Five lead compounds, namely, CNP02386963,4,5-trihydroxy-2-[({7,13,14-trihydroxy-3,10-dioxo-2,9-dioxatetracyclo[6.6.2.04,16.011,15]hexadeca-1(14),4,6,8(16),11(15),12-hexaen-6-yl}oxy)methyl]benzoic acid, CNP0129274 4-(dimethylamino)-1,5,6,10,12,12a-hexahydroxy-6-methyl-3,11-dioxo-3,4,4a,5,5a,6,11,12a-octahydrotetracene-2-carboxamide, CNP0242717 3-Hydroxyasebotin, CNP0361624 3,4,5-trihydroxy-6-[(2,4,5,6,7-pentahydroxy-1-oxooctan-3-yl)oxy]oxane-2-carboxylic acid, and CNP0285058 2-{[4,5-dihydroxy-6-(hydroxymethyl)-3-[(3,4,5-trihydroxy-6-methyloxan-2-yl)oxy]oxan-2-yl]oxy}-2-(4-hydroxyphenyl)acetonitrile were obtained using the two-step virtual screening process. The molecular dynamics study revealed that the CNP0238696 was found to be stable in the binding pocket of AgrA. ADME profiles show that this compound has two Lipinski violations and low bioavailability. Further studies should be performed to assess the anti-biofilm activity of the lead compound in vitro.

Discovery and verification of bitter components in Panax notoginseng based on the integrated strategy of pharmacophore model, system separation and bitter tracing technology

Panax notoginseng is a world-renowned tonic herb, which has been used as a characteristic food in Southwest China for hundreds of years. However, the taste of Panax notoginseng is extremely bitter and serious after tasting, and its bitter components are unknown. This manuscript proposes a new strategy for discovering bitter components of Panax notoginseng based on the integrated analysis of pharmacophore model, system separation and bitter tracing technology. Firstly, 16 potential bitter components were obtained by UPLC-Q-Orbitrap HRMS combined with virtual screening, most of which were saponins.Then, the bitter components were further separated by system component separation and 5 potential bitter components were obtained. Finally, the main contributors of bitterness in Panax notoginseng were verified to be Ginsenoside Rg1, Ginsenoside Rb1 and Ginsenoside Rd by components knock-in and fNIRS. In general, this paper is the first literature report on the relatively systematic study of bitter components in Panax notoginseng.

Identification of Anti-Influenza A Compounds Inhibiting the Viral Non-Structural Protein 1 (NS1) Using a Type I Interferon-Driven Screening Strategy

There is an urgent need to identify efficient antiviral compounds to combat existing and emerging RNA virus infections, particularly those related to seasonal and pandemic influenza outbreaks. While inhibitors of the influenza viral integral membrane proton channel protein (M2), neuraminidase (NA), and cap-dependent endonuclease are available, circulating influenza viruses acquire resistance over time. Thus, the need for the development of additional anti-influenza drugs with novel mechanisms of action exists. In the present study, a cell-based screening assay and a small molecule library were used to screen for activities that antagonized influenza A non-structural protein 1 (NS1), a highly conserved, multifunctional accessory protein that inhibits the type I interferon response against influenza. Two potential anti-influenza agents, compounds 157 and 164, were identified with anti-NS1 activity, resulting in the reduction of A/PR/8/34(H1N1) influenza A virus replication and the restoration of IFN-β expression in human lung epithelial A549 cells. A 3D pharmacophore modeling study of the active compounds provided a glimpse of the structural motifs that may contribute to anti-influenza virus activity. This screening approach is amenable to a broader analysis of small molecule compounds to inhibit other viral targets.

Structure-Based Drug Design of RdRp Inhibitors against SARS-CoV-2

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a worldwide pandemic since 2019, spreading rapidly and posing a significant threat to human health and life. With over 6 billion confirmed cases of the virus, the need for effective therapeutic drugs has become more urgent than ever before. RNA-dependent RNA polymerase (RdRp) is crucial in viral replication and transcription, catalysing viral RNA synthesis and serving as a promising therapeutic target for developing antiviral drugs. In this article, we explore the inhibition of RdRp as a potential treatment for viral diseases, analysing the structural information of RdRp in virus proliferation and summarizing the reported inhibitors' pharmacophore features and structure-activity relationship profiles. We hope that the information provided by this review will aid in structure-based drug design and aid in the global fight against SARS-CoV-2 infection.

Recent Advances in Computer-Aided Structure-Based Drug Design on Ion Channels

Ion channels play important roles in fundamental biological processes, such as electric signaling in cells, muscle contraction, hormone secretion, and regulation of the immune response. Targeting ion channels with drugs represents a treatment option for neurological and cardiovascular diseases, muscular degradation disorders, and pathologies related to disturbed pain sensation. While there are more than 300 different ion channels in the human organism, drugs have been developed only for some of them and currently available drugs lack selectivity. Computational approaches are an indispensable tool for drug discovery and can speed up, especially, the early development stages of lead identification and optimization. The number of molecular structures of ion channels has considerably increased over the last ten years, providing new opportunities for structure-based drug development. This review summarizes important knowledge about ion channel classification, structure, mechanisms, and pathology with the main focus on recent developments in the field of computer-aided, structure-based drug design on ion channels. We highlight studies that link structural data with modeling and chemoinformatic approaches for the identification and characterization of new molecules targeting ion channels. These approaches hold great potential to advance research on ion channel drugs in the future.

Discovery of selective HDAC6 inhibitors based on a multi-layer virtual screening strategy

The abnormal enhancement of histone deacetylase 6 (HDAC6) has been demonstrated to be closely related to the occurrence and development of various malignant tumors, attracting extensive attention as a promising target for cancer therapy. Currently, only limited selective HDAC6 inhibitors have entered clinical trials, making the rapid discovery of selective HDAC6 inhibitors with safety profiles particularly urgent. In this study, a multi-layer virtual screening workflow was established, and the representative compounds screened were biologically evaluated in combination with enzyme inhibitory and anti-tumor cell proliferation experiments. The experimental results showed that the screened compounds L-25, L-32, L-45 and L-81 exhibited nanomolar inhibitory activity against HDAC6, and exerted a certain degree of anti-proliferative activities against tumor cells, especially the cytotoxicity of L-45 to A375 (IC₅₀ = 11.23 ± 1.27 μM) and the cytotoxicity of L-81 against HCT-116 (IC₅₀ = 12.25 ± 1.13 μM). Additionally, the molecular mechanisms underlying the subtype selective inhibitory activities of the selected compounds were further elucidated using computational approaches, and the hotspot residues on HDAC6 contributing to the ligands' binding were identified. In summary, this study established a multi-layer screening scheme to quickly and effectively screen out hit compounds with enzyme inhibitory activity and anti-tumor cell proliferation, providing novel scaffolds for the subsequent anti-tumor drug design based on HDAC6 target.

Developing new PI3Kγ inhibitors by combining pharmacophore modeling, molecular dynamic simulation, molecular docking, fragment-based drug design, and virtual screening

Since dysregulation of the phosphatidylinositol 3-kinase gamma (PI3Kγ) signaling pathway is associated with the pathogenesis of cancer, inflammation, and autoimmunity, PI3Kγ has emerged as an attractive target for drug development. IPI-549 is the only selective PI3Kγ inhibitor that has advanced to clinical trials, thus, IPI-549 could serve as a promising template for designing novel PI3Kγ inhibitors. In this present study, a modeling strategy consisting of common feature pharmacophore modeling, receptor-ligand pharmacophore modeling, and molecular dynamics simulation was utilized to identify the key pharmacodynamic characteristic elements of the target compound and the key residue information of the PI3Kγ interaction with the inhibitors. Then, 10 molecules were designed based on the structure-activity relationships, and some of them exhibited satisfactory predicted binding affinities to PI3Kγ. Finally, a hierarchical multistage virtual screening method, involving the developed common feature and receptor-ligand pharmacophore model and molecular docking, was constructed for screening the potential PI3Kγ inhibitors. Overall, we hope these findings would provide some guidance for the development of novel PI3Kγ inhibitors.