SEARCH · PubMed Health

Results for “molecular docking”

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 253 records · Page 14Linked to original sources

PubMed2026

Identification of potential key genes involved in iron deficiency for sepsis: A retrospective cohort and transcriptomic study.

Iron overload has been associated with sepsis, but the role of iron deficiency and its molecular links remain unclear. We investigated the association between iron deficiency and sepsis and identified candidate genes potentially linking these conditions. MIMIC-IV data were used to assess the association between serum iron and sepsis status. Transcriptomic datasets from dietary iron-deficient mice (GSE10421), LPS-induced septic mice (GSE267388), and a human blood sepsis cohort (GSE137340) were sequentially analyzed to identify and externally evaluate candidate genes. IEU Open GWAS summary statistics were used for exploratory Mendelian randomization (MR). Exploratory drug prediction was performed using L1000FWD, followed by molecular docking analysis. Patients with sepsis had significantly lower serum iron levels, and restricted cubic spline analysis showed a nonlinear association between serum iron and the odds of sepsis. Cross-tissue transcriptomic analysis identified Sqle, Lss, and Rdh11 as candidate genes. In the human blood cohort, SQLE and RDH11 were significantly increased, whereas LSS was not significantly altered. Exploratory MR showed that genetically proxied SQLE expression was associated with higher odds of sepsis (odds ratio [OR] = 1.23, P = 1.67 × 10-3), whereas LSS expression was associated with lower odds (OR = 0.97, P = 8.90 × 10-4); RDH11 showed no significant association (OR = 1.01, P = .90). Drug prediction identified ML106 as the top-ranked candidate drug, and molecular docking predicted potential binding poses with SQLE and LSS. Serum iron showed a nonlinear association with sepsis status. SQLE, LSS, and RDH11 emerged as candidate genes, with concordant expression changes of SQLE and RDH11 observed in human blood. MR findings for SQLE and LSS were exploratory and require further validation. ML106 was identified through exploratory drug prediction and requires experimental validation before its therapeutic relevance can be established.

Sepsis↗

PubMed2026

Whole-genome analysis of Brevibacterium sanguinis AZMABM HM27: a bacterial isolate from the sea anemone Radianthus magnifica and exhibiting promising multi-therapeutic properties.

BACKGROUND: The marine anemone Radianthus magnifica harbors symbiotic microbes with promising biomedical potential, yet their diversity and therapeutic properties remain underexplored. This study aimed to characterize a symbiotic bacterium isolated from R. magnifica collected from Samalona Island, Indonesia, and to evaluate its multi-therapeutic potential. METHODS: Strain AZMABM HM27 was characterized using whole-genome sequencing, functional annotation, biosynthetic gene cluster prediction, molecular docking, and in vitro bioactivity assays. RESULTS: Phylogenetic and genome-based analyses confirmed AZMABM HM27 as Brevibacterium sanguinis, with an OrthoANI value of 97.37% and a dDDH value of 76.50% against the type strain. The genome comprises a 3,834,082 bp chromosome encoding 3,362 protein-coding genes, including 95 genes involved in secondary metabolite biosynthesis. Five biosynthetic gene clusters were predicted, including those associated with ectoine, terpene, and siderophore production. The crude extract demonstrated antioxidant activity (IC₅₀ = 0.87 mg/mL), anti-inflammatory activity (up to 60% inhibition), antidiabetic activity through α-glucosidase inhibition (up to 40% inhibition), and dose-dependent antiproliferative activity against MCF-7 breast cancer cells (74.10% viability at 1 mg/mL). Molecular docking identified a lead compound, 8,9,9,10,10,11-hexafluoro-4,4-dimethyl-3,5-dioxatetracyclo [5.4.1.0(2,6)0.0(8,11)] dodecane, with strong binding affinities to selected therapeutic targets. CONCLUSIONS: B. sanguinis AZMABM HM27 represents a marine symbiotic strain associated with R. magnifica and a promising source of bioactive compounds with antioxidant, anti-inflammatory, antidiabetic, and antiproliferative potential. Further purification, structural elucidation, and in vivo studies are warranted to validate its therapeutic potential.

Animals↗

PubMed2026

Resveratrol Attenuates Gemcitabine Resistance in Hepatocellular Carcinoma Cells by Inhibiting Thymidylate Synthase.

BACKGROUND: Hepatocellular carcinoma (HCC) is a leading cause of cancer death worldwide. Gemcitabine (Gem) is a commonly used drug against HCC, but its efficacy is limited by the development of resistance. Resveratrol (Res), a natural polyphenol with antitumor activity, may reverse Gem resistance in HCC, although the mechanism remains unclear. METHODS: The effects of Res on the proliferation, apoptosis, cell cycle, and invasion of Hep3B and HuH-7 cells were assessed via cell counting kit-8 (CCK-8), clonogenic, flow cytometry, and Transwell assays, respectively. Potential Res targets were predicted by network pharmacology, and markers of HCC prognosis were identified from the cancer genome atlas (TCGA) data. The interaction between Res and thymidylate synthase (TYMS) was validated by molecular docking and dynamics simulation. A Gem-resistant HuH-7 cell line (HuH-7/GR) was established, and when these cells were treated with Res combined with Gem, the effect on Gem sensitivity was detected by CCK-8 assay, clonogenic assay, and flow cytometry. Finally, a subcutaneous nude mouse model of HCC was used to evaluate the in vivo effects of Res combined with Gem. RESULTS: Res inhibited HCC cell proliferation, induced apoptosis and G2/M arrest, and suppressed invasion in a concentration-dependent manner. Network pharmacology and TCGA analysis identified TYMS as an important target gene for Res. TYMS was highly expressed in HCC tissues and correlated with poor prognosis. Res treatment reduced TYMS expression, while molecular docking and simulation showed stable binding of Res to TYMS. TYMS levels were elevated in HuH-7/GR resistant cells. Res combined with Gem was found to reverse drug resistance, inhibit proliferation and colony formation, and induce apoptosis. The Res + Gem combination group showed the smallest tumor volume in the in vivo model. CONCLUSION: By attenuating Gem resistance through TYMS inhibition, Res holds promise as a clinically viable adjunct to Gem-based chemotherapy, offering a potential strategy to improve outcomes in HCC patients.

Resveratrol↗

PubMed2026

Active components and potential mechanisms of Wuzhuyu decoction in the treatment of ethanol-induced acute gastric mucosal injury: a network pharmacology and experimental verification.

OBJECTIVE: To investigate the underlying mechanisms and active components of Wuzhuyu decoction (, WD) in alleviating ethanol-induced acute gastric mucosal injury (GMI) using an integrated approach of network pharmacology and experimental verification. METHODS: Sprague-Dawley rats were randomly divided into six groups: control (Con), model (Mod), bismuth potassium citrate (BPC), WD at low (WD-L), medium (WD-M), and high (WD-H) doses. Following seven days of continuous intragastric administration of the respective treatments, an ethanol-induced gastric mucosal injury model was established in all groups except the control group by oral gavage of anhydrous ethanol. The gastric mucosal injury index was evaluated, and pathological changes were assessed viahematoxylin and eosin (HE) staining. Levels of tumor necrosis factor-alpha (TNF-α), interleukin-1 beta (IL-1β), malondialdehyde (MDA), superoxide dismutase (SOD), and glutathione peroxidase (GSH-Px) were measured by enzyme-linked immunosorbent assay (ELISA). The chemical composition was identified by ultra-performance liquid chromatography-tandem mass spectrometry. Active compounds were screened using the Swiss-absorption, distribution, metabolism, and excretion database, and their potential targets were predicted using the Swiss Target Prediction database and bioinformatics annotation database for molecular mechanism. Simultaneously, disease targets related to GMI were retrieved from the online mendelian inheritance in man and GeneCards databases. A protein-protein interaction (PPI) network was constructed, and functional enrichment analyses of gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were performed using the Metascape database. Key predictions from the network pharmacology analysis were subsequently verified through animal experiments. Protein expression levels of B-cell lymphoma-2 (Bcl-2), Bcl-2-associated X protein (Bax), Cleaved Caspase-3, and Cleaved Caspase-9 were analyzed by Western blot. Finally, molecular docking was performed using AutoDock Vina to investigate the interactions between the active components and core targets. RESULTS: WD treatment significantly reduced the gastric mucosal injury index and the levels of TNF-α, IL-1β, MDA, while it increased the activities of SOD and GSH-Px. Histopathological examination revealed marked improvement in gastric tissue morphology. A total of 145 compounds were identified in WD. Network pharmacology analysis identified 440 overlapping targets between WD and GMI. GO and KEGG enrichment analyses highlighted the apoptosis signaling pathway as a key mechanism for WD's protective effect against ethanol-induced GMI. Experimental validation demonstrated that WD treatment reduced the apoptosis of gastric mucosal epithelial cells, promoted the expression of Bcl-2, and inhibited the expression of Bax, Cleaved Caspase-3 and Cleaved Caspase-9. Molecular docking results indicated that dehydroevodiamine, rutaecarpine, evodiamine, hexahydrocurcumin, and isorhamnetin are potential active components in WD that contribute to the inhibition of apoptosis. CONCLUSIONS: WD alleviates ethanol-induced acute GMI, at least in part, by inhibiting the apoptosis. The primary active components responsible for this effect are dehydroevodiamine, rutaecarpine, evodiamine, hexahydrocurcumin, and isorhamnetin.

Drugs, Chinese Herbal↗

PubMed2025

Coordinated regulation of glutathione S-transferases confers metabolic flexibility in multi-insecticide-resistant Frankliniella occidentalis (Pergande).

INTRODUCTION: The evolution of multi-insecticide resistance in insect pests threatens global food security. Although glutathione S-transferases (GSTs) are implicated in detoxification, the coordinated mechanism by which specific gene subfamilies interact to confer broad-spectrum resistance remains inadequately characterized. OBJECTIVE: To dissect the functional allocation and cooperation of GST subfamilies in multi-insecticide-resistant strains of Frankliniella occidentalis. METHODS: We integrated comparative genomics (20 GST genes cloned), transcriptomics (qRT-PCR), RNAi-mediated silencing, molecular docking (AutoDock Vina), and in vitro metabolism assays (UPLC-MS/MS) across susceptible and resistant thrips strains. RESULTS: The two resistant strains (NIL-R and FS-R) exhibited moderate to high resistance to five insecticides (chlorfenapyr, emamectin benzoate, spinetoram, spinosad, and thiamethoxam), accompanied by significantly elevated GSTs activity. Phylogenetic analysis indicates that GSTs include 10 conserved delta and 7 diverse sigma members. The sigma subfamily has undergone a marked expansion due to gene duplication. Delta (FoGSTd1, d4, and d9) and sigma (FoGSTs1, s2, and s6) genes were significantly up-regulated in the resistant strains. RNAi showed specialized functional allocation among GSTs: delta GSTs mediated resistance to spinosad and chlorfenapyr, sigma GSTs were responsible for thiamethoxam resistance, and notably, cooperation between these subfamilies contributed to resistance against emamectin benzoate and spinetoram. Molecular docking and in vitro metabolism assays of FoGSTd9 and FoGSTs1 proteins further supported the functional allocation and cooperative roles of GST subfamilies. CONCLUSION: Our results indicate that F. occidentalis may coordinate GST subfamilies to achieve metabolic flexibility in response to multi-insecticide pressure. This survival strategy, mediated by mechanistic functional allocation and cooperative interactions among subfamilies, may contribute to energy conservation and reduced adaptive costs. Disruption of this coordinated mechanism represents a potential approach for overcoming resistance in agricultural pest populations.

Animals↗

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↗

PubMed2001

High throughput docking for library design and library prioritization.

The prioritization of the screening of combinatorial libraries is an extremely important task for the rapid identification of tight binding ligands and ultimately pharmaceutical compounds. When structural information for the target is available, molecular docking is an approach that can be used for prioritization. Here, we present the initial validation of a new rapid approach to molecular docking developed for prioritizing combinatorial libraries. The algorithm is tested on 103 individual cases from the protein data bank and in nearly 90% of these cases docks the ligand to within 2.0 A of the observed binding mode. Because the mean CPU time is <5 s/mol, this approach can process hundreds of thousands of compounds per week. Furthermore, if a somewhat less thorough search is performed, the search time drops to 1 s/mol, thus allowing millions of compounds to be docked per week and tested for potential activity. Proteins 2001;43:113-124.

Binding Sites↗

PubMed2000

Synthesis and evaluation of 4-(N,N-diarylamino)piperidines with high selectivity to the delta-opioid receptor: a combined 3D-QSAR and ligand docking study.

A series of 4-(N,N-diarylamino)piperidines are synthesized and evaluated for high affinity binding and selectivity to the delta-opioid receptor using a combination of 3D-QSAR and molecular docking techniques. Based on experimental ligand binding data to both mu- and delta- opioid receptors, CoMFA fields are generated and applied to identify potential ligand modifications to further optimize lead compounds. Molecular docking experiments to the delta-receptor are also reported that explain the CoMFA trends predicted as well as the differential binding and selectivity displayed by various compounds in the series. An analysis of the binding site model proposed indicates the piperidines take advantage of 3 key sites or binding domains within the delta-receptor. These include an aromatic pocket (approximately 1/3 into the receptor cavity), an aspartic acid residue (which serves as a docking point for the piperidinyl cationic amine) and a hydrophobic pocket at the extracellular boundary of the receptor cavity. Links are established between ligand modification and amino acid composition at these sites in mu and delta, providing new insight to the structural basis to binding and selectivity across the series and for related piperazines (i.e. SNC80 and BW373U86). Results are also presented that indicate delta- and mu-selectivity may be determined at alternate sites, suggesting opioid receptors may display multiple binding domains. The model is further supported by comparisons with opiate binding modes and site directed mutagenesis studies and is finally applied to suggest new strategies in ligand design.

Amino Acid Sequence↗

PubMed2006

Prediction of the binding mode between BMS-378806 and HIV-1 gp120 by docking and molecular dynamics simulation.

BMS-378806 is a newly discovered small molecule that effectively blocks the binding of CD4 with gp120. The binding mode of this kind of inhibitor remains unknown. In this paper, AutoDock 3.0 in conjunction with molecular dynamics simulation, accommodating the receptor's flexibility, was used to explore the binding mode between BMS-378806 and gp120. Two structures, Mode I and Mode II, with the lowest docking energy were selected as different representative binding modes. The analysis of the results from the molecular dynamics simulation indicated that the binding of BMS-348806 in Mode II is more stable. The average structure of Mode II was analyzed and compared with the experimental data. The conclusion was that BMS-378806 inserts the azaindole ring deeply into the PHE43 cavity and makes contact with a number of residues in the cavity, on the cavity and near the cavity. This study benefits the understanding of the mechanism of this kind of inhibitor and may provide useful information for rational drug design.

Binding Sites↗

PubMed2000

Automated docking and molecular dynamics simulations of nimesulide in the cyclooxygenase active site of human prostaglandin-endoperoxide synthase-2 (COX-2).

Molecular models of the complex between the selective COX-2 inhibitor nimesulide and the cyclooxygenase active site of human prostaglandin-endoperoxide synthase-2 have been built using a combination of homology modelling, conformational searching and automated docking techniques. The stability of the resulting complexes has been assessed by molecular dynamics simulations and interaction energy decomposition. It is found that nimesulide exploits the extra space made available by the replacement at position 523 of an isoleucine residue in COX-1 by a valine in COX-2 and establishes electrostatic interactions with both Arg-106 and Arg-499 (Arg-120 and Arg-513 in PGHS-1 numbering). Two alternate binding modes are proposed which are compatible with the pharmacological profile of this agent as a COX-2 selective inhibitor.

Amino Acid Sequence↗

PubMed2026

Computational discovery of emodin-based anthraquinones as PARP-1 inhibitors with relevance to ovarian and prostate cancer.

Cancer is a disease characterized by genomic instability and aberrant DNA repair. Poly (ADP-ribose) polymerase-1 (PARP-1) represents a well-established therapeutic target, particularly in ovarian and prostate cancer. However, the currently approved PARP inhibitors face challenges such as resistance, toxicity, and reduced efficacy. The search for alternative scaffolds has therefore become increasingly urgent. In this study, we used an integrated approach combining computer-aided methods to search for potential lead compounds among emodin-based anthraquinone derivatives as PARP-1 inhibitors. Using a PASS-based QSAR approach, drug-likeness prediction, and in silico ADMET assessment, we pre-screened a large set of anthraquinones and identified several potential hits for interaction with PARP-1. These hits were studied using molecular docking with the PARP-1 catalytic domain (PDB ID: 7KK4). The most stable and compact complexes were further explored by 500 ns molecular dynamics (MD) simulations and various dynamic properties (RMSD, RMSF, Rg, SASA, MolSA, hydrogen bonds, PCA, DCCM). The key finding of this study is that several emodin-derived anthraquinones exhibited binding behavior and ADMET profiles comparable to, or better than, the reference PARP-1 inhibitor. Among them, CID-10425624 emerged as the most promising candidate, exhibiting stable binding, reduced conformational fluctuation, compact complex formation, persistent hydrogen-bond interactions, and enhanced dynamic residue correlations within the PARP-1 catalytic domain. These findings suggest that the anthraquinone scaffold can provide a valuable starting point for developing structurally distinct PARP-1 inhibitors. In summary, this study identified several emodin-derived anthraquinones, particularly CID-10425624, as computationally prioritized lead candidates for PARP-1 inhibition, providing a novel anthraquinone-based scaffold for further experimental validation and optimization.

Anthraquinones↗

PubMed2001

A new molecular simulation software package--Peking University Drug Design System (PKUDDS) for structure-based drug design.

We present a comprehensive molecular simulation program package, the Peking University Drug Design System (PKUDDS), which runs on personal computers. PKUDDS has been developed mainly for computer-aided drug design using the methods of two-dimensional quantitative structure-activity relationships, three-dimensional quantitative structure-activity relationships, molecular docking, and database screening. This study presents an overview of its functionality, especially of methods developed in our group. PKUDDS uses genetic algorithms in molecular docking, conformational analysis, and quantitative structure-activity relationships as the most useful optimization technique. A user-friendly graphical interface provides easy access to many functions of PKUDDS. We report some examples of our considerable research using PKUDDS.

Antineoplastic Agents↗

PubMed2005

Virtual screening against metalloenzymes for inhibitors and substrates.

Molecular docking uses the three-dimensional structure of a receptor to screen databases of small molecules for potential ligands, often based on energetic complementarity. For many docking scoring functions, which calculate nonbonded interactions, metalloenzymes are challenging because of the partial covalent nature of metal-ligand interactions. To investigate how well molecular docking can identify potential ligands of metalloenzymes using a "standard" scoring function, we have docked the MDL Drug Data Report (MDDR), a functionally annotated database of 95,000 small molecules, against the X-ray crystal structures of five metalloenzymes. These enzymes included three zinc proteases, the nickel analogue of an iron enzyme, and a molybdenum metalloenzyme. The ability of the docking program to retrospectively enrich the annotated ligands as high-scoring hits for each enzyme and to calculate proper geometries was evaluated. In all five systems, the annotated ligands within the MDDR were enriched at least 20 times over random. To test the approach prospectively, a sixth target, the zinc beta-lactamase from Bacteroides fragilis, was screened against the fragment-like subset of the ZINC database. We purchased and tested 15 compounds from among the top 50 top-ranked ligands from docking, and found 5 inhibitors with apparent K(i) values less than 120 microM, the best of which was 2 microM. A more ambitious test still was predicting actual substrates for a seventh target, a Zn-dependent phosphotriesterase from Pseudomonas diminuta. Screening the Available Chemicals Directory (ACD) identified 25 thiophosphate esters as potential substrates within the top 100 ranked compounds. Eight of these, all previously uncharacterized for this enzyme, were acquired and tested, and all were confirmed experimentally as substrates. These results suggest that a simple, noncovalent scoring function may be used to identify inhibitors of at least some metalloenzymes.

Crystallography, X-Ray↗

PubMedDate not supplied

[Computer drug design based on analysis of a target macromolecule structure. I. Search and description of a ligand binding site in a target molecule].

The applicability of molecular docking method for finding of ligand binding site in a target protein was tested. The basic principle of tested approach consists in generation of hypotheses of protein-ligand complexes by molecular docking of small ligand to all surface of protein. The subsequent scoring of these hypotheses utilizes the values of contact surfaces and complex formation energy. The docking procedure was executed using the original software DockSearch running on PC Pentium or SGI computers. A set of hypotheses of probable ligand positions on the protein surface was created and evaluated with the help of DockSearch. Energy minimization of molecular complexes was done using Sybyl 6.5 (Tripos Inc.) running on SGI server Origin 200. The final set of the best hypotheses of complexes were selected by the values of contact surfaces and complex formation energy. The applicability and the limits of this approach was tested using known 3D structures of different proteins in free state and in complexes with ligands. For most target proteins the method allows to find the ligand binding site correctly. The accuracy of description of ligand binding site is adequate for subsequent searching of lead compounds by database mining. The reasons of some negative results obtained in testing of this approach are also discussed.

Binding Sites↗

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↗

PubMed2006

Molecular modeling, docking and ADMET studies applied to the design of a novel hybrid for treatment of Alzheimer's disease.

Alzheimer's disease (AD) is the most common form of dementia in adults, which is characterized by senile plaquets and cholinergic deficit as the disease progresses. Improvement of cholinergic neurotransmission is the basis of some drugs currently used in the treatment of AD. It is achieved by acetylcholinesterase (AChE) inhibition, the enzyme responsible for acetylcholine hydrolysis. Molecular modeling techniques were of utmost importance to design a new pharmaceutical against Alzheimer's disease, with potential inhibitory activity over AChE, since the inhibition of human plasma butyrylcholinesterase (BChE) may cause side effects. Some of the drugs currently used in the treatment of AD are capable of increasing the cholinergic transmission through the AChE inhibition. In this work we proposed molecular hybrids of tacrine with donepezil (fusion of these structures), in order to suggest new proposals of AChE inhibitors for future treatment of AD. We have analyzed all the structures by docking, density functional studies and drug like properties.

Acetylcholinesterase↗

PubMed1996

Scoring noncovalent protein-ligand interactions: a continuous differentiable function tuned to compute binding affinities.

Exploitation of protein structures for potential drug leads by molecular docking is critically dependent on methods for scoring putative protein-ligand interactions. An ideal function for scoring must exhibit predictive accuracy and high computational speed, and must be tolerant of variations in the relative protein-ligand molecular alignment and conformation. This paper describes the development of an empirically derived scoring function, based on the binding affinities of protein-ligand complexes coupled with their crystallographically determined structures. The function's primary terms involve hydrophobic and polar complementarity, with additional terms for entropic and solvation effects. The issue of alignment/conformation dependence was solved by constructing a continuous differentiable nonlinear function with the requirement that maxima in ligand conformation/alignment space corresponded closely to crystallographically determined structures. The expected error in the predicted affinity based on cross-validation was 1.0 log unit. The function is sufficiently fast and accurate to serve as the objective function of a molecular-docking search engine. The function is particularly well suited to the docking problem, since it has spatially narrow maxima that are broadly accessible via gradient descent.

Bacterial Proteins↗

PubMed1996

Computational methods for biomolecular docking.

With the rapidly increasing amount of molecular biological data available, the computer-based analysis of molecular interactions becomes more and more feasible. Methods for computer-aided molecular docking have to include a reasonably accurate model of energy and must be able to deal with the combinatorial complexity incurred by the molecular flexibility of the docking partners. In both respects, recent years have seen substantial progress.

Databases, Factual↗