Molecular docking study of flavonoid compounds for possible matrix metalloproteinase-13 inhibition

Exploring molecular targets: herbal isolates in cervical cancer therapy

OBJECTIVE

Cervical cancer (CxCa) stands as a significant global health challenge, ranking fourth in cancer-related mortality among the female population. While chemotherapy regimens have demonstrated incremental progress in extending overall survival, the outlook for recurrent CxCa patients remains disheartening. An imperative necessity arises to delve into innovative therapeutic avenues, with molecular targeted therapy emerging as a promising candidate. Previous investigations have shed light on the therapeutic effectiveness of five distinct herbal compounds, epicatechin, curcumin, myricetin, jatrorrhizine, and arborinine, within the context of CxCa.

METHODS

A systems biology approach was employed to discern differentially expressed genes (DEGs) in CxCa tissues relative to healthy cervical epithelial tissues. A protein-protein interaction network (PPIN) was constructed, anchored in the genes related to CxCa. The central genes were discerned within the PPIN, and Kaplan-Meier survival curves explored their prognostic significance. An assessment of the binding affinity of the selected herbal compounds to the master regulator of prognostic markers in CxCa was conducted.

RESULTS

A significant correlation between the overexpression of MYC, IL6, JUN, RRM2, and VEGFA and an adverse prognosis in CxCa was indicated. The regulation of these markers is notably influenced by the transcription factor CEBPD. Molecular docking analysis indicated that the binding affinity between myricetin and the CEBPD DNA binding site was robust.

CONCLUSION

The findings presented herein have unveiled pivotal genes and pathways that play a central role in the malignant transformation of CxCa. CEBPD has emerged as a potential target for harnessing the therapeutic potential of myricetin in this context.

Potential JAK2 Inhibitors from Selected Natural Compounds: A Promising Approach for Complementary Therapy in Cancer Patients

Background

Janus-activated kinase 2 (JAK2) plays a pivotal role in numerous essential biological processes, including proliferation, apoptosis, and metastasis in human cells. Prior studies have indicated that inhibiting JAK2 could be a promising strategy to mitigate cell proliferation and induce apoptosis in tumor cells.

Objectives

This study aimed to estimate the binding affinity of 79 herbal compounds, comprising 46 flavonoids, 21 anthraquinones, and 12 cinnamic acids, to the ATP-binding cleft of JAK2 to identify potential herbal inhibitors of JAK2.

Methods

The binding affinities between ligands and JAK2 were calculated utilizing AutoDock 4.0 software in conjunction with the Cygwin environment. Cross-validation was conducted using the Schrödinger tool. Molecular dynamics simulations were employed to evaluate the stability of docked poses for the most significant JAK2 inhibitors. Furthermore, the Discovery Studio Visualizer tool was utilized to elucidate interactions between the top-ranked JAK2 inhibitors and residues within the JAK2 ATP-binding site.

Results

Twelve flavonoids, two anthraquinones, and three cinnamic acids demonstrated substantial binding affinities to the protein kinase domain of the receptor, with a criterion of ΔGbinding < -10 kcal/mol. Among the studied flavonoids, anthraquinones, and cinnamic acid derivatives, orientin, chlorogenic acid, and pulmatin emerged as the most potent JAK2 inhibitors, exhibiting ΔGbinding scores of -14.49, -11.87, and -10.76 kcal/mol, respectively. Furthermore, the docked poses of orientin, pulmatin, and chlorogenic acid remained stable throughout 60 ns computer simulations. The average root mean square deviation values calculated for JAK2 when complexed with orientin, chlorogenic acid, and pulmatin were 2.04 Å, 2.06 Å, and 1.95 Å, respectively.

Conclusion

This study underscores the robust inhibitory potential of orientin, pulmatin, and chlorogenic acid against JAK2. The findings hold promise for the development of novel and effective drugs for cancer treatment.

BACE1 Inhibition Utilizing Organic Compounds Holds Promise as a Potential Treatment for Alzheimer's and Parkinson's Diseases

Background

Neurological disorders like Alzheimer's disease (AD) and Parkinson's disease (PD) manifest through gradually deteriorating cognitive functions. An encouraging strategy for addressing these disorders involves the inhibition of precursor-cleaving enzyme 1 (BACE1).

Objectives

In the current research, a virtual screening technique was employed to identify potential BACE1 inhibitors among selected herbal isolates.

Methods

This study evaluated 79 flavonoids, anthraquinones (AQs), and cinnamic acid derivatives for their potential blood-brain barrier (BBB) permeability. Using the AutoDock 4.0 tool, molecular docking analysis was conducted to determine the binding affinity of BBB permeable compounds to the BACE1 active site. Molecular dynamics (MD) simulations were performed to assess the stability of the docked poses of the most potent inhibitors. The interactions between the most effective plant-based inhibitors and the residues within the BACE1 catalytic site were examined before and after MD simulations.

Results

Ponciretin, danthron, chrysophanol, and N-p-coumaroyltyramine were among the highest-ranking BACE1 inhibitors, with inhibition constant values calculated in the nanomolar range. Furthermore, during 10 ns simulations, the docked poses of these ligands were observed to be stable.

Conclusion

The findings propose that ponciretin, danthron, chrysophanol, and N-p-coumaroyltyramine might serve as potential choices for the treatment of AD and PD, laying the groundwork for the creation of innovative BACE1 inhibitors.

Mechanisms Involved in Therapeutic Effects of Scutellaria baicalensis Georgi in Oral Squamous Cell Carcinoma Based on Systems Biology and Structural Bioinformatics Approaches

Objective

Oral squamous cell carcinoma (OSCC) is the most frequent oral cancer, constituting more than 90% of all oral carcinomas. The 5-year survival rate of OSCC patients is not satisfactory, and therefore, there is an urgent need for new practical therapeutic approaches besides the current therapies to overcome OSCC. Scutellaria baicalensis Georgi (SBG) is a plant of the family Lamiaceae with several pharmaceutical properties such as antioxidant, anti-inflammatory, and anticancer effects. Previous studies have demonstrated the curative effects of SBG in OSCC.

Methods

A systems biology approach was conducted to identify differentially expressed miRNAs (DEMs) in OSCC patients with a dismal prognosis compared to OSCC patients with a favorable prognosis. A protein interaction map (PIM) was built based on DEMs targets, and the hub genes within the PIM were indicated. Subsequently, the prognostic role of the hubs was studied using Kaplan-Meier curves. Next, the binding affinity of SBG's main components, including baicalein, wogonin, oroxylin-A, salvigenin, and norwogonin, to the prognostic markers in OSCC was evaluated using molecular docking analysis.

Results

Survival analysis showed that overexpression of CAV1, SERPINE1, ACTB, SMAD3, HMGA2, MYC, EIF2S1, HSPA4, HSPA5, and IL6 was significantly related to a poor prognosis in OSCC. Besides, molecular docking analysis demonstrated the ΔGbinding and inhibition constant values between SBG's main components and SERPINE1, ACTB, HMGA2, EIF2S1, HSPA4, and HSPA5 were as <-8.00 kcal/mol and nanomolar concentration, respectively. The most salient binding affinity was observed between wogonin and SERPINE1 with a criterion of ΔGbinding < -10.02 kcal/mol.

Conclusion

The present results unraveled potential mechanisms involved in therapeutic effects of SBG in OSCC based on systems biology and structural bioinformatics analyses.

Herbal Metabolites as Potential Carbonic Anhydrase Inhibitors: Promising Compounds for Cancer and Metabolic Disorders

Background

Human carbonic anhydrases (CAs) play a role in various pathological mechanisms by controlling intracellular and extracellular pH balance. Irregular expression and function of CAs have been associated with multiple human diseases, such as obesity, cancer, glaucoma, and epilepsy. In this work, we identify herbal compounds that are potential inhibitors of CA VI.

Methods

We used the AutoDock tool to evaluate binding affinity between the CA VI active site and 79 metabolites derived from flavonoids, anthraquinones, or cinnamic acids. Compounds ranked at the top were chosen for molecular dynamics (MD) simulations. Interactions between the best CA VI inhibitors and residues within the CA VI active site were examined before and after MD analysis. Additionally, the effects of the most potent CA VI inhibitor on cell viability were ascertained in vitro through the 2,5-diphenyl-2H-tetrazolium bromide (MTT) assay.

Results

Kaempferol 3-rutinoside-4-glucoside, orientin, kaempferol 3-rutinoside-7-sophoroside, cynarin, and chlorogenic acid were estimated to establish binding with the CA VI catalytic domain at the picomolar scale. The range of root mean square deviations for CA VI complexes with kaempferol 3-rutinoside-4-glucoside, aloe-emodin 8-glucoside, and cynarin was 1.37 to 2.05, 1.25 to 1.85, and 1.07 to 1.54 Å, respectively. The MTT assay results demonstrated that cynarin had a substantial effect on HCT-116 cell viability.

Conclusion

This study identified several herbal compounds that could be potential drug candidates for inhibiting CA VI.

Cinnamic acids as promising bioactive compounds for cancer therapy by targeting MAPK3: a computational simulation study

OBJECTIVES

Mitogen-activated protein kinase-3 (MAPK3) is the upstream regulator in the MAPK cascade and is involved in many critical signaling pathways and biological processes, such as cell proliferation, survival, and apoptosis. MAPK3 overexpression is linked to onset, development, metastasis, and drug resistance in several human cancers. Thus, identifying novel and effective MAPK3 inhibitors is highly demanded. Herein, we aimed to discover organic compounds from cinnamic acid derivatives as potential MAPK3 inhibitors.

METHODS

The binding affinity of 20 cinnamic acids to the MAPK3 active site was tested using the AutoDock 4.0 software. Top-ranked cinnamic acids were ranked based on the ΔG binding values between the ligands and the receptor's active site. Interaction modes between top-ranked cinnamic acids and MAPK3 catalytic site were indicated using the Discovery Studio Visualizer tool. Molecular dynamics (MD) simulation was carried out to study the stability of the docked pose for the most potent MAPK3 inhibitor in this study.

RESULTS

Cynarin, chlorogenic acid, rosmarinic acid, caffeic acid 3-glucoside, and cinnamyl caffeate exhibited a salient binding affinity to the MAPK3 active site with the criteria of ΔG binding <-10 k cal/mol. Further, the inhibition constant value for cynarin was calculated at the picomolar concentration. The docked pose of cynarin within the MAPK3 catalytic domain was stable in 100 ns simulation.

CONCLUSIONS

Cynarin, chlorogenic acid, rosmarinic acid, caffeic acid 3-glucoside, and cinnamyl caffeate might be helpful in cancer therapy by inhibiting MAPK3.

Cinnamic acid derivatives as potential matrix metalloproteinase-9 inhibitors: molecular docking and dynamics simulations

Matrix metalloproteinase-9 (MMP-9) is a zinc and calcium-dependent proteolytic enzyme involved in extracellular matrix degradation. Overexpression of MMP-9 has been confirmed in several disorders, including cancers, Alzheimer′s disease, autoimmune diseases, cardiovascular diseases, and dental caries. Therefore, MMP-9 inhibition is recommended as a therapeutic strategy for combating various diseases. Cinnamic acid derivatives have shown therapeutic effects in different cancers, Alzheimer′s disease, cardiovascular diseases, and dental caries. A computational drug discovery approach was performed to evaluate the binding affinity of selected cinnamic acid derivatives to the MMP-9 active site. The stability of docked poses for top-ranked compounds was also examined. Twelve herbal cinnamic acid derivatives were tested for possible MMP-9 inhibition using the AutoDock 4.0 tool. The stability of the docked poses for the most potent MMP-9 inhibitors was assessed by molecular dynamics (MD) in 10 nanosecond simulations. Interactions between the best MMP-9 inhibitors in this study and residues incorporated in the MMP-9 active site were studied before and after MD simulations. Cynarin, chlorogenic acid, and rosmarinic acid revealed a considerable binding affinity to the MMP-9 catalytic domain (ΔGbinding < -10 kcal/mol). The inhibition constant value for cynarin and chlorogenic acid were calculated at the picomolar scale and assigned as the most potent MMP-9 inhibitor from the cinnamic acid derivatives. The root-mean-square deviations for cynarin and chlorogenic acid were below 2 Å in the 10 ns simulation. Cynarin, chlorogenic acid, and rosmarinic acid might be considered drug candidates for MMP-9 inhibition.

Prognostic biomarkers and molecular pathways mediating Helicobacter pylori–induced gastric cancer: a network-biology approach

Cancer of the stomach is the second most frequent cancer-related death worldwide. The survival rate of patients with gastric cancer (GC) remains fragile. There is a requirement to discover biomarkers for prognosis approaches. Helicobacter pylori in the stomach is closely associated with the progression of GC. We identified the genes associated with poor/favorable prognosis in H. pylori–induced GC. Multivariate statistical analysis was applied on the Gene Expression Omnibus (GEO) dataset GSE54397 to identify differentially expressed miRNAs (DEMs) in gastric tissues with H. pylori–induced cancer compared with the H. pylori–positive with non-cancerous tissue. A protein interaction map (PIM) was built and subjected to DEMs targets. The enriched pathways and biological processes within the PIM were identified based on substantial clusters. Thereafter, the most critical genes in the PIM were illustrated, and their prognostic impact in GC was investigated. Considering p-value less than 0.01 and |Log2 fold change| as >1, five microRNAs demonstrated significant changes among the two groups. Gene functional analysis revealed that the ubiquitination system, neddylation pathway, and ciliary process are primarily involved in H. pylori–induced GC. Survival analysis illustrated that the overexpression of DOCK4, GNAS, CTGF, TGF-b1, ESR1, SELE, TIMP3, SMARCE1, and TXNIP was associated with poor prognosis, while increased MRPS5 expression was related to a favorable prognosis in GC patients. DOCK4, GNAS, CTGF, TGF-b1, ESR1, SELE, TIMP3, SMARCE1, TXNIP, and MRPS5 may be considered prognostic biomarkers for H. pylori–induced GC. However, experimental validation is necessary in the future.

MicroRNA-Based Markers of Oral Tongue Squamous Cell Carcinoma and Buccal Squamous Cell Carcinoma: A Systems Biology Approach

Objective

Oral tongue squamous cell carcinoma (OTSCC) and buccal squamous cell carcinoma (BSCC) are the first and second leading causes of oral cancer, respectively. OTSCC and BSCC are associated with poor prognosis in patients with oral cancer. Thus, we aimed to indicate signaling pathways, Gene Ontology terms, and prognostic markers mediating the malignant transformation of the normal oral tissue to OTSCC and BSCC.

Methods

The dataset GSE168227 was downloaded and reanalyzed from the GEO database. Orthogonal partial least square (OPLS) analysis identified common differentially expressed miRNAs (DEMs) in OTSCC and BSCC compared to their adjacent normal mucosa. Next, validated targets of DEMs were identified using the TarBase web server. With the use of the STRING database, a protein interaction map (PIM) was created. Using the Cytoscape program, hub genes and clusters within the PIM were shown. Next, gene-set enrichment analysis was carried out using the g:Profiler tool. Using the GEPIA2 web tool, analyses of gene expression and survival analysis were also performed.

Results

Two DEMs, including has-miR-136 and has-miR-377, were common in OTSCC and BSCC (p value <0.01; |Log2 FC| > 1). A total of 976 targets were indicated for common DEMs. PIM included 96 hubs, and the upregulation of EIF2S1, CAV1, RAN, ANXA5, CYCS, CFL1, MYC, HSP90AA1, PKM, and HSPA5 was significantly associated with a poor prognosis in the head and neck squamous cell carcinoma (HNSCC), while NTRK2, HNRNPH1, DDX17, and WDR82 overexpression was significantly linked to favorable prognosis in the patients with HNSCC. "Clathrin-mediated endocytosis" was considerably dysregulated in OTSCC and BSCC.

Conclusion

The present study suggests that has-miR-136 and has-miR-377 are underexpressed in OTSCC and BSCC than in normal oral mucosa. Moreover, EIF2S1, CAV1, RAN, ANXA5, CYCS, CFL1, MYC, HSP90AA1, PKM, HSPA5, NTRK2, HNRNPH1, DDX17, and WDR82 demonstrated prognostic markers in HNSCC. These findings may benefit the prognosis and management of individuals with OTSCC/BSCC. However, additional experimental verification is required.

Flavonoids as Strong Inhibitors of MAPK3: A Computational Drug Discovery Approach

Background

Mitogen-activated protein kinase 3 (MAPK3) mediates the onset, progression, metastasis, drug resistance, and poor prognosis in various malignancies, including glioma, liver, ovarian, thyroid, lung, breast, gastric, and oral cancers. Negative regulation of MAPK3 expression using miRNAs has led to therapeutic effects in cancer.

Objectives

The present study performed molecular docking and dynamics simulation to identify potential MAPK3 inhibitors from natural flavonoids, possibly leading to drug development in cancer therapy.

Methods

A computational drug discovery approach was performed using the AutoDock tool to identify potential MAPK3 inhibitors from 46 plant-based flavonoids. A cross-validation study was executed using the Schrödinger Maestro docking tool. Molecular dynamics (MD) was executed to evaluate the stability of docked poses between the top-ranked compounds and the MAPK3 catalytic domain. Interactions among the most potent MAPK3 inhibitors and residues within the receptor's active site were studied using the BIOVIA Discovery Studio Visualizer before and after 100 ns MD simulations.

Results

Kaempferol 3-rutinoside-4'-glucoside, kaempferol 3-rutinoside-7-sophoroside, rutin, and vicenin-2 exhibited a magnificent binding affinity to the receptor's active site. In addition, the stability of the docked poses of these compounds seemed to be stable after ∼45 ns computer simulations.

Conclusion

The present study suggests that kaempferol 3-rutinoside-4'-glucoside, kaempferol 3-rutinoside-7-sophoroside, rutin, and vicenin-2 could strongly bind to the MAPK3 catalytic site and could be assigned as a potent inhibitor for MAPK3. These findings may be helpful in the treatment of various cancers. However, further validation experiments are needed.

Anti-Aging, Anti-Acne, and Cytotoxic Activities of Houttuynia cordata Extracts and Phytochemicals Analysis by LC-MS/MS

Although many biological properties of Houttuynia cordata have been found, its anti-aging and anti-acne effects have not yet been investigated. This study was aimed to evaluate the in vitro anti-aging and anti-acne activities of H. cordata extracts and their cytotoxic activities and phytochemicals analyzed with liquid chromatography with tandem mass spectrometry (LC-MS/MS). Dried aerial parts of H. cordata were given different extractions. The aqueous and ethanolic extracts obtained were named HCA and HCE, respectively, and used to screen total phenolic and flavonoid contents. In vitro anti-aging, skin-related antimicrobial, scanning electron microscopy (SEM), in vitro cytotoxic, and LC-MS/MS analyses were performed. The total phenolic contents of the HCA and HCE were 5.11 ± 0.25 and 27.02 ± 1.07 mg gallic acid equivalent (GAE)/g dry extract while their total flavonoid contents were 104.94 ± 5.16 and 571.86 ± 2.86 mg quercetin equivalent (QE)/g dry extract, respectively. The HCA and HCE inhibited the activities of collagenase (28.33–46.00%), elastase (30.00–34.33%), and hyaluronidase (93.87–98.72%). The minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC) of the HCA against Cutibacterium acnes DMST14916 were 5.77 and 5.77 mg/mL while those of the HCE were 2.47 and 2.47 mg/mL, respectively. Cell collapses of C. acnes after treatment with the extracts were observed with SEM. The HCE was not toxic to macrophages, keratinocytes, and fibroblasts up to 400 mg/mL. The HCA showed toxicity against macrophages at 62.5 mg/mL and both skin cells at 250 mg/mL. The main phytochemicals in the extracts were identified with LC-MS/MS. Phenolic compounds, flavonoids, and flavonoid derivatives in H. cordata extracts could be major phytochemicals to possess a broad spectrum of biological activities including antioxidant, antimicrobial, and anti-aging activities. The findings from this study showed that the HCE has potential anti-aging and anti-acne properties while having non-cytotoxic activities on the immune and skin cells. These results indicate that the extract is probably advantageous in the development of skincare cosmeceutics and beauty treatments.

Potential Biomarkers and Signaling Pathways Associated with the Pathogenesis of Primary Ameloblastoma: A Systems Biology Approach

Objective

Ameloblastoma is a benign odontogenic tumor that may lead to ameloblastic carcinoma. This study aimed to determine potential signaling pathways and biological processes, critical genes and their regulating transcription factors (TFs), and miRNAs, as well as protein kinases involved in the etiology of primary ameloblastoma.

Methods

The dataset GSE132472 was obtained from the GEO database, and multivariate statistical analyses were applied to identify differentially expressed genes (DEGs) in primary ameloblastoma tissues compared to the corresponding normal gingiva samples. A protein-protein interaction (PPI) map was built using the STRING database. The Cytoscape software identified significant modules and the hub genes within the PPI network. Gene Ontology annotation and signaling pathway analyses were executed by employing the DAVID and Reactome databases, respectively. Significant TFs and miRNAs acting on the hub genes were identified using the iRegulon plugin and MiRWalk 2.0 database, respectively. A protein kinase enrichment analysis was conducted using the online Kinase Enrichment Analysis 2 (KEA2) web server. The approved drugs acting on the hub genes were also found.

Results

A total of 1,629 genes were differentially expressed in primary ameloblastoma (P value <0.01 and |Log2FC| > 1). HRAS, CDK1, MAPK3, ERBB2, COL1A1, CYCS, and BRCA1 demonstrated high degree and betweenness centralities in the PPI network. E2F4 was the most significant TF acting on the hub genes. BTK was the protein kinase significantly enriched by the TFs. Cholesterol biosynthesis was considerably involved in primary ameloblastoma.

Conclusions

This study provides an intuition into the potential mechanisms involved in the etiology of ameloblastoma.

Molecular Docking, Tyrosinase, Collagenase, and Elastase Inhibition Activities of Argan By-Products

The argan tree (Argania spinosa (L.) Skeels) is one of the most important floristic resources in Morocco. This Moroccan endemic tree is known for its numerous therapeutic and medicinal uses. In addition to some medicinal and cosmetic uses, argan fruit pulp and press cake are traditionally used by the Berber population for heating and feeding livestock. Molecular docking is an in silico approach that predicts the interaction between a ligand and a protein. This approach is mainly used in chemistry and pharmacology of natural products as a prediction tool with the purpose of selecting plant extracts or fractions for in vitro tests. The aim of this research is to study the evaluation of potential tyrosinase, collagenase, and elastase inhibitory activities of argan fruit press-cake and pulp extracts. Extracts were evaluated for their total phenolic content (TPC), and the major polyphenols of both press-cake and pulp extracts were submitted to molecular docking in order to determine the mechanisms of action of these compounds. Obtained results revealed that fruit pulp had the strongest dermocosmetic activities, as well as the highest TPC, with values above 55 mg gallic-acid equivalent per gram of dry matter (mgeq AG/gDM). Moreover, those results were positively correlated with the docking findings, suggesting that the pulp lead compounds have higher affinity with tyrosinase, collagenase, and elastase action sites. The results here presented are very promising and open new perspectives for the exploitation of argan-tree by-products as cosmetic agents towards the development of new anti-aging products.

Computational drug repurposing study of antiviral drugs against main protease, RNA polymerase, and spike proteins of SARS-CoV-2 using molecular docking method

OBJECTIVES

The new Coronavirus (SARS-CoV-2) created a pandemic in the world in late 2019 and early 2020. Unfortunately, despite the increasing prevalence of the disease, there is no effective drug for the treatment. A computational drug repurposing study would be an appropriate and rapid way to provide an effective drug in the treatment of the coronavirus disease of 2019 (COVID-19) pandemic. In this study, the inhibitory potential of more than 50 antiviral drugs on three important proteins of SARS-CoV-2, was investigated using the molecular docking method.

METHODS

By literature review, three important proteins, including main protease, RNA-dependent RNA polymerase (RdRp), and spike, were selected as the drug targets. The three-dimensional (3D) structure of protease, spike, and RdRp proteins was obtained from the Protein Data Bank. Proteins were energy minimized. More than 50 antiviral drugs were considered as candidates for protein inhibition, and their 3D structure was obtained from Drug Bank. Molecular docking settings were defined using Autodock 4.2 software and the algorithm was executed.

RESULTS

Based on the estimated binding energy of docking and hydrogen bond analysis and the position of drug binding, five drugs including, indinavir, lopinavir, saquinavir, nelfinavir, and remdesivir, had the highest inhibitory potential for all three proteins.

CONCLUSIONS

According to the results, among the mentioned drugs, saquinavir and lopinavir showed the highest inhibitory potential for all three proteins compared to the other drugs. This study suggests that saquinavir and lopinavir could be included in the laboratory phase studies as a two-drug treatment for SARS-CoV-2 inhibition.