Pub Date : 2026-01-21DOI: 10.1016/j.jmgm.2026.109307
Shuxin Song, Yusen Su, Qingyang Guo, Taigang Liu
Clathrin is a key structural protein in intracellular vesicle transport, mainly mediating clathrin-mediated endocytosis (CME) through trimeric assembly. Its functional abnormalities are closely associated with various diseases, including neurodegenerative disorders, tumor metastasis, and immune system dysregulation. Traditional experimental methods for identifying the presence of Clathrin have limitations such as high cost and time consumption. Therefore, it is particularly urgent to develop efficient and reliable computational methods to assist in Clathrin recognition. In this study, we propose a model named ClathPLM, which integrates sequence embeddings from three pre-trained protein language models (PPLMs), i.e., ProtT5, ProtBert, and ESM-3, and performs deep representation learning on each feature through an independent branch composed of a convolutional neural network (CNN) and a multi-head Attention (MHA) mechanism, finally fusing the representations of the three views to accomplish the classification task. To validate the effectiveness of this design, we further examined variants of the fusion strategy and attention mechanism. Evaluation results show that ClathPLM demonstrates excellent overall classification performance and robustness, surpassing current state-of-the-art methods. Moreover, the model performs strongly on an additional case-study dataset and shows good scalability on an extra vesicular transport proteins (VTPs) dataset. We anticipate that ClathPLM may contribute to a deeper understanding of the role of Clathrin in cellular regulation and disease mechanisms, and facilitate future biological studies as well as potential clinical applications.
{"title":"ClathPLM: Deep multi-view feature extraction with CNN and attention enhances clathrin protein identification","authors":"Shuxin Song, Yusen Su, Qingyang Guo, Taigang Liu","doi":"10.1016/j.jmgm.2026.109307","DOIUrl":"10.1016/j.jmgm.2026.109307","url":null,"abstract":"<div><div>Clathrin is a key structural protein in intracellular vesicle transport, mainly mediating clathrin-mediated endocytosis (CME) through trimeric assembly. Its functional abnormalities are closely associated with various diseases, including neurodegenerative disorders, tumor metastasis, and immune system dysregulation. Traditional experimental methods for identifying the presence of Clathrin have limitations such as high cost and time consumption. Therefore, it is particularly urgent to develop efficient and reliable computational methods to assist in Clathrin recognition. In this study, we propose a model named ClathPLM, which integrates sequence embeddings from three pre-trained protein language models (PPLMs), i.e., ProtT5, ProtBert, and ESM-3, and performs deep representation learning on each feature through an independent branch composed of a convolutional neural network (CNN) and a multi-head Attention (MHA) mechanism, finally fusing the representations of the three views to accomplish the classification task. To validate the effectiveness of this design, we further examined variants of the fusion strategy and attention mechanism. Evaluation results show that ClathPLM demonstrates excellent overall classification performance and robustness, surpassing current state-of-the-art methods. Moreover, the model performs strongly on an additional case-study dataset and shows good scalability on an extra vesicular transport proteins (VTPs) dataset. We anticipate that ClathPLM may contribute to a deeper understanding of the role of Clathrin in cellular regulation and disease mechanisms, and facilitate future biological studies as well as potential clinical applications.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109307"},"PeriodicalIF":3.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146046802","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-21DOI: 10.1016/j.jmgm.2026.109305
Samir Azizov , Vusala Nabi Jafarova , Khayala Ajdar Hasanova
An integrated experimental and first-principles investigation was carried out to elucidate the dielectric relaxation and optical properties of the chlorobenzene–n-butyl alcohol binary system containing 0.25 vol fraction of alcohol. The real (ε′) and imaginary (ε″) parts of the complex dielectric permittivity were measured over microwave wavelengths λ = 2.14–37.9 cm and radio frequencies of 0.047–15 MHz in a wide temperature range from −40 °C to 100 °C. The dielectric spectra in the liquid phase reveal two well-defined Debye-type relaxation regions, which are assigned to the independent reorientation of chlorobenzene and n-butyl alcohol molecules, as well as an asymmetric dispersion near the liquidus–solidus interval that is satisfactorily described by the Davidson–Cole model and attributed to hydrogen-bonded alcohol clusters. These features indicate pronounced microheterogeneity and cluster formation governed by hydrogen bonding. Complementary density functional theory (DFT) calculations were performed to provide atomistic insight into the polarization mechanisms and electronic structure of the chlorobenzene–n-butyl alcohol complex. Optical properties were systematically investigated using both the GGA-PBE approach and the more advanced LDA-RPA framework, allowing assessment of polarization screening and collective excitation effects. The calculations predict a wide band gap of about 4–4.5 eV and reveal strong optical anisotropy in the dielectric function, refractive index, optical conductivity, reflectivity, and absorption spectra, with the dominant response along the hydrogen-bonded molecular axis. The onset of intense π→π∗ and n→σ∗ electronic transitions above 4 eV in both GGA-PBE and LDA-RPA spectra is consistent with the experimentally observed dielectric relaxation and absorption behavior. The close agreement between experimental dielectric measurements and DFT-based optical responses confirms that hydrogen-bond-assisted dipole alignment governs both dielectric relaxation and optical polarization. As a result, the chlorobenzene–n-butyl alcohol system can be classified as a low-loss, wide-band-gap dielectric material, making it a promising candidate for microwave resonators, optoelectronic coatings, and polarization-sensitive sensing applications.
{"title":"Integrated experimental and DFT study of dielectric relaxation and optical properties in the chlorobenzene–n-butyl alcohol","authors":"Samir Azizov , Vusala Nabi Jafarova , Khayala Ajdar Hasanova","doi":"10.1016/j.jmgm.2026.109305","DOIUrl":"10.1016/j.jmgm.2026.109305","url":null,"abstract":"<div><div>An integrated experimental and first-principles investigation was carried out to elucidate the dielectric relaxation and optical properties of the chlorobenzene–n-butyl alcohol binary system containing 0.25 vol fraction of alcohol. The real (<em>ε</em>′) and imaginary (<em>ε</em>″) parts of the complex dielectric permittivity were measured over microwave wavelengths λ = 2.14–37.9 cm and radio frequencies of 0.047–15 MHz in a wide temperature range from −40 °C to 100 °C. The dielectric spectra in the liquid phase reveal two well-defined Debye-type relaxation regions, which are assigned to the independent reorientation of chlorobenzene and n-butyl alcohol molecules, as well as an asymmetric dispersion near the liquidus–solidus interval that is satisfactorily described by the Davidson–Cole model and attributed to hydrogen-bonded alcohol clusters. These features indicate pronounced microheterogeneity and cluster formation governed by hydrogen bonding. Complementary density functional theory (DFT) calculations were performed to provide atomistic insight into the polarization mechanisms and electronic structure of the chlorobenzene–n-butyl alcohol complex. Optical properties were systematically investigated using both the GGA-PBE approach and the more advanced LDA-RPA framework, allowing assessment of polarization screening and collective excitation effects. The calculations predict a wide band gap of about 4–4.5 eV and reveal strong optical anisotropy in the dielectric function, refractive index, optical conductivity, reflectivity, and absorption spectra, with the dominant response along the hydrogen-bonded molecular axis. The onset of intense π→π∗ and n→σ∗ electronic transitions above 4 eV in both GGA-PBE and LDA-RPA spectra is consistent with the experimentally observed dielectric relaxation and absorption behavior. The close agreement between experimental dielectric measurements and DFT-based optical responses confirms that hydrogen-bond-assisted dipole alignment governs both dielectric relaxation and optical polarization. As a result, the chlorobenzene–n-butyl alcohol system can be classified as a low-loss, wide-band-gap dielectric material, making it a promising candidate for microwave resonators, optoelectronic coatings, and polarization-sensitive sensing applications.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109305"},"PeriodicalIF":3.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146079022","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}
A theoretical study of the reaction forming a cyclic carbonate from the cycloaddition of CO2 with a typical substituted epoxide, was carried out using first principles methodologies. As epoxide, we considered trans 1-(4-methoxyphenyl)-2-methoxy epoxyethane, exhibiting an optimal reactivity for cycloaddition with CO2. We found that this reaction exhibits a high energy barrier in the gas phase, thus requiring an activation method. We also investigated the role of various substituted epoxides (mono- and disubstituted) and solvent types, including nonpolar, polar protic, and polar aprotic solvents. Our findings revealed a preference for polar protic solvents, particularly water. Besides, we used catalysts, including Lewis bases and Lewis acids, to selectively favor one reaction pathway over the other. In sum, this study provides deep insights into the factors influencing the reactivity of this cycloaddition and identifies the most suitable solvents and catalysts to promote the reaction, with potential implications for the development of more efficient CO2 valorization processes.
{"title":"Green chemistry strategies for CO2 valorization: Impact of substituents, solvents and catalysts on the cycloaddition of CO2 with epoxides forming cyclic carbonates","authors":"Imen Ferchichi , Frédéric Guégan , Muneerah Mogren Al-Mogren , Majdi Hochlaf , Youssef Arfaoui","doi":"10.1016/j.jmgm.2026.109306","DOIUrl":"10.1016/j.jmgm.2026.109306","url":null,"abstract":"<div><div>A theoretical study of the reaction forming a cyclic carbonate from the cycloaddition of CO<sub>2</sub> with a typical substituted epoxide, was carried out using first principles methodologies. As epoxide, we considered <em>trans</em> 1-(4-methoxyphenyl)-2-methoxy epoxyethane, exhibiting an optimal reactivity for cycloaddition with CO<sub>2</sub>. We found that this reaction exhibits a high energy barrier in the gas phase, thus requiring an activation method. We also investigated the role of various substituted epoxides (mono- and disubstituted) and solvent types, including nonpolar, polar protic, and polar aprotic solvents. Our findings revealed a preference for polar protic solvents, particularly water. Besides, we used catalysts, including Lewis bases and Lewis acids, to selectively favor one reaction pathway over the other. In sum, this study provides deep insights into the factors influencing the reactivity of this cycloaddition and identifies the most suitable solvents and catalysts to promote the reaction, with potential implications for the development of more efficient CO<sub>2</sub> valorization processes.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109306"},"PeriodicalIF":3.0,"publicationDate":"2026-01-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146046762","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-20DOI: 10.1016/j.jmgm.2026.109299
Mubarak A. Alamri , Mohamed Enneiymy , Yassine Riadi , Ali Altharawi , Taibah Aldakhil , Bharath Kumar Chagaleti , Ali Oubella , Reda A. Haggam
Gamma-hydroxybutyrate (GHB) is a psychoactive compound of high clinical and forensic concern due to its involvement in drug-facilitated intoxications and its rapid metabolic clearance, which complicates reliable detection. The reliance of conventional analytical techniques on centralized laboratories and costly instrumentation highlights the urgent need for fast, sensitive, and on-site sensing platforms. While carbon nanomaterial-based sensors have been widely investigated for narcotics detection, the sensing potential of coronene remains largely unknown. In this work, we present the first systematic DFT and TD-DFT investigation of pristine coronene (Cor) and its boron- and nitrogen-doped derivatives (B-Cor and N-Cor) for GHB detection, revealing a clear dopant-dependent sensing functionality. B-Cor emerges as the most effective GHB adsorbent, exhibiting a strong adsorption energy of −13.52 kcal mol−1, a large increase in dipole moment from 0.24 to 2.07 Debye, and enhanced polarizability (290.04–332.70 a.u.). These effects lead to an exceptional bathochromic shift in the UV–Vis spectrum (λmax from 373 to 594 nm, Δλ = 221 nm) and a decrease in exciton energy from 3.3 to 2.1 eV, establishing B-Cor as a highly promising single-use adsorptive and colorimetric sensor for naked-eye GHB detection. In contrast, N-Cor is identified as an optimal electrochemical sensor, characterized by a dramatically reduced HOMO-LUMO gap (0.57 eV), high chemical softness (1.75 eV-1), elevated electrical conductivity (2.75 × 109 A.m−2), and an ultra-fast recovery time of 5.63 × 10−6 s, enabling rapid and reusable sensing. Moderate increases in dipole moment (0.21–0.72 Debye) and polarizability (283.06–325.85 a.u.) further support its strong electronic responsiveness upon GHB binding. By distinctly identifying B-Cor as a superior adsorbent and colorimetric sensor and N-Cor as an efficient electrochemical sensor, this study introduces coronene as a versatile and tunable sensing scaffold and provides a robust theoretical foundation for the rational design of next-generation GHB sensing platforms for forensic and clinical applications.
{"title":"Boron and nitrogen-doped coronene as high-performance sensors for gamma-hydroxybutyrate drug sensing: A DFT/TD-DFT study","authors":"Mubarak A. Alamri , Mohamed Enneiymy , Yassine Riadi , Ali Altharawi , Taibah Aldakhil , Bharath Kumar Chagaleti , Ali Oubella , Reda A. Haggam","doi":"10.1016/j.jmgm.2026.109299","DOIUrl":"10.1016/j.jmgm.2026.109299","url":null,"abstract":"<div><div>Gamma-hydroxybutyrate (GHB) is a psychoactive compound of high clinical and forensic concern due to its involvement in drug-facilitated intoxications and its rapid metabolic clearance, which complicates reliable detection. The reliance of conventional analytical techniques on centralized laboratories and costly instrumentation highlights the urgent need for fast, sensitive, and on-site sensing platforms. While carbon nanomaterial-based sensors have been widely investigated for narcotics detection, the sensing potential of coronene remains largely unknown. In this work, we present the first systematic DFT and TD-DFT investigation of pristine coronene (Cor) and its boron- and nitrogen-doped derivatives (B-Cor and N-Cor) for GHB detection, revealing a clear dopant-dependent sensing functionality. B-Cor emerges as the most effective GHB adsorbent, exhibiting a strong adsorption energy of −13.52 kcal mol<sup>−1</sup>, a large increase in dipole moment from 0.24 to 2.07 Debye, and enhanced polarizability (290.04–332.70 a.u.). These effects lead to an exceptional bathochromic shift in the UV–Vis spectrum (<em>λ</em><sub>max</sub> from 373 to 594 nm, Δλ = 221 nm) and a decrease in exciton energy from 3.3 to 2.1 eV, establishing B-Cor as a highly promising single-use adsorptive and colorimetric sensor for naked-eye GHB detection. In contrast, N-Cor is identified as an optimal electrochemical sensor, characterized by a dramatically reduced HOMO-LUMO gap (0.57 eV), high chemical softness (1.75 eV<sup>-1</sup>), elevated electrical conductivity (2.75 × 10<sup>9</sup> A.m<sup>−2</sup>), and an ultra-fast recovery time of 5.63 × 10<sup>−6</sup> s, enabling rapid and reusable sensing. Moderate increases in dipole moment (0.21–0.72 Debye) and polarizability (283.06–325.85 a.u.) further support its strong electronic responsiveness upon GHB binding. By distinctly identifying B-Cor as a superior adsorbent and colorimetric sensor and N-Cor as an efficient electrochemical sensor, this study introduces coronene as a versatile and tunable sensing scaffold and provides a robust theoretical foundation for the rational design of next-generation GHB sensing platforms for forensic and clinical applications.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109299"},"PeriodicalIF":3.0,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146035985","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-20DOI: 10.1016/j.jmgm.2026.109304
Sadaf Noreen , Mamduh J. Aljaafreh , Ashour M. Ahmed
To understand the structural foundation of organic compounds is crucial in fields like chemistry and materials science. This study is a machine learning quest for predicting the experimental carbonyl peaks in the infrared (IR) spectrum of organic compounds from Modred and RDKit descriptors. The results show that FractionCSP3 is the most correlating descriptor for both types of descriptors. The Extra Trees (ET) regression yields the best performance with its coefficient of determination (R2) of 0.72–0.78. The analysis of models with SHapley Additive exPlanations (SHAP) shows that BCUT2D_MRLOW (RDKit) and FCSP3 (Modred) are the most influential descriptors. Hyperparameter tuning with 50 estimators optimizes model performance. Additionally, the calculated synthetic accessibility (SA) scores have 0.00–0.15 range to provide insights into the feasibility of synthesis. The current findings demonstrate the power of machine learning in uncovering the structural basis of organic compounds and predicting their experimental IR peaks.
{"title":"Exploring the structural basis of organic compounds by predicting experimental IR peaks: a machine learning analysis","authors":"Sadaf Noreen , Mamduh J. Aljaafreh , Ashour M. Ahmed","doi":"10.1016/j.jmgm.2026.109304","DOIUrl":"10.1016/j.jmgm.2026.109304","url":null,"abstract":"<div><div>To understand the structural foundation of organic compounds is crucial in fields like chemistry and materials science. This study is a machine learning quest for predicting the experimental carbonyl peaks in the infrared (IR) spectrum of organic compounds from Modred and <em>RDKit</em> descriptors. The results show that FractionCSP3 is the most correlating descriptor for both types of descriptors. The Extra Trees (<em>ET</em>) regression yields the best performance with its coefficient of determination (<em>R</em><sup><em>2</em></sup>) of 0.72–0.78. The analysis of models with SHapley Additive exPlanations (<em>SHAP</em>) shows that <em>BCUT2D_MRLOW</em> (<em>RDKit</em>) and <em>FCSP3</em> (Modred) are the most influential descriptors. Hyperparameter tuning with 50 estimators optimizes model performance. Additionally, the calculated synthetic accessibility (<em>SA</em>) scores have 0.00–0.15 range to provide insights into the feasibility of synthesis. The current findings demonstrate the power of machine learning in uncovering the structural basis of organic compounds and predicting their experimental <em>IR</em> peaks.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109304"},"PeriodicalIF":3.0,"publicationDate":"2026-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146029911","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-19DOI: 10.1016/j.jmgm.2026.109301
Bithia R., George Priya Doss C.
The extracellular domain (ECD) of the RET receptor tyrosine kinase depends on its cysteine-rich domain (CRD) for calcium coordination, structural stability, and assembly with GFRα1 and GDNF. Mutations close to the CRD CaLM motif have been associated with disease, but their molecular effects remain understudied. In this study, we analyzed two clinically reported variants, T564N and T564P, using all-atom molecular dynamics simulations of both the isolated CRD and the RET/GFRα1/GDNF ternary complex. Our analysis showed that both the mutations introduced localized structural changes in the CRD monomer. T564N caused increased residue fluctuations at the mutation site and solvent exposure, whereas T564P enhanced flexibility across all calcium-coordinating residues and slightly decreased stabilizing contacts. These effects became more noticeable in the ternary complex. Within the complex, interactions with the neighbouring domains caused the CRD to adopt conformations that compensated for the structural changes observed in the CRD monomer. In this context, each mutation affected calcium-binding energetics differently, resulting in more favourable binding in the mutants than in the wild-type. Although calcium binding was energetically favourable, the overall interaction energy within the complex was still affected. The complex highlighted mutation-specific differences in RET's interactions with GFRα1 and GDNF. The comparison between monomeric and complex simulations indicates that the functional impact of T564 mutations cannot be inferred from the isolated CRD. Together, these results show that the structural and energetic consequences of CRD CaLM mutations depend strongly on the full signaling assembly. This underscores the need to assess RET variants within their native multiprotein environment to understand how disease-associated mutations may alter receptor function.
{"title":"Dynamic consequences of threonine mutations in the CaLM motif of RET/GFRα1/GDNF ternary complex","authors":"Bithia R., George Priya Doss C.","doi":"10.1016/j.jmgm.2026.109301","DOIUrl":"10.1016/j.jmgm.2026.109301","url":null,"abstract":"<div><div>The extracellular domain (ECD) of the RET receptor tyrosine kinase depends on its cysteine-rich domain (CRD) for calcium coordination, structural stability, and assembly with GFRα1 and GDNF. Mutations close to the CRD CaLM motif have been associated with disease, but their molecular effects remain understudied. In this study, we analyzed two clinically reported variants, T564N and T564P, using all-atom molecular dynamics simulations of both the isolated CRD and the RET/GFRα1/GDNF ternary complex. Our analysis showed that both the mutations introduced localized structural changes in the CRD monomer. T564N caused increased residue fluctuations at the mutation site and solvent exposure, whereas T564P enhanced flexibility across all calcium-coordinating residues and slightly decreased stabilizing contacts. These effects became more noticeable in the ternary complex. Within the complex, interactions with the neighbouring domains caused the CRD to adopt conformations that compensated for the structural changes observed in the CRD monomer. In this context, each mutation affected calcium-binding energetics differently, resulting in more favourable binding in the mutants than in the wild-type. Although calcium binding was energetically favourable, the overall interaction energy within the complex was still affected. The complex highlighted mutation-specific differences in RET's interactions with GFRα1 and GDNF. The comparison between monomeric and complex simulations indicates that the functional impact of T564 mutations cannot be inferred from the isolated CRD. Together, these results show that the structural and energetic consequences of CRD CaLM mutations depend strongly on the full signaling assembly. This underscores the need to assess RET variants within their native multiprotein environment to understand how disease-associated mutations may alter receptor function.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109301"},"PeriodicalIF":3.0,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146079025","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-19DOI: 10.1016/j.jmgm.2026.109303
Areen Rasool, Jamshaid Ul Rahman, Qasim Ali
In bioinformatics, deep learning-based methods for Compound-Protein Interaction (CPI) prediction play a vital role in virtual screening, drug discovery, and drug repositioning. Recent improvements in computational methods have shown great possibility to save costs of experiment and speed up target identification. Nevertheless, the current CPI forecasting methods remain severely limited. Many rely on shallow Graph Neural Networks (GNNs) that struggle to capture the global structural context of compounds, while conventional Convolutional Neural Networks (CNNs) focus primarily on local sequence motifs and fail to model long-range dependencies in proteins. Even though a number of recent architectures strive to solve these problems by adding complexity to models, or by adding complex modules, these additions often cause significant computational overhead. To overcome these challenges, we propose DeepHybridCPI, a hybrid deep learning framework designed for accurate and efficient CPI prediction. Our hybrid model integrates a multiscale, densely connected GNN to extract compound features capturing both local substructures and global molecular topology, and employs CNNs with Long Short-Term Memory (LSTM) networks to model both local motifs and extended dependencies in protein sequences. The learned compound and protein representations are fused into a unified latent space to enable effective interaction modeling. Experimental evaluations on benchmark Human and C. elegans datasets demonstrate that DeepHybridCPI consistently outperforms existing state-of-the-art baseline methods in terms of AUC, Precision, and Recall. These findings highlight the importance of combining multiscale compound representations with hybrid sequence encoders within a single unified framework, providing a promising avenue for accelerating computational drug discovery. We release our source code and dataset at: https://github.com/jamshaidwarraich/DeepHybridCPI.
{"title":"DeepHybridCPI: A hybrid deep learning framework for compound–protein interaction prediction","authors":"Areen Rasool, Jamshaid Ul Rahman, Qasim Ali","doi":"10.1016/j.jmgm.2026.109303","DOIUrl":"10.1016/j.jmgm.2026.109303","url":null,"abstract":"<div><div>In bioinformatics, deep learning-based methods for Compound-Protein Interaction (CPI) prediction play a vital role in virtual screening, drug discovery, and drug repositioning. Recent improvements in computational methods have shown great possibility to save costs of experiment and speed up target identification. Nevertheless, the current CPI forecasting methods remain severely limited. Many rely on shallow Graph Neural Networks (GNNs) that struggle to capture the global structural context of compounds, while conventional Convolutional Neural Networks (CNNs) focus primarily on local sequence motifs and fail to model long-range dependencies in proteins. Even though a number of recent architectures strive to solve these problems by adding complexity to models, or by adding complex modules, these additions often cause significant computational overhead. To overcome these challenges, we propose DeepHybridCPI, a hybrid deep learning framework designed for accurate and efficient CPI prediction. Our hybrid model integrates a multiscale, densely connected GNN to extract compound features capturing both local substructures and global molecular topology, and employs CNNs with Long Short-Term Memory (LSTM) networks to model both local motifs and extended dependencies in protein sequences. The learned compound and protein representations are fused into a unified latent space to enable effective interaction modeling. Experimental evaluations on benchmark Human and <em>C. elegans</em> datasets demonstrate that DeepHybridCPI consistently outperforms existing state-of-the-art baseline methods in terms of AUC, Precision, and Recall. These findings highlight the importance of combining multiscale compound representations with hybrid sequence encoders within a single unified framework, providing a promising avenue for accelerating computational drug discovery. We release our source code and dataset at: <span><span>https://github.com/jamshaidwarraich/DeepHybridCPI</span><svg><path></path></svg></span>.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109303"},"PeriodicalIF":3.0,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146018798","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-19DOI: 10.1016/j.jmgm.2026.109302
Mingyu Wei , Zien Yu , Chi Huang , Shuwei Han , Sitian Zhou , Ting Wang , Shibin Chang , Xiaoying Chen , Feisheng Zhong , Shaojie Ma
Aquaporin-1 (AQP1),a key water channel protein, is aberrantly overexpressed in multiple malignancies, rendering it a compelling therapeutic target. The natural products Bacopaside I and Bacopaside II have demonstrated inhibitory activity against AQP1, yet their molecular mechanisms remain elusive. To elucidate the atomic basis of this inhibition, we employed a comprehensive computational approach combining molecular docking, molecular dynamics (MD) simulations, and extensive Gaussian accelerated molecular dynamics (GaMD) simulations with molecular mechanics generalized Born surface area (MM/GBSA) analysis. Our simulations indicate that both compounds exert spatial effects by occupying pore space, physically blocking channels, and forming van der Waals interactions with hydrophobic amino acids. In addition, ligands form hydrogen bonds with amino acids near these regions, resulting in narrower channels compared to other parts of AQP1. MM/GBSA calculations indicate that Bacopaside Ⅱ (ΔGbind = −34.48 kcal/mol) has a higher binding affinity than Bacopaside I (ΔGbind = −31.76 kcal/mol). Energy decomposition analysis identifies key interacting residues Pro171, Ile174, and Ala66 that anchor the inhibitors. Although both ligands induce subtle constriction of the channel pores, Bacopaside II establishes a more persistent hydrogen bonding network, underscoring its unique energetic contribution to the inhibition profile. Overall, these findings provide a detailed mechanistic blueprint for AQP1 inhibition by Bacopasides and offer a structural framework for the rational design of next-generation AQP1-targeted anticancer therapies.
{"title":"Molecular mechanisms of aquaporin 1 inhibition by Bacopaside I and Bacopaside II: Insights from molecular dynamics simulations","authors":"Mingyu Wei , Zien Yu , Chi Huang , Shuwei Han , Sitian Zhou , Ting Wang , Shibin Chang , Xiaoying Chen , Feisheng Zhong , Shaojie Ma","doi":"10.1016/j.jmgm.2026.109302","DOIUrl":"10.1016/j.jmgm.2026.109302","url":null,"abstract":"<div><div>Aquaporin-1 (AQP1),a key water channel protein, is aberrantly overexpressed in multiple malignancies, rendering it a compelling therapeutic target. The natural products Bacopaside I and Bacopaside II have demonstrated inhibitory activity against AQP1, yet their molecular mechanisms remain elusive. To elucidate the atomic basis of this inhibition, we employed a comprehensive computational approach combining molecular docking, molecular dynamics (MD) simulations, and extensive Gaussian accelerated molecular dynamics (GaMD) simulations with molecular mechanics generalized Born surface area (MM/GBSA) analysis. Our simulations indicate that both compounds exert spatial effects by occupying pore space, physically blocking channels, and forming van der Waals interactions with hydrophobic amino acids. In addition, ligands form hydrogen bonds with amino acids near these regions, resulting in narrower channels compared to other parts of AQP1. MM/GBSA calculations indicate that Bacopaside Ⅱ (Δ<em>G</em><sub>bind</sub> = −34.48 kcal/mol) has a higher binding affinity than Bacopaside I (Δ<em>G</em><sub>bind</sub> = −31.76 kcal/mol). Energy decomposition analysis identifies key interacting residues Pro171, Ile174, and Ala66 that anchor the inhibitors. Although both ligands induce subtle constriction of the channel pores, Bacopaside II establishes a more persistent hydrogen bonding network, underscoring its unique energetic contribution to the inhibition profile. Overall, these findings provide a detailed mechanistic blueprint for AQP1 inhibition by Bacopasides and offer a structural framework for the rational design of next-generation AQP1-targeted anticancer therapies.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109302"},"PeriodicalIF":3.0,"publicationDate":"2026-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146029923","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-17DOI: 10.1016/j.jmgm.2026.109298
Noha A. Saleh , Bushra M.D. Mansoor , A. Alkhaldi , Munirah A. Almessiere
Quercetin is a crucial natural bioactive substance present in many plants offering numerous health benefits. This computational study uses DFT calculations at the B3LYP/6-31G∗∗ level to computationally investigate novel quercetin derivatives modified by a tetra-amino acid sequence (Ser-Gly-Lys-Arg) in various structural forms including α-amino acids, β-amino acids, mixed-amino acids, and ketoacid. These modifications aim to improve the physicochemical and biological properties of quercetin providing a helpful framework for producing more potent bioactive compounds. Frontier molecular orbital analysis, molecular electrostatic potential mapping and reactivity indices reveal that β-amino acid and ketoacid derivatives exhibit enhanced electronic reactivity and thermodynamic stability relative to native quercetin. In addition, quantitative structure activity relationship (QSAR) descriptors and ADMET predictions are employed as comparative indicators to assess trends in polarity, solubility, and molecular permeability among the investigated compounds. The results highlight how peptide conjugation modulates the electronic and physicochemical properties of quercetin providing a rational computational basis for the future design of quercetin based peptidomimetic systems for further experimental investigation.
{"title":"Quantum chemical investigations of electronic properties, chemical reactivity, FTIR, thermodynamic behaviors and biological activity for some novel quercetin derivatives","authors":"Noha A. Saleh , Bushra M.D. Mansoor , A. Alkhaldi , Munirah A. Almessiere","doi":"10.1016/j.jmgm.2026.109298","DOIUrl":"10.1016/j.jmgm.2026.109298","url":null,"abstract":"<div><div>Quercetin is a crucial natural bioactive substance present in many plants offering numerous health benefits. This computational study uses DFT calculations at the B3LYP/6-31G∗∗ level to computationally investigate novel quercetin derivatives modified by a tetra-amino acid sequence (Ser-Gly-Lys-Arg) in various structural forms including α-amino acids, β-amino acids, mixed-amino acids, and ketoacid. These modifications aim to improve the physicochemical and biological properties of quercetin providing a helpful framework for producing more potent bioactive compounds. Frontier molecular orbital analysis, molecular electrostatic potential mapping and reactivity indices reveal that β-amino acid and ketoacid derivatives exhibit enhanced electronic reactivity and thermodynamic stability relative to native quercetin. In addition, quantitative structure activity relationship (QSAR) descriptors and ADMET predictions are employed as comparative indicators to assess trends in polarity, solubility, and molecular permeability among the investigated compounds. The results highlight how peptide conjugation modulates the electronic and physicochemical properties of quercetin providing a rational computational basis for the future design of quercetin based peptidomimetic systems for further experimental investigation.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109298"},"PeriodicalIF":3.0,"publicationDate":"2026-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146003762","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-16DOI: 10.1016/j.jmgm.2026.109297
Meiling Liu , Yuan Zhao , Linfan Shi , Zhongyang Ren , Wuyin Weng , Meitian Xiao
The effects of temperature (25–95 °C) and pH (7.0–9.0) on the film-forming mechanisms of soybean β-conglycinin (7S) and glycinin (11S) were investigated using molecular dynamics (MD) simulations and experimental validation. All MD simulations achieved equilibrium within 50 ns. The elevated temperatures and lower pH conditions reduced the center-of-mass distance and the radius-of-gyration (Rg) between 7S and 11S, while increasing hydrogen bond formation and binding free energy. Solvent-accessible surface area decreased with temperature, while root-mean-square fluctuation remained stable at pH 7.0 but increased with temperature at pH 9.0. The 7S-11S films prepared at higher temperatures exhibited enhanced tensile strength and higher proportion of hydrophobic interactions. With increasing temperature of the 7S-11S solution, the elongation at break increased at pH 7.0, but initially increased and then decreased at pH 9.0. Fourier transform infrared spectra revealed that hydrogen bonds and β-sheet structures increased with increasing temperature. In conclusion, heating the film-forming solution at pH 7.0 promoted 7S-11S molecular interactions, thereby improving the mechanical properties.
{"title":"Molecular dynamic simulation elucidates temperature- and pH-dependent film-forming mechanisms of soybean β-conglycinin and glycinin","authors":"Meiling Liu , Yuan Zhao , Linfan Shi , Zhongyang Ren , Wuyin Weng , Meitian Xiao","doi":"10.1016/j.jmgm.2026.109297","DOIUrl":"10.1016/j.jmgm.2026.109297","url":null,"abstract":"<div><div>The effects of temperature (25–95 °C) and pH (7.0–9.0) on the film-forming mechanisms of soybean <em>β</em>-conglycinin (7S) and glycinin (11S) were investigated using molecular dynamics (MD) simulations and experimental validation. All MD simulations achieved equilibrium within 50 ns. The elevated temperatures and lower pH conditions reduced the center-of-mass distance and the radius-of-gyration (Rg) between 7S and 11S, while increasing hydrogen bond formation and binding free energy. Solvent-accessible surface area decreased with temperature, while root-mean-square fluctuation remained stable at pH 7.0 but increased with temperature at pH 9.0. The 7S-11S films prepared at higher temperatures exhibited enhanced tensile strength and higher proportion of hydrophobic interactions. With increasing temperature of the 7S-11S solution, the elongation at break increased at pH 7.0, but initially increased and then decreased at pH 9.0. Fourier transform infrared spectra revealed that hydrogen bonds and <em>β</em>-sheet structures increased with increasing temperature. In conclusion, heating the film-forming solution at pH 7.0 promoted 7S-11S molecular interactions, thereby improving the mechanical properties.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109297"},"PeriodicalIF":3.0,"publicationDate":"2026-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146003781","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}