SEARCH · PubMed Health

Results for “Molecular Docking 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 217 records · Page 12Linked to original sources

PubMed2026

Domesticated Argania spinosa in Eastern Morocco: HPLC-DAD/GC-MS Chemical Profiling, Antioxidant and Antidiabetic Activities, and Network Pharmacology-Guided Molecular Docking.

The argan tree (Argania spinosa) is an endemic Moroccan species known for its primary product, argan oil, which possesses exceptional nutritional and medicinal properties. The current study aimed to evaluate and compare the antidiabetic and antioxidant activities of argan oil obtained from the introduced and native argan tree in eastern Morocco, to analyze its chemical composition using HPLC-DAD and GC-MS, and to investigate the molecular mechanisms behind the obtained pharmacological activities through an in silico pharmacological networking and molecular docking study. The results revealed that argan oil from all three regions of Morocco (Oujda, Agadir, and Chouihya) is rich in oleic and linoleic acids as major constituents, along with the presence of significant tocopherols. Regarding the antioxidant assays, including DPPH radical scavenging and iron-reducing power tests, argan oil from Oujda exhibited the highest activity, with the lowest IC50 values of 15.25 ± 0.022 mg/mL and 28.5 ± 1.7 mg/mL, respectively. Concerning the antidiabetic activity, we found that oil from Chaouihya showed the strongest α-amylase inhibition, while Oujda oil had the highest antiglycation activity, indicating that even introduced argan trees retain potent bioactivity. The results of the in silico investigation suggested that tocopherols may contribute to the antioxidant and antidiabetic potential of argan oil, showing predicted antioxidant activity (Pa = 0.843-0.967) and favorable binding affinities toward iNOS (ΔG = -9.3 kcal mol-1) and α-glucosidase (ΔG = -8.2 kcal mol-1). The identified fatty acids also showed predicted insulin-promoting activity (Pa = 0.59-0.75) and moderate enzyme-binding potential. Pharmacological network analysis identified 51 shared genes associated with antioxidant, antidiabetic, and argan-related targets, with enrichment of the AGE-RAGE signaling pathway. These computational findings provide possible molecular associations that may help explain the observed biological activities, although they remain predictive and require experimental validation. Overall, the in silico analysis suggests that tocopherols could be among the contributors to the multi-target profile of Argania spinosa oil, while fatty acids may provide complementary effects related to glycemic regulation.

Sapotaceae↗

PubMed2026

A Network Pharmacology and Molecular Docking Study of TongBi Formula for Osteoarthritis.

This study applied network pharmacology combined with molecular docking to predict the potential therapeutic targets and molecular mechanisms of TongBi Formula (TBF) in osteoarthritis (OA). Active components and corresponding targets of TBF were retrieved from the traditional Chinese medicine Systems Pharmacology Database and Analysis Platform, while OA-related targets were collected from Online Mendelian Inheritance in Man, GeneCards, DrugBank, and Therapeutic Target Database. A network visualization and analysis software was used to construct compound-target and protein-protein interaction (PPI) networks. Gene Ontology functional annotation and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed using the Database for Annotation, Visualization and Integrated Discovery platform. Molecular docking analysis was conducted using a molecular docking software to evaluate the predicted binding affinity between key active compounds and core target proteins. A total of 47 overlapping targets between TBF and OA were identified. PPI network analysis highlighted JUN, RELA, IL6, MAPK1, and IL10 as potential hub targets. Enrichment analysis suggested that TBF may regulate inflammation, lipid metabolism, and multiple intracellular signaling pathways associated with OA progression. Molecular docking results demonstrated favorable predicted binding affinities between core active compounds and key OA-related protein targets. These findings provide a computational framework for understanding the potential mechanisms of TBF against OA and support further experimental validation.

Molecular Docking Simulation↗

PubMed2026

Integrated Network Pharmacology and Molecular Docking Analysis of Sishen Decoction Identifies Potential Targets and Pathways in Gout.

Gout is a disease characterized by hyperuricemia and the deposition of urate crystals in joints and soft tissues, leading to recurrent acute arthritis. Its increasing prevalence imposes substantial clinical and socioeconomic burdens. Sishen Decoction (SSD) has been used in the treatment of gout, but its potential molecular mechanisms remain unclear. This study applied an integrated network pharmacology and molecular docking approach to identify potential targets and signaling pathways associated with SSD in gout. Active compounds and corresponding targets of SSD were retrieved from the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP), while gout-related targets were collected from the GeneCards and Online Mendelian Inheritance in Man (OMIM) databases. Overlapping targets were identified and used to construct a drug-component-target-disease network. A protein-protein interaction (PPI) network was established using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed, followed by molecular docking using the docking server analysis module. A total of 37 bioactive compounds were associated with 116 overlapping gout-related targets. The top hub targets included TP53, IL6, IL1B, TNF, AKT1, EGFR, CASP3, JUN, BCL2, and MMP9. GO analysis suggested that these targets are involved in gene expression regulation and signal transduction. KEGG enrichment analysis indicated significant associations with the mitogen-activated protein kinase (MAPK), phosphoinositide 3-kinase/protein kinase B (PI3K-Akt), interleukin-17 (IL-17), and tumor necrosis factor (TNF) signaling pathways. Molecular docking predicted favorable interactions between key compounds and hub targets, with all binding energies of ≤-5 kcal/mol. These computational findings provide potential mechanistic hypotheses for the action of SSD in gout and may support future experimental validation.

Molecular Docking Simulation↗

PubMed2026

[Effects and mechanisms of Jiawei Yigong San on the Th17/Treg balance in food allergy].

Objective To explore the effects and mechanisms of Jiawei Yigong San (JWYGS) on the T helper type 17 (Th17)/regulatory T (Treg) cell balance in food allergy (FA). Methods Active components, action targets of JWYGS, and FA-related targets, were screened via network pharmacology. Overlapping targets between JWYGS and FA were used to construct a protein-protein interaction (PPI) network. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were performed to predict key signaling pathways. Molecular docking was conducted to validate the binding affinity between the main active components and the predicted targets. Mice were randomly divided into control group, model group, JWYGS low-dose, medium-dose, and high-dose groups, and dexamethasone (DXM) group. An ovalbumin (OVA)-induced FA model was established. During the OVA challenge period, mice received daily intragastric administration, after which allergy and diarrhea scores were assessed. Small intestinal pathology was evaluated by HE staining. Serum ovalbumin-specific immunoglobulin E (OVA-sIgE), interleukin 6 (IL-6), IL-17, IL-2, and IL-10 were measured by ELISA. Small intestinal IL-6, IL-17, and IL-10 protein expression was detected by immunohistochemistry. Splenic Th17 and Treg cells were quantified by flow cytometry, and the Th17/Treg ratio was calculated. The mRNA expression of IL-6, retinoic acid receptor-related orphan receptor γt (RORγt), and forkhead box protein P3 (FOXP3) in the small intestine was detected by qPCR. Results Network pharmacology identified 235 active components of JWYGS and 136 common targets. GO and KEGG enrichment analyses highlighted cytokine response and Th17 differentiation. Molecular docking confirmed stable binding between core components and targets. Compared with the control group, the model group exhibited aggravated allergy and diarrhea scores, marked small intestinal inflammation and mucosal damage, elevated serum levels of OVA-sIgE, IL-6, IL-17 and IL-2, along with increased splenic Th17 cell frequency and Th17/Treg ratio. Intestinal IL-6 and IL-17 protein levels as well as IL-6 and RORγt mRNA expression were upregulated, whereas serum IL-10 levels were decreased, and intestinal expression of IL-10 protein and FOXP3 mRNA was downregulated. After JWYGS treatment, allergy and diarrhea scores were significantly reduced. Small intestinal inflammation and mucosal damage were alleviated. Serum levels of OVA-sIgE, IL-17, IL-6 and IL-2, Th17 cell frequency and the Th17/Treg ratio, intestinal IL-6 and IL-17 protein levels were decreased. IL-6 and RORγt mRNA expression was downregulated. Serum IL-10 levels were increased and intestinal expression of IL-10 protein and FOXP3 mRNA was upregulated. Moreover, the JWYGS high-dose group demonstrated comparable efficacy to the DXM group. Conclusion JWYGS can ameliorate symptoms and reverse the Th17/Treg imbalance in FA mice, potentially by inhibiting IL-6 transcription and regulating RORγt/FOXP3 expression.

Animals↗

PubMed2006

Molecular dynamics simulation of the P2Y14 receptor. Ligand docking and identification of a putative binding site of the distal hexose moiety.

A rhodopsin-based homology model of the P2Y14 receptor was inserted into a phospholipid bilayer and refined by molecular dynamics (MD) simulation. The binding modes of several known agonists, namely UDP-glucose and its analogues, were proposed using automatic molecular docking combined with Monte Carlo Multiple Minimum calculations. Compared to other P2Y receptors, the P2Y14 receptor has an atypical binding mode of the nucleobase, ribose, and phosphate moieties. The diphosphate moiety interacts with only one cationic residue, namely Lys171 of EL2, while in other P2Y receptor subtypes three Arg or Lys residues interact with the phosphate chain. Two other conserved cationic residues, namely Arg253 (6.55) and Lys277 (7.35) of the P2Y14 receptor together with two anionic residues (Glu166 and Glu174, located in EL2), are likely involved in interactions with the distal hexose moiety.

Binding Sites↗

PubMed1999

Docking of flexible ligands to flexible receptors in solution by molecular dynamics simulation.

In this paper, a method of simulating the docking of small flexible ligands to flexible receptors in water is reported. The method is based on molecular dynamics simulations and is an extension of an algorithm previously reported by Di Nola et al. (Di Nola et al., Proteins 1994;19:174-182). The method allows a fast exploration of the receptor surface, using a high temperature of the center of mass translational motion, while the ligand internal motions, the solvent, and the receptor are simulated at room temperature. In addition, the method allows a fast center of mass motion of the ligand, even in solution. The dampening effect of the solvent can be overcome by applying different weights to the interactions between system subsets (solvent, receptor, and ligand). Specific ligand-receptor distances have been used to compare the results of the simulations with the crystal structure. The method is applied, as a test system, to the docking of the phosphocholine to the immunoglobulin McPC603. The results show the similarity of structure between the complex in solution and in the crystal.

Algorithms↗

PubMed2001

Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA.

In this work, a new ansatz is presented that combines molecular dynamics simulations with MM-PBSA (Molecular Mechanics Poisson-Boltzmann/surface area) to rank the binding affinities of 12 TIBO-like HIV-1 RT inhibitors. Encouraging results have been obtained not only for the relative binding free energies, but also for the absolute ones, which have a root-mean-square deviation of 1.0 kcal/mol (the maximum error is 1.89 kcal/mol). Since the root-mean-square error is rather small, this approach can be reliably applied in ranking the ligands from the databases for this important target. Encouraged by the results, we decided to apply MM-PBSA combined with molecular docking to determine the binding mode of efavirenz SUSTIVA(TM) another promising HIV-1 RT inhibitor for which no ligand-protein crystal structure had been published at the time of this work. To proceed, we define the following ansatz: Five hundred picosecond molecular dynamics simulations were first performed for the five binding modes suggested by DOCK 4.0, and then MM-PBSA was carried out for the collected snapshots. MM-PBSA successfully identified the correct binding mode, which has a binding free energy about 7 kcal/mol more favorable than the second best mode. Moreover, the calculated binding free energy (-13.2 kcal/mol) is in reasonable agreement with experiment (-11.6 kcal/mol). In addition, this procedure was also quite successful in modeling the complex and the structure of the last snapshot was quite close to that of the measured 2,3 A resolution crystal (structure the root-mean-square deviation of the 54 C(alpha) around the binding site and the inhibitor is 1.1 A). We want to point out that this result was achieved without prior knowledge of the structure of the efavirenz/RT complex. Therefore, molecular docking combined with MD simulations followed by MM-PBSA analysis is an attractive approach for modeling protein complexes a priori.

Alkynes↗

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

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↗

PubMedDate not supplied

Structure based prediction of binding affinity of human immunodeficiency virus-1 protease inhibitors.

A series of computations were performed to derive a strategy for the prediction of binding affinities of non-peptidic human immunodeficiency virus-1 (HIV-1) protease inhibitors. This paper describes the development of a 3D quantitative structure-activity relationship (3D-QSAR) methodology by using receptor information of HIV-1 protease. The docking and molecular dynamics simulations were performed on a model ligand/enzyme complex to optimize the variables involved in the generation of ligand/enzyme models. The protonation scheme of the active site aspartic acid residues of HIV-1 protease was derived from a computational study. The active site aspartate is monoprotonated with a proton placed on the OD1 atom of the ASP B25. This protocol of docking and molecular dynamics (MD) simulation was then used to derive the ligand-enzyme complexes of the molecules used in the present study. The molecular mechanics interaction descriptors were calculated from these ligand/enzyme models. A partial least squares (PLS) method was used to derive a linear correlation between the interaction descriptors and the biological activity. A good correlation was observed when the change in the energy of the ligand upon complex formation and the electrostatic contributions to the solvation energy of the ligand were included in the QSAR analysis. A highest cross-validated q2 value of 0.649 was observed. This model had a conventional r2 of 0.725, and when this model was used to predict the activity of the external test set, it produced a r2pred of 0.761. The total interaction energy was partitioned into interactions in different subsites and interactions with each of the amino acid residues of the enzyme. The PLS analysis using these descriptors helped to identify the important interactions which can be exploited for the design of HIV-1 protease inhibitors.

HIV Protease↗

PubMed2026

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↗

PubMed2003

Targeting aurora2 kinase in oncogenesis: a structural bioinformatics approach to target validation and rational drug design.

The aurora kinases are a novel oncogenic family of mitotic serine/threonine kinases (S/T kinases) that are overexpressed in a number of solid tumors, including pancreas and colorectal cancer. A PSI-BLAST search [National Center for Biotechnology Information (NCBI)] with the sequence of the S/T kinase domain of human aurora1 kinase [also known as AUR1, ARK2, AIk2, AIM-1, and STK12] and human aurora2 kinase (also known as AUR2, ARK1, AIK, BTAK, and STK15) showed a high sequence similarity to the three-dimensional structures of bovine cAMP-dependent kinase [Brookhaven Protein Data Bank code 1CDK], murine cAMP-dependent kinase (1APM), and Caenorhabditis elegans twitchin kinase (1KOA). When the aurora1 or aurora2 sequence was input into the tertiary structure prediction programs THREADER and 3D-PSSM (three-dimensional position-sensitive scoring matrix), the top structural matches were 1CDK, 1APM, and 1KOA, confirming that these domains are structurally conserved. The structural models of aurora1 and aurora2 were built using 1CDK as the template structure. Molecular dynamics and docking simulations, targeting the ATP binding site of aurora2 with adenylyl imidodiphosphate (AMP-PNP), staurosporine, and six small molecular S/T kinase inhibitors, identified active-site residues that interact with these inhibitors differentially. The docked structures of the aurora2-AMP-PNP and aurora2-staurosporine complexes indicated that the adenine ring of AMP-PNP and the indolocarbazole moiety of staurosporine have similar positions and orientations and provided the basis for the docking of the other S/T kinase inhibitors. Inhibitors with isoquinoline and quinazoline moieties were recognized by aurora2 in which H-89 and 6,7-dimethoxyquinazoline compounds exhibited high binding energies compared with that of staurosporine. The calculated binding energies for the docked small-molecule inhibitors were qualitatively consistent with the IC(50) values generated using an in vitro kinase assay. The aurora2 structural model provides a rational basis for site-directed mutagenesis of the active site; design of novel H-89, staurosporine, and quinazoline analogues; and the screening of the available chemical database for the identification of other novel, small-molecular entities.

Amino Acid Sequence↗

PubMed2002

Molecular modeling of protein tyrosine phosphatase 1B (PTP 1B) inhibitors.

Binding modes of a series of aryloxymethylphosphonates and monoanionic biosteres of phosphate group from a series of benzylic alpha,alpha-diflluoro phosphate and its biosteres as protein tyrosine phosphatase 1B (PTP 1B) inhibitors have been identified by molecular modeling techniques. We have performed docking and molecular dynamics simulations of these inhibitors with PTP 1B enzyme. The initial conformation of the inhibitors for docking was obtained from simulated annealing technique. Solvent accessible surface area calculations suggested that active site of PTP 1B is highly hydrophobic. The results indicate that for aryloxymethylphosphonates, in addition to hydrogen bonding interactions, Tyr46, Arg47, Asp48, Val49, Glu115, Lys116, Lys120 amino acid residues of PTP 1B are responsible for governing inhibitor potency of the compounds. The sulfonate and tetrazole functional groups have been identified as effective monoanionic biosteres of phosphate group and biphenyl ring system due to its favorable interactions with Glu115, Lys116, Lys120 residues of PTP 1B found to be more suitable aromatic functionality than naphthalene ring system for benzylic alpha,alpha-diflluoro phosphate and its biosteres. The information generated from the present study should be useful in the design of more potent PTP 1B inhibitors as anti diabetic agents.

Binding Sites↗

PubMed2006

Subtype selectivity and flexibility of ionotropic glutamate receptors upon antagonist ligand binding.

The binding modes of a set of known ionotropic glutamate receptor antagonist-ligands have been studied using homology modeling, molecular docking, molecular dynamics (MD) simulations and ab initio quantum mechanical calculations. The core structure of the studied ligands is the decahydroisoquinoline ring, which has a carboxylic acid group at position three and different negatively-charged substituents (R) at position six. The binding affinities of these molecules have been reported earlier. From the current study, the carboxylate group of the decahydroisoquinoline ring hydrogen bonds with Arg485, the amino group with Pro478 and Thr480, and the negatively charged substituent R interacts with the positively charged N-terminus of helix-F. The subtype selectivity of these ligands seems to be strongly dependent on the amino acid at position 650 (GluR2: leucine, GluR5: valine), which affects the conformation of the ligand and ligand-receptor interactions, but depends considerably on the size of the R-group of the ligand. In addition, the MD simulations also revealed that the relative positions of the S1 and S2 domains can alter significantly showing different "closure" and "rotational movements" depending on the antagonist-ligand that is bound. Accordingly, molecular docking of antagonist ligands into static crystal structures cannot sufficiently explain ligand binding and subtype selectivity.

Amino Acid Sequence↗

PubMed2003

Putative caveolin-binding sites in SARS-CoV proteins.

AIM: To obtain the information of protein-protein interaction between the SARS-CoV proteins and caveolin-1, identify the possible caveolin-binding sites in SARS-CoV proteins. METHODS: On the basis of three related caveolin-binding motifs, amino acid motif search was employed to predict the possible caveolin-1 related interaction domains in the SARS-CoV proteins. The molecular modeling and docking simulation methods were used to confirm the interaction between caveolin-1 and SARS-CoV proteins. RESULTS: Thirty six caveolin-binding motifs in the SARS-CoV proteins have been mapped out using bioinformatics analysis tools. Molecular modeling and simulation have confirmed 8 caveolin-binding sites. These caveolin-binding sites located in replicase 1AB, spike protein, orf3 protein, and M protein, respectively. CONCLUSION: Caveolin-1 may serve as a possible receptor of the SARS-CoV proteins, which may be associated with the SARS-CoV infection, replication, assembly, and budding.

Amino Acid Motifs↗

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↗

PubMed2006

Possible dynamic anchor points in a benzoxazinone derivative-human oxytocin receptor system--a molecular docking and dynamics calculation.

In this study, we performed a molecular docking and dynamics simulation for a benzoxazinone-human oxytocin receptor system to determine the possible hydrophobic and electrostatic interaction points in the dynamic complex. After the homology modeling, the ligand was docked into the putative active using AutoDock 3.05. After the application of energetic and structural filters, the complexes obtained were further refined with a simulated annealing protocol (AMBER8) to remove steric clashes. Three complexes were selected for subjection to the molecular dynamics simulation (5 ns), and the results on the occurrence of average anchor points showed a stable complex between the benzoxazinone derivative and the receptor. The complex could be used as a good starting point for further analysis with site-directed mutagenesis, or further computational research.

Animals↗

PubMed2004

Homology modeling, force field design, and free energy simulation studies to optimize the activities of histone deacetylase inhibitors.

As an effort to develop therapeutics for cancer treatments, a number of effective histone deacetylase inhibitors with structural diversity have been discovered. To gain insight into optimizing the activity of an identified lead compound, a computational protocol sequentially involving homology modeling, docking experiments, molecular dynamics simulation, and free energy perturbation calculations was applied for rationalizing the relative activities of known histone deacetylase inhibitors. With the newly developed force field parameters for the coordination environment of the catalytic zinc ion in hand, the computational strategy proved to be successful in predicting the rank orders for 12 derivatives of three hydroxamate-based inhibitor scaffolds with indole amide, pyrrole, and sulfonamide moieties. The results showed that the free energy of an inhibitor in aqueous solution should be an important factor in determining the binding free energy. Hence, in order to enhance the inhibitory activity by adding or substituting a chemical group, the increased stabilization in solution due to the structural changes must be overcome by a stronger enzyme-inhibitor interaction. It was also found that to optimize inhibitor potency, the hydrophobic head of an inhibitor should be elongated or enlarged so that it can interact with Pro29 and His28 that are components of the flexible loop at the top of the active site.

Amino Acid Sequence↗