Pub Date : 2026-01-13DOI: 10.1016/j.jmgm.2026.109280
Yingxuan Huang , W. Eltayeb Ahmed , Muhammad Farhan Hanif , Saba Hanif , Muhammad Imran , Muhammad Kamran Siddiqui
Quantitative structure property relationship(QSPR) has emerged as an indispensable tool for the estimation of physicochemical properties in drug molecules using mathematical and computational methods. Here, we introduce novel reverse degree based topological indices to see their applicability in case of selected antibiotic compounds property prediction. Reliable models to predict properties such as the boiling point, molar refractivity and enthalpy of vaporization exist to correlate molecular structure with experimentally reported physicochemical parameters. We have analyzed structurally different antibiotics with regression models developed in Python and SPSS in order to guarantee the robustness and reproducibility. We note here that predictive measures of cubic regression models seem to perform better, as observed through generally greater correlation coefficients. The results show that the reverse topological indices are efficient for recording structural differences in antibiotic molecules and they can be excellent descriptors for predicting their physical and chemical properties. It also stresses that, the use of reverse degree based descriptors on antibiotic compounds is new, providing a basis for further QSPR modeling for more general drug families. This work is part of a growing trend to study the interfaces between graph theory and cheminformatics where new indices help to improve our understanding over molecular properties with importance for drug design.
{"title":"On QSPR analysis for predicting efficacy of physicochemical properties of antibiotics drugs via topological indices and regression models","authors":"Yingxuan Huang , W. Eltayeb Ahmed , Muhammad Farhan Hanif , Saba Hanif , Muhammad Imran , Muhammad Kamran Siddiqui","doi":"10.1016/j.jmgm.2026.109280","DOIUrl":"10.1016/j.jmgm.2026.109280","url":null,"abstract":"<div><div>Quantitative structure property relationship(QSPR) has emerged as an indispensable tool for the estimation of physicochemical properties in drug molecules using mathematical and computational methods. Here, we introduce novel reverse degree based topological indices to see their applicability in case of selected antibiotic compounds property prediction. Reliable models to predict properties such as the boiling point, molar refractivity and enthalpy of vaporization exist to correlate molecular structure with experimentally reported physicochemical parameters. We have analyzed structurally different antibiotics with regression models developed in Python and SPSS in order to guarantee the robustness and reproducibility. We note here that predictive measures of cubic regression models seem to perform better, as observed through generally greater correlation coefficients. The results show that the reverse topological indices are efficient for recording structural differences in antibiotic molecules and they can be excellent descriptors for predicting their physical and chemical properties. It also stresses that, the use of reverse degree based descriptors on antibiotic compounds is new, providing a basis for further QSPR modeling for more general drug families. This work is part of a growing trend to study the interfaces between graph theory and cheminformatics where new indices help to improve our understanding over molecular properties with importance for drug design.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109280"},"PeriodicalIF":3.0,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980543","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-13DOI: 10.1016/j.jmgm.2026.109288
Pınar Siyah , Firat Baris Barlas
Cancer is the second leading cause of death globally and remains a priority due to its impact on life quality, treatment complexity, and high costs. To expedite drug development, researchers are increasingly repurposing FDA-approved drugs and clinical candidates, reducing time and costs through in silico methods. In this study, 3235 FDA-approved and clinical molecules were screened for EGFR inhibition, a significant target due to its role in cancer progression and treatment resistance. A pharmacophore model was generated based on erlotinib's co-crystallized structure and quantitative structure-activity relationships. Molecules meeting the pharmacophoric criteria underwent SP and XP docking, with thresholds of −6.00 kcal/mol and −7.00 kcal/mol, respectively, followed by anti-cancer potential analysis via MetaCore/MetaDrug and MD simulations at 1, 10, and 100 ns to assess EGFR-binding stability. For the molecule Ticagrelor, which demonstrated particularly promising results, and Erlotinib cell culture viability assays were conducted across three cell lines—cancerous A549, U87, and healthy BEAS-2B— (IC50) of, 8.2576 μM, 9.4058 μM, and 15.893 μM, respectively for Ticagrelor and 11.708 μM, 12.747 μM and 14.6709 μM, respectively for Erlotinib. In silico results highlight Ticagrelor's significant EGFR-inhibiting potential with enhanced binding stability compared to the reference.
{"title":"Repurposing drugs for EGFR-targeted cancer therapy: An in silico and in vitro study with pharmacophore-based insights","authors":"Pınar Siyah , Firat Baris Barlas","doi":"10.1016/j.jmgm.2026.109288","DOIUrl":"10.1016/j.jmgm.2026.109288","url":null,"abstract":"<div><div>Cancer is the second leading cause of death globally and remains a priority due to its impact on life quality, treatment complexity, and high costs. To expedite drug development, researchers are increasingly repurposing FDA-approved drugs and clinical candidates, reducing time and costs through in silico methods. In this study, 3235 FDA-approved and clinical molecules were screened for EGFR inhibition, a significant target due to its role in cancer progression and treatment resistance. A pharmacophore model was generated based on erlotinib's co-crystallized structure and quantitative structure-activity relationships. Molecules meeting the pharmacophoric criteria underwent SP and XP docking, with thresholds of −6.00 kcal/mol and −7.00 kcal/mol, respectively, followed by anti-cancer potential analysis via MetaCore/MetaDrug and MD simulations at 1, 10, and 100 ns to assess EGFR-binding stability. For the molecule Ticagrelor, which demonstrated particularly promising results, and Erlotinib cell culture viability assays were conducted across three cell lines—cancerous A549, U87, and healthy BEAS-2B— (IC50) of, 8.2576 μM, 9.4058 μM, and 15.893 μM, respectively for Ticagrelor and 11.708 μM, 12.747 μM and 14.6709 μM, respectively for Erlotinib. In silico results highlight Ticagrelor's significant EGFR-inhibiting potential with enhanced binding stability compared to the reference.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109288"},"PeriodicalIF":3.0,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146010770","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Alzheimer's disease is an unavoidable neurological disorder in which the death of brain cells brings on memory loss, cognitive decline, and eventual dementia. There is no recognized treatment for Alzheimer's illness. By the year 2050, it is expected that the global population will witness approximately 100 million cases of Alzheimer's disease (AD). Despite recognizing AD as a formidable illness for over a century, no effective cure has been discovered thus far. Synaptic dysfunction could result from disturbed synaptic calcium handling caused by excessive activation of glutamate receptors, particularly the N-methyl-D-aspartate receptors (NMDARs). Glutamate serves as the brain's primary excitatory neurotransmitter, acting on ionotropic and metabotropic glutamate receptors. In recent years, several pharmacologically active substances derived from plants, animals, and microbes have shown promise in treating AD by focusing on various pathogenic processes. Initially, we used virtual screening to assess natural product-like compounds against NMDA receptors. In this research study, we have screened a natural compound database derived from zinc15. The best candidate was then validated through molecular dynamics simulation (MDS). The results revealed that out of 4221 compounds tested, only 165 showed superior binding interactions compared to native ligands, making them inhibitors for protein. Further analysis using ADMET indicates favorable drug-like properties, particularly for CNS drug-likeness. The MDS results, including RMSD, RMSF, Rg, and residue interactions, indicated a strong and stable association between top molecules and target protein. This confirms that top molecules can effectively remain within the binding pockets of the target proteins, forming stable protein-ligand complexes.
阿尔茨海默病是一种不可避免的神经系统疾病,脑细胞死亡会导致记忆丧失、认知能力下降,最终导致痴呆。阿尔茨海默病目前还没有公认的治疗方法。到2050年,预计全球人口将见证约1亿阿尔茨海默病(AD)病例。尽管一个多世纪以来人们就认识到阿尔茨海默病是一种可怕的疾病,但迄今为止还没有发现有效的治疗方法。突触功能障碍可能是由于谷氨酸受体,特别是n -甲基- d -天冬氨酸受体(NMDARs)的过度激活引起的突触钙处理紊乱。谷氨酸是大脑的主要兴奋性神经递质,作用于嗜离子性和代谢性谷氨酸受体。近年来,从植物、动物和微生物中提取的一些药理活性物质通过关注各种致病过程,在治疗AD方面显示出了希望。最初,我们使用虚拟筛选来评估抗NMDA受体的天然产物样化合物。在本研究中,我们筛选了一个来源于zinc15的天然化合物数据库。然后通过分子动力学模拟(MDS)验证最佳候选物。结果显示,在测试的4221种化合物中,只有165种与天然配体相比表现出更好的结合相互作用,使其成为蛋白质的抑制剂。进一步的ADMET分析显示了良好的药物样性质,特别是对中枢神经系统药物相似。MDS结果包括RMSD、RMSF、Rg和残基相互作用,表明顶分子与靶蛋白之间存在强而稳定的关联。这证实了顶部分子可以有效地留在靶蛋白的结合口袋内,形成稳定的蛋白质-配体复合物。
{"title":"Investigating the anti-Alzheimer potential of biogenic compounds from Zinc15 database as NMDA antagonist: An in-silico approach","authors":"Somdatta Chaudhari , Asavari Shinde , Mukund Salunke , Shriram Bairagi , Azad Dhage , Pinkal Patel , Vivek Rathod , Sandeep Pathare , Nojood Altwaijry , Mohd Shahnawaz Khan","doi":"10.1016/j.jmgm.2026.109277","DOIUrl":"10.1016/j.jmgm.2026.109277","url":null,"abstract":"<div><div>Alzheimer's disease is an unavoidable neurological disorder in which the death of brain cells brings on memory loss, cognitive decline, and eventual dementia. There is no recognized treatment for Alzheimer's illness. By the year 2050, it is expected that the global population will witness approximately 100 million cases of Alzheimer's disease (AD). Despite recognizing AD as a formidable illness for over a century, no effective cure has been discovered thus far. Synaptic dysfunction could result from disturbed synaptic calcium handling caused by excessive activation of glutamate receptors, particularly the N-methyl-D-aspartate receptors (NMDARs). Glutamate serves as the brain's primary excitatory neurotransmitter, acting on ionotropic and metabotropic glutamate receptors. In recent years, several pharmacologically active substances derived from plants, animals, and microbes have shown promise in treating AD by focusing on various pathogenic processes. Initially, we used virtual screening to assess natural product-like compounds against NMDA receptors. In this research study, we have screened a natural compound database derived from zinc15. The best candidate was then validated through molecular dynamics simulation (MDS). The results revealed that out of 4221 compounds tested, only 165 showed superior binding interactions compared to native ligands, making them inhibitors for protein. Further analysis using ADMET indicates favorable drug-like properties, particularly for CNS drug-likeness. The MDS results, including RMSD, RMSF, Rg, and residue interactions, indicated a strong and stable association between top molecules and target protein. This confirms that top molecules can effectively remain within the binding pockets of the target proteins, forming stable protein-ligand complexes.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109277"},"PeriodicalIF":3.0,"publicationDate":"2026-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146157389","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-12DOI: 10.1016/j.jmgm.2026.109283
Mohamed Bakhit, Sina Karimzadeh, Tien-Chien Jen
This study investigates the potential of Mo2C MXene as a hydrogen storage material using density functional theory (DFT) and molecular dynamics (MD) simulations to examine its structural stability, electronic properties, and hydrogen adsorption behavior. The optimized Mo2C structure exhibits a hexagonal lattice with favorable adsorption sites over Mo atoms and shows a surface area expansion of approximately 4 % after hydrogen loading while maintaining lattice symmetry. Thermodynamic stability is confirmed through adsorption energy calculations, which reveal a clear relationship between energy levels and hydrogen concentration. The results indicate that H2 adsorption on Mo2C is a thermodynamically favorable and exothermic process, with adsorption energies ranging from −0.184 to −0.528 eV, satisfying the criteria for practical hydrogen storage applications. Charge transfer analysis identifies Mo atoms as electron acceptors. Density of States (DOS) calculations reveal a near-zero band gap, confirming the metallic nature of Mo2C, while Projected DOS (PDOS) and orbital maps show significant hybridization and electronic polarization among H, Mo, and C atoms. Charge density difference maps highlight effective charge redistribution with strong electric fields around Mo atoms. MD simulations further confirm the structural stability of the Mo2C–H2 system, showing minimal deformation during a 100 ps simulation and supporting efficient hydrogen adsorption. Overall, these findings establish Mo2C MXene as a promising candidate for hydrogen storage applications and provide valuable insights for experimental validation and further development of sustainable energy storage technologies.
{"title":"A first-principles study of hydrogen storage on MXene Mo2C monolayer","authors":"Mohamed Bakhit, Sina Karimzadeh, Tien-Chien Jen","doi":"10.1016/j.jmgm.2026.109283","DOIUrl":"10.1016/j.jmgm.2026.109283","url":null,"abstract":"<div><div>This study investigates the potential of Mo<sub>2</sub>C MXene as a hydrogen storage material using density functional theory (DFT) and molecular dynamics (MD) simulations to examine its structural stability, electronic properties, and hydrogen adsorption behavior. The optimized Mo<sub>2</sub>C structure exhibits a hexagonal lattice with favorable adsorption sites over Mo atoms and shows a surface area expansion of approximately 4 % after hydrogen loading while maintaining lattice symmetry. Thermodynamic stability is confirmed through adsorption energy calculations, which reveal a clear relationship between energy levels and hydrogen concentration. The results indicate that H<sub>2</sub> adsorption on Mo<sub>2</sub>C is a thermodynamically favorable and exothermic process, with adsorption energies ranging from −0.184 to −0.528 eV, satisfying the criteria for practical hydrogen storage applications. Charge transfer analysis identifies Mo atoms as electron acceptors. Density of States (DOS) calculations reveal a near-zero band gap, confirming the metallic nature of Mo<sub>2</sub>C, while Projected DOS (PDOS) and orbital maps show significant hybridization and electronic polarization among H, Mo, and C atoms. Charge density difference maps highlight effective charge redistribution with strong electric fields around Mo atoms. MD simulations further confirm the structural stability of the Mo<sub>2</sub>C–H<sub>2</sub> system, showing minimal deformation during a 100 ps simulation and supporting efficient hydrogen adsorption. Overall, these findings establish Mo<sub>2</sub>C MXene as a promising candidate for hydrogen storage applications and provide valuable insights for experimental validation and further development of sustainable energy storage technologies.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109283"},"PeriodicalIF":3.0,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980546","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-12DOI: 10.1016/j.jmgm.2026.109285
Sidra Manzoor , Nadeem Raza , Faheem Abbas
Organic solar cells (OSCs) offer lightweight, flexible, and cost−effective energy solutions. However, fullerene−based systems face limitations in stability, tunability, and absorption, prompting the exploration of non−fullerene alternatives to enhance efficiency and scalability. In this work, five newly designed molecules (2TT−1A to 2TT−5A) were systematically studied using Density Functional Theory (DFT) and Time−Dependent DFT (TD−DFT) at the B3LYP/6−311G (d, p) level in both gas and solvent (chloroform) phases. Key optoelectronic properties, including HOMO−LUMO gaps, absorption spectra, dipole moments, and excitation energies, were analyzed to evaluate their photovoltaic performance. All compounds demonstrated strong light−harvesting abilities, with a notable redshift in the absorption spectra observed in the solvent phase. Among them, 2TT−5A stood out with the narrowest energy gap (1.35 eV), the longest absorption wavelength (861 nm), the highest dipole moment (12.17 D), and the lowest excitation energy (1.43 eV), indicating efficient charge transfer and exciton dissociation. Open−circuit voltages (Voc) ranging from 0.54 to 1.38 V also suggest good photovoltaic potential. Additionally, the nonlinear optical (NLO) and organic light−emitting diodes (OLED) properties of 2TT−5A were explored, revealing significant hyperpolarizability and a favorable emission profile. These results suggested that 2TT−5A is an exceptional multifunctional candidate, encouraging experimental synthesis and validating this material's stability, potentially accelerating the development of multifunctional organic optoelectronic devices.
{"title":"Unveiling the potential of π-conjugated 2TT-R non-fullerene alternatives for multifunctional optoelectronic applications: A first-principles study","authors":"Sidra Manzoor , Nadeem Raza , Faheem Abbas","doi":"10.1016/j.jmgm.2026.109285","DOIUrl":"10.1016/j.jmgm.2026.109285","url":null,"abstract":"<div><div>Organic solar cells (OSCs) offer lightweight, flexible, and cost−effective energy solutions. However, fullerene−based systems face limitations in stability, tunability, and absorption, prompting the exploration of non−fullerene alternatives to enhance efficiency and scalability. In this work, five newly designed molecules <strong>(2</strong><strong>TT−1A to 2</strong><strong>TT−5A)</strong> were systematically studied using Density Functional Theory (DFT) and Time−Dependent DFT (TD−DFT) at the B3LYP/6−311G (d, p) level in both gas and solvent (chloroform) phases. Key optoelectronic properties, including HOMO−LUMO gaps, absorption spectra, dipole moments, and excitation energies, were analyzed to evaluate their photovoltaic performance. All compounds demonstrated strong light−harvesting abilities, with a notable redshift in the absorption spectra observed in the solvent phase. Among them, <strong>2</strong><strong>TT−5A</strong> stood out with the narrowest energy gap (1.35 eV), the longest absorption wavelength (861 nm), the highest dipole moment (12.17 D), and the lowest excitation energy (1.43 eV), indicating efficient charge transfer and exciton dissociation. Open−circuit voltages (<em>V</em><sub>oc</sub>) ranging from 0.54 to 1.38 V also suggest good photovoltaic potential. Additionally, the nonlinear optical (NLO) and organic light−emitting diodes (OLED) properties of <strong>2</strong><strong>TT−5A</strong> were explored, revealing significant hyperpolarizability and a favorable emission profile. These results suggested that <strong>2</strong><strong>TT−5A</strong> is an exceptional multifunctional candidate, encouraging experimental synthesis and validating this material's stability, potentially accelerating the development of multifunctional organic optoelectronic devices.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109285"},"PeriodicalIF":3.0,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146018703","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Human metapneumovirus (HMPV) remains a major respiratory pathogen without approved antivirals, highlighting the urgent need for novel therapeutics. This study implemented an integrative computational pipeline combining virtual screening, molecular docking, 2 μs molecular dynamics (MD) simulations, density functional theory (DFT), pharmacophore modeling, and ADMET profiling to identify potent HMPV inhibitors from Traditional Chinese Medicine. Among 180 screened phytoconstituents, glycyrrhizin (–9.3 kcal mol−1), hesperidin (–9.1 kcal mol−1), and saikosaponins (–9.0 kcal mol−1) exhibited strong binding affinities toward the HMPV matrix protein (PDB ID: 5WB0). Extended MD simulations confirmed complex stability with RMSD 0.17–0.22 nm, average of 3–5 persistent H-bonds, and DCCM correlation coefficient = 0.86 for glycyrrhizin. MM-PBSA binding free energies (ΔG_bind) of –46.2 ± 2.5, –44.7 ± 2.8, and –43.9 ± 2.2 kJ mol−1 for glycyrrhizin, hesperidin, and oseltamivir respectively, validated strong and stable interactions. DFT results indicated favorable electronic reactivity (HOMO–LUMO gap = 3.86 eV; electrophilicity = 2.74 eV), enhancing ligand-target complementarity. ADMET analysis predicted low systemic toxicity (LD50= 380–530 mg kg−1) but revealed moderate CYP3A4/CYP2C9 inhibition, suggesting the need for metabolic stability evaluation. Compared with reported fusion inhibitors such as EGCG and rutin, this matrix-targeted strategy introduces a distinct therapeutic mechanism. Overall, these findings establish a robust computational foundation for developing and experimentally validating potent natural inhibitors against HMPV.
{"title":"Exploring traditional Chinese medicine for antiviral drug discovery: A computational approach to combat human metapneumovirus (HMPV)","authors":"Amit Dubey , Manish Kumar , Aisha Tufail , Vivek Dhar Dwivedi","doi":"10.1016/j.jmgm.2026.109290","DOIUrl":"10.1016/j.jmgm.2026.109290","url":null,"abstract":"<div><div>Human metapneumovirus (HMPV) remains a major respiratory pathogen without approved antivirals, highlighting the urgent need for novel therapeutics. This study implemented an integrative computational pipeline combining virtual screening, molecular docking, 2 μs molecular dynamics (MD) simulations, density functional theory (DFT), pharmacophore modeling, and ADMET profiling to identify potent HMPV inhibitors from Traditional Chinese Medicine. Among 180 screened phytoconstituents, <strong>glycyrrhizin (–9.3 kcal mol<sup>−1</sup>)</strong>, <strong>hesperidin (–9.1 kcal mol<sup>−1</sup>)</strong>, and <strong>saikosaponins (–9.0 kcal mol<sup>−1</sup>)</strong> exhibited strong binding affinities toward the HMPV matrix protein (<strong>PDB ID: 5WB0</strong>). Extended MD simulations confirmed complex stability with <strong>RMSD 0.17</strong>–<strong>0.22 nm</strong>, average of <strong>3</strong>–<strong>5 persistent H-bonds</strong>, and <strong>DCCM correlation coefficient = 0.86</strong> for glycyrrhizin. <strong>MM-PBSA</strong> binding free energies (ΔG_bind) of <strong>–46.2 ± 2.5</strong>, <strong>–44.7 ± 2.8</strong>, and <strong>–43.9 ± 2.2 kJ mol<sup>−1</sup></strong> for glycyrrhizin, hesperidin, and oseltamivir respectively, validated strong and stable interactions. <strong>DFT</strong> results indicated favorable electronic reactivity (HOMO–LUMO gap = 3.86 eV; electrophilicity = 2.74 eV), enhancing ligand-target complementarity. <strong>ADMET</strong> analysis predicted low systemic toxicity (<strong>LD<sub>50</sub></strong> <strong>= 380</strong>–<strong>530 mg kg<sup>−1</sup></strong>) but revealed moderate <strong>CYP3A4/CYP2C9 inhibition</strong>, suggesting the need for metabolic stability evaluation. Compared with reported fusion inhibitors such as EGCG and rutin, this <strong>matrix-targeted strategy</strong> introduces a distinct therapeutic mechanism. Overall, these findings establish a robust computational foundation for developing and experimentally validating potent natural inhibitors against HMPV.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109290"},"PeriodicalIF":3.0,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980434","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This study employs DFT calculations to investigate the structural stability and electronic properties of pristine and transition-metal-doped boron nitride (BN) nanosheets, using yttrium (Y) and zirconium (Zr) as dopants, as well as their gas-sensing response toward formamide (FO). The findings show that introducing Y or Zr atoms leads to notable modifications in the electronic structure of the BN nanosheet, substantially improving its chemical reactivity and adsorption performance. In the aqueous phase, the interaction between FO and Y/Zr-doped BN nanosheets becomes moderately weaker, with adsorption energies decreasing to – 4.23 to – 24.97 kcal mol−1; however, the most stable complexes still exhibit comparatively strong binding. Solvation also alters the electronic structure of the nanosheets, leading to noticeable variations in their energy gaps. Despite this reduction in interaction strength, both doped materials retain high sensitivity toward FO in water, with ZrBN reaching 99.9 %/1.43 × 103 % and YBN achieving 55.9 %/86.5 %. Moreover, the nanosheets exhibit extremely short recovery times in the liquid phase, with values of 1.27 × 10−15 s for ZrBN and 2.06 s for YBN, enabling rapid FO desorption and efficient restoration of active metal sites. These combined features confirm the strong potential of Y- and Zr-doped BN nanosheets as reusable and high-performance sensors for formamide detection in aqueous environments.
{"title":"Y- and Zr-modified boron nitride nanosheets as efficient sensors for formamide: A first-principles approach","authors":"Meryem Derdare, Abdel-Ghani Boudjahem, Nedjoua Cheghib","doi":"10.1016/j.jmgm.2026.109278","DOIUrl":"10.1016/j.jmgm.2026.109278","url":null,"abstract":"<div><div>This study employs DFT calculations to investigate the structural stability and electronic properties of pristine and transition-metal-doped boron nitride (BN) nanosheets, using yttrium (Y) and zirconium (Zr) as dopants, as well as their gas-sensing response toward formamide (FO). The findings show that introducing Y or Zr atoms leads to notable modifications in the electronic structure of the BN nanosheet, substantially improving its chemical reactivity and adsorption performance. In the aqueous phase, the interaction between FO and Y/Zr-doped BN nanosheets becomes moderately weaker, with adsorption energies decreasing to – 4.23 to – 24.97 kcal mol<sup>−1</sup>; however, the most stable complexes still exhibit comparatively strong binding. Solvation also alters the electronic structure of the nanosheets, leading to noticeable variations in their energy gaps. Despite this reduction in interaction strength, both doped materials retain high sensitivity toward FO in water, with ZrBN reaching 99.9 %/1.43 × 10<sub>3</sub> % and YBN achieving 55.9 %/86.5 %. Moreover, the nanosheets exhibit extremely short recovery times in the liquid phase, with values of 1.27 × 10<sup>−15</sup> s for ZrBN and 2.06 s for YBN, enabling rapid FO desorption and efficient restoration of active metal sites. These combined features confirm the strong potential of Y- and Zr-doped BN nanosheets as reusable and high-performance sensors for formamide detection in aqueous environments.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109278"},"PeriodicalIF":3.0,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980542","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2026-01-12DOI: 10.1016/j.jmgm.2026.109286
Ali M. Alaseem , Glowi Alasiri , Mohamed M. El-Wekil , Al-Montaser Bellah H. Ali , Ahmed K. Hamdy
Cancer persists as a leading cause of global mortality, and the mitogen-activated protein kinase (MAPK) pathway plays a pivotal role in tumor progression and drug resistance. Among MAPK regulators, TAOK3 has emerged as a promising therapeutic target due to its oncogenic role in various cancers. Despite its significance, no clinically approved TAOK3 inhibitors exist. In this study we implemented a structure-based virtual screening approach to identify potential TAOK3 inhibitors from a library of 10,000 lead-like compounds. Molecular docking identified ten top-ranked candidates, with compound Z1 (ZINC ID: 77585305) demonstrating the strongest binding affinity (ΔG = −8.42 kcal/mol), outperforming reported inhibitors NCGC00188382 and SBI-581. ADMET profiling confirmed Z1's favorable drug-like properties, including high gastrointestinal absorption and minimal toxicity risks. Molecular dynamics simulations (100 ns) confirmed stable binding of Z1 to TAOK3, as indicated by low RMSD (<0.25 nm), consistent RMSF profiles, and compact radius of gyration. End-state free energy calculations using MM/GBSA also supported favorable binding, with Z1 showing excellent van der Waals interactions (−39.82 kcal/mol). Dynamic cross-correlation matrices and free energy landscape analysis further validated the stability of the TAOK3-Z1 complex. Collectively, these findings highlight Z1 as a promising TAOK3 inhibitor and a potential lead compound for further experimental validation in anticancer drug development.
{"title":"Structure-based discovery of novel TAOK3 inhibitor via virtual screening, molecular dynamics simulations, and MM/GBSA analysis","authors":"Ali M. Alaseem , Glowi Alasiri , Mohamed M. El-Wekil , Al-Montaser Bellah H. Ali , Ahmed K. Hamdy","doi":"10.1016/j.jmgm.2026.109286","DOIUrl":"10.1016/j.jmgm.2026.109286","url":null,"abstract":"<div><div>Cancer persists as a leading cause of global mortality, and the mitogen-activated protein kinase (MAPK) pathway plays a pivotal role in tumor progression and drug resistance. Among MAPK regulators, TAOK3 has emerged as a promising therapeutic target due to its oncogenic role in various cancers. Despite its significance, no clinically approved TAOK3 inhibitors exist. In this study we implemented a structure-based virtual screening approach to identify potential TAOK3 inhibitors from a library of 10,000 lead-like compounds. Molecular docking identified ten top-ranked candidates, with compound Z1 (ZINC ID: 77585305) demonstrating the strongest binding affinity (ΔG = −8.42 kcal/mol), outperforming reported inhibitors NCGC00188382 and SBI-581. ADMET profiling confirmed Z1's favorable drug-like properties, including high gastrointestinal absorption and minimal toxicity risks. Molecular dynamics simulations (100 ns) confirmed stable binding of Z1 to TAOK3, as indicated by low RMSD (<0.25 nm), consistent RMSF profiles, and compact radius of gyration. End-state free energy calculations using MM/GBSA also supported favorable binding, with Z1 showing excellent van der Waals interactions (−39.82 kcal/mol). Dynamic cross-correlation matrices and free energy landscape analysis further validated the stability of the TAOK3-Z1 complex. Collectively, these findings highlight Z1 as a promising TAOK3 inhibitor and a potential lead compound for further experimental validation in anticancer drug development.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109286"},"PeriodicalIF":3.0,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980544","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The prediction of the binding affinity of proteins and ligands in computational drug discovery with high accuracy is critical when evaluating the effectiveness of potential therapeutic compounds. This research work introduces RLBindDeep, a novel deep learning architecture based on the amalgamation of the ResNet and LSTM architectures, for improved accuracy in predicting protein–ligand binding affinities. Most traditional methodologies utilizing conventional molecular docking techniques suffer from poor accuracy owing to semi-flexible modeling approaches and limited considerations of complex interactions. On the other hand, RLBindDeep, which is formulated as a pose-independent binding affinity regression model that directly predicts experimental protein–ligand binding affinities from fixed complex structures, without performing docking or rescoring multiple poses, has performed well in extracting important features of the protein–ligand interaction. Specifically, the extracted features encompass ligand physicochemical descriptors (e.g., molecular weight, LogP, TPSA), protein-level features such as amino acid composition, and detailed interaction features including van der Waals, electrostatic, and hydrogen-bond energies. The model has been tested rigorously over the CASF-2016 benchmark dataset and has returned Pearson’s coefficient , Spearman’s coefficient , and Root Mean Square Error . This significantly outperforms existing state-of-the-art models, such as HAC-Net and AutoDock Vina. Improved accuracy and robustness in RLBindDeep further highlight the possibility of deep learning to revolutionize computational drug discovery processes, making strategies for drug development more efficient and targeted.
{"title":"RLBindDeep: A ResNet-LSTM based novel framework for protein–ligand binding affinity prediction","authors":"Ekarsi Lodh , Shalini Majumder , Tapan Chowdhury , Manashi De","doi":"10.1016/j.jmgm.2026.109282","DOIUrl":"10.1016/j.jmgm.2026.109282","url":null,"abstract":"<div><div>The prediction of the binding affinity of proteins and ligands in computational drug discovery with high accuracy is critical when evaluating the effectiveness of potential therapeutic compounds. This research work introduces RLBindDeep, a novel deep learning architecture based on the amalgamation of the ResNet and LSTM architectures, for improved accuracy in predicting protein–ligand binding affinities. Most traditional methodologies utilizing conventional molecular docking techniques suffer from poor accuracy owing to semi-flexible modeling approaches and limited considerations of complex interactions. On the other hand, RLBindDeep, which is formulated as a pose-independent binding affinity regression model that directly predicts experimental protein–ligand binding affinities from fixed complex structures, without performing docking or rescoring multiple poses, has performed well in extracting important features of the protein–ligand interaction. Specifically, the extracted features encompass ligand physicochemical descriptors (e.g., molecular weight, LogP, TPSA), protein-level features such as amino acid composition, and detailed interaction features including van der Waals, electrostatic, and hydrogen-bond energies. The model has been tested rigorously over the CASF-2016 benchmark dataset and has returned Pearson’s coefficient <span><math><mrow><mi>R</mi><mo>=</mo><mn>0</mn><mo>.</mo><mn>875</mn></mrow></math></span>, Spearman’s coefficient <span><math><mrow><mi>ρ</mi><mo>=</mo><mn>0</mn><mo>.</mo><mn>864</mn></mrow></math></span>, and Root Mean Square Error <span><math><mrow><mi>RMSE</mi><mo>=</mo><mn>0</mn><mo>.</mo><mn>993</mn></mrow></math></span>. This significantly outperforms existing state-of-the-art models, such as HAC-Net and AutoDock Vina. Improved accuracy and robustness in RLBindDeep further highlight the possibility of deep learning to revolutionize computational drug discovery processes, making strategies for drug development more efficient and targeted.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109282"},"PeriodicalIF":3.0,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980545","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Phosphatidylinositol-4,5-bisphosphate 3-kinase alpha (PI3Kα) is a central signaling enzyme driving cell proliferation and growth in cancers including breast cancer. Selective inhibition of PI3Kα isoform has become a promising therapeutic approach. In this work, 2000 in-house natural compounds were virtually screened against the ATP-binding site of PI3Kα. Of these, 618 compounds were predicted to have acceptable drug-likeness, pharmacokinetic, and toxicity properties based on in silico ADMET screening. Docking analysis highlighted four candidates forming stable hydrogen bonds with key residues V851, S854, and Q859 in the PI3Kα binding pocket. Molecular dynamics simulations were then used to assess their structural features and dynamic stability. Hit 2 was found to form strong hydrogen bonds with E849 and V851 of the PI3Kα protein. MM/GBSA-based binding free energy analysis supported that Hit 2 possessed the most favorable binding affinity to PI3Kα among the identified candidates. In vitro cytotoxicity assays were then performed in MCF-7 and MDA-MB-231 breast cancer cell lines, with alpelisib as a reference compound. Hit 2 reduced cell viability in both cell lines, but its effect was particularly pronounced in MDA-MB-231 cells, a model of triple-negative breast cancer (TNBC). These results suggest that Hit 2 represents a promising natural scaffold for further design and development in breast cancer therapy, with particular relevance for aggressive TNBC.
{"title":"Discovery of a novel PI3Kα inhibitor for breast cancer therapy via virtual screening method, molecular dynamics simulation and biological evaluation","authors":"Thitiya Boonma , Bodee Nutho , Phongthon Kanjanasirirat , Chananya Rajchakom , Nadtanet Nunthaboot","doi":"10.1016/j.jmgm.2026.109289","DOIUrl":"10.1016/j.jmgm.2026.109289","url":null,"abstract":"<div><div>Phosphatidylinositol-4,5-bisphosphate 3-kinase alpha (PI3Kα) is a central signaling enzyme driving cell proliferation and growth in cancers including breast cancer. Selective inhibition of PI3Kα isoform has become a promising therapeutic approach. In this work, 2000 in-house natural compounds were virtually screened against the ATP-binding site of PI3Kα. Of these, 618 compounds were predicted to have acceptable drug-likeness, pharmacokinetic, and toxicity properties based on <em>in silico</em> ADMET screening. Docking analysis highlighted four candidates forming stable hydrogen bonds with key residues V851, S854, and Q859 in the PI3Kα binding pocket. Molecular dynamics simulations were then used to assess their structural features and dynamic stability. Hit 2 was found to form strong hydrogen bonds with E849 and V851 of the PI3Kα protein. MM/GBSA-based binding free energy analysis supported that Hit 2 possessed the most favorable binding affinity to PI3Kα among the identified candidates. <em>In vitro</em> cytotoxicity assays were then performed in MCF-7 and MDA-MB-231 breast cancer cell lines, with alpelisib as a reference compound. Hit 2 reduced cell viability in both cell lines, but its effect was particularly pronounced in MDA-MB-231 cells, a model of triple-negative breast cancer (TNBC). These results suggest that Hit 2 represents a promising natural scaffold for further design and development in breast cancer therapy, with particular relevance for aggressive TNBC.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109289"},"PeriodicalIF":3.0,"publicationDate":"2026-01-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145980548","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}