Expected distributions of root-mean-square positional deviations in proteins

Understanding the interactions between repurposed drugs sertindole and temoporfin with receptor for advanced glycation endproducts: Therapeutic implications in cancer and metabolic diseases

CONTEXT

In the pursuit of novel therapeutic possibilities, repurposing existing drugs has gained prominence as an efficient strategy. The findings from our study highlight the potential of repurposed drugs as promising candidates against receptor for advanced glycation endproducts (RAGE) that offer therapeutic implications in cancer, neurodegenerative conditions and metabolic syndromes. Through careful analyses of binding affinities and interaction patterns, we identified a few promising candidates, ultimately focusing on sertindole and temoporfin. These candidates exhibited exceptional binding affinities, efficacy, and specificity within the RAGE binding pocket. Notably, they displayed a pronounced propensity to interact with the active site of RAGE. Our investigation further revealed that sertindole and temoporfin possess desirable pharmacological properties that highlighted them as attractive candidates for targeted drug development. Overall, our integrated computational approach provides a comprehensive understanding of the interactions between repurposed drugs, sertindole and temoporfin and RAGE that pave the way for future experimental validation and drug development endeavors.

METHODS

We present an integrated approach utilizing molecular docking and extensive molecular dynamics (MD) simulations to evaluate the potential of FDA-approved drugs, sourced from DrugBank, against RAGE. To gain deeper insights into the binding mechanisms of the elucidated candidate repurposed drugs, sertindole and temoporfin with RAGE, we conducted extensive all-atom MD simulations, spanning 500 nanoseconds (ns). These simulations elucidated the conformational dynamics and stability of the RAGE-sertindole and RAGE-temoporfin complexes.

Potential mechanisms of traditional Chinese medicine in treating insomnia: A network pharmacology, GEO validation, and molecular-docking study

The purpose of this study is to investigate the potential mechanisms of Chinese herbs for the treatment of insomnia using a combination of data mining, network pharmacology, and molecular-docking validation. All the prescriptions for insomnia treated by the academician Qi Wang from 2020 to 2022 were collected. The Ancient and Modern Medical Case Cloud Platform v2.3 was used to identify high-frequency Chinese medicinal herbs and the core prescription. The Traditional Chinese Medicine Systems Pharmacology and UniProt databases were utilized to predict the effective active components and targets of the core herbs. Insomnia-related targets were collected from 4 databases. The intersecting targets were utilized to build a protein-protein interaction network and conduct gene ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis using the STRING database, Cytoscape software, and clusterProfiler package. Gene chip data (GSE208668) were obtained from the Gene Expression Omnibus database. The limma package was applied to identify differentially expressed genes (DEGs) between insomnia patients and healthy controls. To create a "transcription factor (TF)-miRNA-mRNA" network, the differentially expressed miRNAs were entered into the TransmiR, FunRich, Targetscan, and miRDB databases. Subsequently, the overlapping targets were validated using the DEGs, and further validations were conducted through molecular docking and molecular dynamics simulations. Among the 117 prescriptions, 65 herbs and a core prescription were identified. Network pharmacology and bioinformatics analysis revealed that active components such as β-sitosterol, stigmasterol, and canadine acted on hub targets, including interleukin-6, caspase-3, and hypoxia-inducible factor-1α. In GSE208668, 6417 DEGs and 7 differentially expressed miRNAs were identified. A "TF-miRNA-mRNA" network was constructed by 4 "TF-miRNA" interaction pairs and 66 "miRNA-mRNA" interaction pairs. Downstream mRNAs exert therapeutic effects on insomnia by regulating circadian rhythm. Molecular-docking analyses demonstrated good docking between core components and hub targets. Molecular dynamics simulation displayed the strong stability of the complex formed by small molecule and target. The core prescription by the academician Qi Wang for treating insomnia, which involves multiple components, targets, and pathways, showed the potential to improve sleep, providing a basis for clinical treatment of insomnia.

Repurposing Drugs to Modulate Sortilin: Structure-Guided Strategies Against Atherogenesis, Coronary Artery Disease, and Neurological Disorders

Sortilin (SORT1) is a multifunctional protein intricately involved in atherogenesis, coronary artery disease (CAD), and various neurological disorders. It has materialized as a potential pharmacological target for therapeutic development due to its diverse biological roles in pathological processes. Despite its central role under these conditions, effective therapeutic strategies targeting SORT1 remain challenging. In this study, we introduce a drug repurposing strategy guided by structural insights to identify potent SORT1 inhibitors with broad therapeutic potential. Our approach combines molecular docking, virtual screening, and molecular dynamics (MD) simulations, enabling the systematic evaluation of 3648 FDA-approved drugs for their potential to modulate SORT1. The investigation reveals a subset of repurposed drugs exhibiting highly favorable binding profiles and stable interactions within the binding site of SORT1. Notably, two hits, ergotamine and digitoxin, were carefully chosen based on their drug profiles and subjected to analyze their interactions with SORT1 and stability assessment via all-atom MD simulations spanning 300 ns (ns). The structural analyses uncover the complex binding interactions between these identified compounds and SORT1, offering essential mechanistic insights. Additionally, we explore the clinical implications of repurposing these compounds as potential therapeutic agents, emphasizing their significance in addressing atherogenesis, CAD, and neurological disorders. Overall, this study highlights the efficacy of structure-guided drug repurposing and provides a solid foundation for future research endeavors aimed at the development of effective therapies targeting SORT1 under diverse pathological conditions.

Evaluation of the Binding Mechanism of Dietary Phytochemical, Ellagic Acid, with Human Transferrin: Spectroscopic, Calorimetric, and Computational Approaches Targeting Neurodegenerative Diseases

Human transferrin (Htf) is vital in maintaining iron within the brain cells; any disruption results in the development of neurodegenerative diseases (NDs) and other related pathologies, especially Alzheimer's disease (AD). Ellagic acid (EA), a naturally occurring phenolic antioxidant, possesses neuroprotective potential and is present in a broad variety of fruits and vegetables. The current work explores the binding mechanism of dietary polyphenol, EA, with Htf by a combination of experimental and computational approaches. Molecular docking studies unveiled the binding of EA to Htf with good affinity. Molecular dynamic (MD) simulation further provided atomistic details of the binding process, demonstrating a stable Htf-EA complex formation without causing substantial alterations to the protein's conformation. Furthermore, fluorescence binding measurements indicated that EA forms a high-affinity interaction with Htf. Isothermal titration calorimetric measurements advocated the spontaneous nature of binding and also revealed the binding process to be exothermic. In conclusion, the study deciphered the binding mechanism of EA with Htf. The results demonstrated that EA binds with Htf with an excellent affinity spontaneously, thereby laying the groundwork for potential applications of EA in the realm of therapeutics for NDs in the context of iron homeostasis.

Exploring the potential of phytochemicals as inhibitors of 3'-phosphoadenosine 5'-phosphosulfate synthase 1 targeting cancer therapy

3'-phosphoadenosine 5'-phosphosulfate synthase 1 (PAPSS1) is an enzyme that critically synthesises the biologically active form of sulfate (PAPS) for all sulfation reactions. The discovery of PAPSS1 as a possible drug target for cancer therapy, specifically in non-small cell lung cancer, has prompted us to investigate potential small-molecule inhibitors of PAPSS1. Here, a structure-based virtual screening method was used to search for phytochemicals in the IMPPAT database to find potential inhibitors of PAPSS1. The primary hits were selected based on their physicochemical, ADMET, and drug-like properties. Then, the binding affinities were calculated and analyzed the interactions to identify safer and more effective hits. The research identified two phytochemicals, Guggulsterone and Corylin, that exhibited significant affinity and specific interaction with the ATP-binding pocket of PAPSS1. Structural observations made by molecular docking were further accompanied by molecular dynamics (MD) simulations and principal component analysis (PCA) to examine the conformational changes and stability of PAPSS1 with the elucidated compounds Guggulsterone and Corylin. MD simulation results suggested that the binding of Guggulsterone and Corylin stabilizes the PAPSS1 structure, leading to fewer conformational changes. This implies that these compounds may be useful in developing PAPSS1 inhibitors for the therapeutic development against non-small cell lung cancer (NSCLC). This study highlights the potential of phytochemicals as PAPSS1 inhibitors and the utility of computational approaches in drug discovery.Communicated by Ramaswamy H. Sarma.

Identification of β-cycloidal-derived mono-carbonyl curcumin analogs as potential interleukin-6 inhibitor to treat wound healing through QSAR, molecular docking, MD simulation, MM-GBSA calculation

Interleukin-6 (IL-6) is a cytokine that involved in the different phases of wound healing. It is responsible for promoting inflammation, regulating tissue repair scar formation, stimulating the production of extracellular matrix components and recruiting immune cells to the wound site. Therefore, suppressing IL-6 is beneficial for wound healing. However, no small molecules are currently available in the market against the IL-6. As a result, this research gap motivates us to find a potential inhibitor. This study aimed to investigate the wound healing potential of novel β-cycloidal-derived mono-carbonyl curcumin analogs reported in the literature through screening a series of computational studies. The calculated pIC50 value of 18 compounds (below 10) showed that all compounds may have potential therapeutic efficacy. Molecular docking studies revealed that compound C12 (-45.6044 kcal/mol) bound most strongly in the active site of IL-6 compared to the FDA-approved drug clindamycin (-42.3223). The Molecular Dynamic (MD) simulation displayed that lead compound C12 had the highest stability in the active site of IL-6 compared to the reference drug clindamycin. Furthermore, MMGBSA results indicated that C12 (-20.28 kcal/mol) had the highest binding energy compared to clindamycin (-8.36 kcal/mol). The ADMET analysis predicted that C12 are favourable for drug candidates. This study recommended compound C12 as a lead IL-6 inhibitor for future testing and development as therapeutics for wound healing.Communicated by Ramaswamy H. Sarma.

The development of prognostic gene markers associated with disulfidptosis in gastric cancer and their application in predicting drug response

Background

Gastric cancer (GC) is a malignancy known for its high fatality rate. Disulfidptosis, a potentially innovative therapeutic strategy for cancer treatment, has been proposed. Nevertheless, the specific involvement of disulfidptosis in the context of GC remains uncertain.

Methods

The mRNA expression profiles were obtained from the TCGA and GEO databases. Univariate and LASSO Cox regression analyses were employed to identify differentially expressed genes and develop a risk model for disulfidptosis-related genes. The performance of the model was evaluated using Kaplan-Meier curve, ROC curve, and nomogram. Univariate and multivariate Cox regression analyses were conducted to determine if the risk model could serve as an independent prognostic factor. The biological function of the identified genes was assessed through GO, KEGG, and GSEA analyses. The prediction of drug response was conducted employing the package "pRRophetic". Furthermore, gene expression was determined using qRT-PCR.

Results

An eight-gene signature were identified and utilized to categorize patients into low- and high-risk groups. Survival, receiver operating characteristic (ROC) curve, and Cox analyses provided clarification that these eight hub genes served as a favorable independent prognostic factor for patients with GC. A nomogram was constructed by integrating clinical parameters with the risk signatures, demonstrating high precision in predicting 1-, 3-, and 5-year survival rates. Additionally, drug sensitivity was different in the high-risk and low-risk groups, and the expression of three genes was verified by qRT-PCR.

Conclusion

The prognostic risk model developed in this study demonstrates the potential to accurately forecast the prognosis of patients with GC.

Identification of Phytochemicals from Arabian Peninsula Medicinal Plants as Strong Binders to SARS-CoV-2 Proteases (3CLPro and PLPro) by Molecular Docking and Dynamic Simulation Studies

We provide promising computational (in silico) data on phytochemicals (compounds 1-10) from Arabian Peninsula medicinal plants as strong binders, targeting 3-chymotrypsin-like protease (3CLPro) and papain-like proteases (PLPro) of SARS-CoV-2. Compounds 1-10 followed the Lipinski rules of five (RO5) and ADMET analysis, exhibiting drug-like characters. Non-covalent (reversible) docking of compounds 1-10 demonstrated their binding with the catalytic dyad (CYS145 and HIS41) of 3CLPro and catalytic triad (CYS111, HIS272, and ASP286) of PLPro. Moreover, the implementation of the covalent (irreversible) docking protocol revealed that only compounds 7, 8, and 9 possess covalent warheads, which allowed the formation of the covalent bond with the catalytic dyad (CYS145) in 3CLPro and the catalytic triad (CYS111) in PLPro. Root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), and radius of gyration (Rg) analysis from molecular dynamic (MD) simulations revealed that complexation between ligands (compounds 7, 8, and 9) and 3CLPro and PLPro was stable, and there was less deviation of ligands. Overall, the in silico data on the inherent properties of the above phytochemicals unravel the fact that they can act as reversible inhibitors for 3CLPro and PLPro. Moreover, compounds 7, 8, and 9 also showed their novel properties to inhibit dual targets by irreversible inhibition, indicating their effectiveness for possibly developing future drugs against SARS-CoV-2. Nonetheless, to confirm the theoretical findings here, the effectiveness of the above compounds as inhibitors of 3CLPro and PLPro warrants future investigations using suitable in vitro and in vivo tests.

Identifying repurposed drugs as potential inhibitors of Apolipoprotein E: A bioinformatics approach to target complex diseases associated with lipid metabolism and neurodegeneration

Apolipoprotein E (ApoE), a pivotal contributor to lipid metabolism and neurodegenerative disorders, emerges as an attractive target for therapeutic intervention. Within this study, we deployed an integrated in-silico strategy, harnessing structure-based virtual screening, to identify potential compounds from DrugBank database. Employing molecular docking, we unveil initial hits by evaluating their binding efficiency with ApoE. This first tier of screening narrows our focus to compounds that exhibit a strong propensity to bind with ApoE. Further, a detailed interaction analysis was carried out to explore the binding patterns of the selected hits towards the ApoE binding site. The selected compounds were then evaluated for the biological properties in PASS analysis, which showed anti-neurodegenerative properties. Building upon this foundation, we delve deeper, employing all-atom molecular dynamics (MD) simulations extending over an extensive 500 ns. In particular, Ergotamine and Dihydroergocristine emerge as noteworthy candidates, binding to ApoE in a competitive mode. This intriguing binding behavior positions these compounds as potential candidates warranting further analysis in the pursuit of novel therapeutics targeting complex diseases associated with lipid metabolism and neurodegeneration. This approach holds the promise of catalyzing advancements in therapeutic intervention for complex disorders, thereby reporting a meaningful pace towards improved healthcare outcomes.

Identification of mitogen-activated protein kinase 7 inhibitors from natural products: Combined virtual screening and dynamic simulation studies

Mitogen-activated protein kinase 7 (MAPK7) is a serine/threonine protein kinase that belongs to the MAPK family and plays a vital role in various cellular processes such as cell proliferation, differentiation, gene transcription, apoptosis, metabolism, and cell survival. The elevated expression of MAPK7 has been associated with the onset and progression of multiple aggressive tumors in humans, underscoring the potential of targeting MAPK7 pathways in therapeutic research. This pursuit holds promise for the advancement of anticancer drug development by developing potential MAPK7 inhibitors. To look for potential MAPK7 inhibitors, we exploited structure-based virtual screening of natural products from the ZINC database. First, the Lipinski rule of five criteria was used to filter a large library of ~90,000 natural compounds, followed by ADMET and pan-assay interference compounds (PAINS) filters. Then, top hits were chosen based on their strong binding affinity as determined by molecular docking. Further, interaction analysis was performed to find effective and specific compounds that can precisely bind to the binding pocket of MAPK7. Consequently, two compounds, ZINC12296700 and ZINC02123081, exhibited significant binding affinity and demonstrated excellent drug-like properties. All-atom molecular dynamics simulations for 200 ns confirmed the stability of MAPK7-ZINC12296700 and MAPK7-ZINC02123081 docked complexes. According to the molecular mechanics Poisson-Boltzmann surface area investigation, the binding affinities of both complexes were considerable. Overall, the result suggests that ZINC12296700 and ZINC02123081 might be used as promising leads to develop novel MAPK7 inhibitors. Since these compounds would interfere with the kinase activity of MAPK7, therefore, may be implemented to control cell growth and proliferation in cancer after required validations.

Unveiling the Molecular Interactions Between Human Transferrin and Limonene: Natural Compounds in Alzheimer's Disease Therapeutics

Background

Neurodegeneration is a term describing an irreversible process of neuronal damage. In recent decades, research efforts have been directed towards deepening our knowledge of numerous neurodegenerative disorders, with a particular focus on conditions such as Alzheimer's disease (AD). Human transferrin (htf) is a key player in maintaining iron homeostasis within brain cells. Any disturbance in this equilibrium gives rise to the emergence of neurodegenerative diseases and associated pathologies, particularly AD. Limonene, a natural compound found in citrus fruits and various plants, has shown potential neuroprotective properties.

Objective

In this study, our goal was to unravel the binding of limonene with htf, with the intention of comprehending the interaction mechanism of limonene with htf.

Methods

Binding was scrutinized using fluorescence quenching and UV-Vis spectroscopic analyses. The binding mechanism of limonene was further investigated at the atomic level through molecular docking and extensive 200 ns molecular dynamic simulation (MD) studies.

Results

Molecular docking uncovered that limonene interacted extensively with the deep cavity located within the htf binding pocket. MD results indicated that binding of limonene to htf did not induce substantial structural alterations, ultimately forming stable complex. The findings from fluorescence binding indicated a pronounced interaction between limonene and htf, limonene binds to htf with a binding constant (K) of 0.1×105 M-1. UV spectroscopy also advocated stable htf-limonene complex formation.

Conclusions

The study deciphered the binding mechanism of limonene with htf, providing a platform to use limonene in AD therapeutics in context of iron homeostasis.

Prediction of molecular interaction of Phosphodiesterase 10A inhibition by natural compounds: insights from structure-based screening and molecular dynamics simulations

Phosphodiesterase 10 A (PDE10A) is an enzyme that regulates cyclic adenosine monophosphate (cAMP) and cyclic guanosine monophosphate (cGMP) levels in the brain, particularly in the striatum, which plays a critical role in movement control and reward processing. Inhibition of PDE10A can increase cAMP and cGMP levels, improving neuronal signaling and reducing symptoms of neuropsychiatric disorders such as schizophrenia, Huntington's disease, and Parkinson's disease. In this study, a structure-based virtual screening was conducted to identify potential anti-neuropsychiatric disorders compounds from phytoconstituents in the IMPPAT database. The ligands were docked against PDE10A, resulting in 40 compounds with appreciable docking scores. These 40 compounds underwent further ADMET predictions and drug likeliness, resulting in five potential compounds. Finally, based on the specific interactions, two compounds (Colladonin and Isopongachromene), were subjected to molecular dynamics (MD) simulation and MM-PBSA studies. The MM-PBSA analysis validated and captured the intermolecular interactions, indicating that Colladonin and Isopongachromene had appreciable binding affinities of -155.60 kJ.mol-1 and -108.28 kJ.mol-1, respectively and were promising candidates against neuropsychiatric disorders, targeting PDE10A. Overall, this study provides insight into the potential of PDE10A inhibitors as therapeutic agents for treating neuropsychiatric disorders, and Colladonin and Isopongachromene are promising compounds for further development.Communicated by Ramaswamy H. Sarma.

Exploring natural compound, Panicutine as leucine-rich repeat kinase 2 inhibitor against Parkinson's disease: a structure-guided approach

Leucine-rich repeat kinase 2 (LRRK2) is a promising drug target for the therapeutic management of Parkinson's disease (PD) and other neurodegenerative disorders. LRRK2 inhibitors have the potential to modulate neuroinflammation, reduce alpha-synuclein aggregation and improve motor symptoms in PD patients. Although LRRK2 inhibitors are still in the early stages of clinical development, the identification of potent and selective inhibitors through structure-guided approaches provides a promising avenue for the development of effective therapies for PD and other neurodegenerative disorders. In this study, natural compounds from the IMPPAT database were screened using a state-of-the-art computational virtual screening approach to identify potential inhibitors of LRRK2. We carried out a docking screening on a library of natural compounds and identified a few compounds with strong binding affinity, docking score and specificity towards LRRK2 as the top hits. These hits were then subjected to further analysis based on multiple parameters for the Pan-assay interference compounds and their physicochemical and pharmacokinetics evaluation followed by a detailed interaction analysis. After careful evaluation, one natural compound, Panicutine, was identified as a promising candidate for LRRK2 due to its significant affinity and specificity towards the LRRK2 binding pocket. Additionally, it exhibited drug-like properties with blood-brain barrier permeability as determined by ADMET properties. To gain a deeper understanding of the stability and conformational changes of the LRRK2-ligand complex, MD simulations were conducted for 100 nanoseconds under explicit solvent conditions followed by principal component analysis and free energy dynamics. The simulation results demonstrated that the LRRK2-Panicutine complex remained stable throughout the simulation trajectories. Based on these findings, it is concluded that Panicutine has the potential to act as a LRRK2 inhibitor against PD and other neurodegenerative disorders.Communicated by Ramaswamy H. Sarma.

Phytoconstituents as potential therapeutic agents against COVID-19: a computational study on inhibition of SARS-CoV-2 main protease

The Coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) has become a global health crisis, and the urgent need for effective treatments is evident. One potential target for COVID-19 therapeutics is the main protease (Mpro) of SARS‑CoV‑2, an essential enzyme for viral replication. Natural compounds have been explored as a source of potential inhibitors for Mpro due to their safety and availability. In this study, we employed a computational approach to screen a library of phytoconstituents and identified potential Mpro inhibitors based on their binding affinities and molecular interactions. The top-ranking compounds were further validated through molecular dynamics simulations (MDS) and free energy calculations. As a result of the above procedures, we identified two phytoconstituents, Khelmarin B and Neogitogenin, with appreciable binding affinity and specificity towards the Mpro binding pocket. Our results suggest that Khelmarin B and Neogitogenin could potentially serve as Mpro inhibitors and have the potential to be developed as COVID-19 therapeutics. Further experimental studies are required to confirm the efficacy and safety of these compounds.Communicated by Ramaswamy H. Sarma.

Micro-simulation insights into the functional and mechanistic understanding of glycyrrhizin against asthma

Asthma is a common chronic respiratory disease, which causes inflammation and airway stenosis, leading to dyspnea, wheezing and chest tightness. Using transgelin-2 as a target, we virtually screened the lead compound glycyrrhizin from the self-built database of anti-asthma compounds by molecular docking technology, and found that it had anti-inflammatory, anti-oxidative and anti-asthma pharmacological effects. Then, molecular dynamics simulations were used to confirm the stability of the glycyrrhizin-transgelin-2 complex from a dynamic perspective, and the hydrophilic domains of glycyrrhizin was found to have the effect of targeting transgelin-2. Due to the self-assembly properties of glycyrrhizin, we explored the formation process and mechanism of the self-assembly system using self-assembly simulations, and found that hydrogen bonding and hydrophobic interactions were the main driving forces. Because of the synergistic effect of glycyrrhizin and salbutamol in improving asthma, we revealed the mechanism through simulation, and believed that salbutamol adhered to the surface of the glycyrrhizin nano-drug delivery system through hydrogen bonding and hydrophobic interactions, using the targeting effect of the hydrophilic domains of glycyrrhizin to reach the pathological parts and play a synergistic anti-asthmatic role. Finally, we used network pharmacology to predict the molecular mechanisms of glycyrrhizin against asthma, which indicated the direction for its clinical transformation.

Identification of potential Escherichia coli DNA gyrase B inhibitors targeting antibacterial therapy: an integrated docking and molecular dynamics simulation study

The alarming rise in the rate of antibiotic resistance is a matter of significant concern. DNA gyrase B (GyrB), a critical bacterial enzyme involved in DNA replication, transcription, and recombination, has emerged as a promising target for antibacterial agents. Inhibition of GyrB disrupts bacterial DNA replication, leading to cell death, making it an attractive candidate for antibiotic development. Although several classes of antibiotics targeting GyrB are currently in clinical use, the emergence of antibiotic resistance necessitates the exploration of novel inhibitors. In this study, we aimed to identify potential Escherichia coli GyrB inhibitors from a database of phytoconstituents sourced from Indian medicinal plants. Utilizing virtual screening, we performed a rigorous search to identify compounds with the most promising inhibitory properties against GyrB. Two compounds, namely Zizogenin and Cucurbitacin S, were identified based on their favorable drug likeliness and pharmacokinetic profiles. Employing advanced computational techniques, we analyzed the binding interactions of Zizogenin and Cucurbitacin S with the ATP-binding site of GyrB through molecular docking simulations. Both compounds exhibited robust binding interactions, evidenced by their high docking energy scores. To assess the stability of these interactions, we conducted extensive 100 ns molecular dynamics (MD) simulations, which confirmed the stability of Zizogenin and Cucurbitacin S when bound to GyrB. In conclusion, our study highlights Zizogenin and Cucurbitacin S as promising candidates for potential antibacterial agents targeting GyrB. Experimental validation of these compounds is warranted to further explore their efficacy and potential as novel antibiotics to combat antibiotic-resistant bacteria.Communicated by Ramaswamy H. Sarma.

Genome-wide analysis of the Populus trichocarpa laccase gene family and functional identification of PtrLAC23

INTRODUCTION

Biofuel is a kind of sustainable, renewable and environment friendly energy. Lignocellulose from the stems of woody plants is the main raw material for "second generation biofuels". Lignin content limits fermentation yield and is therefore a major obstacle in biofuel production. Plant laccase plays an important role in the final step of lignin formation, which provides a new strategy for us to obtain ideal biofuels by regulating the expression of laccase genes to directly gain the desired lignin content or change the composition of lignin.

METHODS

Multiple sequence alignment and phylogenetic analysis were used to classify PtrLAC genes; sequence features of PtrLACs were revealed by gene structure and motif composition analysis; gene duplication, interspecific collinearity and Ka/Ks analysis were conducted to identify ancient PtrLACs; expression levels of PtrLAC genes were measured by RNA-Seq data and qRT-PCR; domain analysis combine with cis-acting elements prediction together showed the potential function of PtrLACs. Furthermore, Alphafold2 was used to simulate laccase 3D structures, proLAC23::LAC23-eGFP transgenic Populus stem transects were applied to fluorescence observation.

RESULTS

A comprehensive analysis of the P. trichocarpa laccase gene (PtLAC) family was performed. Some ancient PtrLAC genes such as PtrLAC25, PtrLAC19 and PtrLAC41 were identified. Gene structure and distribution of conserved motifs clearly showed sequence characteristics of each PtrLAC. Combining published RNA-Seq data and qRT-PCR analysis, we revealed the expression pattern of PtrLAC gene family. Prediction results of cis-acting elements show that PtrLAC gene regulation was closely related to light. Through above analyses, we selected 5 laccases and used Alphafold2 to simulate protein 3D structures, results showed that PtrLAC23 may be closely related to the lignification. Fluorescence observation of proLAC23::LAC23-eGFP transgenic Populus stem transects and qRT-PCR results confirmed our hypothesis again.

DISCUSSION

In this study, we fully analyzed the Populus trichocarpa laccase gene family and identified key laccase genes related to lignification. These findings not only provide new insights into the characteristics and functions of Populus laccase, but also give a new understanding of the broad prospects of plant laccase in lignocellulosic biofuel production.

Unveiling Phytoconstituents with Inhibitory Potential Against Tyrosine-Protein Kinase Fyn: A Comprehensive Virtual Screening Approach Targeting Alzheimer's Disease

BACKGROUND

Tyrosine-protein kinase Fyn (Fyn) is a critical signaling molecule involved in various cellular processes, including neuronal development, synaptic plasticity, and disease pathogenesis. Dysregulation of Fyn kinase has been implicated in various complex diseases, including neurodegenerative disorders such as Alzheimer's and Parkinson's diseases, as well as different cancer types. Therefore, identifying small molecule inhibitors that can inhibit Fyn activity holds substantial significance in drug discovery.

OBJECTIVE

The aim of this study was to identify potential small-molecule inhibitors among bioactive phytoconstituents against tyrosine-protein kinase Fyn.

METHODS

Through a comprehensive approach involving molecular docking, drug likeliness filters, and molecular dynamics (MD) simulations, we performed a virtual screening of a natural compounds library. This methodology aimed to pinpoint compounds potentially interacting with Fyn kinase and inhibiting its activity.

RESULTS

This study finds two potential natural compounds: Dehydromillettone and Tanshinone B. These compoundsdemonstrated substantial affinity and specific interactions towards the Fyn binding pocket. Their conformations exhibitedcompatibility and stability, indicating the formation of robust protein-ligand complexes. A significant array of non-covalentinteractions supported the structural integrity of these complexes.

CONCLUSION

Dehydromillettone and Tanshinone B emerge as promising candidates, poised for further optimization as Fynkinase inhibitors with therapeutic applications. In a broader context, this study demonstrates the potential of computationaldrug discovery, underscoring its utility in identifying compounds with clinical significance. The identified inhibitors holdpromise in addressing a spectrum of cancer and neurodegenerative disorders. However, their efficacy and safety necessitatevalidation through subsequent experimental studies.

Nature-Derived Compounds as Potential Bioactive Leads against CDK9-Induced Cancer: Computational and Network Pharmacology Approaches

Given the importance of cyclin-dependent kinases (CDKs) in the maintenance of cell development, gene transcription, and other essential biological operations, CDK blockers have been generated to manage a variety of disorders resulting from CDK irregularities. Furthermore, CDK9 has a crucial role in transcription by regulating short-lived anti-apoptotic genes necessary for cancer cell persistence. Addressing CDK9 with blockers has consequently emerged as a promising treatment for cancer. This study scrutinizes the effectiveness of nature-derived compounds (geniposidic acid, quercetin, geniposide, curcumin, and withanolide C) against CDK9 through computational approaches. A molecular docking study was performed after preparing the protein and the ligands. The selected blockers of the CDK9 exerted reliable binding affinities (−8.114 kcal/mol to −13.908 kcal/mol) against the selected protein, resulting in promising candidates compared to the co-crystallized ligand (LCI). The binding affinity of geniposidic acid (−13.908 kcal/mol) to CDK9 is higher than quercetin (−10.775 kcal/mol), geniposide (−9.969 kcal/mol), curcumin (−9.898 kcal/mol), withanolide C (−8.114 kcal/mol), and the co-crystallized ligand LCI (−11.425 kcal/mol). Therefore, geniposidic acid is a promising inhibitor of CDK9. Moreover, the molecular dynamics studies assessed the structure–function relationships and protein–ligand interactions. The network pharmacology study for the selected ligands demonstrated the auspicious compound–target–pathway signaling pathways vital in developing tumor, tumor cell growth, differentiation, and promoting tumor cell progression. Moreover, this study concluded by analyzing the computational approaches the natural-derived compounds that have potential interacting activities against CDK9 and, therefore, can be considered promising candidates for CKD9-induced cancer. To substantiate this study’s outcomes, in vivo research is recommended.

Battling BTK mutants with noncovalent inhibitors that overcome Cys481 and Thr474 mutations in Waldenström macroglobulinemia therapy: structural mechanistic insights on the role of fenebrutinib

Recently, the non-covalent Bruton tyrosine kinase (BTK) inhibitor fenebrutinib was presented as a therapeutic option with strong inhibitory efficacy against a single (C481S) and double (T474S/C481S) BTK variant in the treatment of Waldenström macroglobulinemia (WM). However, the molecular events surrounding its inhibition mechanism towards this variant remain unresolved. Herein, we employed in silico methods such as molecular dynamic simulation coupled with binding free energy estimations to explore the mechanistic activity of the fenebrutinib on (C481S) and (T474S/C481S) BTK variant, at a molecular level. Our investigations reveal that amino acid arginine contributed immensely to the total binding energy, this establishing the cruciality of amino acid residues, Arg132 and Arg156 in (C481S) and Arg99, Arg137, and Arg132 in (T474S/C481S) in the binding of fenebrutinib towards both BTK variants. The structural orientations of fenebrutinib within the respective hydrophobic pockets allowed favorable interactions with binding site residues, accounting for its superior binding affinity by 24.5% and relative high hydrogen bond formation towards (T474S/C481S) when compared with (C481S) BTK variants. Structurally, fenebrutinib impacted the stability, flexibility, and solvent accessible surface area of both BTK variants, characterized by various alterations observed in the bound and unbound structures, which proved enough to disrupt their biological function. Findings from this study, therefore, provide insights into the inhibitory mechanism of fenebrutinib at the atomistic level and reveal its high selectivity towards BTK variants. These insights could be key in designing and developing BTK mutants' inhibitors to treat Waldenström macroglobulinemia (WM).

Network pharmacology and bioinformatics analyses identify the intersection genes and mechanism of Huang Bai for recurrent aphthous stomatitis

Recurrent aphthous stomatitis (RAS) are complex inflammatory diseases caused by multi-factors, which severely impact patient quality of life. However, there is still no effective treatment method for RAS without side effects. Traditionally, Cortex Phellodendri known as "Huang Bai" was used to treat RAS for antibacterial and anti-inflammatory properties in China. Network pharmacology methods and bioinformatics analysis were utilized to search and fish incorporating target. Network analysis and silico validation were used to discover the pharmacological mechanisms of "Huang Bai" for the treatment of RAS. A total of 25 active ingredients in HB, 200 drug targets, and 578 differentially expressed genes (DEGs) between Recurrent aphthous stomatitis and normal samples were obtained. The Gene Ontology enrichment analysis revealed that the immune response was the most significantly enriched term within the DEGs. The KEGG pathway analysis identified 60 significant pathways, most of which involved in the inhibition of inflammation and regulation of immunological response. The functions are dependent on a multi-pathway, particularly the TNF signaling pathway and the HIF-1 signaling pathway. We identified six hub genes in the PPI network, most of which were validated as highly expressed in oral ulcers by DiseaseMeth databases. In addition, molecular docking displayed that the primary molecule combined well with the key targets. "Huang Bai" contains potential anti-RAS active compounds. This study reflects the multi-component multi-target multi-pathway action characteristics of "Huang Bai." Our study provides potential biomarkers or treatment targets for further research.

A Dentin Biomimetic Remineralization Material with an Ability to Stabilize Collagen

The integrity of collagen matrix structure is a prerequisite for effectively inducing biomimetic remineralization. Repeated low pH stimulation activates matrix metalloproteinases (MMPs) in dental caries. Activated MMPs cause the breakdown of collagen fibrils. Collagen stabilization is a major obstacle to the clinical application of remineralization templates. Here, galardin-loaded poly(amido amine) (PAMAM)-NGV (PAMAM-NGV@galardin, PNG) is constructed to induce collagen stabilization and dentin biomimetic remineralization simultaneously, in order to combat early caries in dentin. PAMAM acts in the role of nucleation template for dentin remineralization, while galardin acts as the role of MMPs inhibitor. NGV peptides modified on the surface of dendrimer core can form small clusters with synergistic movement in short range, and those short-range clusters can form domain areas with different properties on the surface of PAMAM core and restrict the movement of collagen, favoring collagen crosslinking, which can be explained through the computational simulation analysis results. NGV peptides and galardin show a dual collagen-protective effect, laying the foundation for the dentin remineralization effect induced by PAMAM. PNG induces dentin remineralization in an environment with collagenase, meanwhile showsing anti-dentin caries efficacy in vivo. These findings indicate that PNG has great potential to combat early dentin caries for future clinical application.

Repurposing antiviral phytochemicals from the leaf extracts of Spondias mombin (Linn) towards the identification of potential SARSCOV-2 inhibitors

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), a pneumonia-like disease with a pattern of acute respiratory symptoms, currently remains a significant public health concern causing tremendous human suffering. Although several approved vaccines exist, vaccine hesitancy, limited vaccine availability, high rate of viral mutation, and the absence of approved drugs account for the persistence of SARS-CoV-2 infections. The investigation of possibly repurposing of phytochemical compounds as therapeutic alternatives has gained momentum due to their reported affordability and minimal toxicity. This study investigated anti-viral phytochemical compounds from ethanolic leaf extracts of Spondias mombin L as potential inhibitor candidates against SARS-CoV-2. We identified Geraniin and 2-O-Caffeoyl-(+)-allohydroxycitric acid as potential SARS-CoV-2 inhibitor candidates targeting the SARS-CoV-2 RNA-dependent polymerase receptor-binding domain (RBD) of SARS-CoV-2 viral S-protein and the 3C-like main protease (3CLpro). Geraniin exhibited binding free energy (ΔGbind) of - 25.87 kcal/mol and - 21.74 kcal/mol towards SARS-CoV-2 RNA-dependent polymerase and receptor-binding domain (RBD) of SARS-CoV-2 viral S-protein respectively, whereas 2-O-Caffeoyl-(+)-allohydroxycitric acid exhibited a ΔGbind of - 32 kcal/mol towards 3CLpro. Molecular Dynamics simulations indicated a possible interference to the functioning of SARS-CoV-2 targets by the two identified inhibitors. However, further in vitro and in vivo evaluation of these potential SARS-CoV-2 therapeutic inhibitor candidates is needed.

"Mango in all her majesty"-the potential of mangiferin and its analogues in the inhibition of Eimeria tenella hexokinase-a computational approach

The potential of natural products in mitigating infections and diseases are being considered lately. Herein, via in silico methods, we report the possible molecular mechanism of mangiferin (isolated from the fruit, peel, bark and leaves of mango tree) and its derivatives in inhibiting Eimeria tenella hexokinase. We evaluated the binding affinity of these inhibitors to the glucose binding site of EtHK and thereafter proceeded to molecular dynamics simulation. The Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) reveals that three of the derivatives (CPAMM, MxPAMM and NAMM) had better total binding free energy than mangiferin. The ADMET and physicochemical properties assessed shows that inhibitors also hold a potential to be drug-likely. Finally, in mediating their inhibitory potentials, the ligands stabilize both the global and local structures of the protein. This study provides a theoretical premise on which the anti-coccidial propensities of mangiferin most especially its derivatives can be investigate in vitro and in vivo.Communicated by Ramaswamy H. Sarma.

Unravelling the Structural Mechanism of Action of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione in Dual-Targeting Tankyrase 1 and 2: A Novel Avenue in Cancer Therapy

Tankyrase (TNKS) belonging to the poly(ADPribose) polymerase family, are known for their multi-functioning capabilities, and play an essential role in the Wnt β-catenin pathway and various other cellular processes. Although showing inhibitory potential at a nanomolar level, the structural dual-inhibitory mechanism of the novel TNKS inhibitor, 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione, remains unexplored. By employing advanced molecular modeling, this study provides these insights. Results of sequence alignments of binding site residues identified conserved residues; GLY1185 and ILE1224 in TNKS-1 and PHE1035 and PRO1034 in TNKS-2 as crucial mediators of the dual binding mechanism of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione, corroborated by high per-residue energy contributions and consistent high-affinity interactions of these residues. Estimation of the binding free energy of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione showed estimated total energy of -43.88 kcal/mol and -30.79 kcal/mol towards TNKS-1 and 2, respectively, indicating favorable analogous dual binding as previously reported. Assessment of the conformational dynamics of TNKS-1 and 2 upon the binding of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione revealed similar structural changes characterized by increased flexibility and solvent assessible surface area of the residues inferring an analogous structural binding mechanism. Insights from this study show that peculiar, conserved residues are the driving force behind the dual inhibitory mechanism of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione and could aid in the design of novel dual inhibitors of TNKS-1 and 2 with improved therapeutic properties.

Screening of novel and selective inhibitors for neuronal nitric oxide synthase (nNOS) via structure-based drug design techniques

NO, or nitric oxide, is produced by a family of enzymes called nitric oxide synthase (NOS) from L-arginine. NO regulates many physiological functions such as smooth muscle relaxation, immune defense, and memory function. The overproduction of NO by the neuronal isoform of nitric oxide synthase (nNOS) is implicated in neurodegeneration and neuropathic pain, making nNOS inhibition a promising therapeutic approach. Many developed nNOS inhibitors, generally L-arginine mimetics, have some issues in selectivity and bioavailability. According to earlier studies, targeting nNOS has the advantage of decreasing excess NO in the brain while avoiding the negative consequences of inhibiting the two isozymes: endothelial NOS (eNOS) and inducible NOS (iNOS). This study applied structure-based virtual screening, molecular docking, and molecular dynamics simulations to design potent and selective inhibitors against nNOS over related isoforms (eNOS and iNOS) using human X-ray crystal structures of the NOS isoforms. It was discovered that some compounds displayed a very good inhibitory potency for hnNOS and moderate selectivity for the other isozymes, eNOS and iNOS, in addition to good solubility and desirable physiochemical properties. The compounds which showed good stability and selectivity with nNOS, such as ZINC000013485422, can be interesting and informative guidance for designing more potent human nNOS inhibitors.Communicated by Ramaswamy H. Sarma.

Conjugated Dienones from Differently Substituted Cinnamaldehyde as Highly Potent Monoamine Oxidase-B Inhibitors: Synthesis, Biochemistry, and Computational Chemistry

Fifteen multiconjugated dienones (MK1-MK15) were synthesized and evaluated to determine their inhibitory activities against monoamine oxidases (MAOs) A and B. All derivatives were found to be potent and highly selective MAO-B inhibitors. Compound MK6, with an IC₅₀ value of 2.82 nM, most effectively inhibited MAO-B, like MK12 (IC₅₀ = 3.22 nM), followed by MK5, MK13, and MK14 (IC₅₀ = 4.02, 4.24, and 4.89 nM, respectively). The selectivity index values of MK6 and MK12 for MAO-B over MAO-A were 7361.5 and 1780.5, respectively. Compounds MK6 and MK12 were competitive reversible inhibitors of MAO-B, with K _i values of 1.10 ± 0.20 and 3.0 ± 0.27 nM, respectively. Cytotoxic studies showed that MK5, MK6, MK12, and MK14 exhibited low toxicities on Vero cells, with IC₅₀ values of 218.4, 149.1, 99.96, and 162.3 μg/mL, respectively, which were much higher than those for their effective nanomolar-level concentrations. Also, MK5, MK6, MK12, and MK14 decreased cell damage in H₂O₂-induced cells via a significant scavenging effect of reactive oxygen species. Molecular modeling was performed to rationalize the potential inhibitory activities of MK5, MK6, MK12, and MK14 toward MAO-B and their possible binding mechanisms, showing high-affinity binding pocket interactions and conformation perturbations of the compounds with MAO-B, which were interpreted as the conformational dynamics of MAO-B. This study concluded that all the compounds tested were more potent MAO-B inhibitors than the reference drugs, and leading compounds could be further explored for their effectiveness in various kinds of neurodegenerative disorders.

An Allosteric Inhibitory Potential of Triterpenes from Combretum racemosum on the Structural and Functional Dynamics of Plasmodium falciparum Lactate Dehydrogenase Binding Landscape

Multidrug resistance is a significant drawback in malaria treatment, and mutations in the active sites of the many critical antimalarial drug targets have remained challenging. Therefore, this has necessitated the global search for new drugs with new mechanisms of action. Plasmodium falciparum lactate dehydrogenase (pfLHD), a glycolytic enzyme, has emerged as a potential target for developing new drugs due to the parasite reliance on glycolysis for energy. Strong substrate-binding is required in pfLDH enzymatic catalysis; however, there is a lack of information on small molecules' inhibitory mechanism bound to the substrate-binding pocket. Therefore, this study investigated a potential allosteric inhibition of pfLDH by targeting the substrate-binding site. The structural and functional behaviour of madecassic acid (MA), the most promising among the six triterpenes bound to pfLDH, were unravelled using molecular dynamic simulations at 300 ns to gain insights into its mechanism of binding and inhibition and chloroquine as a standard drug. The docking studies identified that the substrate site has the preferred position for the compounds even in the absence of a co-factor. The bound ligands showed comparably higher binding affinity at the substrate site than at the co-factor site. Mechanistically, a characteristic loop implicated in the enzyme catalytic activity was identified at the substrate site. This loop accommodates key interacting residues (LYS174, MET175, LEU177 and LYS179) pivotal in the MA binding and inhibitory action. The MA-bound pfLHD average RMSD (1.60 Å) relative to chloroquine-bound pfLHD RMSD (2.00 Å) showed higher stability for the substrate pocket, explaining the higher binding affinity (-33.40 kcal/mol) observed in the energy calculations, indicating that MA exhibited profound inhibitory activity. The significant pfLDH loop conformational changes and the allostery substrate-binding landscape suggested inhibiting the enzyme function, which provides an avenue for designing antimalarial compounds in the future studies of pfLDH protein as a target.

In vitro and in silico analysis of a E559K mutation on cartilage oligomeric matrix protein

Pseudoachondroplasia (PSACH) is known as an autosomal dominant disorder associated with mutations in the gene of cartilage oligomeric matrix protein (COMP). The pathomolecular mechanisms of PSACH as a result of C-terminal globular region (CTD) mutations remain unclear. A heterozygous mutation (E559 K) in a Chinese family diagnosed with PSACH was reported in this study. To understand the pathogenesis of this mutation, we studied chondrogenic differentiation of patient menstrual blood-derived stem cells (MenSCs), and the impact of the mutation on structural changes of COMP was investigated using all-atom molecular dynamics simulation. The results suggested that the interactions with calcium and other molecules in the mutant structure were affected resulting in misfolding of the protein, which leads to ER stress and finally affects the survival of chondrocytes. The findings may promote the understanding of the pathomolecular mechanisms of PSACH, and possibly the development of drugs to treat the disease.

Determination of tyrosinase-cyanidin-3-O-glucoside and (-/+)-catechin binding modes reveal mechanistic differences in tyrosinase inhibition

Tyrosinase, exquisitely catalyzes the phenolic compounds into brown or black pigment, inhibition is used as a treatment for dermatological or neurodegenerative disorders. Natural products, such as cyanidin-3-O-glucoside and (-/+)-catechin, are considered safe and non-toxic food additives in tyrosinase inhibition but their ambiguous inhibitory mechanism against tyrosinase is still elusive. Thus, we presented the mechanistic insights into tyrosinase with cyanidin-3-O-glucoside and (-/+)-catechin using computational simulations and in vitro assessment. Initial molecular docking results predicted ideal docked poses (- 9.346 to - 5.795 kcal/mol) for tyrosinase with selected flavonoids. Furthermore, 100 ns molecular dynamics simulations and post-simulation analysis of docked poses established their stability and oxidation of flavonoids as substrate by tyrosinase. Particularly, metal chelation via catechol group linked with the free 3-OH group on the unconjugated dihydropyran heterocycle chain was elucidated to contribute to tyrosinase inhibition by (-/+)-catechin against cyanidin-3-O-glucoside. Also, predicted binding free energy using molecular mechanics/generalized Born surface area for each docked pose was consistent with in vitro enzyme inhibition for both mushroom and murine tyrosinases. Conclusively, (-/+)-catechin was observed for substantial tyrosinase inhibition and advocated for further investigation for drug development against tyrosinase-associated diseases.

Mechanistic Insights into the Selective Dual BET and PLK1 Inhibitory Activity of a Novel Benzamide Compound in Castration-Resistant Prostrate Cancer

Though multifactorial, BET and PLK1 proteins have been found to be key players in the oncogenic process leading to castration-resistant prostate cancer through regulation of AR and MYC-mediated transcription. Hence, dual inhibition of these proteins appears to be an auspicious approach for CRPC therapy. WNY0824 has been reported to exhibit nanomolar range inhibition as well as significant anti-proliferative activity on AR-positive CRPC cells in vitro. However, structural, and mechanistic events associated with its dual inhibitory and anti-proliferative mechanisms remain unclear. Utilizing integrative computer-assisted atomistic techniques, analyses revealed that the dual-inhibitory activity of WNY0824 against BRD4 and PLK1 proteins is mediated by conserved residues present in the binding cavities of both proteins which are shown to elicit various strong intermolecular interactions and thus favour binding affinity. Also, binding orientation of the ligand at the protein binding cavities allowed for important hydrophobic interactions which resulted in high binding free energy of -42.50 kcal/mol and -51.64 kcal/mol towards BRD4 and PLK1, respectively. While van der Waals interactions are very important to ligand binding in BRD4-WNY complex, electrostatic interactions are pertinent to PLK1-WNY complex. Intriguingly, WNY0824 triggered conformational alterations in both proteins through increased structural instability, decreased structural compactness and mitigation in exposure of residues to solvent surface area. Consequently, critical interactions peculiar to the oncogenic activities of BRD4 and PLK1 were inhibited, a phenomenon that results in an antagonism of CRPC progression. The mechanistic insights presented in this report would further assist in the structure-based design of improved inhibitors useful in CRPC therapy.

Leveraging on Active Site Similarities; Identification of Potential Inhibitors of Zinc-Finger and UFSP domain Protein (ZUFSP)

BACKGROUND

ZUFSP (Zinc-finger and UFSP domain protein) is a novel representative member of the recently characterized seventh class of deubiquitinating enzymes (DUBs). Due to the roles DUBs play in genetic instability, they have become a major drug target in cancer and neurodegenerative diseases. ZUFSP, being a DUB enzyme has also been implicated in genetic stability. However, no lead compound has been developed to target ZUFSP.

OBJECTIVE/METHODS

Therefore, in this study, we used a combined drug repurposing, virtual screening and per-Residue Energy Decomposition (PRED) to identify ZUFSP inhibitors with therapeutic potential. 3-bromo-6-{[4-hydroxy-1-3(3-phenylbutanoyl)piperidin-4-yl]methyl}-4H,5H,6H,7H-thieno[2,3- C]pyridine-7-one (BHPTP) which is an inhibitor of USP7 was repurposed to target ZUFSP. The rationale behind this is based on the similarity of the active between USP7 and ZUFSP.

RESULTS

PRED of the binding between BHPTP and ZUFSP revealed Cys223, Arg408, Met410, Asn460, and Tyr465 as the crucial residues responsible for this interaction. The pharmacophoric moieties of BHPTP responsible for this binding along with other physiochemical properties were used as a filter to retrieve potential ligands. 799 compounds were retrieved, ZINC083241427, ZINC063648749, and ZINC063648753 were selected due to the binding energy they exhibited. Cheminformatics analysis revealed that the compounds possess high membrane permeability, however, BHPTP had a low membrane permeability. Furthermore, the compounds are drug like, having obeyed Lipinski's rule of five.

CONCLUSION

Taken together, findings from this study put ZINC083241427, ZINC063648749, and ZINC063648753 as potential ZUFSP inhibitor, however, more experimental validation is required to unravel the mechanism of actions of these compounds.

Molecular Dynamics Simulations Reveal the Modulated Mechanism of STING Conformation

Stimulator of interferon genes (STING), which is an integral ER-membrane protein, could induce an antiviral state and boost antitumor immunity. Recent experiments reported that different small molecules could modulate the conformation of the STING. However, the mechanism of small molecules modulating the conformation of STING is still unknown. To illustrate the conformational modulated mechanism of STING by small molecules at atomic level, we investigated the interactions between STING and the small molecules: cGAMP and diABZI with molecular dynamics (MD) simulations method. Interestingly, we found that the residues of STING in the binding pocket are more flexible in the monomers of STING than that in the dimer of STING. We also demonstrated that cGAMP and diABZI have a similar binding mode to STING monomers/dimer, and π-π stacking interactions play important roles for the agonists and STING. Our study proposed mechanistic insights into the STING conformation modulated by small molecules and we suggested that the special molecule (e. g. diABZI) could induce the conformational transition of STING from the "open" monomers to the "closed" dimer state. Our research may provide a clue for the development of cancer immunotherapy.

Development of Halogenated Pyrazolines as Selective Monoamine Oxidase-B Inhibitors: Deciphering via Molecular Dynamics Approach

Halogens have been reported to play a major role in the inhibition of monoamine oxidase (MAO), relating to diverse cognitive functions of the central nervous system. Pyrazoline/halogenated pyrazolines were investigated for their inhibitory activities against human monoamine oxidase-A and -B. Halogen substitutions on the phenyl ring located at the fifth position of pyrazoline showed potent MAO-B inhibition. Compound 3-(4-ethoxyphenyl)-5-(4-fluorophenyl)-4,5-dihydro-1H-pyrazole (EH7) showed the highest potency against MAO-B with an IC₅₀ value of 0.063 µM. The potencies against MAO-B were increased in the order of –F (in EH7) > –Cl (EH6) > –Br (EH8) > –H (EH1). The residual activities of most compounds for MAO-A were > 50% at 10 µM, except for EH7 and EH8 (IC₅₀ = 8.38 and 4.31 µM, respectively). EH7 showed the highest selectivity index (SI) value of 133.0 for MAO-B, followed by EH6 at > 55.8. EH7 was a reversible and competitive inhibitor of MAO-B in kinetic and reversibility experiments with a Ki value of 0.034 ± 0.0067 µM. The molecular dynamics study documented that EH7 had a good binding affinity and motional movement within the active site with high stability. It was observed by MM-PBSA that the chirality had little effect on the overall binding of EH7 to MAO-B. Thus, EH7 can be employed for the development of lead molecules for the treatment of various neurodegenerative disorders.

Dual-Knockout of Mutant Isocitrate Dehydrogenase 1 and 2 Subtypes Towards Glioma Therapy: Structural Mechanistic Insights on the Role of Vorasidenib

Recently, Vorasidenib (AG-881) has been reported as a therapeutic alternative that exerts potent dual inhibitory activity against mIDH1/2 towards the treatment of low-grade glioma. However, structural and dynamic events associated with its dual inhibition mechanism remain unclear. As such, we employ integrative computer-assisted atomistic techniques to provide thorough structural and dynamic insights. Our analysis proved that the dual-targeting ability of AG-881 is mediated by Val255/Val294 within the binding pockets of both mIDH1 and mIDH2 which are shown to elicit a strong intermolecular interaction, thus favoring binding affinity. The structural orientations of AG-881 within the respective hydrophobic pockets allowed favorable interactions with binding site residues which accounted for its high binding free energy of -28.69 kcal/mol and -19.89 kcal/mol towards mIDH1 and mIDH2, respectively. Interestingly, upon binding, AG-881 was found to trigger systemic alterations of mIDH1 and mIDH2 characterized by restricted residue flexibility and a reduction in exposure of residues to the solvent surface area. As a result of these structural alterations, crucial interactions of the mutant enzymes were inhibited, a phenomenon that results in a suppression of the production of oncogenic stimulator 2-HG. Findings therefore provide thorough structural and dynamic insights associated with the dual inhibitory activity of AG-881 towards glioma therapy.

Configurational Entropy of Folded Proteins and Its Importance for Intrinsically Disordered Proteins

Many pairwise additive force fields are in active use for intrinsically disordered proteins (IDPs) and regions (IDRs), some of which modify energetic terms to improve the description of IDPs/IDRs but are largely in disagreement with solution experiments for the disordered states. This work considers a new direction-the connection to configurational entropy-and how it might change the nature of our understanding of protein force field development to equally well encompass globular proteins, IDRs/IDPs, and disorder-to-order transitions. We have evaluated representative pairwise and many-body protein and water force fields against experimental data on representative IDPs and IDRs, a peptide that undergoes a disorder-to-order transition, for seven globular proteins ranging in size from 130 to 266 amino acids. We find that force fields with the largest statistical fluctuations consistent with the radius of gyration and universal Lindemann values for folded states simultaneously better describe IDPs and IDRs and disorder-to-order transitions. Hence, the crux of what a force field should exhibit to well describe IDRs/IDPs is not just the balance between protein and water energetics but the balance between energetic effects and configurational entropy of folded states of globular proteins.

Probing Protein-Protein Interactions and Druggable Site Identification: Mechanistic Binding Events between Ubiquitin and Zinc Finger with UFM1-Specific Peptidase Domain Protein (ZUFSP)

BACKGROUND

Deubiquitinating enzymes (DUBs) protein family have been implicated in some deregulated pathways involved in carcinogenesis such as cell cycle, gene expression, and DNA damage response (DDR). Zinc finger with UFM1-specific peptidase domain protein (ZUFSP) is one of the recently discovered members of the DUBs Objectives: To identify and cross validate the ZUFSP binding site using the bioinformatic tools including SiteMap & Metapocket respectively. To understand the molecular basis of complementary ZUFSP-Ub interaction and associated structural events using MD Simulation Methods: In this study, four binding pockets were predicted, characterized, and cross-validated based on physiochemical features such as site score, druggability score, site volume, and site size. Also, Molecular dynamics simulation technique was employed to determine the impact of ubiquitin-binding on ZUFSP Results: Site 1 with a site score 1.065, Size 102, D scores 1.00, and size volume 261 was predicted to be the most druggable site. Structural studies revealed that upon ubiquitin-binding, the motional movement of ZUFSP was reduced when compared to the unbound ZUFSP. Also, the ZUFSP helical arm (ZHA) domain orient in such a way that it moves closer to the Ub, this orientation enables the formation of a UBD which is very peculiar to ZUFSP.

CONCLUSION

The impact of ubiquitin on ZUFSP movement and the characterization of its predicted druggable site can be targeted in the development of therapeutics.

A probable means to an end: exploring P131 pharmacophoric scaffold to identify potential inhibitors of Cryptosporidium parvum inosine monophosphate dehydrogenase

Compound P131 has been established to inhibit Cryptosporidium parvum's inosine monophosphate dehydrogenase (CpIMPDH). Its inhibitory activity supersedes that of paromomycin, which is extensively used in treating cryptosporidiosis. Through the per-residue energy decomposition approach, crucial moieties of P131 were identified and subsequently adopted to create a pharmacophore model for virtual screening in the ZINC database. This search generated eight ADMET-compliant hits that were examined thoroughly to fit into the active site of CpIMPDH via molecular docking. Three compounds ZINC46542062, ZINC58646829, and ZINC89780094, with favorable docking scores of - 8.3 kcal/mol, - 8.2 kcal/mol, and - 7.5 kcal/mol, were selected. The potential inhibitory mechanism of these compounds was probed using molecular dynamics simulation and Molecular Mechanics Generalized Poisson Boltzmann Surface Area (MM/PBSA) analyses. Results revealed that one of the hits (ZINC46542062) exhibited a lower binding free energy of - 39.52 kcal/mol than P131, which had - 34.6 kcal/mol. Conformational perturbation induced by the binding of the identified hits to CpIMPDH was similar to P131, suggesting a similarity in inhibitory mechanisms. Also, in silico investigation of the properties of the hit compounds implied superior physicochemical properties with regards to their synthetic accessibility, lipophilicity, and number of hydrogen bond donors and acceptors in comparison with P131. ZINC46542062 was identified as a promising hit compound with the highest binding affinity to the target protein and favorable physicochemical and pharmacokinetic properties relative to P131. The identified compounds can serve as a basis for conducting further experimental investigations toward the development of anticryptosporidials, which can overcome the challenges of existing therapeutic options. Graphical abstract P131 and the identified compounds docked in the NAD⁺ binding site of Cryptosporidium parvum IMPDH.

Identification of potential SARS-CoV-2 inhibitors from South African medicinal plant extracts using molecular modelling approaches

The coronavirus is a group of viruses found in animals as well as humans and have been detected since the 1960s. However, a newly identified form, SARS-CoV-2, has triggered a recent pandemic of respiratory disease now called COVID-19. There is currently no specific antiviral drug for the treatment of this pandemic, with most treatment strategies focused on symptomatic management and supportive therapy. As such, several drug discovery efforts are ongoing for potent treatment agents, with medicinal plants gradually gaining prominence. Approximately 80% of the South African population use traditional medicines to meet their primary health care needs. The current study aimed to identify potential COVID-19 therapeutic agents from a list of 29 bioactive compounds isolated from commonly used South African medicinal plants using molecular docking and molecular dynamics. Molecular docking identified arabic acid from Acacia senegal and L-canavanine found in Sutherlandia frutescens as a potential inhibitor of SARS-CoV-2 3C-like main protease. Similarly, hypoxoside isolated from Hypoxis hemerocallidea and uzarin from Xysmalobium undulatum, were identified as a potential inhibitor of SARS-CoV-2 receptor binding domain and SARS-CoV-2 RNA-dependent polymerase. These four bioactive compounds exhibited favourable binding orientations characterized by strong molecular interactions within respective inhibitors binding pockets of the target enzymes. Molecular dynamics simulations revealed that the binding of the identified inhibitors are characterized by structural perturbations which favour the inhibitory potency of these bioactive compounds. Additionally, in silico pharmacokinetic assessment of the compounds demonstrated favourable anti-SARS-CoV-2 properties. Although not conclusive, further experimental exploration of these compounds could serve as a starting point for the discovery of novel SARS-CoV-2 therapeutic.

'Polymorphism-aided' Selective Targeting and Inhibition of Caspase-6 by a Novel Allosteric Inhibitor Towards Efficient Alzheimer's Disease Treatment

The predominance of Alzheimer's disease (AD) among the aged remains a global challenge. As such, the search for alternative and effective therapeutic options continuous unabated. Among the therapeutic targets explored over the years toward impeding the progression of AD is caspase-6 (Casp6), although selectively targeting Casp6 remains a challenge due to high homology with other members of the caspase family. Methyl 3-[(2,3-dihydro-1-benzofuran-2-yl formamido) methyl]-5-(furan-2-amido) benzoate (C13), a novel allosteric inhibitor, is reportedly shown to exhibit selective inhibition against mutant human Casp6 variants (E35K). However, structural and atomistic insights accounting for the reported inhibitory prowess of C13 remains unresolved. In this study, we seek to unravel the mechanistic selectivity of C13 coupled with the complementary effects of E35K single-nucleotide polymorphism (SNP) relative to Casp6 inhibition. Analyses of binding dynamics revealed that the variant Lysine-35 mediated consistent high-affinity interactions with C13 at the allosteric site, possibly forming the molecular basis of the selectivity of C13 as well as its high binding free energy as estimated. Analysis of residue interaction network around Glu35 and Lys35 revealed prominent residue network distortions in the mutant Casp6 conformation evidenced by a decrease in node degree, reduced number of edges and an increase short in path length relative to a more compact conformation in the wild system. The relatively higher binding free energy of C13 coupled with the stronger intermolecular interactions elicited in the mutant conformation further suggests that the mutation E35K probably favours the inhibitory activity of C13. Further analysis of atomistic changes showed increased C-α atom deviations consistent with structural disorientations in the mutant Casp6. Structural Insights provided could open up a novel paradigm of structure-based design of selective allosteric inhibition of Casp6 towards the treatment of neurodegenerative diseases.

In Silico Repurposing of J147 for Neonatal Encephalopathy Treatment: Exploring Molecular Mechanisms of Mutant Mitochondrial ATP Synthase

BACKGROUND

Neonatal Encephalopathy (NE) is a mitochondrial ATP synthase (mATPase) disease, which results in the death of infants. The case presented here is reportedly caused by complex V deficiency as a result of mutation of Arginine to Cysteine at residue 329 in the mATPase. A recent breakthrough was the discovery of J147, which targets mATPase in the treatment of Alzheimer's disease. Based on the concepts of computational target-based drug design, this study investigated the possibility of employing J147 as a viable candidate in the treatment of NE.

OBJECTIVE/METHODS

The structural dynamic implications of this drug on the mutated enzyme are yet to be elucidated. Hence, integrative molecular dynamics simulations and thermodynamic calculations were employed to investigate the activity of J147 on the mutated enzyme in comparison to its already established inhibitory activity on the wild-type enzyme.

RESULTS

A correlated structural trend occurred between the wild-type and mutant systems whereby all the systems exhibited an overall conformational transition. Equal observations in favorable free binding energies further substantiated uniformity in the mobility, and residual fluctuation of the wild-type and mutant systems. The similarity in the binding landscape suggests that J147 could as well modulate mutant mATPase activity in addition to causing structural modifications in the wild-type enzyme.

CONCLUSION

Findings suggest that J147 can stabilize the mutant protein and restore it to a similar structural state as the wild-type which depicts functionality. These details could be employed in drug design for potential drug resistance cases due to mATPase mutations that may present in the future.

Identifying new liver X receptor alpha modulators and distinguishing between agonists and antagonists by crystal ligand pocket screening

Background: Modulators of LXRα are of high pharmacological interest as LXRα regulates fatty acid metabolism, inflammatory processes and cancer. We aim to identify new LXRα modulators and to recognize a distinguishable feature of agonists. Results&methodology: The ligand self-dock and largest-cavity-size searching purposely located two appropriate ligand-binding sites to reach the two aims. One is identifying the new modulators from Maybridge library. 20 new compounds are confirmed by the in vitro reporter gene assay. The other is denoting an agonist by at least one best docking pose having one hydrogen bond to LXRα Helix12 His421. Conclusion: Based on the quality x-ray binding pocket, we can identify new LXRα modulators and distinguish between agonists and antagonists by molecular docking.

The Dual-Targeting Activity of the Metabolite Substrate of Para-amino Salicyclic Acid in the Mycobacterial Folate Pathway: Atomistic and Structural Perspectives

Therapeutic targeting of folate biosynthetic pathway has recently been explored as a viable strategy in the treatment of tuberculosis. The bioactive metabolite substrate of Para-amino salicyclic acid (PAS-M) reportedly dual-targets dihydrofolate reductase (DHFR) and flavin-dependent thymidylate synthase (FDTS), two essential enzymes in folate biosynthetic pathway. However, the molecular mechanisms and structural dynamics of this dual inhibitory activity of the PAS-M remain elusive. Molecular dynamics simulations revealed that binding of PAS-M towards DHFR is characterized by a recurrence of strong conventional hydrogen bond interactions between a peculiar DHFR binding site residue (Asp27) and the 2-amino-decahydropteridin-4-ol group of PAS-M. Similarly, the binding of PAS-M towards FDTS also involved consistent strong conventional hydrogen bond interactions between some specific residues (Tyr101, Arg172, Thr4, Gln103, Arg87 and Gln106) and, the 2-amino-decahydropteridin-4-ol group, thus establishing the cruciality of the group. Structural dynamics of the bound complexes of both enzymes revealed that, upon binding, PAS-M is anchored at the entrance of hydrophobic pockets by strong hydrogen bond interactions while the rest of the structure gains access to deeper hydrophobic residues to engage in favorable interactions. Further analysis of atomistic changes of both enzymes showed increased C-α atom deviations as well as an increase C-α atoms radius of gyration consistent with structural disorientations. These conformational changes possibly interfered with the biological functions of the enzymes and hence their inhibition as experimentally reported. Structural Insights provided could open up a novel paradigm of structure-based design of multi-targeting inhibitors of biological targets in the folate biosynthetic pathway toward tuberculosis therapy.

Structural responses of kiwifruit allergen Act d 2 to thermal and electric field stresses based on molecular dynamics simulations and experiments

Kiwifruit is considered to be the most common plant-based food causing allergic reactions, after peanuts, soybeans, and wheat.

Tweaking α -Galactoceramides: Probing the Dynamical Mechanisms of Improved Recognition for Invariant Natural Killer T-cell Receptor in Cancer Immunotherapeutics

BACKGROUND

The last few decades have witnessed groundbreaking research geared towards immune surveillance mechanisms and have yielded significant improvements in the field of cancer immunotherapy. This approach narrows down on the development of therapeutic agents that either activate or enhance the recognitive function of the immune system to facilitate the destruction of malignant cells. The α -galactosylceramide derivative, KRN7000, is an immunotherapeutic agent that has gained attention due to its pharmacological ability to activate CD1d-restricted invariant natural killer T(iNKT) cells with notable potency against cancer cells in mouse models; a therapeutic success was not well replicated in human models. Dual structural modification of KRN7000 entailing the incorporation of hydrocinnamoyl ester on C6" and C4-OH truncation of the sphingoid base led to the development of AH10-7 which, interestingly, exhibited high potency in human cells.

OBJECTIVE/METHODS

Therefore, to gain molecular insights into the structural dynamics and selective mechanisms of AH10-7 for human variants, we employed integrative molecular dynamics simulations and thermodynamic calculations to investigate the inhibitory activities of KRN7000 andAH10-7 on hTCR-CD1d towards activating iNKT.

RESULTS

Interestingly, our findings revealed that AH10-7 exhibited higher affinity binding and structural effects on hTCR-CD1d, as mediated by the incorporated hydrocinnamoyl ester moiety which accounted for stronger intermolecular interactions with 'non-common' binding site residues.

CONCLUSION

Findings extracted from this study further reveal important molecular and structural perspectives that could aid in the design of novel α-GalCer derivatives for cancer immunotherapeutics.

Broadening the horizon: Integrative pharmacophore-based and cheminformatics screening of novel chemical modulators of mitochondria ATP synthase towards interventive Alzheimer's disease therapy

The proven efficacy of J147 in the treatment of Alzheimer's disease (AD) has been emphatic, particularly since its selective modulatory roles towards mitochondrial ATP synthase (mATPase) were defined. This prospect, if methodically probed, could further pave way for the discovery of novel anti-AD drugs with improved pharmacokinetics and therapeutic potential. To this effect, for the first time, we employed a four-step paradigm that integrated our in-house per-residue energy decomposition (PRED) protocol coupled with molecular dynamics, cheminformatics and analytical binding free energy methods. This was geared towards the screening and identification of new leads that exhibit modulatory potentials towards mATPase in a J147-similar pattern. Interestingly, from a large-scale library of compounds, we funnelled down on three potential hits that demonstrated selective and high-affinity binding activities towards α-F1-ATP synthase (ATP5A) relative to J147. Moreover, these compounds exhibited higher binding propensity with a differential ΔGs greater than -1 kcal/mol comparative to J147, and also elicited distinct modulatory effects on ATP5A domain structures. More interestingly, per-residue pharmacophore modeling of these lead compounds revealed similar interactive patterns with crucial residues at the α-site region of ATP5A characterized by high energy contributions based on binding complementarity. Recurrent target residues involved in high-affinity interactions with the hit molecules relative to J147 include Arg1112 and Gln426. Furthermore, assessments of pharmacokinetics revealed that the lead compounds were highly drug-like with minimal violations of the Lipinski's rule of five. As developed in this study, the most extrapolative pharmacophore model of the selected hits encompassed three electron donors and one electron acceptor. We speculate that these findings will be fundamental to the reformation of anti-AD drug discovery procedures.