PubMed HealthSearch

SEARCH · PubMed Health

Results for “Molecular docking and molecular dynamics simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Validating the potential mechanism and therapeutic effect of Qinlian Jiangxia decoction in the treatment of type 2 diabetes mellitus complicated with hyperlipidemia through network pharmacology, molecular docking, molecular dynamics simulation, andexperiments.

OBJECTIVE: To investigate the mechanism of action of Qinlian Jiangxia decoction (, QLJXD) in the treatment of type 2 diabetes mellitus (T2DM) complicated by hyperlipidemia using network pharmacology, molecular docking, molecular dynamics simulation and in vivo experiments. METHODS: Drug components, targets and disease targets were identified using databases such as TCM systems pharmacology database and analysis platform and GeneCards. The intersecting targets were subjected to protein-protein interaction analysis using the search tool for the retrieval of interacting genes/proteins database. Subsequently, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis of the intersecting targets were conducted using the Metascape platform to identify core components and targets. The results were validated using molecular docking, molecular dynamics simulations and in vivo experiments. RESULTS: QLJXD contains 76 active ingredients and 136 disease targets. The core ingredients are quercetin, β-sitosterol, wogonin and baicalein, while the core targets are fatty acid binding protein 4 (FABP4) and peroxisome proliferative activated receptor gamma (PPARG). Molecular docking and molecular dynamics simulations revealed that the core ingredients bound well to the core targets. Animal experiments demonstrated that QLJXD effectively inhibited the expression of FABP4 and increased the expression of PPARG, thereby enhancing disorders of glycolipid metabolism. CONCLUSION: The putative therapeutic efficacy of QLJXD in the management of T2DM complicated with hyperlipidemia may be ascribed to the synergistic actions of multiple components, such as quercetin, β-sitosterol, wogonin, and baicalein, which collectively modulate FABP4 and PPARG molecular targets.

Molecular Docking Simulation

Molecular docking, molecular dynamics simulation, and enzyme inhibitory studies of vitamin K family members on aldose reductase.

Aldose reductase (AR) is a key enzyme in the polyol pathway and plays a major role in the progression of secondary complications of diabetes. Despite extensive efforts to develop natural and synthetic aldose reductase inhibitors (ARIs), most candidates have shown limited clinical efficacy, highlighting the need for more potent and selective inhibitors. In this study, we have systematically evaluated the inhibitory potential of vitamin K family members (vitamin K1, vitamin K2, and vitamin K3) using molecular docking, protein-ligand interaction analysis, molecular dynamics simulations, and enzyme kinetics. Docking analysis predicted that vitamin K2 has the highest binding affinity for AR. Subsequent molecular dynamics simulations revealed that both vitamin K1 and vitamin K2 formed stable complexes with the protein, exhibiting comparable RMSD (∼0.5 Å difference), similar RMSF profiles, and reduced radius of gyration, indicating compact and stable binding. Interaction analysis demonstrated that ligand binding is predominantly driven by hydrophobic interactions, with vitamin K2 forming a higher number of hydrophobic contacts, while vitamin K1 exhibited slightly more hydrogen bonding. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) results further supports stronger binding of vitamin K2 (-56 kcal/mol) compared to vitaminK1 (-51 kcal/mol). Consistent with these findings, enzyme kinetics showed a slightly lower Ki value for vitamin K2 than vitamin K1. In contrast, vitamin K3 failed to maintain stable binding and moved out of the active site during simulation. Overall, the study highlights that hydrophobic interaction-driven stabilization plays a key role in ligand binding, and identifies vitamin K1 and vitamin K2 as promising inhibitors against AR, with vitamin K2 exhibiting more favourable hydrophobic interactions and binding stability.

Aldose Reductase

Elucidating the Mechanism of Xiaoqinglong Decoction in Chronic Urticaria Treatment: An Integrated Approach of Network Pharmacology, Bioinformatics Analysis, Molecular Docking, and Molecular Dynamics Simulations.

INTRODUCTION: Xiaoqinglong Decoction (XQLD) is a traditional Chinese medicinal formula commonly used to treat chronic urticaria (CU). However, its underlying therapeutic mechanisms remain incompletely characterized. This study employed an integrated approach combining network pharmacology, bioinformatics, molecular docking, and molecular dynamics simulations to identify the active components, potential targets, and related signaling pathways involved in XQLD's therapeutic action against CU, thereby providing a mechanistic foundation for its clinical application. METHODS: The active components of XQLD and their corresponding targets were identified using the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. CU-related targets were retrieved from the OMIM and GeneCards databases. Subsequently, core components and targets were determined via protein-protein interaction (PPI) network analysis and component-target-pathway network construction. Topological analyses were performed using Cytoscape software to prioritize core nodes within these networks. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted via the DAVID database to identify enriched biological processes and signaling pathways. Molecular docking was performed to evaluate binding interactions between key components and core targets, while molecular dynamics (MD) simulations were employed to assess the stability of the component-target complexes with the lowest binding energy. Finally, CU-related targets of XQLD were validated using datasets from the Gene Expression Omnibus (GEO) database. RESULTS: A total of 135 active components and 249 potential targets of XQLD were identified, alongside 1,711 CU-related targets. Core components, such as quercetin, kaempferol, beta-sitosterol, naringenin, stigmasterol, and luteolin, exhibited high degree values in the constructed networks. The core targets identified included AKT1, TNF, IL6, TP53, PTGS2, CASP3, BCL2, ESR1, PPARG, and MAPK3. GO and KEGG pathway enrichment analyses revealed the PI3K-Akt signaling pathway as a central regulatory mechanism. Molecular docking studies demonstrated strong binding affinities between active components and core targets, with the stigmasterol-AKT1 complex exhibiting the lowest binding energy (-11.4 kcal/mol) and high stability in MD simulations. Validation using GEO datasets identified 12 core genes shared between CU-related targets and XQLD-associated targets, including PTGS2 and IL6, which were also prioritized as core targets in the network pharmacology analyses. DISCUSSION: This study comprehensively integrates multidisciplinary approaches to clarify the potential molecular mechanisms of XQLD in treating CU, highlighting its multitarget and multipathway synergistic effects. Molecular docking and dynamics simulations confirm the stable interaction between stigmasterol and the core target AKT1. Additionally, GEO dataset analysis verifies the pathogenic relevance of targets such as PTGS2 and IL6, significantly enhancing the credibility of our findings. These results provide a modern scientific basis for the traditional therapeutic effects of XQLD on CU and have important implications for developing multitarget treatments for this condition. However, this study mainly relies on database mining and computational simulations. Further in vitro and in vivo experimental validations are needed to confirm the predicted component-target-pathway interactions. CONCLUSION: This study identifies the active components, potential targets, and pathways through which XQLD exerts therapeutic effects on CU. These findings provide a theoretical foundation for further mechanistic studies and support their clinical application in the treatment of CU.

Molecular Docking Simulation

Zhiling Jiangya decoction treats hypertension in rats: An integrative study of network pharmacology, immune infiltration, molecular simulation, and 16S rDNA sequencing.

OBJECTIVE: This study integrated network pharmacology, immune infiltration analysis, molecular docking, molecular dynamics simulation, ADMET prediction, 16S rDNA sequencing, and rat experiments to elucidate the potential mechanisms underlying the antihypertensive effects of Zhiling Jiangya Decoction (ZLJYD). METHODS: Active compounds and their potential targets were screened from the PubChem, TCMSP, NovoPro, and SwissTargetPrediction databases. Hypertension-related targets were retrieved from the OMIM and GeneCards databases, and overlapping targets were identified. The STRING database and Cytoscape 3.10.1 software were used to construct a protein-protein interaction network and a herb-component-target-disease network. Gene Ontology functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were performed to identify the key biological processes and signaling pathways involved. Using the CIBERSORT algorithm combined with correlation analysis, we investigated the association between key targets and immune cell infiltration. Molecular docking, molecular dynamics simulations, and ADMET predictions were performed to assess the binding stability and pharmacokinetic properties of the main compounds with their corresponding targets. Finally, the antihypertensive efficacy of ZLJYD was validated using a spontaneously hypertensive rat model, and alterations in gut microbiota were analyzed using 16S rDNA sequencing. RESULTS: A total of 123 active compounds and 267 hypertension-related targets of ZLJYD were identified. Enrichment analysis revealed that these targets were primarily associated with the PI3K-Akt signaling pathway and lipid and atherosclerosis pathways. Immune infiltration analysis suggested that the therapeutic effects of ZLJYD may involve the regulation of follicular helper T cells, naïve B cells, and naïve CD4⁺ T cells. Molecular docking and dynamics simulations supported the stable binding of key compounds to their target proteins, while ADMET predictions indicated favorable pharmacokinetic properties and safety profiles. Rat experiments demonstrated that ZLJYD significantly reduced blood pressure in spontaneously hypertensive rats, partially alleviated gut microbiota dysbiosis, and altered microbial community structure and phylogenetic diversity. CONCLUSION: This study systematically elucidates the potential mechanisms underlying the antihypertensive effects of ZLJYD through multiple components, targets, and pathways, particularly immune regulation and gut microbiota remodeling. These findings provide mechanistic insights into its potential therapeutic application.

16S rDNA sequencing

In silico screening of anti-atherosclerotic compounds from Morus alba leaves by machine learning and network pharmacology.

OBJECTIVE: This study integrates machine learning with network pharmacology, molecular docking, and molecular dynamics simulations to screen bioactive compounds from Mulberry leaves and elucidate their potential mechanisms against atherosclerosis (AS). METHODS: A training dataset of anti-AS active compounds was compiled and encoded as Morgan fingerprints. Three machine learning classifiers, specifically Random Forest (RF), Support Vector Machine (SVM), and Extreme Gradient Boosting (XG-Boost), were constructed and evaluated using multiple performance metrics. Potential active components from Mulberry leaves and AS-related targets were retrieved, followed by protein-protein interaction network construction and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Molecular docking was then performed to evaluate binding affinities between core targets and candidate compounds, and the most stable complex was subjected to molecular dynamics simulations using GROMACS (2025). RESULTS: The RF model achieved superior performance (accuracy= 0.8354, F1 = 0.8408, AUC = 0.9119) with 100% external validation accuracy. Thirteen anti-AS candidates were prioritized from mulberry leaves, four of which have been previously documented. Network pharmacology revealed AKT1 and IL6 as core targets, enriched in pathways such as endocrine resistance. Molecular docking and dynamics simulations confirmed strong binding between oxysanguinarine and AKT1, with the complex exhibiting high stability. CONCLUSION: The RF model provides a reliable computational tool for prioritizing anti-AS compounds from Mulberry leaves. The integrated analysis reveals that Mulberry leaves exert anti-atherosclerotic effects through multi-target (e.g., AKT1, IL6) and multi-pathway (e.g., PI3K-Akt) mechanisms, offering a framework for further experimental validation.

Morus

Structure-based virtual screening, multi-score docking, and molecular dynamics simulation of novel small molecules targeting the epidermal growth factor receptor for potential management of oral squamous cell carcinoma.

UNLABELLED: Oral squamous cell carcinoma (OSCC) is a major global health burden, with epidermal growth factor receptor (EGFR) serving as an important therapeutic target. However, resistance to currently available EGFR inhibitors limits the efficacy of long-term treatment. In this study, a structure-based virtual screening approach was employed using the Mcule database to identify novel small molecules with potential EGFR-inhibitory activity. The top-ranked compounds were subjected to consensus docking using multiple docking platforms and compared with established EGFR inhibitors. The most promising complexes were further evaluated using 500 ns molecular dynamics simulations to investigate their structural stability, conformational flexibility, and binding persistence. ADMET and pharmacokinetic analyses were performed to assess the drug-like and safety profiles. Five lead compounds (C1-C5) demonstrated significant binding affinities toward EGFR, ranging from - 9.9 to - 9.2 kcal/mol, while satisfying the major drug-likeness criteria. Molecular dynamics simulations suggested that C1 and C4 may form relatively stable EGFR-ligand complexes, as supported by stable RMSD convergence and persistent interactions with key active-site residues throughout the simulation period. Trajectory-based interaction analyses further indicated a sustained binding behavior. ADMET profiling predicted favorable oral bioavailability and low predicted toxicity for most compounds, particularly C3 and C5, although a potential risk of cytochrome P450-mediated drug-drug interactions was observed. Overall, the shortlisted compounds exhibited docking and dynamic stability profiles comparable to those of the reference inhibitor Lapatinib. These findings suggest the potential therapeutic relevance of structurally novel EGFR-targeting scaffolds in OSCC and provide a foundation for future experimental validation through in vitro and in vivo studies. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at https://doi.org/10.1007/s40203-026-00722-4.

ADMET

Matrine Alleviates Sepsis-Induced Acute Lung Injury by Reinforcing NQO1/SLC7A11/GPX4-Associated Anti-Ferroptotic Defenses and Attenuating NF-κB-Driven Inflammation.

BACKGROUND: Sepsis triggers dysregulated systemic inflammation and multiple-organ dysfunction, with the lungs being particularly susceptible to injury. Sepsis-induced acute respiratory distress syndrome (ARDS) is mainly driven by TLR4/NF-κB-mediated hyperinflammation and alveolar macrophage activation. Matrine, a bioactive alkaloid derived from Sophora flavescens, has been reported to modulate redox homeostasis and ferroptosis-associated lipid peroxidation. However, the target-specific mechanisms underlying its effects on ferroptosis and inflammatory signaling in sepsis-induced acute lung injury (SALI) remain incompletely understood. PURPOSE: This study aimed to evaluate the therapeutic effects of matrine in a cecal ligation and puncture (CLP)-induced SALI model and to determine whether its protective effects involve reinforcement of NQO1/SLC7A11/GPX4-associated anti-ferroptotic defenses and suppression of NF-κB-driven inflammation. METHODS: We analyzed the single-cell RNA-sequencing (scRNA-seq) dataset GSE273924 to characterize CD45-enriched pulmonary immune-cell subsets in sham mice and mice with intratracheal Escherichia coli-induced pneumonia. Network pharmacology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to predict Kushen (KS)-related targets and pathways associated with SALI. Differential expression analysis and weighted gene co-expression network analysis (WGCNA) of GSE245013 were used to identify candidate targets. Matrine-NQO1 binding and intracellular target engagement were evaluated using molecular docking, molecular dynamics simulations, surface plasmon resonance (SPR), and the cellular thermal shift assay (CETSA). The therapeutic effects of matrine were assessed in mice with CLP-induced SALI and in lipopolysaccharide (LPS)-stimulated MH-S cells. Lung histopathology, inflammatory cytokine production, target protein expression, ferroptosis-associated indicators, and NF-κB activation were evaluated using molecular, biochemical, and histological assays. The functional contribution of NQO1 was further examined using the NQO1 inhibitor ES936. RESULTS: scRNA-seq analysis of GSE273924 revealed substantial remodeling of the CD45-enriched pulmonary immune-cell landscape in mice with intratracheal E. coli-induced pneumonia, including macrophage transcriptional programs associated with ferroptosis and inflammatory signaling. Integrated network pharmacology and bioinformatics analyses prioritized NQO1 as a candidate target of matrine and identified NF-κB signaling as a potentially relevant pathway. Molecular docking, molecular dynamics simulations, SPR, and CETSA supported matrine-NQO1 binding and intracellular target engagement. Functionally, matrine improved survival, attenuated lung injury, reinforced NQO1/SLC7A11/GPX4-associated anti-ferroptotic defenses, and suppressed NF-κB activation in CLP mice. Similar protective effects were observed in LPS-stimulated MH-S cells. ES936 partially attenuated the matrine-mediated improvements in cell viability, redox homeostasis, ferroptosis-associated indicators, and NF-κB p65 phosphorylation, supporting a functional contribution of NQO1 to the protective effects of matrine. CONCLUSION: Matrine alleviates SALI by reinforcing NQO1/SLC7A11/GPX4-associated anti-ferroptotic defenses and attenuating NF-κB-driven inflammation.

Animals

Integrated immunoinformatics for the design of novel multi-epitope vaccine and identification of new drug targets against Stenotrophomonas maltophilia, a multidrug-resistant superbug.

BACKGROUND: Stenotrophomonas maltophilia is a multidrug-resistant opportunistic pathogen causing severe hospital-acquired infections, especially in immunocompromised patients. The absence of an effective vaccine and rising antibiotic resistance underscore the need for novel interventions. This study employed an integrated reverse vaccinology and computational analyses to identify new immunogenic targets, design a multi-epitope vaccine (MEV), and propose potential drug targets. METHODS: A comprehensive immunoinformatics pipeline was employed to assess antigenicity, allergenicity, human similarity, and physicochemical properties of S. maltophilia proteins. Both B- and T-cell epitopes were screened; however, only the top B-cell epitopes were selected for MEV construction, given the extracellular nature of S. maltophilia. MEV-TLR interactions were analyzed through molecular docking and dynamics simulations. In parallel, cytoplasmic proteins were screened via a subtractive genomics approach to identify essential, non-human homologous, and non-microbiome-similar proteins, which were further evaluated for druggability and interaction networks to propose novel therapeutic targets. RESULTS: From a total of 4111 proteins, seven potential immunogenic targets were identified: GspD (WP_108270537.1), FhuE (WP_049451370.1), fimbrial protein (WP_012479122.1), TonB-dependent receptor (WP_169448402.1), TolC family protein (WP_108270106.1), autotransporter beta-barrel OMP (WP_169448945.1), and a hypothetical protein (WP_005407892.1). Subsequently, an MEV was designed using five immunogenic epitopes derived from four of these targets: WP_005407892.1 (ADQDSSNM), WP_049451370.1 (SGKAEQ and GEESKTPS), WP_108270537.1 (GVTSTQSDSERT), and WP_169448945.1 (RELGGDRNE). Molecular docking and molecular dynamics simulations demonstrated strong, stable, and feasible interactions between the MEV and TLR-2 and TLR-4 receptors. Moreover, nine novel drug targets were predicted for S. maltophilia, providing new therapeutic insights. CONCLUSION: The designed MEV and identified immunogenic targets represent promising vaccine candidates against S. maltophilia. Further in vitro and in vivo studies are essential to confirm their safety, immunogenicity, and protective efficacy. Additionally, subtractive genomics analysis revealed nine novel, non-homologous drug targets, offering safer and more specific therapeutic avenues.

Drug targets

LARS promotes hepatocellular carcinoma progression via the PI3K/AKT/mTOR pathway and interaction with RPS5, and serves as a prognostic biomarker.

BACKGROUND: Hepatocellular carcinoma (HCC) caused many cancer deaths around the world. Its progression involves complex mechanisms, creating an urgent need to identify new therapeutic targets. Leucine-tRNA synthetase (LARS) is a key enzyme for protein synthesis, but its specific role and mechanism in HCC are not well understood. PURPOSE: This research aims to investigate the biological function, molecular mechanism, and clinical relevance of the LARS gene in HCC progression, to assess its potential as a treatment target. METHODS: LARS expression was assessed in HCC cell lines (PLC-PRF-5, HCC-LM3) and in mouse subcutaneous tumor models using siRNA and adeno-associated virus (AAV). Techniques including Cell Counting Kit-8(CCK-8), colony formation, EdU, Transwell, wound healing, and flow cytometry were used to measure cell proliferation, migration, invasion, and apoptosis. RNA-seq, proteomics (TMT), western blot, co-immunoprecipitation (Co-IP) with mass spectrometry, molecular docking, and molecular dynamics simulation were employed to study the affected signaling pathway (PI3K/AKT/mTOR) and interacting protein (RPS5). The TCGA (The Cancer Genome Atlas) database and UALCAN platform were used to analyze links between LARS expression and clinicopathological features or prognosis in HCC patients. RESULTS: Reducing LARS expression significantly inhibited the proliferation, colony formation, migration, and invasion of HCC cells, while promoting apoptosis. In mice, LARS knockdown markedly slowed tumor growth. Mechanistic studies showed that reducing LARS expression levels affected the PI3K/AKT/mTOR signaling pathway and led to decreased levels of the key interacting protein RPS5. Overexpressing RPS5 partly reversed the proliferation inhibition caused by LARS depletion. Molecular docking and dynamics simulations suggested that the environmental contaminant triphenyl phosphate (TPP) might bind to the LARS protein. Clinical data analysis revealed that LARS expression is higher in HCC tissues. High LARS expression was significantly associated with shorter overall survival (OS) in patients and correlated positively with various clinical features like tumor stage, grade, and TP53 mutation status. CONCLUSION: LARS helped HCC become worse by affecting the PI3K/AKT/mTOR pathway and working with RPS5. High LARS meant a worse outcome for patients. This suggested LARS could be used to predict disease or as a treatment target in HCC.

Carcinoma, Hepatocellular

A Computational Workflow for Prioritizing Microbial Metabolite-Associated Host Genes in Constipation-Predominant Irritable Bowel Syndrome.

No standardized computational pipeline exists for systematically prioritizing microbial metabolite-associated host genes and protein-ligand complexes from publicly available chemical, genomic, and structural databases. This article describes an eight-stage workflow that accepts a user-defined set of gut microbiota-derived metabolites and produces a ranked shortlist of candidate metabolite-associated host genes, enriched biological pathways, and structurally prioritized protein-ligand complexes for experimental follow-up. The pipeline integrates (i) chemoinformatic metabolite profiling; (ii) multi-database candidate target prediction using protein-chemical interaction and ligand-based target-prediction tool and a molecular docking program; (iii) differential gene expression analysis of publicly available transcriptomic data; (iv) target-differentially expressed gene overlap; (v) protein-protein interaction network construction and pathway enrichment; (vi) molecular docking with a molecular docking program; (vii) 200 ns molecular dynamics simulation using a molecular dynamics engine with a protein force field used for molecular dynamics simulations; and (viii) MM-PBSA binding free-energy estimation. As a worked example, nine gut microbiota-derived or microbiota-modified metabolites representing short-chain fatty acids, bile acids, tryptophan-derived metabolites, and urolithin A were processed using the public IBS-C rectal mucosal transcriptomic dataset GSE36701. The workflow ranked 17 unique predicted metabolite-associated genes that were differentially expressed in this dataset. Docking, molecular dynamics simulation, and MM-PBSA analyses structurally prioritized five metabolite-protein complexes: lithocholic acid-VDR, lithocholic acid-NR1H4/FXR, ursodeoxycholic acid-NR1H4/FXR, tryptamine-HTR2A (simulated in an explicit 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer), and urolithin A-CASP3. The protocol is designed to be adaptable to other metabolite sets, disease transcriptomic datasets, and target classes; all outputs are hypothesis-generating computational predictions that require independent transcriptomic replication, protein-level validation, and functional ligand-response assays before causal or therapeutic conclusions can be drawn.

Irritable Bowel Syndrome

Novel environmental contaminant 6PPD-quinone promotes malignant phenotypes in colorectal cancer cells and identifies candidate response-associated genes.

6PPD-quinone (6PPDQ), an oxidative transformation product of the widely used tire antioxidant 6PPD, is a ubiquitous environmental contaminant with bioaccumulation potential and widespread human exposure. Recent epidemiological evidence indicates a positive association between urinary 6PPDQ levels and colorectal cancer (CRC) risk; however, its biological effects on CRC-related phenotypes and associated molecular responses remain unclear. We integrated bioinformatics analysis, prognostic modeling, molecular docking and dynamics simulations, and in vitro experiments to investigate cellular and molecular responses to 6PPDQ in CRC models. Predicted 6PPDQ targets were intersected with CRC prognosis-related genes from The Cancer Genome Atlas, followed by functional enrichment and LASSO regression to construct a prognostic risk model, with 1-, 3-, and 5-year AUC values of 0.727, 0.754, and 0.778, respectively. Molecular docking and 100-ns molecular dynamics simulations suggested interactions between 6PPDQ and candidate proteins, including CPT2, SHC2, SRMS, and STK35. Functional assays showed that 6PPDQ exposure altered proliferation, wound-closure capacity, and invasion in Caco-2 and HCT116 cells across the nanomolar concentration range, with non-monotonic and cell-line-dependent responses. In contrast, NCM460 cells showed no increase in EdU incorporation at 10 or 100 nM, whereas reduced proliferation at higher concentrations was accompanied by increased LDH release. 6PPDQ also altered the expression of several prognosis-associated candidate genes. These findings identify cellular phenotypes and candidate molecular responses associated with 6PPDQ exposure under the tested in vitro conditions, but do not establish their causal roles or in vivo relevance. Further mechanistic and in vivo studies are required.

Humans

The Mechanism of Celosiae Semen in the Treatment of Diabetic Cataract: Based on Network Pharmacology.

INTRODUCTION: Diabetes mellitus can be complicated by a variety of ocular diseases, among which the postoperative complications of diabetic cataract (DC) are significantly higher than those of non-DC patients. Therefore, finding drugs with natural active ingredients is an urgent challenge in the prevention and treatment of DC. Discovering the potential molecular mechanism of celosiae semen (CS) for the treatment of DC and providing new ideas and programs for the treatment and prevention of DC. METHODS: In this study, network pharmacology, molecular docking, and molecular dynamics simulations were utilized to predict the binding and functional enrichment of the main active ingredients of CS with DC-related targets, and to explore the potential pathways and mechanisms of CS for the treatment of DC. RESULTS: Through database searching and screening, a total of 45 potential targets of CS for the treatment of DC were identified, functionally enriched, and a protein-protein interaction network was constructed, and the key target, SRC, was finally found. The results of molecular docking and molecular dynamics simulation showed that the main active ingredient of CS, stigmasterol, could bind stably to the key target SRC protein. DISCUSSION: This study not only elucidates the phyto-pharmacological basis of CS in DC management but also provides a framework for developing natural product-derived targeted therapies against diabetic ocular complications. The integration of modern genomics and computational chemistry to deconstruct the therapeutic effects of traditional Chinese herbal medicines has great clinical significance in expanding the scope of traditional Chinese medicines for the treatment of DC and promoting precision targeting. However, this requires verification through basic experiments. CONCLUSION: These computational findings suggest that CS may exert its anti-cataract effects through the multi-target modulation of diabetic metabolic pathways and SRC-mediated signaling cascades.

Humans

Exploring the Potential Molecular Targets of Cyanidin-3-O-glucoside for Type 2 Diabetes Mellitus Treatment.

INTRODUCTION: This study aims to elucidate the multi-target molecular mechanism of cyanidin-3-O-glucoside (C3G) in treating Type 2 Diabetes (T2DM) through network pharmacology methods. METHODS: The study was designed to predict the targets of C3G through public databases and to screen for T2DM-related targets. Protein-protein interaction (PPI) network analysis, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the common targets. Core targets were further validated through molecular docking and molecular dynamics (MD) simulations. RESULTS: This research identified a total of 57 potential targets of C3G in the treatment of T2DM. Subsequent PPI analysis identified ALB (Degree=43), AKT1 (Degree=41), and TNF (Degree=41) as the top three hub proteins. Pathway analysis indicated significant involvement in the insulin signaling pathway (P = 4.205×10-9), AMPK signaling pathway (P = 9.582×10-7), and FoxO signaling pathway (P = 1.315×10-6). Molecular docking revealed strong binding affinities between C3G and NOS3 (-9.5 kcal/mol), PPARG (-9.0 kcal/mol), TNF (-8.5 kcal/mol), and INSR (-8.4 kcal/mol). MD simulations further confirmed that the C3G-target complex has excellent binding stability. DISCUSSION: C3G may intervene in the pathological progression of T2DM by regulating key pathways such as insulin sensitivity, inflammatory responses, and oxidative stress. Further studies suggest that INSR and NOS3 may be new targets through which C3G exerts its effects, but their specific mechanisms and in vivo biological functions still need to be elucidated by subsequent experiments. CONCLUSION: C3G may intervene in the progression of T2DM in a multi-pathway synergistic manner by targeting key molecules such as INSR and NOS3.

Anthocyanins

Leveraging bioinformatics approaches for drug repositioning in space radiation protection.

The health effects of space radiation, primarily Galactic Cosmic Rays (GCRs), on humans remain largely unknown, with potential cardiovascular consequences posing a significant threat to astronauts on long-duration spaceflight missions. Currently, there are no established pharmacological countermeasures for GCR exposure. Drug repositioning offers a promising strategy to accelerate pharmaceutical research in space medicine. This study leverages existing bioinformatics techniques to identify and prioritize potential drug candidates associated with proteomic perturbations following simulated GCR exposure using previously published murine cardiac proteomic data. A protein-protein interaction (PPI) network was constructed using the top differentially expressed proteins (DEPs) from murine heart tissue following exposure to 5-ion GCRs as seed nodes, focusing on experimentally supported interactions. Network topology, Markov clustering, and functional enrichment analyses were used to characterize biologically relevant proteins and pathways. Drug-protein interactions were predicted using Drugst.One and mapped to PPI clusters of interest to identify candidate drugs. Selected drug-macromolecule interactions were further explored using CB-Dock2 molecular docking and short-duration molecular dynamics simulations as hypothesis-generating structural assessments. Analysis of a key PPI network cluster consisting of several ATP synthase proteins identified 23 unique drug candidates. These analyses demonstrate a systematic approach for leveraging bioinformatics techniques to identify candidate molecular targets and generate pharmacological hypotheses in the context of space radiation countermeasures. Ultimately, this strategy introduces a hypothesis-generating framework for the prioritization of potential drug candidates for future computational characterization and experimental investigation against spaceflight stressors.

Animals

Exploring the mechanism of Shengmai San in treating lung adenocarcinoma based on bioinformatics and molecular dynamics simulation.

To investigate the mechanism of Shengmai San (SMS) in the treatment of lung adenocarcinoma (LUAD) based on an integrated strategy combining "network pharmacology, bioinformatics, molecular docking, and molecular dynamics simulation," aiming to provide a precise combination therapy strategy and identify potential bioactive compounds. Differentially expressed genes in LUAD were identified from the Gene Expression Omnibus database using R (originally developed at Bell Laboratories and currently managed by Lucent Technologies). SMS components (ginseng, Ophiopogon japonicus, and Schisandra chinensis) were retrieved from encyclopaedia of traditional Chinese medicine, with Lipinski-compliant compounds selected. Compound targets were predicted via SwissTargetPrediction and Similarity Ensemble Approach. Intersecting targets between differentially expressed genes and compound targets were identified for "herbs-compounds-targets-disease" network construction. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed. Hub targets were identified by analyzing the protein-protein interaction network. High-prognostic relevance targets were screened from The Cancer Genome Atlas. Compounds targeting these were identified through the herbs-compounds-targets-disease network, and absorption, distribution, metabolism, excretion, and toxicity-compliant compounds were selected using SwissADME (a web-based tool provided by the Molecular Modeling Group of the Swiss Institute of Bioinformatics). Core regulatory targets were identified through molecular docking, with complex stability assessed by molecular dynamics simulations. The key bioactive compounds of SMS for treating LUAD were identified as 7-hydroxy-2,5-dimethyl-4H-1-benzopyran-4-one, N-trans-feruloyltyramine, paprazine, and (E)-N-[(2S)-2-hydroxy-2-(4-hydroxyphenyl)ethyl]-3-(4-hydroxyphenyl)prop-2-enamide. Hub targets included AURKA, CCNA2, CCNB1, CDK1, CHEK1, KIF11, NEK2, PLK1, TTK, and TYMS. Among these, CDK1, CHEK1, and PLK1 demonstrated both high-prognostic relevance and strong binding affinity with SMS, emerging as core regulatory targets for SMS in LUAD treatment. Mechanistically, SMS exerts its anticancer effects primarily by modulating the tumor necrosis factor, interleukin-17, cell cycle, and Lipid and atherosclerosis signaling pathways. The active components of SMS, such as paprazine, may exert antitumor effects partly through downregulating CDK1, CHEK1, and PLK1 expression. Although the present study did not examine drug-resistance models or combination regimens, our findings raise the possibility that, in patients with high expression of these genes, combining SMS with standard chemotherapy or targeted therapy could potentially enhance chemosensitivity and mitigate the development of resistance. This hypothesis, however, requires formal testing in appropriate preclinical models and functional validation studies.

Molecular Dynamics Simulation

Integrated dual transcriptome sequencing and experimental validation reveal potential mechanisms of baicalin against pneumocystis pneumonia in immunosuppressed rats.

BACKGROUND: Pneumocystis pneumonia (PCP) remains a major cause of morbidity and mortality in immunocompromised individuals. Although baicalin (Ba), a natural bioactive flavonoid, has demonstrated protective and therapeutic effects against PCP, its molecular mechanisms remain undefined. We employed dual RNA sequencing (dual RNA-seq) to characterize host and pathogen transcriptional responses to Ba treatment in an immunosuppressed rat model of PCP. METHODS: Comparative transcriptomic analyses identified differentially expressed genes in both the host and Pneumocystis, followed by Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and gene set enrichment analyses. Candidate targets were further investigated using network pharmacology, protein-protein interaction analysis, molecular docking, and molecular dynamics simulations. Key findings were validated by immunohistochemistry, enzyme-linked immunosorbent assay, and quantitative PCR. RESULTS: Ba markedly remodeled host and pathogen transcriptomes. Host transcriptomic analyses showed that Ba attenuated inflammatory and oxidative stress responses by modulating immune-related pathways, including Toll-like receptor, NF-κB, cytokine-cytokine receptor interaction, chemokine signaling, Th17 cell differentiation, and antigen processing and presentation. Experimental validation demonstrated that Ba reduced pulmonary expression of indoleamine 2,3-dioxygenase 1 (IDO1), Toll-like receptor 2 (TLR2), and TLR4 while increasing nuclear factor erythroid 2-related factor 2 (Nrf2) and its downstream antioxidant enzyme heme oxygenase-1 (HO-1). Pathogen transcriptomic analysis identified Pneumocystis Rtt109 (PcRtt109), a fungal histone acetyltransferase, as a potential pathogen-specific target that was significantly downregulated after Ba treatment. Molecular docking and molecular dynamics simulations supported stable interactions between Ba and IDO1, Nrf2, TLR2, TLR4, and PcRtt109, with the strongest predicted binding observed for PcRtt109. CONCLUSION: Dual RNA-seq revealed that Ba exerts anti-PCP activity through coordinated modulation of host and pathogen molecular networks. Its therapeutic effects are associated with suppression of inflammatory signaling, enhancement of antioxidant defenses, and inhibition of a fungal virulence-associated target. These findings provide mechanistic insights into host-pathogen interactions during PCP and support Ba as a potential therapeutic candidate for PCP.

Nrf2

Targeting EGFR in cancer using Terminalia arjuna: An integrated In Silico, molecular dynamics, experimental validation, and network pharmacology study.

The Epidermal Growth Factor Receptor (EGFR) plays a pivotal role in 20-60% of cancer cases, including glioblastoma, lung adenocarcinoma, and head and neck squamous cell carcinoma, as reported in The Cancer Genome Atlas (TCGA) dataset. The present study employed an integrated in silico and experimental workflow to evaluate EGFR-targeted compounds from Terminalia arjuna. Drug-likeness and ADMET screening were performed, followed by molecular docking and 1000 ns molecular dynamics simulations. In vitro validation was conducted using cancer cell-based assays and network pharmacology to explore the molecular mechanisms associated with the identified compound. Screening shortlisted eight compounds from T. arjuna. Molecular docking identified Arjunaside C (-8.2 kcal/mol), Arjunapthanoloside (-7.7 kcal/mol), and Beta-sitosterol (-7.4 kcal/mol) as potential EGFR inhibitors compared to Erlotinib (-6.6 kcal/mol). Arjunapthanoloside formed more H-bonds and exhibited most stable interactions with EGFR. MD simulations at 1000 ns revealed lower RMSD, RMSF, SASA, and Rg values for the Arjunapthanoloside-EGFR complex, indicating enhanced stability. Direct binding validation was limited by the unavailability of purified Arjunapthanoloside; therefore, Arjuna extract was evaluated, which demonstrated potent cytotoxicity with an IC₅₀ of 9 µg/mL in H357 oral cancer cells. Flow cytometry confirmed apoptosis-mediated cell death by increased early- and late-apoptotic cell populations. Network pharmacology analysis further identified additional targets (MMP3, MMP7, MMP9, and HRAS) that are directly involved in various cancers. Overall, the findings provide new insights into the therapeutic potential of Arjunapthanoloside as a stable compound that interacts with EGFR from T. arjuna, highlighting its significance in EGFR-targeted anticancer research.

ErbB Receptors

Investigating the mechanisms of PhIP-induced colorectal cancer through network toxicology, machine learning, and molecular dynamics simulation.

BACKGROUND: Over the past few years, 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP)- a compound from grilled or processed meats-has emerged as a major player in cancer development, especially colorectal cancer (CRC). This work dives into its potential links to CRC and uncovers the key genes that bridge this connection. METHODS: We tapped into various databases to pinpoint target genes tied to PhIP and CRC, then ran protein-protein interaction (PPI) analyses for visualization. Next, we explored underlying mechanisms through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. To nail down predictions, we tested 107 machine learning pipelines and picked the best one, validating its accuracy and the core genes' prognostic value across datasets. Next, molecular docking and dynamics simulations probed the interactions between these genes and PhIP. Finally, cell proliferation was assessed using Cell Counting Kit-8 (CCK-8) and 5-ethynyl-2'-deoxyuridine (EdU) assays, and polymerase chain reaction (PCR) was performed to validate the expression levels of the hub genes. RESULTS: Our analysis identified 39 overlapping genes, from which a machine learning model (glmBoost + Enet) identified six candidate targets: CDK4, CEBPB, COMT, SOX9, TIMP1, and TOP2A. To prioritize these, a hierarchical screening framework was applied. Molecular docking and dynamics simulations identified CDK4, COMT, and TIMP1 as the most stable interactors with PhIP. Functional assays confirmed that PhIP treatment significantly enhanced the proliferation of CRC cells. Crucially, quantitative PCR (qPCR) validation in multiple CRC cell lines identified TIMP1 as the primary target, showing the most consistent and significant upregulation upon PhIP exposure. CONCLUSIONS: In essence, these genes drive PhIP is role in CRC, offering novel insights into its molecular pathways. This could reshape how we tackle food-related pollutants, paving the way for better prevention and targeted therapies.

Colorectal cancer (CRC)