Nocardiosis, an infectious disease most commonly caused by the nosocomial pathogen Nocardia asteroides, predominantly affects immunosuppressed individuals and often presents as a pulmonary infection transmitted through inhalation. Although initial drugs like sulfadiazine and sulfisoxazole were used, their poor tissue penetration and the requirement for prolonged therapy have led to their exclusion from current treatment regimens. Consequently, there is a growing demand for novel therapeutics to suppress the pathogenic activity of N. asteroides. The pathogenic strain NCTC11293 of N. asteroides demonstrated a high completeness level of 98.19% in the CheckM analysis, with PGAP annotation identifying 6485 proteins. Given its clinical relevance and resistance profile, the NOCAP analysis revealed glycosyltransferase as a potential drug target associated with drug resistance and pathogenicity. The predicted and validated structure of the glycosyltransferase protein (accession WP_019050588.1) was subjected to structure-based virtual screening using the Schrödinger Suite against marine and phytochemical databases, including CMNPD, SWMD and IMPPAT. This multi-tier screening identified four promising hit molecules: CMNPD1518, CMNPD23345, SWMDGA001, and IMPY466077. The docking scores and MM-GBSA binding energies ranged from -7.20 to -8.90 kcal/mol and -49.81 to -59.90 kcal/mol, respectively. These hit molecules adhered to Lipinski's Rule of Five and exhibited accepted pharmacokinetic properties, along with balanced electronic structural characteristics, illustrating their higher reactivity and stability. Molecular dynamics (MD), essential dynamics (ED) analyses, including PCA-based FEL, and DSSP assessments reinforced the structural stability of the glycosyltransferase-phytochemical complexes. These computational findings suggest that the identified molecules represent promising therapeutic candidates for the treatment of nocardiosis.
{"title":"Structural and functional insights into glycosyltransferase from Nocardia asteroides NCTC11293 and structure-guided discovery of marine and phytochemical leads through pharmacokinetic screening and molecular dynamics studies.","authors":"Maharaja Muthuvairam Subbulakshmi, Manikandan Jayaraman, Prabhu Dhamodharan, Jeyakanthan Jeyaraman","doi":"10.1016/j.jmgm.2026.109317","DOIUrl":"https://doi.org/10.1016/j.jmgm.2026.109317","url":null,"abstract":"<p><p>Nocardiosis, an infectious disease most commonly caused by the nosocomial pathogen Nocardia asteroides, predominantly affects immunosuppressed individuals and often presents as a pulmonary infection transmitted through inhalation. Although initial drugs like sulfadiazine and sulfisoxazole were used, their poor tissue penetration and the requirement for prolonged therapy have led to their exclusion from current treatment regimens. Consequently, there is a growing demand for novel therapeutics to suppress the pathogenic activity of N. asteroides. The pathogenic strain NCTC11293 of N. asteroides demonstrated a high completeness level of 98.19% in the CheckM analysis, with PGAP annotation identifying 6485 proteins. Given its clinical relevance and resistance profile, the NOCAP analysis revealed glycosyltransferase as a potential drug target associated with drug resistance and pathogenicity. The predicted and validated structure of the glycosyltransferase protein (accession WP_019050588.1) was subjected to structure-based virtual screening using the Schrödinger Suite against marine and phytochemical databases, including CMNPD, SWMD and IMPPAT. This multi-tier screening identified four promising hit molecules: CMNPD1518, CMNPD23345, SWMDGA001, and IMPY466077. The docking scores and MM-GBSA binding energies ranged from -7.20 to -8.90 kcal/mol and -49.81 to -59.90 kcal/mol, respectively. These hit molecules adhered to Lipinski's Rule of Five and exhibited accepted pharmacokinetic properties, along with balanced electronic structural characteristics, illustrating their higher reactivity and stability. Molecular dynamics (MD), essential dynamics (ED) analyses, including PCA-based FEL, and DSSP assessments reinforced the structural stability of the glycosyltransferase-phytochemical complexes. These computational findings suggest that the identified molecules represent promising therapeutic candidates for the treatment of nocardiosis.</p>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"109317"},"PeriodicalIF":3.0,"publicationDate":"2026-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146149857","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}
Molecular dynamics (MD) simulations have carried out to investigate the solvation scenario of ions (Na+, K+, Cl-, I-, NH4+, CH3COO-) in the binary water-taurine and ternary water-taurine-urea mixtures at 308 K. The hydrogen bond (H-bond) properties and dynamics of water‒water, taurine‒water, and ammonium/acetate ion‒water are also explored. Our study revealed that monovalent cations have a higher tendency towards oxygen of taurine than water. NH4+ preferably interacts with the oxygen of the sulfonate group, whereas acetate ions are attributed towards the amino-hydrogens of taurine compared to the oxygen and hydrogens of water. It is found that taurine forms stronger H-bonds with water molecules through its amino group indicating the higher donating tendency of taurine. The translational and orientational relaxation dynamics of water molecules are slowed down in both the binary and ternary mixtures. The average lifetime of the water's H-bonds is higher in the ternary mixtures than the binary mixtures. Addition of urea enhances the overall lifetime of the ammonium ion-water (HAMM…OWAT) and acetate ion-water (OACE…HWAT) H-bonds. We have also observed that taurine counteracts the denaturation effect of urea on human islet amyloid polypeptide (hIAPP), as evidenced by the secondary structure. The α-helix content of hIAPP is restored with the addition of 0.45 M taurine in 5.60 M urea concentrations.
{"title":"Effects of taurine on the ion solvation scenario and its counteraction against the urea-induced denaturation of protein.","authors":"Chinmay Parida, Pragin Chettiyankandy, Deepak Kumar Mahanta, Snehasis Chowdhuri","doi":"10.1016/j.jmgm.2026.109315","DOIUrl":"https://doi.org/10.1016/j.jmgm.2026.109315","url":null,"abstract":"<p><p>Molecular dynamics (MD) simulations have carried out to investigate the solvation scenario of ions (Na<sup>+</sup>, K<sup>+</sup>, Cl<sup>-</sup>, I<sup>-</sup>, NH<sub>4</sub><sup>+</sup>, CH<sub>3</sub>COO<sup>-</sup>) in the binary water-taurine and ternary water-taurine-urea mixtures at 308 K. The hydrogen bond (H-bond) properties and dynamics of water‒water, taurine‒water, and ammonium/acetate ion‒water are also explored. Our study revealed that monovalent cations have a higher tendency towards oxygen of taurine than water. NH<sub>4</sub><sup>+</sup> preferably interacts with the oxygen of the sulfonate group, whereas acetate ions are attributed towards the amino-hydrogens of taurine compared to the oxygen and hydrogens of water. It is found that taurine forms stronger H-bonds with water molecules through its amino group indicating the higher donating tendency of taurine. The translational and orientational relaxation dynamics of water molecules are slowed down in both the binary and ternary mixtures. The average lifetime of the water's H-bonds is higher in the ternary mixtures than the binary mixtures. Addition of urea enhances the overall lifetime of the ammonium ion-water (H<sub>AMM</sub>…O<sub>WAT</sub>) and acetate ion-water (O<sub>ACE</sub>…H<sub>WAT</sub>) H-bonds. We have also observed that taurine counteracts the denaturation effect of urea on human islet amyloid polypeptide (hIAPP), as evidenced by the secondary structure. The α-helix content of hIAPP is restored with the addition of 0.45 M taurine in 5.60 M urea concentrations.</p>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"109315"},"PeriodicalIF":3.0,"publicationDate":"2026-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146142674","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-02-05DOI: 10.1016/j.jmgm.2026.109314
Sha Ding, Yong Xia
Spiropyran (SP) derivatives, as stimuli-responsive fluorescent materials, undergo reversible isomerization between closed-ring (SP) and open-ring (merocyanine, MC) forms upon light irradiation, holding promise for applications in optical switches and sensing. However, their excited-state dynamics and the roles of substituents/solvents in photoinduced ring-opening remain elusive. Herein, we use TD-CAM-B3LYP/6-31G(d) to clarify the photoinduced ring-opening mechanism of carboxyl-functionalized SP derivatives (SP-R, R = NO2, H, OCH3). Ground-state (S0) SP-R is stable with high C2-O1 cleavage barriers (15.5-31.4 kcal/mol), impeding spontaneous ring-opening; the electron-withdrawing NO2 reduces barriers, while the electron-donating OCH3 increases them. In the excited state (S1), these barriers are drastically lowered or even eliminated, facilitating SP→MC conversion, with "bright states" identified as S2 for SP-NO2 and S1 for SP-H/SP-OCH3. Structural changes promote ring-opening, with SP-OCH3 exhibiting the most pronounced effects. AIM analysis reveals weakened C2-O1 bond covalency in S1, particularly for SP-H and SP-OCH3. Solvent polarity modulates stability, dipole moments, and emission properties: polar solvents enhance trans-MC stability, and Trans-MC-OCH3 emits the longest wavelengths (602.2-639.3 nm) due to extended π-conjugation. This study clarifies SP→MC mechanisms, highlights the regulatory roles of substituents and solvents, and provides guidance for designing photoresponsive materials.
{"title":"A theoretical study on the mechanism of photoinduced ring-opening reaction of spiropyran derivatives.","authors":"Sha Ding, Yong Xia","doi":"10.1016/j.jmgm.2026.109314","DOIUrl":"https://doi.org/10.1016/j.jmgm.2026.109314","url":null,"abstract":"<p><p>Spiropyran (SP) derivatives, as stimuli-responsive fluorescent materials, undergo reversible isomerization between closed-ring (SP) and open-ring (merocyanine, MC) forms upon light irradiation, holding promise for applications in optical switches and sensing. However, their excited-state dynamics and the roles of substituents/solvents in photoinduced ring-opening remain elusive. Herein, we use TD-CAM-B3LYP/6-31G(d) to clarify the photoinduced ring-opening mechanism of carboxyl-functionalized SP derivatives (SP-R, R = NO<sub>2</sub>, H, OCH<sub>3</sub>). Ground-state (S<sub>0</sub>) SP-R is stable with high C2-O1 cleavage barriers (15.5-31.4 kcal/mol), impeding spontaneous ring-opening; the electron-withdrawing NO<sub>2</sub> reduces barriers, while the electron-donating OCH<sub>3</sub> increases them. In the excited state (S<sub>1</sub>), these barriers are drastically lowered or even eliminated, facilitating SP→MC conversion, with \"bright states\" identified as S<sub>2</sub> for SP-NO<sub>2</sub> and S<sub>1</sub> for SP-H/SP-OCH<sub>3</sub>. Structural changes promote ring-opening, with SP-OCH<sub>3</sub> exhibiting the most pronounced effects. AIM analysis reveals weakened C2-O1 bond covalency in S<sub>1</sub>, particularly for SP-H and SP-OCH<sub>3</sub>. Solvent polarity modulates stability, dipole moments, and emission properties: polar solvents enhance trans-MC stability, and Trans-MC-OCH<sub>3</sub> emits the longest wavelengths (602.2-639.3 nm) due to extended π-conjugation. This study clarifies SP→MC mechanisms, highlights the regulatory roles of substituents and solvents, and provides guidance for designing photoresponsive materials.</p>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"109314"},"PeriodicalIF":3.0,"publicationDate":"2026-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146142630","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-31DOI: 10.1016/j.jmgm.2026.109310
Nemat Ali, Ali M Alaseem, M Arockia Babu, Thakur Gurjeet Singh, Glowi Alasiri, Mohammad Fareed, Yogita Tyagi, Nisha Bansal
EGFR is a validated drug target in anticancer drug discovery. However, the existing issues of drug resistance, toxicities, and evolving mutations in EGFR have led to the global decline in the putative and potent anti-EGFR therapeutics. To address this gap, drug repurposing or repositioning is a valid technique for drug identification in a quick time and in an economical way. We herein employed the drug repurposing techniques from the library of 3500 approved USFDA small molecules. The initial HTVS (high-throughput virtual screening), followed by XP analysis and molecular mechanics (MMGBSA) and ADME profiling, led to the identification of Cycloguanil (antimalarial) and Metformin (antidiabetic) as the putative lead molecules with the potential to inhibit the EGFR. The top-scoring ligand was subjected to a molecular dynamics simulation study that stabilized the stability of cycloguanil within the ATP pocket of EGFR. The biological investigations further corroborated the in-silico studies. Cycloguanil inhibited the EGFR with a half maximal inhibitory concentration (IC50) of IC50 of 490 nM compared to erlotinib with an IC50 of 222 nM. Besides this, cycloguanil was able to halt the cell cycle progression at the G1 phase (46.54%), a peculiar feature of the kinase inhibitors that affect the CDKs and cyclins required for the passage of cancer cells through G1. Annexin V assay revealed that cycloguanil induced profound apoptosis in A549 cells. The lead molecules were also found to possess cytotoxicity profiles in MCF-7, A549, and HCT-116, which were reported to harbor the expression of EGFR. From the analysis, it was deduced that cycloguanil exhibited the most potent cytotoxicity towards the A549 cell with an IC50 of 6.83 μM, followed by HCT-116 with an IC50 of 9.32 μM, while in MCF-7, it exhibited an IC50 of 14.82 μM. The lead molecule, cycloguanil, may plausibly serve as an important template that may be optimized by performing bioisosteric replacements, leading to a putative kinase inhibitor with a potent anticancer profile.
{"title":"Drug repurposing strategy to identify the putative leads against the Epidermal growth factor receptor (EGFR) from the USFDA-approved drug pool: Investigating the utility as an anticancer agent.","authors":"Nemat Ali, Ali M Alaseem, M Arockia Babu, Thakur Gurjeet Singh, Glowi Alasiri, Mohammad Fareed, Yogita Tyagi, Nisha Bansal","doi":"10.1016/j.jmgm.2026.109310","DOIUrl":"https://doi.org/10.1016/j.jmgm.2026.109310","url":null,"abstract":"<p><p>EGFR is a validated drug target in anticancer drug discovery. However, the existing issues of drug resistance, toxicities, and evolving mutations in EGFR have led to the global decline in the putative and potent anti-EGFR therapeutics. To address this gap, drug repurposing or repositioning is a valid technique for drug identification in a quick time and in an economical way. We herein employed the drug repurposing techniques from the library of 3500 approved USFDA small molecules. The initial HTVS (high-throughput virtual screening), followed by XP analysis and molecular mechanics (MMGBSA) and ADME profiling, led to the identification of Cycloguanil (antimalarial) and Metformin (antidiabetic) as the putative lead molecules with the potential to inhibit the EGFR. The top-scoring ligand was subjected to a molecular dynamics simulation study that stabilized the stability of cycloguanil within the ATP pocket of EGFR. The biological investigations further corroborated the in-silico studies. Cycloguanil inhibited the EGFR with a half maximal inhibitory concentration (IC50) of IC<sub>50</sub> of 490 nM compared to erlotinib with an IC<sub>50</sub> of 222 nM. Besides this, cycloguanil was able to halt the cell cycle progression at the G1 phase (46.54%), a peculiar feature of the kinase inhibitors that affect the CDKs and cyclins required for the passage of cancer cells through G1. Annexin V assay revealed that cycloguanil induced profound apoptosis in A549 cells. The lead molecules were also found to possess cytotoxicity profiles in MCF-7, A549, and HCT-116, which were reported to harbor the expression of EGFR. From the analysis, it was deduced that cycloguanil exhibited the most potent cytotoxicity towards the A549 cell with an IC50 of 6.83 μM, followed by HCT-116 with an IC50 of 9.32 μM, while in MCF-7, it exhibited an IC<sub>50</sub> of 14.82 μM. The lead molecule, cycloguanil, may plausibly serve as an important template that may be optimized by performing bioisosteric replacements, leading to a putative kinase inhibitor with a potent anticancer profile.</p>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"109310"},"PeriodicalIF":3.0,"publicationDate":"2026-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146125163","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-31DOI: 10.1016/j.jmgm.2026.109312
Juhaer Anjum, Sabiha Enam Spriha, Mir Sabit Hossain, S M Abdur Rahman
Tinospora crispa, a plant with widespread use in ethnomedicine, has several extracts as well as phytoconstituents that have been reported to possess anticancer properties. Clerodane-type furano-diterpenoids from this genus of plants show promise as anticancer agents, although only a few of these have been assayed in vitro. 38 clerodane-type furano-diterpenoids from T. crispa were assayed computationally against 13 common anticancer macromolecular targets. The assayed ligand library displayed higher binding affinities than the respective controls in six of the thirteen anticancer macromolecular targets. Binding affinity values that were higher than the respective controls were obtained for two ligands against BRCA1, for four against CDK9/Cyclin T1, for two against C-ROS1, for thirty-seven against EGFR, for three against MEK 1-kinase, and for one against STAT3. Many of these ligands displayed greater binding affinity values than the respective controls against multiple macromolecules. Notable among these are tinocrispide (against BRCA1, C-ROS1, EGFR and MEK1), tinosporol C (against CDK9/Cyclin T1 and MEK1-kinase), crispene E (against BRCA1 and CDK9/Cyclin T1) and (5R,6R,8S,9R,10R,12S)-15,16-Epoxy-2-oxo-6-O-(β-D-glucopyranosyl)-cleroda-3,13(16),14-trien-17,12-olid-18-oic acid methyl ester (against BRCA1 and EGFR). The top scoring compounds of each macromolecule were subjected to molecular dynamics simulations, and the obtained protein-ligand complexes were observed to be stable for all macromolecules except EGFR. The ligand library was mostly predicted to be soluble, but GI absorption was predicted to vary. A few were predicted to be carcinogenic, immunotoxic, and/or cytotoxic. Overall, the clerodane-type furano-diterpenoids from T. crispa show encouraging results in in-silico studies and may be considered for further modification and lead development.
{"title":"Computational assessment of clerodane-type furano-diterpenoids from Tinospora crispa: A potential source for anticancer lead compounds.","authors":"Juhaer Anjum, Sabiha Enam Spriha, Mir Sabit Hossain, S M Abdur Rahman","doi":"10.1016/j.jmgm.2026.109312","DOIUrl":"https://doi.org/10.1016/j.jmgm.2026.109312","url":null,"abstract":"<p><p>Tinospora crispa, a plant with widespread use in ethnomedicine, has several extracts as well as phytoconstituents that have been reported to possess anticancer properties. Clerodane-type furano-diterpenoids from this genus of plants show promise as anticancer agents, although only a few of these have been assayed in vitro. 38 clerodane-type furano-diterpenoids from T. crispa were assayed computationally against 13 common anticancer macromolecular targets. The assayed ligand library displayed higher binding affinities than the respective controls in six of the thirteen anticancer macromolecular targets. Binding affinity values that were higher than the respective controls were obtained for two ligands against BRCA1, for four against CDK9/Cyclin T1, for two against C-ROS1, for thirty-seven against EGFR, for three against MEK 1-kinase, and for one against STAT3. Many of these ligands displayed greater binding affinity values than the respective controls against multiple macromolecules. Notable among these are tinocrispide (against BRCA1, C-ROS1, EGFR and MEK1), tinosporol C (against CDK9/Cyclin T1 and MEK1-kinase), crispene E (against BRCA1 and CDK9/Cyclin T1) and (5R,6R,8S,9R,10R,12S)-15,16-Epoxy-2-oxo-6-O-(β-D-glucopyranosyl)-cleroda-3,13(16),14-trien-17,12-olid-18-oic acid methyl ester (against BRCA1 and EGFR). The top scoring compounds of each macromolecule were subjected to molecular dynamics simulations, and the obtained protein-ligand complexes were observed to be stable for all macromolecules except EGFR. The ligand library was mostly predicted to be soluble, but GI absorption was predicted to vary. A few were predicted to be carcinogenic, immunotoxic, and/or cytotoxic. Overall, the clerodane-type furano-diterpenoids from T. crispa show encouraging results in in-silico studies and may be considered for further modification and lead development.</p>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"109312"},"PeriodicalIF":3.0,"publicationDate":"2026-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146119257","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-27DOI: 10.1016/j.jmgm.2026.109308
Mohammad Asad , Mohamed Usman , Anisha Ashokan, Naveen Kulkarni
Angiotensin-converting enzyme (ACE) inhibitors are widely prescribed for cardiovascular disorders, yet their conventional chemical synthesis involves multiple steps, high energy consumption, and poor stereoselectivity. In this work, we present a fully insilico framework for the computational design and prediction of transaminase (TA) variants capable of catalyzing the asymmetric amination of ethyl 2-oxo-4-phenylbutanoate (OPBE) to yield L-homophenylalanine ethyl ester (L-HPE), a key chiral intermediate in ACE inhibitor synthesis. A homology model of the Silicibacter transaminase (62.7 % identity to 5KR6) was constructed in its dimeric form. Eleven active-site variants were designed and screened through molecular docking, followed by 100 ns molecular dynamics simulations. The top variants SbTA10, SbTA01, and SbTA11 exhibited reactive distances below 6 Å, binding energies between −17.7 and −20.8 kcal/mol, and substrate RMSD values under 2.0 Å, indicating stable enzyme–substrate complexes. A composite QZ-score integrating productive conformations, binding energy, pocket contacts, and structural stability ranked SbTA10 highest (0.89), followed by SbTA01 (0.88) and SbTA11 (0.85). Free-energy profiles derived from umbrella sampling revealed binding minima of −10.2 kcal/mol (SbTA10) and −8.1 kcal/mol (SbTA11), suggesting distinct substrate retention characteristics that may influence catalytic turnover. Collectively, these results identify plausible transaminase variants with favorable structural and energetic features for the proposed OPBE to L-HPE transformation. This study presents a computational framework for predicting transaminase variants, providing a basis for rational biocatalyst design that warrants future experimental validation to confirm catalytic efficiency and stereoselectivity.
{"title":"In silico engineering of transaminase variants for enhanced biocatalytic conversion of an ACE inhibitor precursor","authors":"Mohammad Asad , Mohamed Usman , Anisha Ashokan, Naveen Kulkarni","doi":"10.1016/j.jmgm.2026.109308","DOIUrl":"10.1016/j.jmgm.2026.109308","url":null,"abstract":"<div><div>Angiotensin-converting enzyme (ACE) inhibitors are widely prescribed for cardiovascular disorders, yet their conventional chemical synthesis involves multiple steps, high energy consumption, and poor stereoselectivity. In this work, we present a fully <em>in</em> <em>silico</em> framework for the computational design and prediction of transaminase (TA) variants capable of catalyzing the asymmetric amination of ethyl 2-oxo-4-phenylbutanoate (OPBE) to yield L-homophenylalanine ethyl ester (L-HPE), a key chiral intermediate in ACE inhibitor synthesis. A homology model of the <em>Silicibacter</em> transaminase (62.7 % identity to 5KR6) was constructed in its dimeric form. Eleven active-site variants were designed and screened through molecular docking, followed by 100 ns molecular dynamics simulations. The top variants SbTA10, SbTA01, and SbTA11 exhibited reactive distances below 6 Å, binding energies between −17.7 and −20.8 kcal/mol, and substrate RMSD values under 2.0 Å, indicating stable enzyme–substrate complexes. A composite QZ-score integrating productive conformations, binding energy, pocket contacts, and structural stability ranked SbTA10 highest (0.89), followed by SbTA01 (0.88) and SbTA11 (0.85). Free-energy profiles derived from umbrella sampling revealed binding minima of −10.2 kcal/mol (SbTA10) and −8.1 kcal/mol (SbTA11), suggesting distinct substrate retention characteristics that may influence catalytic turnover. Collectively, these results identify plausible transaminase variants with favorable structural and energetic features for the proposed OPBE to L-HPE transformation. This study presents a computational framework for predicting transaminase variants, providing a basis for rational biocatalyst design that warrants future experimental validation to confirm catalytic efficiency and stereoselectivity.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109308"},"PeriodicalIF":3.0,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146079023","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-27DOI: 10.1016/j.jmgm.2026.109309
Uttara Bakshi, Ashis Datta
Epilepsy remains a debilitating neurological disorder affecting approximately 50 million individuals worldwide. Current Anti-seizure medications (ASM) targeting Synaptic Vesicle Protein 2A (SV2A), such as levetiracetam, exhibit variable efficacy and adverse effect profiles, necessitating exploration of novel therapeutic approaches. This study investigates Bacopa monnieri-derived compounds as potential dual inhibitors of SV2A and Carbonic Anhydrase II (CA-II), two mechanistically distinct targets implicated in seizure pathophysiology. An integrated computational workflow comprising ADMET profiling, molecular docking, molecular dynamics (MD) simulations, and MM-GBSA binding free energy calculations was implemented to evaluate selected phytoconstituents. ADMET analysis revealed favourable pharmacokinetic parameters for several compounds, particularly Ebelin Lactone and Jujubogenin, with superior blood-brain barrier permeability (BBB >0.9) and human intestinal absorption (>95 %) compared to reference drugs. Molecular docking identified compelling binding affinities, with Ebelin Lactone (−11.2 kcal/mol) and Jujubogenin (−10.9 kcal/mol) exhibiting stronger interactions with SV2A than brivaracetam (−6.8 kcal/mol). Similarly, these compounds demonstrated robust binding to CA-II (−8.4 and −8.7 kcal/mol, respectively). Protein-ligand interaction profiler analysis elucidated key stabilizing interactions, including hydrogen bonds with Thr199 and His64 in CA-II and hydrophobic contacts within the SV2A binding pocket. MD simulations confirmed structural stability of the protein-ligand complexes, evidenced by equilibrated RMSD trajectories (0.2–0.3 nm). MM-GBSA calculations substantiated the thermodynamic favourability of these interactions, with binding free energies for top candidates significantly exceeding reference compounds. These findings establish an in silico framework for prioritizing Bacopa monnieri phytoconstituents as computationally predicted dual-target candidates, highlighting their potential relevance for seizure management, with the possibility of improved efficacy through simultaneous modulation of distinct pathophysiological mechanisms. All findings are based solely on in silico analyses and require experimental validation.
{"title":"Structure-based computational assessment of Bacopa monnieri-derived compounds as potential dual target anti-seizure medications: An integrated docking and molecular dynamics simulation approach","authors":"Uttara Bakshi, Ashis Datta","doi":"10.1016/j.jmgm.2026.109309","DOIUrl":"10.1016/j.jmgm.2026.109309","url":null,"abstract":"<div><div>Epilepsy remains a debilitating neurological disorder affecting approximately 50 million individuals worldwide. Current Anti-seizure medications (ASM) targeting Synaptic Vesicle Protein 2A (SV2A), such as levetiracetam, exhibit variable efficacy and adverse effect profiles, necessitating exploration of novel therapeutic approaches. This study investigates <em>Bacopa monnieri</em>-derived compounds as potential dual inhibitors of SV2A and Carbonic Anhydrase II (CA-II), two mechanistically distinct targets implicated in seizure pathophysiology. An integrated computational workflow comprising ADMET profiling, molecular docking, molecular dynamics (MD) simulations, and MM-GBSA binding free energy calculations was implemented to evaluate selected phytoconstituents. ADMET analysis revealed favourable pharmacokinetic parameters for several compounds, particularly Ebelin Lactone and Jujubogenin, with superior blood-brain barrier permeability (BBB >0.9) and human intestinal absorption (>95 %) compared to reference drugs. Molecular docking identified compelling binding affinities, with Ebelin Lactone (−11.2 kcal/mol) and Jujubogenin (−10.9 kcal/mol) exhibiting stronger interactions with SV2A than brivaracetam (−6.8 kcal/mol). Similarly, these compounds demonstrated robust binding to CA-II (−8.4 and −8.7 kcal/mol, respectively). Protein-ligand interaction profiler analysis elucidated key stabilizing interactions, including hydrogen bonds with Thr199 and His64 in CA-II and hydrophobic contacts within the SV2A binding pocket. MD simulations confirmed structural stability of the protein-ligand complexes, evidenced by equilibrated RMSD trajectories (0.2–0.3 nm). MM-GBSA calculations substantiated the thermodynamic favourability of these interactions, with binding free energies for top candidates significantly exceeding reference compounds. These findings establish an in silico framework for prioritizing <em>Bacopa monnieri</em> phytoconstituents as computationally predicted dual-target candidates, highlighting their potential relevance for seizure management, with the possibility of improved efficacy through simultaneous modulation of distinct pathophysiological mechanisms. All findings are based solely on in silico analyses and require experimental validation.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109309"},"PeriodicalIF":3.0,"publicationDate":"2026-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146079021","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}
Incorporating hole-transporting materials (HTMs) with optimal hole mobility and solution-processability is crucial for modifying effective materials of solar cells. In this investigation, we designed eight molecules with a D-A-type arrangement. The modified hole-transporting materials were studied using a quantum computation approach using density functional theory to found structural properties related to the electrochemical, charge transfer, quantum physical, solubility, and photovoltaic properties. The outcomes reveal that accepting fragments manifested hole-transport materials appropriate band alignment with deeper levels (ranging from −6.50 to −6.76 eV), higher absorption coefficients, remarkable solution processibility, and hole mobility with low exciton binding energy. These features revealed a higher photocurrent-generating ability, as estimated from transition density calculations across the molecular frameworks, a low charge-coupling estimated by the lower reorganization energy, and robust exciton dissociation. These notable outcomes unveiled that modified molecules are comparatively better than WWC-103 as HTMs for fabricating efficient material in the photovoltaic industry.
{"title":"Significantly improved optoelectronic properties of WWC-103 engineered for efficient perovskite solar cells: A DFT approach","authors":"Umer Yaqoob , Sidra Rafiq , Muhammad Zohaib Sabir , Ali Raza Ayub , Salba , Makhvela Anwer , Asma Parveen","doi":"10.1016/j.jmgm.2026.109300","DOIUrl":"10.1016/j.jmgm.2026.109300","url":null,"abstract":"<div><div>Incorporating hole-transporting materials (HTMs) with optimal hole mobility and solution-processability is crucial for modifying effective materials of solar cells. In this investigation, we designed eight molecules with a D-A-type arrangement. The modified hole-transporting materials were studied using a quantum computation approach using density functional theory to found structural properties related to the electrochemical, charge transfer, quantum physical, solubility, and photovoltaic properties. The outcomes reveal that accepting fragments manifested hole-transport materials appropriate band alignment with deeper <span><math><mrow><msub><mi>E</mi><mrow><mi>H</mi><mi>O</mi><mi>M</mi><mi>O</mi></mrow></msub></mrow></math></span> levels (ranging from −6.50 to −6.76 eV), higher absorption coefficients, remarkable solution processibility, and hole mobility with low exciton binding energy. These features revealed a higher photocurrent-generating ability, as estimated from transition density calculations across the molecular frameworks, a low charge-coupling estimated by the lower reorganization energy, and robust exciton dissociation. These notable outcomes unveiled that modified molecules are comparatively better than WWC-103 as HTMs for fabricating efficient material in the photovoltaic industry.</div></div>","PeriodicalId":16361,"journal":{"name":"Journal of molecular graphics & modelling","volume":"144 ","pages":"Article 109300"},"PeriodicalIF":3.0,"publicationDate":"2026-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146035986","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.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}