SEARCH · Search PubMed

Results for “Molecular Dynamics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

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

At least 19 recordsLinked to original sources

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

Sanguinarine as a multi-target therapeutic candidate for laryngeal cancer: insights from network pharmacology, molecular dynamics and in vitro validation.

OBJECTIVE: To identify the core targets and elucidate the potential molecular mechanisms of sanguinarine (SA) against laryngeal squamous cell carcinoma (LSCC), and to validate its antitumor effects in vitro. METHODS: Potential targets of SA were predicted using SwissTargetPrediction, TargetNet, and SuperPred and intersected with LSCC-related targets obtained from the GeneCards, OMIM, and DISEASES databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed. A protein-protein interaction (PPI) network was constructed using the STRING database (combined score > 0.900), and topological parameters including degree centrality (DC), betweenness centrality (BC), closeness centrality (CC), eigenvector centrality (EC), and local average connectivity (LAC) were calculated in Cytoscape to identify core genes based on median thresholds. Molecular docking and 100-ns molecular dynamics (MD) simulations were conducted for epidermal growth factor receptor (EGFR), Phosphatidylinositide-3-kinase catalytic subunit alpha (PIK3CA), phosphatidylinositol-4,5-biphosphate 3-kinase catalytic subunit β (PIK3CB), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit delta (PIK3CD), and Non-Receptor Tyrosine Kinase (SRC). The effects of SA on LSCC were evaluated using CCK-8, colony formation, Transwell migration, and wound-healing assays in TU177 cells and TU212. RESULTS: A total of 213 common targets were identified, which were significantly enriched in PI3K-Akt signaling, EGFR tyrosine kinase inhibitor resistance, and adhesion- and migration-related pathways. The PPI network comprised 259 nodes and 259 edges, from which five core genes-PIK3CA, PIK3CB, PIK3CD, EGFR, and SRC-were identified. Molecular docking revealed strong binding affinities between SA and the PI3K family proteins (- 9.79 to - 10.96 kcal/mol), as well as EGFR (- 8.58 kcal/mol) and SRC (- 6.77 kcal/mol). MD simulations indicated greater stability of SA complexes with EGFR and PI3K family members compared with SRC. In vitro assays demonstrated that SA significantly inhibited TU177 cell and TU212 cell proliferation, colony formation, and migration. CONCLUSION: SA may exert anti-laryngeal cancer effects through synergistic multi-target inhibition centered on the EGFR/SRC/PI3K signaling axis, highlighting its potential as a promising therapeutic candidate for LSCC.

Humans↗

PubMed2023

A new horizon in the phosphorylated sites of AGA: the structural impact of C163S mutation in aspartylglucosaminuria through molecular dynamics simulation.

Aspartylglucosaminuria (AGU) is a lysosomal storage disorder caused by insufficient aspartylglucosaminidase (AGA) activity leading to chronic neurodegeneration. We utilized the PhosphoSitePlus tool to identify the AGA protein's phosphorylation sites. The phosphorylation was induced on the specific residue of the three-dimensional AGA protein, and the structural changes upon phosphorylation were studied via molecular dynamics simulation. Furthermore, the structural behaviour of C163S mutation and C163S mutation with adjacent phosphorylation was investigated. We have examined the structural impact of phosphorylated forms and C163S mutation in AGA. Molecular dynamics simulations (200 ns) exposed patterns of deviation, fluctuation, and change in compactness of Y178 phosphorylated AGA protein (Y178-p), T215 phosphorylated AGA protein (T215-p), T324 phosphorylated AGA protein (T324-p), C163S mutant AGA protein (C163S), and C163S mutation with Y178 phosphorylated AGA protein (C163S-Y178-p). Y178-p, T215-p, and C163S mutation demonstrated an increase in intramolecular hydrogen bonds, leading to greater compactness of the AGA forms. Principle component analysis (PCA) and Gibbs free energy of the phosphorylated/C163S mutation structures exhibit transition in motion/orientation than Wild type (WT). T215-p may be more dominant among these than the other studied phosphorylated forms. It might contribute to hydrolyzing L-asparagine functioning as an asparaginase, thereby regulating neurotransmitter activity. This study revealed structural insights into the phosphorylation of Y178, T215, and T324 in AGA protein. Additionally, it exposed the structural changes of the C163S mutation and C163S-Y178-p of AGA protein. This research will shed light on a better understanding of AGA's phosphorylated mechanism.Communicated by Ramaswamy H. Sarma.

Molecular Dynamics Simulation↗

PubMed2026

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↗

PubMed2026

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↗

PubMed2026

Integrative computational analysis combining network pharmacology, regulatory network modeling, and molecular dynamics reveals the mechanisms of Quanshen compound in ITP.

UNLABELLED: Immune thrombocytopenia (ITP) is a hemorrhagic disorder caused by immune dysfunction. Quanshen Compound (QSC) is an in-house preparation developed by the Uyghur Hospital in Hotan Prefecture. This study primarily investigates and validates the potential pharmacological basis and mechanism of action of QSC in modulating immune thrombopoiesis. Based on the multi-database screening of the QSC and the related targets of ITP, the intersection was obtained to construct a protein-protein interaction (PPI) network and screen the core targets; the intersection targets were analyzed for gene ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using R packages; a component-target-pathway network was constructed to screen the key active components and their mechanisms of action. At the same time, the TF-mRNA-miRNA regulatory network of the core targets was constructed, and chromosome localization and subcellular localization analysis were performed; further, the binding stability of key components and core targets was verified through molecular docking and molecular dynamics simulation. A total of 227 potential target sites were screened out, among which TNF, IL6, AKT1, TP53 and IL1B were the core targets. The enrichment results indicated that these intersecting target sites mainly participated in inflammatory responses, immune regulation and hemostasis-related biological processes, and were significantly enriched in the PI3K-Akt signaling pathway, Toll-like receptor signaling pathway, Th17 cell differentiation and PD-1/PD-L1 signaling pathway. The core target TF-mRNA-miRNA regulatory network contained 184 nodes and 200 edges, suggesting that the core targets were subject to multi-level regulation. Molecular docking results showed that the main active components had good binding activity with the core targets, and molecular dynamics simulation further verified the stability of the complex. QSC may improve ITP through a multi-component, multi-target, and multi-pathway synergistic mechanism involving key targets such as TNF, IL6, AKT1, TP53, and IL1B, as well as the PI3K-Akt signaling pathway. These findings provide new insights into the potential therapeutic mechanisms of QSC against ITP and warrant further experimental validation. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at https://doi.org/10.1007/s40203-026-00718-0.

Immune thrombocytopenia↗

PubMed2026

Molecular dynamics simulations of positively selected codons in FcγRI reveal novel biochemical binding properties.

FcγRI is a high-affinity receptor for IgG, associated with autoimmune disease pathology and determines clinical responses to antibody-based immunotherapies. FcγRI has a complex evolutionary history that is not fully understood, and to address this we explored signatures of positive selection in the receptor's functional gene, FCGR1A, using codon-based selection tests on aligned 1-1 orthologous sequences from placental mammals (n = 32). Signatures of positive selection have occurred at several locations within the gene, with two sites (H148 (M2a ω 0.997 & M8 ω = 0.993)) and (W149 (M2a ω = 0.999 & M8 ω = 1.000)) exhibiting highest posterior probabilities, suggesting strong evidence of positive selection; these positions are known to form one of the FcγRI-IgG binding interfaces. We employed ancestral reconstruction to statistically infer prior codon sequences at these sites and identified ancestral H148P and W149R codons at different nodes in the phylogeny. Employing molecular dynamics simulations, we determined how evolutionary changes at these sites may have influenced the binding of FcγRI-IgG of modern-day Homo sapiens. Measuring RMSD, free energy, radius of gyration, hydrogen bond formation, and analyzing free energy landscapes, we demonstrate that structural instability between mutant structures vs the WT counterpart; however, overall binding potential increases at position 148, yet decreases at 149 in potential. H148P protonation at physiological pH remains similar, yet during acidotic calculations, protonation is likely reduced, with predicted reduction in affinity for IgG. While ancestral W149R substitutions demonstrate an implication for electron conjugation. Examining key sites at this binding FcγRI-IgG interface, our data demonstrate that these two codons have evolved in humans to be relatively insensitive to shifts in pH promoting a more stable interaction with the Fc portion of IgG during diseases that promote acidosis.

Receptors, IgG↗

PubMed2026

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↗

PubMed2025

Exploring potential targets and molecular mechanisms of traumatic brain injury exacerbated by Benzo(a)pyrene via network toxicology and molecular dynamics simulation.

Benzo(a)pyrene (BaP) is a common environmental pollutant from combustion sources that promotes oxidative stress, neuroinflammation and disruption of blood-brain barrier (BBB). However, its contribution to worsening traumatic brain injury (TBI) remains unclear. In this study, we aimed to assess the contribution of BaP to secondary injury in TBI. By integrating data from e.g., the Comparative Toxicogenomics Database, GeneCards, and Online Mendelian Inheritance in Man, 121 overlapping core targets were identified between BaP and TBI. Enrichment analyses via Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, combined with protein-protein interaction networks and topological algorithms (degree, closeness centrality, betweenness centrality, average shortest path length, topological coefficient and partner of multi-edged node pairs), highlighted five hub genes (TP53, EGFR, AKT1, ACTB, and TNF) implicated in mitogen-activated protein kinase signaling, oxidative stress, and neuroinflammation. Molecular docking showed strong binding affinities of BaP to these hub proteins, with energies from -9.3 to -12.1 kcal/mol, tighter than co-crystal ligands and existing protein-binding drugs. Molecular dynamics simulations confirmed interaction stability through low root-mean-square deviation (< 0.5 nm), fluctuation, and radius of gyration values. Calculation of binding free energies using MM-PBSA validated the strong binding affinity between BaP and binding pockets of each hub genes. Toxicity prediction analysis revealed an oral LD50 of 316 mg/kg for BaP, with high probabilities for neurotoxicity, BBB permeability, carcinogenicity, and mutagenicity, associated with aryl hydrocarbon receptor activation. These findings reveal a "neurovascular homeostasis disruption" network underlying BaP-exacerbated TBI pathology and highlight potential targets to reduce pollution-related risks in TBI management.

Benzo(a)pyrene↗

PubMed2025

Molecular dynamics simulations reveal subtle consequences of H3K9 and H3K27 tri-methylation on chromatin constituents.

Epigenetic modifications of histone tails are key mechanisms of genome regulation. In particular, tri-methylation of lysines (K) 9 and K27 of the histone H3 tail is important for genome silencing. In this work, we explore, using all-atom molecular dynamics simulations, the effect of these two epigenetic marks on the structure and interactions of the H3 tail in several contexts: isolated tails, nucleosomes, chromatosomes, and stacked nucleosomes. Overall, we find that although the isolated tails do not show significant conformational changes upon methylation, a more flexible and extended H3 tail compared to the native tail results in the nucleosome systems, with K9 methylation effects more pronounced. This change could facilitate the interaction of the tail with protein readers like heterochromatin protein 1 or Polycomb group. We also observe that both methylations increase the interactions of the H3 tail with the linker DNA in the context of the chromatosome, producing a chromatosome with tighter linker DNA, which could favor chromatin compaction. For stacked nucleosomes mimicking i±2 zigzag interactions, methylation of either K9 or K27 reduces the interactions of one of the H3 tails with its parental nucleosome and increases its interactions with the nonparental nucleosome, which could also help compact the chromatin fiber. In the three nucleosome-containing systems, we observe an asymmetry between the two tails, especially in the chromatosome, where one tail extends to interact with the linker DNA. This asymmetry modulates the effect that methylation has on each tail. Thus, overall, methylations of K9 and K27 have a subtle but notable impact on the H3 tail structure and its interactions within the chromatin fiber. These results help explain how this epigenetic modification compacts chromatin fibers and promotes longer-range interactions; these changes also guide how to approximate these effects in coarse-grained chromatin models.

Histones↗

PubMed2026

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)↗

PubMed2025

Dynamic Molecular Changes in Brain, Lung, and Heart of Hamsters Infected With SARS-CoV-2: Insights From a Severe and Recovery Phase Model.

The Global pandemic of coronavirus disease 2019 was initiated by the emergence of severe acute respiratory syndrome coronavirus 2. In addition to conventional pulmonary lesions, a range of neurological injury symptoms have been identified in clinical practice, but the aetiology of neurological disorders linked to SARS-CoV-2 infection remains poorly understood. Syrian hamsters, which are highly susceptible to SARS-CoV-2 infection, exhibit a disease phenotype similar to that observed in human COVID-19 patients. In this study, a hamster model of COVID-19 infection was used to analyze molecular changes in different tissues at various time points post infection with distinct strains using proteomic and phosphoproteomic approaches. Multi-omics analysis showed that SARS-COV-2 infection triggers sustained downregulation of the abundance and phosphorylation levels of neuronal and synapse-associated proteins in the brain, suggesting that neuronal damage persists even during the recovery period. Additionally, infections with SARS-CoV-2 may contribute to the onset of long-term symptoms of COVID-19 by impacting energy metabolism, neurotransmitter release, and synaptic transmission pathways. This study provides a comprehensive molecular profile of hamsters infected with different SARS-CoV-2 strains in different tissues, offering foundational insights into the pathogenic mechanisms of COVID-19.

Animals↗

PubMed2026

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↗

PubMed2026

Genomic and structural analysis of dacB variants associated with cephalosporin resistance in Pseudomonas aeruginosa.

The rise of resistance to fourth-generation cephalosporin in Pseudomonas aeruginosa (P. aeruginosa) is a global concern. The resistance is largely driven by variants of chromosomally encoded AmpC β-lactamase, known as Pseudomonas-derived cephalosporinase (PDC), which arise from the mutations in the ampC gene. In addition, alteration in dacB, which encode the penicillin-binding protein 4 (PBP4), can lead to the overexpression of ampC, thereby contributing to β-lactam resistance. Present work analyzed 208 clinical isolates of P. aeruginosa using whole-genome sequencing (WGS) and detected multiple nonsynonymous single nucleotide polymorphisms (nsSNPs), such as Y264C, G444D, and a double mutation (A394P-T428P). All nsSNPs were predicted to be deleterious by several prediction program. Molecular dynamics (MD) simulations suggested that these substitutions destabilize PBP4, increase structural flexibility, and contribute to the resistance mechanism, which favored their selection. To determine the effective therapeutics against these mutations, molecular docking was conducted with various antibiotics. Cefoperazone exhibited the highest binding affinity (-7.3 kcal/mol) among multiple PBP4 variants. The Molecular dynamics (MD) simulations and Molecular Mechanics Poisson Boltzmann Surface Area calculations (MMPBSA) further confirmed the favorable interactions between cefoperazone and PBP4 variants. In vitro MIC analyses supported these findings, indicating that cefoperazone displayed significant activity against clinical dacB mutants of P. aeruginosa. The study offers structural insight of dacB variants leading to antibiotic resistance and emphasizes the need to prioritize specific antibiotics to address the challenges arising from protein mutations.

Pseudomonas aeruginosa↗

PubMed2026

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↗

PubMed2026

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↗