Abstract
Objective:
To promote the development and therapeutic application of new medications, it is crucial to conduct a thorough investigation into the mechanism by which the traditional Chinese herb pair of Haizao-Kunbu (HK) treats Graves’ disease (GD).
Materials and methods:
Chemical ingredients of HK, putative target genes, and GD-associated genes were retrieved from online public databases. Using Cytoscape 3.9.1, a compound-gene target network was established to explore the association between prosperous ingredients and targets. STRING, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes pathway analyses visualized core targets and disease pathways. Additionally, we conducted a refined analysis of the binding interactions between active ingredients and their respective targets. To visualize these findings, we employed precise molecular docking techniques. Furthermore, we carried out molecular dynamics simulations to gain insights into the formation of more tightly bound complexes.
Results:
We found that there were nine key active ingredients in HK, which mainly acted on 21 targets. These targets primarily regulated several biological processes such as cell population proliferation, protein phosphorylation, and regulation of kinase activity, and acted on PI3K-AKT and MAPK pathways to treat GD. Analysis of the molecular interaction simulation under computer technology revealed that the key targets exhibited strong binding activity to active ingredients, and Fucosterol-AKT1 and Isofucosterol-AKT1 complexes were highly stable in humans.
Conclusion:
This study demonstrates that HK exerts therapeutic effects on GD in a multi-component, multi-target, and multi-pathway manner by regulating cell proliferation, differentiation, inflammation, and immunomodulatory-related targets. This study provides a theoretical foundation for further investigation into GD.
Introduction
Graves’ Disease (GD) is an autoimmune illness characterized by an enlarged and hyperactive thyroid gland (Graves’ hyperthyroidism), ocular abnormalities (Graves’ orbitopathy), and regional dermopathy (pretibial myxoedema). After two decades of Universal Salt Iodization (USI), GD prevalence in China has gradually decreased, with the latest epidemiological survey results at 0.53% (). However, this implies that about seven million people still have GD. The development of GD in patients is initiated with a genetic predisposition, further influenced by some environmental factors. The complex interplay manifests through thyroid follicular cell proliferation and hypertrophy and abnormal immune cell proliferation and differentiation, which is essentially the imbalance of cell value-added apoptosis driving GD hyperthyroidism occurrence. Meanwhile, recent studies have shown that elevated levels of Th1 chemokines, such as CXCL10, are found in the serum of patients with relapsed and newly diagnosed hyperthyroidism. Th1 immune response predominates in the immunopathogenesis of GD, causing the emergence and persistence of autoimmune inflammation in the thyroid gland. The onset and treatment of GD are also closely related to cytokines such as TNF-α and IL-6 (, ).
Controlling hyperthyroidism by restoring standard thyroid hormone concentrations is the principal objective of GD treatment. Some therapeutic options for Graves’ hyperthyroidism patients include Antithyroid Drugs (ATDs), radioiodine, and surgery. However, these treatment options have some drawbacks. Antithyroid medicines produce granulocyte count and liver enzyme abnormalities, and recurrence occurs in > 50% of cases after terminating ATDs (). On the other hand, radioiodine or surgery may trigger hypothyroidism, necessitating lifetime levothyroxine implementation as well as clinical and laboratory supervision (). Furthermore, treating special populations of GD patients, such as pregnant and lactating women or patients with malignant tumors, is a difficulty that cannot be overlooked.
There are records of Yingbing (goiter, 瘿病) in Traditional Chinese Medicine (TCM). The records are derived from a Chinese medicine text, Zhu Bing Yuan Hou Lun, which focusses on diagnosis and management of a class of diseases manifested as goiter, and GD is one of them. According to Waike Zhengzong, Yingbing is “not swelling of yin and yang, but stasis of blood, turbidity, and phlegm in the five organs.” Consequently, doctors mostly use Chinese medications that reduce hardness and disperse stagnation for goiter treatment. From inception to the present time, the herb pair Haizao-Kunbu (HK), first appearing in Zhou Hou Bei Ji Fang, has been a common and effective combination therapy for goiter. According to Shennong′s Classic of the Materia Medica, Haizao (Sargassum) can “master goiter, tumor qi, and the nucleus of the neck.” Additionally, Mingyi Bielu reported that Kunbu (Laminaria japonica) could “treat twelve kinds of edema and goiter tumor coalescing gas.” Some TCM formulas containing HK for goiter treatment include the Sihai Shuyu pill from Yangyi Daquan, the Haizao Yuhu decoction from Waike Zhengzong, and the Huaying Micro pill from Rumen Shiqin. Previous pharmacological research has demonstrated that besides reducing autoimmune antibody levels in rats with thyroid disease (–), Haizao has properties that can suppress cell proliferation and apoptosis induction (). Kunbu, on the other hand, has enormous clinical benefits, including anti-cell proliferation, angiogenesis inhibition, apoptosis blocking, anti-inflammatory, and antioxidant properties (, ). Furthermore, Haizao and Kunbu are considered natural medicines with minimal side effects, given their role in protecting liver and kidney function (, ). However, the precise mechanism of HK in GD treatment remains unclear. For the first time, this study will elucidate the molecular mechanisms underlying this well-known combination treatment for GD.
In TCM, there is a focus on adjusting the integrity of the human body based on the balance-regulation theory. However, due to the complex character of TCM, research on its pharmacological mechanism is challenging. Network pharmacology highlights the multi-directional signaling pathway modulation, leading to enhanced therapeutic drug benefits and decreased toxic and adverse effects, increasing the potential efficacy of new pharmaceutical clinical trials, and lowering drug discovery costs. Molecular docking, on the other hand, is a simulated testing approach that models the geometry and interactions between molecules and proteins, allowing for investigations of molecular behavior at target protein binding sites (). Finally, Molecular Dynamics (MD) modeling is a sophisticated in silico tool for exploring biological processes and molecular frameworks underlying linkages among macromolecules and ligands ().
Here, we used network pharmacological analysis, molecular docking, and MD simulation technology to create a “compound-target-pathway” network. The HK chemical compounds were collected and screened for Oral Bioavailability (OB) and Drug-Likeness (DL). Subsequently, public databases were reviewed for target and GD-related genes. To predict core compositions and targets that probably participated in GD treatment, we created a network associating medicinal ingredients with target genes in HK. Additionally, Protein-Protein Interaction (PPI), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were utilized to identify potential targets and pathways. The molecular docking and MD modeling techniques provide a foundation for further research on the HK molecular processes in GD treatment. Figure 1 depicts the study design.
Figure 1
Methods
Screening of active compounds
Potential compounds in Sargassum and Laminaria japonica were detected after searching the Traditional Chinese Medicine Systems Pharmacology (TCMSP) Database (http://tcmspw.com/) () and reassessing pertinent literature. The screening criteria included OB ≥ 30% and DL ≥ 0.18. The SwissADME online platform (www.swissadme.ch) was used to screen all ingredients ().
Prediction of targets for the HK herb pair
The TCMSP database and Swiss Target Prediction platform (www.swisstargetprediction.ch) () were used to predict the gene targets, and the corresponding Sargassum and Laminaria japonica targets were collected. The PubChem Compound Identifier (CID) number of each active ingredient was obtained from TCMSP, and potential targets were updated via Swiss Target Prediction with reference to the SMILE codes and 2D images of compounds obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov) () in SDF format. After officially annotating the potential targets on the UniProt database (https://www.uniprot.org) (), we eventually built the “herb-compound-target” network using Cytoscape 3.9.1 software ().
Acquisition of GD targets
To identify the target genes appertained to GD, we searched the OMIM (https://www.omim.org) (), GeneCards (https://www.genecards.org) (), Drugbank (https://go.drugbank.com) (), and Disgenet (https://www.disgenet.org) () databases using “Graves’ disease” as the passphrase. The online Venn mapping website () was eventually used to map the common HK and GD targets (). Obvolute proteins were then regarded as immanent therapeutic targets for GD intervention.
Protein-protein interaction network building and core target identification
After importing the common targets into STRING, we built a PPI network to further explore the role of crucial targets. In our investigation, we selected interactions with the highest confidence score (≥ 0.900) for ‘Homo sapiens’, among which we focused on the least demanding interrelation score. The unconnected nodes were concealed and the connection network was illustrated using Cytoscape 3.9.1. Additionally, the CytoNCA (Cytoscape software add-in) () was installed to examine the topological characteristics of the targets. The PPI network was used to identify key targets based on Degree value (Degree), Closeness Centrality (CC), and Betweenness Centrality (BC).
Gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment analyses
To further clarify the HK gene occupation and the role of intrinsic signaling pathways in GD, Metascape (https://metascape.org) () and David (https://david.ncifcrf.gov) () were used to assess GO and KEGG pathways of 189 target genes. The Bioinformatics platform (https://www.bioinformatics.com.cn) was used to plot the bar and bubble charts with color gradients for data analysis and visualization. The FDR error control technique was used to establish whether biological processes differed significantly. After correcting the p-value, a significance threshold of p <0.05 was used.
The molecular docking process between active ingredients and key targets
Here, we used computer-assisted technology to further confirm the intensity of the interaction between the targets and the core compound. Molecular forms of key protein targets and mol2 files of the structures of ingredients were obtained from the Protein Data Bank (PDB) (https://www.rcsb.org) () and TCMSP databases, respectively. The PyMOL software was used to eliminate the initial ligands and water molecules of proteins (). The AutoDock 4.2 program () was used to phosphorylate and store the receptor in the PDBQT file. AutoDock Vina () was utilized to dock and determine a superior model. Using binding free energy, all molecules and disease targets were ranked based on their interaction strength after docking simulations. The docking was considered valid when the binding free energy was < 5.0 kcal/mol. Finally, each target’s highest binding energy component was visualized using Discovery Studio 2019 Client, and GraphPad Prism 8.0.2 (GraphPad, CA, USA) was used to draw the binding energy heatmap.
Molecular dynamics simulation
Ligand-receptor docked complex MD simulation was performed using GROMACS (version 2021.2) (). Whereas the ligand topology file was produced by the ACPYPE script using the AMBER forefield, the protein topology file was created using the AMBER99SB-ILDN force field. For MD simulation, TIP3 water molecules were applied in a triclinic box, and periodic boundary conditions were employed. The system was neutralized using NaCl counter ions. Prior to MD simulation, the complex was reduced for 1000 steps and equilibrated by running NVT and NPT for 100 ps. The MD simulation for each system was run for 100 ns under periodic boundaries at 310 K and 1.0 bar of pressure. Finally, the free binding energy of a simulated target-ligand complex was computed using the gmx_mmpbsa tool from GROMACS.
Results
Screening of ingredients and selection of gene targets
We selected 11 components by searching the literature and online platforms and qualified nine ingredients through the screening threshold. Table 1 details the active ingredients. We identified 441 compound targets after deleting the duplicate and unreviewed genes. The HK “Chemical composition-target” network showing the connection between the nine components and 441 target genes was built using Cytoscape 3.9.1. The general characteristic of the network analysis was estimated to be 452 nodes and 751 edges (Figure 2).
Table 1
| Serial number | The scientific name of Chinese medicine | Family | Mol ID | Molecule Name | Molecule Formula | MW | AlogP | Hdon | Hacc | OB (%) | DL |
|---|---|---|---|---|---|---|---|---|---|---|---|
| KB1 | Thalluslaminariae | Laminariaceae | MOL010616 | eckol | C18H12O9 | 372.3 | 3.08 | 6 | 9 | 87.06 | 0.63 |
| KB2 | Thalluslaminariae | Laminariaceae | MOL010617 | Eicosapentaenoic Acid | C20H30O2 | 302.5 | 5.97 | 1 | 2 | 45.66 | 0.21 |
| KB3 | Thalluslaminariae | Laminariaceae | MOL001439 | arachidonic acid | C20H32O2 | 304.52 | 6.41 | 1 | 2 | 45.57 | 0.2 |
| KB4 | Thalluslaminariae | Laminariaceae | MOL000953 | CLR | C27H46O | 386.73 | 7.38 | 1 | 1 | 37.87 | 0.68 |
| KB5 | Thalluslaminariae | Laminariaceae | MOL009622 | Fucosterol | C29H48O | 412.77 | 7.83 | 1 | 1 | 43.78 | 0.76 |
| HZ1 | Sargassum | Sargassaceae | MOL010578 | Aurantiamide | C25H26N2O3 | 402.53 | 3.64 | 3 | 5 | 45.76 | 0.43 |
| HZ2 | Sargassum | Sargassaceae | MOL010580 | Diglycol dibenzoate | C18H18O5 | 314.36 | 3.06 | 0 | 5 | 59.22 | 0.27 |
| HZ3 | Sargassum | Sargassaceae | MOL005440 | Isofucosterol | C29H48O | 412.77 | 7.83 | 1 | 1 | 43.78 | 0.76 |
| HZ4 | Sargassum | Sargassaceae | MOL000098 | quercetin | C15H10O7 | 302.25 | 1.5 | 5 | 7 | 46.43 | 0.28 |
The details of the active compounds in HK herb pair.
Figure 2
GD targets searching
The OMIM, DrugBank, GeneCards, and Disgenet databases yielded 570, 12, 2470, and 585 GD targets, respectively. After collecting all the genes and removing duplicate data, 2,010 targets remained for further research (Supplementary Table S1). Finally, the Venn diagram displayed 189 genes as the herb pair’s implicit aim for GD therapies (Figure 3).
Figure 3
PPI analysis
We imported 189 potential targets into STRING and obtained 165 potential targets through the screening threshold. Subsequently, we created a PPI network (165 nodes and 771 edges) after importing target genes into Cytoscape 3.9.1 software (Figure 3). Increasing quantified values were associated with the improved significance of the node. The key target topological analysis was based on Degree, CC, and BC > one-fold of the median. We first sorted 21 key gene targets using the Degree > 7, CC > 0.393, and BC > 0.003 criteria and then sorted 19 key gene targets using the Degree > 12, CC > 0.433, and BC > 0.013 criteria (Supplementary Table S2).
GO terms and KEGG pathways
We completed GO and KEGG enrichment studies to further demonstrate the proposed targets’ organic functions and prospective mechanisms. The analyses involved 189 HK potential target genes underlying GD. The results were obtained after analysis of the GO and KEGG enrichment results analyzed by Metascape and David databases (Figure 4). We discovered that the common GO biological processes in the two databases were protein phosphorylation (GO: 0006468) and positive cell migration regulation (GO: 0030335), implying that phosphorylation and cell migration may be significantly involved in GD treatment by HK. Additionally, we discovered that protein homodimerization activity (GO: 0042803), protein kinase binding (GO: 0019901), protein kinase activity (GO: 0004672), kinase activity (GO: 0016301), and protein serine/threonine kinase activity (GO: 0004674) were simultaneous entries among the GO molecular functions. This finding indicates that protein kinase activity was crucial throughout the GD treatment process. Furthermore, terms such as cytoplasm perinuclear region (GO: 0048471) and membrane raft (GO: 0045121) were found to be crucial cellular components depending on the cellular composition. Cell membranes and cytoplasm are vital organelles that regulate cell signaling.
Figure 4
A KEGG pathway enrichment assessment was performed to further investigate the possible functions of the targets (P<0.05). The following 19 pathways appeared simultaneously in the top 20 critical pathways analyzed by the two databases: Cancer Pathways (hsa05200), the PI3K-AKT Signaling Pathway (hsa04151), Lipid and Atherosclerosis (hsa05417), Proteoglycans in Cancer (hsa05205), Human Cytomegalovirus (HCMV) Infection (hsa05163), Kaposi Sarcoma-associated Herpesvirus (KSHV) Infection (hsa05167), Fluid Shear Stress and Atherosclerosis (hsa05418), the AGE-RAGE Signaling Pathway in Diabetic Complications (hsa04933), Human Papillomavirus (HPV) Infection (hsa05165), Neurodegeneration-multiple Disease Pathways (hsa05022), the MAPK Signaling Pathway (hsa04010), Hepatitis B (hsa05161), Tuberculosis (hsa05152), Epstein-Barr Virus (EBV) Infection (hsa05169), Chemical Carcinogenesis-Reactive Oxygen Species (hsa05208), MicroRNAs in Cancer (hsa05206), Focal Adhesion (hsa04510), Alzheimer’s Disease (hsa05010), and Endocrine Resistance (hsa01522). Based on this outcome, we finally constructed a “KEGG pathway-gene target” network (Figure 5). Detailed information is provided in Supplementary Table S3.
Figure 5
Molecular docking verification
We used molecular docking to identify chemicals from the PPI system based on GD-related targets. The relationships among essential active chemicals and crucial targets were examined using AutoDock4.2, Discovery Studio 2019 Client, and PyMOL software. Molecular structures of significant targets were obtained from RCSB PDB (Table 2). We did not find the PDB ID of FOS. Since the docking of CAV1, JAK2, and ITGB3 receptors to compounds was not successfully implemented in Autodock Vina, we redocked using the Genetic Algorithm in AutoDock 4.2., and the docking analysis yielded the binding scores (Figure 6). A lower value indicated a greater binding capacity. We discovered that van der Waals forces, hydrogen bonds, and aromatic stacking (Pi-Cation, Pi-Anion, Pi-Sulfur, Pi-alkyl, and alkyl interactions) were involved in the interactions between the active site residues of crucial targets and potential active substances. All active compounds exhibited an excellent binding affinity to specific core targets (only binding energies < -5 kcal/mol were demonstrated). Fucosterol-AKT1, Isofucosterol-AKT1, Isofucosterol-MAPK8, Fucosterol-MAPK8, and Aurantiamide-AKT1 were the top five binding modes (Figure 7). The Fucosterol-AKT1 complex was stabilized by 15 van der Waals forces and 5 Pi-alkyl and alkyl interactions with CYS 77, CYS 60, TRP 80, LEU 264, and VAL 270, respectively (Figure 7A). The Isofucosterol-AKT1 complex was stabilized by one hydrogen bond (1H-bond) with residue ALA 58, and four Pi-alkyl and alkyl interactions with LEU 210, TRP 80, LYS 268, and VAL 270 (Figure 7B). On the other hand, the Isofucosterol-MAPK8 complex was stabilized by 1H-bond with residue ASP 112 and four alkyl interactions with ILE 32, VAL 158, VAL 40, and LEU 168 (Figure 7C). Similarly, a 1H-bond with residue ASP 112 and six alkyl interactions with ILE 32, MET 108, VAL 158, VAL 40, LEU 168, and ILE 86 stabilized the Fucosterol-MAPK8 complex (Figure 7D). Finally, the Aurantiamide-AKT1 complex was stabilized by three Pi-cation and Pi-anion interactions with ARG 273, ARG 86, GLU 298, two H-bonds with residue GLU 85 and GLU 17, one Pi-alkyl bond with ILE 84, and one Pi-sulfur bond with CYS 310 (Figure 7E).
Table 2
| Targets | PDB ID | UniProt ID | Center grid box | ||
|---|---|---|---|---|---|
| X center | Y center | Z center | |||
| CAV1 | 7SC0 | Q03135 | 0 | 0 | 0 |
| MAPK14 | 1OVE | Q16539 | 29.747 | 15.683 | 27.293 |
| JAK2 | 4IVA | O60674 | 0.161 | -11.654 | -8.084 |
| PTPN11 | 3B7O | Q06124 | 28.81 | 8.938 | 63.391 |
| ITGB3 | 4G1M | P05106 | -36.16 | 46.505 | 55.032 |
| PRKCD | 3UFF | Q05655 | 0 | 0 | 0 |
| AKT1 | 4EJN | Q01314 | 30.889 | 52.185 | 19.493 |
| VEGFA | 4QAF | P15692 | 13.309 | 63.036 | -0.942 |
| ITGB1 | 7NXD | P05556 | 130.433 | 182.225 | 111.984 |
| ESR1 | 5AAV | P03372 | 31.01 | 14.202 | 10.787 |
| PIK3R1 | 3I5S | P27986 | 0 | 0 | 0 |
| MAPK1 | 2OJG | P28482 | -13.772 | 13.979 | 41.667 |
| PTK2 | 1MP8 | Q05397 | 36.299 | -3.761 | 24.196 |
| SRC | 1O43 | P12931 | 18.953 | 20.632 | 21.188 |
| IL6 | 1ALU | P05231 | -7.7 | -12.7 | 0 |
| TP53 | 4AGP | P04637 | 91.098 | 96.917 | -46.275 |
| STAT1 | 1YVL | P42224 | -30.304 | -13.959 | 146.805 |
| MAPK8 | 4QTD | P45983 | 14.189 | 15.864 | 19.659 |
| EGFR | 5UG9 | P00533 | -8.371 | 17.712 | -12.846 |
| TNF | 5UUI | P01375 | 41.438 | 43.125 | 1.22 |
| FOS | – | P01100 | – | – | – |
Grid docking parameters in molecular docking.
Figure 6
Figure 7
Molecular dynamics simulation of structural stability and interaction energy
We selected and implemented the top two compound-target dockings (Fucosterol-AKT1 and Isofucosterol-AKT1) for MD simulations. After 100 ns of MD simulations, the dynamic variations of the Fucosterol-AKT1 and Isofucosterol-AKT1 complexes were assessed. Subsequently, we computed the Root-Mean-Square Deviation (RMSD) to understand the complexes’ molecular configurations and the system’s stability during simulation. The RMSD curve represents the location deviations in the protein. Fucosterol-AKT1 and Isofucosterol-AKT1 had average RMSD values of 2.5 and 2.6 Å, respectively. The RMSD of the Fucosterol-AKT1 complex was less than that of the Isofucosterol-AKT1 complex between the 45-85 ns time interval, and the RMSD curves relatively stabilized after 85 ns (Figure 8A). The Root-Mean-Square Fluctuation (RMSF) graph shows the protein amino acid residue variations. According to the findings, most simulations had small alterations in amino acid structure. Furthermore, most of the residues exhibited minor structural modifications. The RMSF values of residue numbers 50-200 in AKT1 after Fucosterol binding showed greater flexibility than the same regions in AKT1 after Isofucosterol binding (Figure 8B). The gyration (Rg) radius curve represents the tightness of the protein’s general configuration. The Fucosterol-AKT1 and the Isofucosterol-AKT1 complexes had stable rotation radii, although the former folded with greater force (Figure 8C). The binding free energy can be used to assess the change in the binding pattern and stability of ligands and proteins. Compared to the isofucosterol-AKT1 complex, which showed an average interaction energy and energy fluctuation of 144.536 kcal/mol and 17.43 kcal/mol, respectively, the fucosterol-AKT1 complex had an average interaction energy and energy fluctuation of 141.412 kcal/mol and13.63 kcal/mol, respectively (Figure 8D). The number of Hydrogen bonds (H-bonds) in a complex might reveal information about its binding strength. The ligands and residues of all five protein compartments created one or several hydrogen bonding connections. During the 100 ns simulations, the Fucosterol-AKT1 complex had a higher H-bond density and size (Figures 8E, F).
Figure 8
Discussion
The GD etiology and pathogenesis in relation to heredity, mental stimulation, environment, infection, and other factors have not been fully elucidated. Severe GD cases might result in life-threatening complications such as liver damage, heart failure, and thyroid storm. Besides being ineffective, current GD therapies have also been linked with severe limitations, such as toxic pharmacological side effects, permanent hypothyroidism, vocal cord paralysis, and other destructive damages, that all cause clinical challenges. There has been a growing demand for novel strategies that utilize traditional approaches to treat GD. In this regard, as a complement, traditional medicine has been shown to generate synergistic results with minimal toxicity. Guided by TCM theory, HK was the most common drug combination used to treat GD. However, the molecular mechanisms of HK in GD treatment have not been fully explored. Herein, we searched and selected QFP compounds and targets from multiple databases and created a “compound-gene target-disease” network. The CytoNCA software was used to perform topological analysis to estimate key protein targets involved in GD treatment and 15 essential targets including AKT1, PIK3R1, MAPK1, MAPK8, MAPK14, VEGFA, EGFR, IL6, TNF, TP53, JAK2, STAT1, ESR1, FOS, and CAV1 were identified. Subsequently, we excluded a list of irrelevant pathways in the top 20 most enriched pathways to elucidate the pathogenesis of GD and establish a prospective network model. The PI3K-AKT and MAPK pathways were identified as the most prominent penetration points in GD treatment. Furthermore, the fundamental biological processes, molecular functions, and cellular components in GO enrichment analysis helped illustrate the multi-dimensional and multi-target therapeutic pathways. Finally, the computational and validation molecular docking and MD simulation processes enabled us to further discuss the biological function of HK in the GD treatment from a micro-level perspective and to present a theoretical framework for GD therapy with TCM.
In the TCM context, HK treats GD by softening firmness and dispersing stagnation. However, Haizao and Kunbu, known as iodine-rich Chinese medicine, have been controversial in GD treatment due to their higher iodine content. Wild Laminaria japonica (0.1-1.4 g) contains approximately 150 μg of iodine, while the iodine concentration of wild Sargassum is 58-629 μg/g (). However, in China, HK refers to specially processed TCM pieces. According to current research, the iodine content range of Haizao and Kunbu pieces are 297.67-814.59 μg/g and 855.33-5481.33 μg/g, respectively (, ). Consistent with an Italian study, dried pieces had higher iodine content than fresh samples (), which can also be attributed to the different plant origins. Here, we studied processed Chinese herbal pieces. Chinese herbal pieces should be decocted before use, and studies have shown that boiling in fresh water can reduce the seaweed iodine content (). Additionally, following gastrointestinal digestion, only 49-82% of seaweed iodine appears to be available for human absorption (). Despite its substantial iodine amounts, ingesting considerable quantities of HK does not necessarily indicate a risk of excessive iodine intake (). There is a clinical belief that GD hyperthyroidism patients should strictly limit iodine intake. However, some new perspectives suggest that limiting iodine intake might not benefit GD patients in areas with adequate or excessive iodine intake (). The considerable Chinese population that participated in the USI program for two decades has reached an iodine adequacy status, and the once-high hyperthyroidism prevalence has consequently been reduced to relatively stable levels (). The latest study proves that the cumulative post-USI hyperthyroidism incidence in different iodine trophic status areas in China has not increased (). Studies have also shown that limiting dietary iodine intake does not improve ATD effectiveness in GD treatment, nor does it increase the GD recurrence rate after ATD withdrawal, and even severe iodine intake restriction can negatively affect GD (–). Furthermore, some GD patients who received an acceptable excess of iodide during treatment experienced reduced hyperthyroidism and thiourea-related side effects (). Moreover, this treatment did not affect the efficacy of radioactive iodine therapy (). Furthermore, potassium iodide has been shown to be effective and safe in specific GD patients, such as pregnant and breastfeeding women, and patients with malignant tumors, or those undergoing radiation and chemotherapy (, ). An in vivo study from Japan recently suggested that the chronic anti-thyroid action of iodine in GD involves hormone secretion inhibition (). In this context, researchers believe that iodine-rich TCMs such as HK can be rationally used in treating some GD patients, such as those who are intolerant to ATD or refuse surgical treatment, and patients with mild and moderate GD whose serum FT4 and Thyroid Stimulating Hormone (TSH) Receptor (TSHR) autoantibody (TRAb) levels are less than the upper limit of the detectable range in a laboratory (–). In our previous study, we demonstrated that an iodine-rich Chinese medicine formula could improve thyroid function and morphology in hyperthyroid rats (). Furthermore, some previous clinical studies have proved that iodine-rich TCM for treating GD has the rapid onset, minor adverse reactions, and reduced serum TRAb advantages (, ). Some researchers in China have recently been conducting large-sample, multi-center, and strictly designed clinical studies on iodine-rich TCM for GD treatment (Registration numbers: ChiCTR2000032706, ChiCTR1900021572), and their findings will provide additional insights and inspiration for our subsequent research. Therefore, in future studies, we recommend a reasonable and comprehensive evaluation of the efficacy and safety of HK and further verification of its clinical effect in GD patients.
Based on our predictions, most HK active ingredients can influence biological processes, including proliferation, apoptosis, and inflammation. Aurantimide is a critical active small-molecule compound, and in vitro experiments have demonstrated that it could act on the PI3K/AKT signaling pathway () and exert anti-inflammatory effects via inhibiting the phosphorylation of the MAPK pathway (). Eckol, a novel natural phizolian derived from marine brown algae, has been shown to downregulate EGFR, p-EGFR, JAK2, and STAT3 expression in tumor cells, indicating pro-apoptotic and anti-proliferative activities (, ). Eicosapentaenoic Acid (EPA), a polyunsaturated omega-3 fatty acid, has previously been reported to affect cell proliferation and inflammatory responses by inhibiting the phosphorylation of AKT (, ). Besides improving hormonal status in hyperthyroid rats, including T3 and TSH levels, EPA can also reduce the levels of pro-inflammatory cytokines, such as serum TNF-α (). Fucosterol, an algae-derived unique plant sterol with various medicinal properties, has been predicted by studies to suppress the phosphorylation of Phosphatidylinositol 3-kinase/protein Kinase B (PI3K/Akt) signaling, reduce Mitogen-Activated Protein Kinase (MAPK) expression, and reduce IL-6 and TNF-α overexpression (–), influencing cell proliferation, apoptosis, and inflammation. On the other hand, quercetin is an excellent antioxidant that exerts sound anti-inflammatory effects and has been reported to inhibit the PI3k/Akt pathway by effectively binding to PIK3R1 (, ). Arachidonic acid metabolites are differentially affected by thyroid hormone status, and elevated levels of AA metabolites have been observed in the serum of patients with hyperthyroidism (). We speculated that it is related to the Inhibitory effects of iodinated derivatives of arachidonic acid on iodine metabolism. Although the mechanism of the above active compounds in GD pathogenesis has not been fully explained, we could predict the optimal mechanism between compounds and target proteins by combining molecular docking results. Furthermore, verifying their biological feeatures in GD models is one of our future research directions.
The PPI network topology analysis revealed 21 core HK targets for treating GD. These protein targets are primarily involved in cell proliferation, apoptosis, and inflammatory processes in GD pathogenesis. The Vascular Endothelial Growth Factor (VEGF) family is required for the proliferation of blood vessels (), and GD patients have elevated serum VEGF levels (). According to research, iodide can decrease the expression of VEGFAs (VEGFs involved in promoting angiogenesis) upregulated by TSH (). Besides eliminating vascular remodeling and inhibiting hair follicle hypertrophy (), blocking VEGFA inhibits hyperthyroidism by increasing lymphatic flow in the Graves thyroid gland (), decreasing thyroid weight during goiter development. Caveolin-1 (CAV1) is a member of the thyrosomal polyprotein complex required for thyroid hormone synthesis and thyroid cell homeostasis (). Low CAV1 expression was observed in fat cells of GD patients (). Pro-Epidermal Growth Factor (EGF) is an essential growth factor in thyroid tissue, and nuclear EGFR expression is elevated in GD tissue samples, implying that the EGFR-dependent modulation of thyroid cell proliferation under physiological conditions may be associated with hyperthyroidism (). The Tumor Protein p53 (TP53) gene is important for inducing apoptosis or cell cycle interruption, and the succession of an insufficiently effective TP53 gene substantially raises the risk of developing GD. Given that the autoimmune thyroid illness may be accompanied by DNA damage and apoptosis, the insufficiently effective TP53 gene may initiate and sustain the autoimmune GD process (). Furthermore, some well-known cytokines (IL6 and TNF-α) are an essential part of the autoimmune response in GD patients, which can promote inflammatory cell proliferation and infiltration into thyroid tissue and affect thyroid follicular cell growth and differentiation. Furthermore, these cytokines have been proven to be related to the recalcitrant nature of the disease and the severity of clinical symptoms (–). Other genes have also been found to be partially responsible for the regulated proliferation and Thyroid Hormone (TH) levels in GD (–). These gene targets imply that we can further investigate the mechanism of HK in GD treatment in terms of cell proliferation, apoptosis, and inflammatory processes.
Additionally, KEGG enrichment analysis demonstrated that the PI3K-AKT and MAPK signaling pathways were involved in GD onset. Furthermore, some of the BP and MF items examined by GO, including protein phosphorylation, positive cell death regulation, cell activation, kinase activity regulation, MAPK cascade regulation, cell population proliferation, Transcription Factor (TF) binding, protein serine/threonine kinase activity, cytokine receptor binding, and Cell Adhesion Molecule (CAM) binding, revealed the significance of cell proliferation, differentiation, protein phosphorylation, and protein kinases in GD pathogenesis. Notably, GD patients have thicker, hypertrophied follicular cells in their thyroid glands that produce active thyroglobulin. The gland exhibits classic lymphocytic infiltrates considered to be principally connected with TSHR autoantibody secretion. Additionally, a histological examination showed occasional apoptotic cells and partial follicular wreckage (86). The presence of immune inflammation, cell proliferation, and apoptosis were all essential components in the pathological mechanism of GD. Although studies have reported that GD patients have all three types of TSHR autoantibodies, stimulating antibodies constitute the distinguishing characteristic of Grave’s hyperthyroidism. When stimulated, TSHR autoantibodies can induce complex signaling cascades, mainly activating GΑs and inducing cAMP/PKA pathways. The production of cAMP activates the cAMP Response Element-Binding Protein (CREB) and protein kinase A, which are directly or indirectly involved in inflammatory mediator generation, thyroid hormone synthesis, and thyroid cell proliferation (, 87). The bioactive exertion of cAMP partially depends on the phosphorylation of the MAPK and PI3K/AKT pathways (88). Additionally, the binding of Β-arrestin attracted by TSHR to the receptor can activate the MAPK pathway (89), inducing protein synthesis, cell differentiation, and angiogenesis via hemodynamic effects (, 90). Protein Kinase B (PKB), also known as AKT, is an intracellular signaling pathway phosphorylated and activated on the plasma membrane. Once activated, AKT regulates cell survival/apoptosis, cell proliferation, and protein synthesis. Recent findings have confirmed the involvement of the PI3K-AKT pathway in GD pathogenesis (91–93). Phosphorylated AKT stimulates the activation of TFs CREB and NF-kB, as well as inflammatory gene expression (such as IL-6) (94). Previous research showed that the cAMP/PKA pathway increased IL-6 production in thyroid cells via processes influencing the stability of IL-6 mRNA, IL-6 gene promoter, and c-Fos expression (95). Increased inflammatory chemical levels further worsen thyroid follicular cell stimulation and destruction, causing the secretion of more thyroid hormones and amplifying the body’s inflammatory response (96). Moreover, Janus Kinase (JAK)/Signal Transducer and Activator of Transcription (STAT) and their downstream effectors are the primary signaling cascades in TSHR activation. Our results and previous findings suggest that JAK2 and STAT1 phosphorylation may occur during GD treatment (97, 98). All the above-mentioned cascades have relevance to cellular development, survival, differentiation, cytokine and chemokine release, and apoptosis induction (99), indicating that thyroid cell activation and proliferation can be partially influenced via the regulation of protein phosphorylation within these central signaling cascades by HK (Figure 9).
Figure 9
MD simulation is an effective technique for examining molecular mechanisms, especially in evaluating the binding stability (100) and the selectivity of specific ligands for their target proteins (101). In this study, we conducted MD simulations to gain deeper insights into the stability of the protein-ligand complexes, specifically focusing on the Fucosterol-AKT1 and Isofucosterol-AKT1 complexes. We assessed the RMSD profiles of these complexes, which exhibited consistent patterns with RMSD values hovering around 0.27 nm. This pattern suggests that the system reached equilibrium during the simulation. We extended the MD simulation duration for these two complexes to generate a stable RMSD profile. In addition to RMSD, we evaluated several other factors, including RMSF, Rg, H-bonds, and binding free energy. These assessments provided further insights into the dynamic properties of the protein-ligand interactions. Overall, these results suggest that during MD simulation, fucosterol-AKT1 and isofucosterol-AKT1 complexes were stable and equilibrated.
Although network pharmacology, molecular docking, and molecular dynamics simulation methods were implemented to characterize the potential chemicals and targets of HK, there are some several limitations that need to be acknowledged. First, to accurately comprehend the behavior of chemical elements functioning on disease targets, network pharmacology can be utilized to predict the up-regulation and down-regulation of targets. Second, in-vivo and in-vitvo experiments should be conducted to accurately explore the therapeutic effects and associated mechanisms of the potential compound and target pairs. It is also necessary to identify the active ingredients with therapeutic effects in HK through basic experiments. These will be the focus of our future research.
Conclusion
GD is an autoimmune disease caused by genetic and environmental factors. The multifaceted nature of GD calls for a multifaceted approach to its treatment, often involving multiple targets. There is a growing interest in innovative treatment strategies that draw from the wisdom of traditional medicine and apply it to the clinical management of GD. These strategies have demonstrated favorable efficacy and safety profiles. Nine compounds including Fucosterol, Isofucosterol, Aurantiamide, Eicosapentaenoic Acid, quercetin and eckol of HK were combined to 21 core targets including AKT1, PIK3R1, MAPK1, MAPK8, MAPK14, VEGFA, EGFR, IL6, TNF, TP53, CAV1, JAK2, and STAT1, and by participating in biological processes such as cell population proliferation, protein phosphorylation, regulation of kinase activity and protein serine/threonine kinase activity and PI3K-AKT, MAPK, and other pathways play a crucial role in the treatment of GD. Molecular docking and molecular dynamics simulation demonstrated that the Fucosterol-AKT1 and Isofucosterol-AKT1 complexes exhibited the highest binding energy, indicating that HK contains key compounds and targets. While future research will necessitate further biological experiments, including additional in vitro and in vivo studies, our current findings suggest that HK holds significant promise as an effective herbal remedy for treating GD.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author/s.
Author contributions
Methodology: MY, YL, XH. Software: MY, XH. Validation: MY, XH. Formal analysis: MY, XH. Visualization: MY, YL, QL, YA. Conceptualization: YL, TG. Investigation: YL, DG, QL, YW, YA. Supervision: YL, TG. Resources: DG, QL, YW, YA. Data Curation: DG, QL, YW. Writing - Original Draft: MY. Writing - Review & Editing: YL, XH, TG. Project administration: YL, TG. Funding acquisition: YL. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Foundation of Shenyang Bureau of Science and Technology [grant number 22-321-34-12].
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2023.1236549/full#supplementary-material
References
1
WangCLiYTengDShiXBaJChenBet al. Hyperthyroidism prevalence in China after universal salt iodization. Front Endocrinol (Lausanne) (2021) 12:651534. doi: 10.3389/fendo.2021.651534
2
AntonelliAFallahiPEliaGRagusaFPaparoSRRuffilliIet al. Graves' Disease: clinical manifestations, immune pathogenesis (Cytokines and chemokines) and therapy. Best Pract Res Clin Endocrinol Metab (2020) 34(1):101388. doi: 10.1016/j.beem.2020.101388
3
AntonelliAFerrariSMRagusaFEliaGPaparoSRRuffilliIet al. Graves' Disease: Epidemiology, genetic and environmental risk factors and viruses. Best Pract Res Clin Endocrinol Metab (2020) 34(1):101387. doi: 10.1016/j.beem.2020.101387
4
DaviesTFAndersenSLatifRNagayamaYBarbesinoGBritoMet al. Graves' Disease. Nat Rev Dis Primers (2020) 6(1):52. doi: 10.1038/s41572-020-0184-y
5
SmithTJHegedüsL. Graves' Disease. N Engl J Med (2016) 375(16):1552–65. doi: 10.1056/NEJMra1510030
6
KokubuSNishiharaGNWatanabeYTsuchiyaYAmamoYTeradaR. The effect of irradiance and temperature on the photosynthesis of a native alga sargassum fusiforme (Fucales) from kagoshima, Japan. Phycologia (2015) 54(3):235–47. doi: 10.2216/15-007.1
7
LiuLHeinrichMMyersSDworjanynSA. Towards a better understanding of medicinal uses of the brown seaweed sargassum in traditional chinese medicine: A phytochemical and pharmacological review. J Ethnopharmacol (2012) 142(3):591–619. doi: 10.1016/j.jep.2012.05.046
8
SongX-hZanR-zYuC-hWangF. Effects of modified haizao yuhu decoction in experimental autoimmune thyroiditis rats. J Ethnopharmacol (2011) 135(2):321–4. doi: 10.1016/j.jep.2011.03.017
9
HeoS-JKimK-NYoonW-JOhCChoiY-UAffanAet al. Chromene induces apoptosis via caspase-3 activation in human leukemia hl-60 cells. Food Chem Toxicol (2011) 49(9):1998–2004. doi: 10.1016/j.fct.2011.05.011
10
AhmedSAMendoncaPElhagRSolimanKFA. Anticancer effects of fucoxanthin through cell cycle arrest, apoptosis induction, angiogenesis inhibition, and autophagy modulation. Int J Mol Sci (2022) 23(24):16091. doi: 10.3390/ijms232416091
11
LuanFZouJRaoZJiYLeiZPengLet al. Polysaccharides from laminaria japonica: An insight into the current research on structural features and biological properties. Food Funct (2021) 12(10):4254–83. doi: 10.1039/d1fo00311a
12
RamuSMuraliANarasimhaiahGJayaramanA. Toxicological evaluation of sargassum wightii greville derived fucoidan in wistar rats: Haematological, biochemical and histopathological evidences. Toxicol Rep (2020) 7:874–82. doi: 10.1016/j.toxrep.2020.07.009
13
Abdel-DaimMMAbushoukAIBahbahEIBungăuSGAlyousifMSAleyaLet al. Fucoidan protects against subacute diazinon-induced oxidative damage in cardiac, hepatic, and renal tissues. Environ Sci pollut Res Int (2020) 27(11):11554–64. doi: 10.1007/s11356-020-07711-w
14
AbdelsattarASDawoudAHelalMA. Interaction of nanoparticles with biological macromolecules: A review of molecular docking studies. Nanotoxicology (2021) 15(1):66–95. doi: 10.1080/17435390.2020.1842537
15
DecherchiSCavalliA. Thermodynamics and kinetics of drug-target binding by molecular simulation. Chem Rev (2020) 120(23):12788–833. doi: 10.1021/acs.chemrev.0c00534
16
RuJLiPWangJZhouWLiBHuangCet al. Tcmsp: A database of systems pharmacology for drug discovery from herbal medicines. J Cheminform (2014) 6:13. doi: 10.1186/1758-2946-6-13
17
DainaAMichielinOZoeteV. Swissadme: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci Rep (2017) 7:42717. doi: 10.1038/srep42717
18
DainaAMichielinOZoeteV. Swisstargetprediction: Updated data and new features for efficient prediction of protein targets of small molecules. Nucleic Acids Res (2019) 47(W1):W357–W64. doi: 10.1093/nar/gkz382
19
KimSChenJChengTGindulyteAHeJHeSet al. Pubchem 2019 update: Improved access to chemical data. Nucleic Acids Res (2019) 47(D1):D1102–D9. doi: 10.1093/nar/gky1033
20
ConsortiumU. Uniprot: The universal protein knowledgebase in 2023. Nucleic Acids Res (2023) 51(D1):D523–D31. doi: 10.1093/nar/gkac1052
21
ShannonPMarkielAOzierOBaligaNSWangJTRamageDet al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res (2003) 13(11):2498–504. doi: 10.1101/gr.1239303
22
KuiFGuWGaoFNiuYLiWZhangYet al. Research on effect and mechanism of xuefu zhuyu decoction on chd based on meta-analysis and network pharmacology. Evid Based Complement Alternat Med (2021) 2021:9473531. doi: 10.1155/2021/9473531
23
StelzerGRosenNPlaschkesIZimmermanSTwikMFishilevichSet al. The genecards suite: From gene data mining to disease genome sequence analyses. Curr Protoc Bioinf (2016) 54:1.30.1–1.30.33. doi: 10.1002/cpbi.5
24
WishartDSFeunangYDGuoACLoEJMarcuAGrantJRet al. Drugbank 5.0: A major update to the drugbank database for 2018. Nucleic Acids Res (2018) 46(D1):D1074–D82. doi: 10.1093/nar/gkx1037
25
PiñeroJRamírez-AnguitaJMSaüch-PitarchJRonzanoFCentenoESanzFet al. The disgenet knowledge platform for disease genomics: 2019 update. Nucleic Acids Res (2020) 48(D1):D845–D55. doi: 10.1093/nar/gkz1021
26
BardouPMarietteJEscudiéFDjemielCKloppC. Jvenn: An interactive venn diagram viewer. BMC Bioinf (2014) 15(1):293. doi: 10.1186/1471-2105-15-293
27
TangYLiMWangJPanYWuF-X. Cytonca: A cytoscape plugin for centrality analysis and evaluation of protein interaction networks. Biosystems (2015) 127:67–72. doi: 10.1016/j.biosystems.2014.11.005
28
ZhouYZhouBPacheLChangMKhodabakhshiAHTanaseichukOet al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun (2019) 10(1):1523. doi: 10.1038/s41467-019-09234-6
29
ShermanBTHaoMQiuJJiaoXBaselerMWLaneHCet al. David: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res (2022) 50(W1):W216–W21. doi: 10.1093/nar/gkac194
30
BurleySKBhikadiyaCBiCBittrichSChenLCrichlowGVet al. Rcsb protein data bank: Powerful new tools for exploring 3d structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences. Nucleic Acids Res (2021) 49(D1):D437–D51. doi: 10.1093/nar/gkaa1038
31
SchrodingerLLC. The pymol molecular graphics system, version 1.8. (2015).
32
MorrisGMHueyRLindstromWSannerMFBelewRKGoodsellDSet al. Autodock4 and autodocktools4: Automated docking with selective receptor flexibility. J Comput Chem (2009) 30(16):2785–91. doi: 10.1002/jcc.21256
33
EberhardtJSantos-MartinsDTillackAFForliS. Autodock vina 1.2.0: New docking methods, expanded force field, and python bindings. J Chem Inf Model (2021) 61(8):3891–8. doi: 10.1021/acs.jcim.1c00203
34
AbrahamMJMurtolaTSchulzRPállSSmithJCHessBet al. Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX (2015) 1-2:19–25. doi: 10.1016/j.softx.2015.06.001
35
BlikraMJHenjumSAakreI. Iodine from brown algae in human nutrition, with an emphasis on bioaccessibility, bioavailability, chemistry, and effects of processing: A systematic review. Compr Rev Food Sci Food Saf (2022) 21(2):1517–36. doi: 10.1111/1541-4337.12918
36
GanDGaoCWangYGaoT. Determination of iodine content in herbal pieces and dispensing granules of kunbu( Thallus laminariae thallus eckloniae) and haizao( Sargassum) from different sources. Chin Arch Traditional Chin Med (2022) 40(12):218–21. doi: 10.13193/j.issn.1673-7717.2022.12.045
37
XiangPZhaoYChenGLiuKWangXLiuC. Determination of iodine content in chinese decoction and herbal for thyroid diseases. Guiding J Traditional Chin Med Pharm (2019) 25(13):94–6+120. doi: 10.13862/j.cnki.cn43-1446/r.2019.13.028
38
FilippiniMBaldisserottoAMenottaSFedrizziGRubiniSGigliottiDet al. Heavy metals and potential risks in edible seaweed on the market in Italy. Chemosphere (2021) 263:127983. doi: 10.1016/j.chemosphere.2020.127983
39
Domínguez-GonzálezMRChiocchettiGMHerbello-HermeloPVélezDDevesaVBermejo-BarreraP. Evaluation of iodine bioavailability in seaweed using in vitro methods. J Agric Food Chem (2017) 65(38):8435–42. doi: 10.1021/acs.jafc.7b02151
40
RoledaMYSkjermoJMarfaingHJónsdóttirRReboursCGietlAet al. Iodine content in bulk biomass of wild-harvested and cultivated edible seaweeds: Inherent variations determine species-specific daily allowable consumption. Food Chem (2018) 254:333–9. doi: 10.1016/j.foodchem.2018.02.024
41
HuangYXuYXuMZhaoXChenM. Application of oral inorganic iodine in the treatment of graves' Disease. Front Endocrinol (Lausanne) (2023) 14:1150036. doi: 10.3389/fendo.2023.1150036
42
ShanZLiYLiYWangHTengDTengXet al. Changing iodine status and the incidence of thyroid disease in mainland China: A prospective 20-year follow-up study. Thyroid (2023) 33(7):858–66. doi: 10.1089/thy.2022.0505
43
HuangHShiYLiangBCaiHCaiQLinR. Optimal iodine supplementation during antithyroid drug therapy for graves' Disease is associated with lower recurrence rates than iodine restriction. Clin Endocrinol (Oxf) (2018) 88(3):473–8. doi: 10.1111/cen.13543
44
ParkSMChoYYJoungJYSohnSYKimSWChungJH. Excessive iodine intake does not increase the recurrence rate of graves' Disease after withdrawal of the antithyroid drug in an iodine-replete area. Eur Thyroid J (2015) 4(1):36–42. doi: 10.1159/000375261
45
SantarosaVAOrlandiDMFiorinLBKasamatsuTSFuruzawaGKKuniiISet al. Low iodine diet does not improve the efficacy of radioiodine for the treatment of graves' Disease. Arch Endocrinol Metab (2015) 59(6):501–6. doi: 10.1590/2359-3997000000082
46
OkamuraKSatoKFujikawaMBandaiSIkenoueHKitazonoT. Remission after potassium iodide therapy in patients with graves' Hyperthyroidism exhibiting thionamide-associated side effects. J Clin Endocrinol Metab (2014) 99(11):3995–4002. doi: 10.1210/jc.2013-4466
47
NishioRUchidaTSuzukiLOnoseHGotoHYamadaEet al. Influence of short-term dietary and therapeutic iodine restriction on the therapeutic effects of radioactive iodine therapy in patients with graves' Disease. Thyroid (2021) 31(3):439–45. doi: 10.1089/thy.2020.0126
48
YoshiharaANohJYWatanabeNMukasaKOhyeHSuzukiMet al. Substituting potassium iodide for methimazole as the treatment for graves' Disease during the first trimester may reduce the incidence of congenital anomalies: A retrospective study at a single medical institution in Japan. Thyroid (2015) 25(10):1155–61. doi: 10.1089/thy.2014.0581
49
UchidaTShimamuraMTakaHKagaNMiuraYNishidaYet al. The effect of long-term inorganic iodine on intrathyroidal iodothyronine content and gene expression in mice with graves' Hyperthyroidism. Thyroid (2023) 33(3):330–7. doi: 10.1089/thy.2022.0496
50
SatoSNohJYSatoSSuzukiMYasudaSMatsumotoMet al. Comparison of efficacy and adverse effects between methimazole 15 mg+Inorganic iodine 38 mg/day and methimazole 30 mg/day as initial therapy for graves' Disease patients with moderate to severe hyperthyroidism. Thyroid (2015) 25(1):43–50. doi: 10.1089/thy.2014.0084
51
SuzukiNYoshimura NohJSugisawaCHoshiyamaAHirumaMKawaguchiAet al. Therapeutic efficacy and limitations of potassium iodide for patients newly diagnosed with graves' Disease. Endocr J (2020) 67(6):631–8. doi: 10.1507/endocrj.EJ19-0379
52
GaoTSWeiJPLiuCLiXZZuoXHHengXPet al. Medicine Professional Committee of Endocrinology of Chinese Association of Integrative Medicine. Consensus on the treatment of graves' Disease with iodine-rich chinese medicine. Chin J Integrated Traditional Western Med (2021) 41(06):663–7.
53
ShiYGaoTYangL. Effect of iodine-rich chinese herbal compound on the thyroid function and morphology in hyperthyroidism wistar rats. J Liaoning Univ Traditional Chin Med (2009) 11(09):186–8. doi: 10.13194/j.jlunivtcm.2009.09.188.shiy.093
54
HongY. The efficacy and safety of iodine-rich nourishing-yin removing-fireand dissipating-stagnation decoction for the treatment of graves disease [硕士]: Liaoning university of traditional chinese medicine. (2019).
55
ZhangYShenJ. Clinical study on hyperthyroidism with chinese herbs for oral administration. World J Integrated Traditional Western Med (2017) 12(05):703–6. doi: 10.13935/j.cnki.sjzx.170526
56
ZhouFHeKGuanYYangXChenYSunMet al. Network pharmacology-based strategy to investigate pharmacological mechanisms of tinospora sinensis for treatment of alzheimer's disease. J Ethnopharmacol (2020) 259:112940. doi: 10.1016/j.jep.2020.112940
57
YoonC-SKimD-CLeeD-SKimK-SKoWSohnJHet al. Anti-neuroinflammatory effect of aurantiamide acetate from the marine fungus aspergillus sp. Sf-5921: Inhibition of nf-Κb and mapk pathways in lipopolysaccharide-induced mouse bv2 microglial cells. Int Immunopharmacol (2014) 23(2):568–74. doi: 10.1016/j.intimp.2014.10.006
58
ZhangM-YGuoJHuX-MZhaoS-QLiS-LWangJ. An in vivo anti-tumor effect of eckol from marine brown algae by improving the immune response. Food Funct (2019) 10(7):4361–71. doi: 10.1039/c9fo00865a
59
ZhangMZhouWZhaoSLiSYanDWangJ. Eckol inhibits reg3a-induced proliferation of human sw1990 pancreatic cancer cells. Exp Ther Med (2019) 18(4):2825–32. doi: 10.3892/etm.2019.7889
60
TahaASharifpanahFWartenbergMSauerH. Omega-3 and omega-6 polyunsaturated fatty acids stimulate vascular differentiation of mouse embryonic stem cells. J Cell Physiol (2020) 235(10):7094–106. doi: 10.1002/jcp.29606
61
DingYMullapudiBTorresCMascariñasEMancinelliGDiazAMet al. Omega-3 fatty acids prevent early pancreatic carcinogenesis via repression of the akt pathway. Nutrients (2018) 10(9):1289. doi: 10.3390/nu10091289
62
GomaaAMSAbd El-AzizEA. Omega-3 fatty acids decreases oxidative stress, tumor necrosis factor-alpha, and interleukin-1 beta in hyperthyroidism-induced hepatic dysfunction rat model. Pathophysiology (2016) 23(4):295–301. doi: 10.1016/j.pathophys.2016.10.001
63
JiangHLiJChenALiYXiaMGuoPet al. Fucosterol exhibits selective antitumor anticancer activity against hela human cervical cell line by inducing mitochondrial mediated apoptosis, cell cycle migration inhibition and downregulation of M-tor/pi3k/akt signalling pathway. Oncol Lett (2023) 25(1):32. doi: 10.3892/ol.2022.13618
64
LeeJ-HJungHAKangMJChoiJSKimG-D. Fucosterol, isolated from ecklonia stolonifera, inhibits adipogenesis through modulation of foxo1 pathway in 3t3-L1 adipocytes. J Pharm Pharmacol (2017) 69(3):325–33. doi: 10.1111/jphp.12684
65
SunZMohamedMAAParkSYYiTH. Fucosterol protects cobalt chloride induced inflammation by the inhibition of hypoxia-inducible factor through pi3k/akt pathway. Int Immunopharmacol (2015) 29(2):642–7. doi: 10.1016/j.intimp.2015.09.016
66
RaoASYadavSSSinghPNandalASinghNGanaieSAet al. A comprehensive review on ethnomedicine, phytochemistry, pharmacology, and toxicity of tephrosia purpurea (L.) pers. Phytother Res (2020) 34(8):1902–25. doi: 10.1002/ptr.6657
67
TuHMaDLuoYTangSLiYChenGet al. Quercetin alleviates chronic renal failure by targeting the pi3k/akt pathway. Bioengineered (2021) 12(1):6538–58. doi: 10.1080/21655979.2021.1973877
68
YaoXSaRYeCZhangDZhangSXiaHet al. Effects of thyroid hormone status on metabolic pathways of arachidonic acid in mice and humans: A targeted metabolomic approach. Prostaglandins Other Lipid Mediat (2015) 118-119:11–8. doi: 10.1016/j.prostaglandins.2015.03.005
69
WongLLLeeNGAmarnaniDChoiCJBielenbergDRFreitagSKet al. Orbital angiogenesis and lymphangiogenesis in thyroid eye disease: An analysis of vascular growth factors with clinical correlation. Ophthalmology (2016) 123(9):2028–36. doi: 10.1016/j.ophtha.2016.05.052
70
OkamotoMWatanabeMInoueNOgawaKHidakaYIwataniY. Gene polymorphisms of vegf and vegfr2 are associated with the severity of hashimoto's disease and the intractability of graves' Disease, respectively. Endocr J (2020) 67(5):545–59. doi: 10.1507/endocrj.EJ19-0480
71
YamadaEYamazakiKTakanoKObaraTSatoK. Iodide inhibits vascular endothelial growth factor-a expression in cultured human thyroid follicles: A microarray search for effects of thyrotropin and iodide on angiogenesis factors. Thyroid (2006) 16(6):545–54. doi: 10.1089/thy.2006.16.545
72
JangJYChoiSYParkIParkDYChoeKKimPet al. Vegfr2 but not vegfr3 governs integrity and remodeling of thyroid angiofollicular unit in normal state and during goitrogenesis. EMBO Mol Med (2017) 9(6):750–69. doi: 10.15252/emmm.201607341
73
RobertsWGPaladeGE. Increased microvascular permeability and endothelial fenestration induced by vascular endothelial growth factor. J Cell Sci (1995) 108(Pt 6):2369–79. doi: 10.1242/jcs.108.6.2369
74
SenouMCostaMJMassartCThimmeschMKhalifaCPoncinSet al. Role of caveolin-1 in thyroid phenotype, cell homeostasis, and hormone synthesis: In vivo study of caveolin-1 knockout mice. Am J Physiol Endocrinol Metab (2009) 297(2):E438–E51. doi: 10.1152/ajpendo.90784.2008
75
Van RegemorterEJorisVVan RegemorterVMariqueLBehetsCLengeléBet al. Downregulation of caveolin-1 and upregulation of deiodinase 3, associated with hypoxia-inducible factor-1Α Increase, are involved in the oxidative stress of graves' Orbital adipocytes. Thyroid (2021) 31(4):627–37. doi: 10.1089/thy.2020.0238
76
GuoKGaoRYuYZhangWYangYYangA. Quantitative mrna expression analysis of selected genes in patients with early-stage hypothyroidism induced by treatment with iodine-131. Mol Med Rep (2015) 12(5):7673–80. doi: 10.3892/mmr.2015.4350
77
LeiteJLBufaloNESantosRBRomaldiniJHWardLS. Herpesvirus type 7 infection may play an important role in individuals with a genetic profile of susceptibility to graves' Disease. Eur J Endocrinol (2010) 162(2):315–21. doi: 10.1530/EJE-09-0719
78
OmmaTYücelÇSertoğluEFıratSNÇulhaCÖzgürtaşT. The role of il-6 and osteoprotegerin in bone metabolism in patients with graves' Disease. Turk J Med Sci (2022) 52(2):338–45. doi: 10.55730/1300-0144.5320
79
ZhuPWuXZhouJWuKLuY. Gene polymorphisms of pro-inflammatory cytokines may affect the risk of graves' Disease: A meta-analysis. J Endocrinol Invest (2021) 44(2):311–9. doi: 10.1007/s40618-020-01300-x
80
LuYXingCZhangCLvXLiuGChenFet al. Promotion of il−17/nf−Κb signaling in autoimmune thyroid diseases. Exp Ther Med (2023) 25(1):51. doi: 10.3892/etm.2022.11750
81
FerrariSMRagusaFPaparoSRNasiniFNardiMFranceschiniSSet al. Differential modulation of cxcl8 versus cxcl10, by cytokines, ppar-gamma, or ppar-alpha agonists, in primary cells from graves' Disease and ophthalmopathy. Autoimmun Rev (2019) 18(7):673–8. doi: 10.1016/j.autrev.2019.05.004
82
ZhuQSuJWangXTangMGaoYZhangD. Serum concentrations of TNF-α and its soluble receptors in graves' Disease. Endocr Connect (2020) 9(7):736–46. doi: 10.1530/EC-20-0162
83
ChenGGVlantisACZengQvan HasseltCA. Regulation of cell growth by estrogen signaling and potential targets in thyroid cancer. Curr Cancer Drug Targets (2008) 8(5):367–77. doi: 10.2174/156800908785133150
84
MatsumotoCItoMYamadaHYamakawaNYoshidaHDateAet al. Genes that characterize T3-predominant graves' Thyroid tissues. Eur J Endocrinol (2013) 168(2):137–44. doi: 10.1530/EJE-12-0507
85
XuSChenGPengWRenkoKDerwahlM. Oestrogen action on thyroid progenitor cells: Relevant for the pathogenesis of thyroid nodules? J Endocrinol (2013) 218(1):125–33. doi: 10.1530/JOE-13-0029
86
MorshedSAMaRLatifRDaviesTF. Cleavage region thyrotropin receptor antibodies influence thyroid cell survival in vivo. Thyroid (2019) 29(7):993–1002. doi: 10.1089/thy.2018.0633
87
WoellerCFRoztocilEHammondCFeldonSE. Tshr signaling stimulates proliferation through pi3k/akt and induction of mir-146a and mir-155 in thyroid eye disease orbital fibroblasts. Invest Ophthalmol Vis Sci (2019) 60(13):4336–45. doi: 10.1167/iovs.19-27865
88
SuhJMSongJHKimDWKimHChungHKHwangJHet al. Regulation of the phosphatidylinositol 3-kinase, akt/protein kinase B, frap/mammalian target of rapamycin, and ribosomal S6 kinase 1 signaling pathways by thyroid-stimulating hormone (Tsh) and stimulating type tsh receptor antibodies in the thyroid gland. J Biol Chem (2003) 278(24):21960–71. doi: 10.1074/jbc.M300805200
89
BoutinAEliseevaEGershengornMCNeumannS. Β-arrestin-1 mediates thyrotropin-enhanced osteoblast differentiation. FASEB J (2014) 28(8):3446–55. doi: 10.1096/fj.14-251124
90
MasoodABenabdelkamelHEkhzaimyAAAlfaddaAA. Plasma-based proteomics profiling of patients with hyperthyroidism after antithyroid treatment. Molecules (2020) 25(12):2831. doi: 10.3390/molecules25122831
91
HersIVincentEETavaréJM. Akt signalling in health and disease. Cell Signal (2011) 23(10):1515–27. doi: 10.1016/j.cellsig.2011.05.004
92
IyerSBahnR. Immunopathogenesis of graves' Ophthalmopathy: The role of the tsh receptor. Best Pract Res Clin Endocrinol Metab (2012) 26(3):281–9. doi: 10.1016/j.beem.2011.10.003
93
KoJKimJ-YLeeEJYoonJS. Inhibitory effect of idelalisib, a selective phosphatidylinositol 3-kinase Δ Inhibitor, on adipogenesis in an in vitro model of graves' Orbitopathy. Invest Ophthalmol Vis Sci (2018) 59(11):4477–85. doi: 10.1167/iovs.18-24509
94
RaychaudhuriNFernandoRSmithTJ. Thyrotropin regulates il-6 expression in cd34+ Fibrocytes: Clear delineation of its camp-independent actions. PLoS One (2013) 8(9):e75100. doi: 10.1371/journal.pone.0075100
95
Szabo-FresnaisNBlondeauJ-PPoméranceM. Activation of the camp pathway synergistically increases il-1-induced il-6 gene expression in frtl-5 thyroid cells: Involvement of ap-1 transcription factors. Mol Cell Endocrinol (2008) 284(1-2):28–37. doi: 10.1016/j.mce.2007.12.017
96
HuYSunYHuangYLiuQRenF. Serum levels of cxcl-13, rbp-4, and il-6, and correlation analysis of patients with graves' Disease. Emerg Med Int (2022) 2022:5131846. doi: 10.1155/2022/5131846
97
LiKLiHXuWLiuWDuYHeJ-Fet al. Research on the potential mechanism of gypenosides on treating thyroid-associated ophthalmopathy based on network pharmacology. Med Sci Monit (2019) 25:4923–32. doi: 10.12659/MSM.917299
98
KimHLeeTHHwangYSBangMAKimKHSuhJMet al. Methimazole as an antioxidant and immunomodulator in thyroid cells: Mechanisms involving interferon-gamma signaling and H(2)O(2) scavenging. Mol Pharmacol (2001) 60(5):972–80. doi: 10.1124/mol.60.5.972
99
MorshedSAAndoTLatifRDaviesTF. Neutral antibodies to the tsh receptor are present in graves' Disease and regulate selective signaling cascades. Endocrinology (2010) 151(11):5537–49. doi: 10.1210/en.2010-0424
100
XieRLinZZhongCLiSChenBWuYet al. Deciphering the potential anti-covid-19 active ingredients in andrographis paniculata (Burm. F.) nees by combination of network pharmacology, molecular docking, and molecular dynamics. RSC Adv (2021) 11(58):36511–7. doi: 10.1039/d1ra06487h
101
FuT-TTuGPingMZhengG-XYangF-YYangJ-Yet al. Subtype-selective mechanisms of negative allosteric modulators binding to group I metabotropic glutamate receptors. Acta Pharmacol Sin (2021) 42(8):1354–67. doi: 10.1038/s41401-020-00541-z
Summary
Keywords
Haizao-Kunbu, Graves’ disease, network pharmacology, molecular docking, molecular dynamic analysis statements and declarations
Citation
Yang M, Lai Y, Gan D, Liu Q, Wang Y, He X, An Y and Gao T (2023) Possible molecular exploration of herbal pair Haizao-Kunbu in the treatment of Graves’ disease by network pharmacology, molecular docking, and molecular dynamic analysis. Front. Endocrinol. 14:1236549. doi: 10.3389/fendo.2023.1236549
Received
07 June 2023
Accepted
22 September 2023
Published
04 October 2023
Volume
14 - 2023
Edited by
Malgorzata Gabriela Wasniewska, University of Messina, Italy
Reviewed by
Silvia Martina Ferrari, University of Pisa, Italy; Guofang Chen, Nanjing University of Chinese Medicine, China
Updates
Copyright
© 2023 Yang, Lai, Gan, Liu, Wang, He, An and Gao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yiwen Lai, hermie1989@126.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.