- Open Access
Clonorchis sinensis MF6p/HDM (CsMF6p/HDM) induces pro-inflammatory immune response in RAW 264.7 macrophage cells via NF-κB-dependent MAPK pathways
Parasites & Vectors volume 13, Article number: 20 (2020)
MF6p/host defense molecules (HDMs) are a broad family of small proteins secreted by helminth parasites. Although the physiological role of MF6p/HDMs in trematode parasites is not fully understood, their potential biological function in maintaining heme homeostasis and modulating host immune response has been proposed.
A gene encoding the MF6p/HDM of Clonorchis sinensis (CsMF6p/HDM) was cloned. Recombinant CsMF6p/HDM (rCsMF6p/HDM) was expressed in Escherichia coli. The biochemical and immunological properties of rCsMF6/HDM were analyzed. CsMF6p/HDM induced pro-inflammatory response in RAW 264.7 cells was analyzed by cytokine array assay, reverse transcription polymerase chain reaction, and enzyme-linked immunosorbent assay. The structural feature of CsMF6p/HDM was analyzed by three-dimensional modeling and molecular docking simulations.
The CsMF6p/HDM shares a high level of amino acid sequence similarity with orthologs from other trematodes and is expressed in diverse developmental stages of the parasite. The rCsMF6p/HDM bound to bacteria-derived lipopolysaccharide (LPS), without effectively neutralizing LPS-induced inflammatory response in RAW 264.7 macrophage cells. Rather, the rCsMF6p/HDM induced pro-inflammatory immune response, which is characterized by the expression of TNF-α and IL-6, in RAW 264.7 cells. The rCsMF6p/HDM-induced pro-inflammatory immune response was regulated by JNK and p38 MAPKs, and was effectively down-regulated via inhibition of NF-κB. The structural analysis of CsMF6p/HDM and the docking simulation with LPS suggested insufficient capture of LPS by CsMF6p/HDM, which suggested that rCsMF6p/HDM could not effectively neutralize LPS-induced inflammatory response in RAW 264.7 cells.
Although rCsMF6p/HDM binds to LPS, the binding affinity may not be sufficient to maintain a stable complex of rCsMF6p/HDM and LPS. Moreover, the rCsMF6p/HDM-induced pro-inflammatory response is characterized by the release of IL-6 and TNF-α in RAW 264.7 macrophage cells. The pro-inflammatory response induced by rCsMF6p/HDM is mediated via NF-κB-dependent MAPK signaling pathway. These results collectively suggest that CsMF6p/HDM mediates C. sinensis-induced inflammation cascades that eventually lead to hepatobiliary diseases.
Clonorchiasis is a zoonotic food-borne parasitic disease with a serious impact on public health . It has been estimated that more than 15 million people worldwide suffer from clonorchiasis and over 200 million are at risk of infection . The disease is caused by a species of carcinogenic liver fluke, Clonorchis sinensis, which is widely distributed in Korea, China and northern Vietnam [1, 3, 4]. Humans are usually infected with the parasite by eating raw or undercooked freshwater fish containing metacercariae. Inside the bile ducts, the parasite induces mechanical damage to biliary epithelium, and generates an inflammatory responses and diverse clinical manifestations such as epithelial hyperplasia, periductal fibrosis and cholangiocarcinoma . Thus, C. sinensis has been regarded as a group I biological carcinogen by the International Agency for Research on Cancer of the World Health Organization .
The MF6p/host defense molecule (HDM) of C. sinensis (CsMF6p/HDM) is a small protein consisting of 90 amino acids and was initially reported as a 7 kDa antigen of unknown function [7, 8]. However, it has recently been classified under a new family of heme-binding proteins (MF6p/HDMs) with homologues from the other flukes including Fasciola hepatica (FhMF6p/HDM), Opisthorchis viverrini (OvMF6p/HDM) and Paragonimus westermani (PwMF6p/HDM) . The physiological role of MF6p/HDMs in trematode parasites has yet to be elucidated. However, several studies investigating FhMF6p/HDM suggested a role in host immune modulation [10,11,12,13]. Moreover, the increased interest in this molecule has been attributed to its heme-binding ability, and its role as a heme scavenger and transporter to maintain heme homeostasis in trematode parasites [9, 14]. The heme-scavenging ability is essential for the survival of blood-feeding trematodes since large amount of heme, which is toxic as a free form, are released from the catabolism of host erythrocytes . These findings collectively suggest that trematode MF6p/HDMs are attractive targets for the development of vaccine and therapeutic drugs.
In this study, the biochemical and immune modulation properties of CsMF6p/HDM were characterized. The recombinant CsMF6p/HDM (rCsMF6p/HDM) bound to lipopolysaccharide (LPS) but did not effectively neutralize LPS-induced pro-inflammatory responses in RAW 264.7 macrophage cells. Further, rCsMF6p/HDM induced pro-inflammatory immune responses in RAW 264.7 macrophage cells via NF-κB-dependent MAPK pathways.
Parasites and sera
Metacercariae of Clonorchis sinensis were collected from the naturally infected intermediate host, Pseudorasbora parva, collected in Korea. Sprague-Dawley rats were experimentally infected with 100 metacercariae by oral feeding. The rats were sacrificed at 2, 4, 6 and 9 weeks post-infection . Juvenile and adult worms were collected from the livers and washed five times with cold phosphate-buffered saline (PBS, pH 7.4) to remove any host contamination. The collected parasites were stored at − 70 °C until use, or used directly for RNA preparation to synthesize cDNA. The sera of infected rats were also collected at the time of sacrifice .
Expression and purification of recombinant CsMF6p/HDM (rCsMF6p/HDM)
The gene encoding CsMF6p/HDM (GenBank: AF281362.1) [7, 8] was amplified by polymerase chain reaction (PCR) from the cDNA of adult C. sinensis and cloned into T&A cloning vector (Real Biotech Corporation, Banqiao City, Taiwan). The nucleotide sequence of the cloned CsMF6p/HDM gene was confirmed by sequencing. To produce the recombinant CsMF6p/HDM (rCsMF6p/HDM), a partial CsMF6p/HDM without the N-terminal signal peptide region was amplified using the following primers; forward (5′-GGA TCC CGT CCC AGT GAG GAG ACC CGT-3′) and reverse (5′-AAG CTT TCA CTC CCC AAC GTA AGT CTC-3′) (restriction sites underlined). The purified PCR product was ligated into the T&A cloning vector (Real Biotech Corporation) and the ligate was transformed into Escherichia coli DH5α. The resulting plasmid DNA was digested with BamHI and HindIII, ligated into the pQE-9 expression vector (Qiagen, Hilden, Germany), and then transformed into E. coli M15 [pREP4] cells (Qiagen). Selected clone was grown and induced with 1 mM isopropyl-1-thio-β-d-galactopyranoside (IPTG). The bacteria were suspended in native lysis buffer (50 mM NaH2PO4, 300 mM NaCl, 10 mM imidazole, pH 8.0), sonicated on ice and centrifuged at 4 °C for 20 min at 12,000× rpm. The rCsMF6p/HDM was purified from the supernatant by nickel-nitrilotriacetic acid (Ni-NTA) resin (Qiagen) according to the manufacturer’s protocols. The purification and purity of the rCsMF6p/HDM was confirmed by 15% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE).
Production of polyclonal antibody for rCsMF6p/HDM (anti-CsMF6p/HDM)
The anti-CsMF6p/HDM was produced by immunizing BALB/c mice with the purified rCsMF6p/HDM. The protein (100 μg) was mixed with Freund’s complete or incomplete adjuvant (Sigma-Aldrich, St. Louis, MO, USA) and injected into mice intraperitoneally three times at 2-week intervals. Two weeks after the final booster, the mice were sacrificed and the sera were collected. The immunoglobulin G (IgG) fraction was further purified with the Protein G agarose (Pierce, Rockford, lL, USA) according to the manufacturer’s instructions. The specificity of the anti-CsMF6p/HDM was confirmed by immunoblot analysis.
Expression profile of CsMF6p/HDM in different developmental stages
The expression level of CsMF6p/HDM was investigated via semi-quantitative reverse transcription PCR (RT-PCR) and immunoblot analysis across different developmental stages of C. sinensis; metacercariae, 2-week-old juveniles, and 4-, 6- and 9-week-old adults. The total RNA was extracted from each developmental stage of worms using RNAiso (Takara, Ostu, Japan) according to the manufacturer’s instructions. The total RNA was quantified by spectrophotometry and equalized, and its purity was verified on a 1% agarose gel. Reverse transcription and subsequent amplification of CsMF6p/HDM transcripts were performed with the same amounts of total RNA (1 μg each) and specific primers using the ProtoScript® II RT-PCR Kit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instruction. PCR was performed using a primer set specific for CsMF6p/HDM (forward: 5′-CGT CCC AGT GAG GAG ACC CGT-3′, reverse: 5′-CAC TCC CCA ACG TAA GTC TC-3′). As an internal control, the C. sinensis actin (Cs actin) gene (GenBank: EU109284) was used as previously described . The amplicons were analyzed on a 2% agarose gel, visualized with RedSafeTM nucleic acid staining solution (iNtRON Biotechnology Inc., Seongnam, Korea) and observed under ultraviolet light. For immunoblot, worms of each developmental stage of C. sinensis were homogenized in PBS containing the CompleteTM protease inhibitor cocktails (Roche, Mannheim, Germany) and the supernatants were collected. The protein concentration in each worm extract was measured with Quick Start™ Bradford Protein Assay kit (Bio-Rad, Hercules, CA, USA) according to the manufacturer’s instructions. The same amount of worm extract (10 μg) from each developmental stage was separated on 15% SDS-PAGE and transferred to nitrocellulose membrane (Bio-Rad). The membrane was blocked with PBS supplemented with 0.05% Tween 20 (PBST, pH 7.4) and 3% skimmed milk for 1 h at room temperature. The membrane was then incubated with anti-CsMF6p/HDM diluted 1:1000 in PBST at room temperature for 2 h. After several washes with PBST, the membrane was incubated with 1:1000 diluted horseradish peroxidase (HRP)-conjugated anti-mouse IgG (Sigma-Aldrich). The immuno-reactive bands were visualized with 4-chloro-1-naphthol (Sigma-Aldrich) and the reaction was terminated by washing the membrane with distilled water.
Lipopolysaccharide (LPS) binding assay
A dot blot was used to evaluate the binding profile of rCsMF6p/HDM with LPS. Nitrocellulose membrane (Bio-Rad, pore size: 0.25 μm) was immersed in PBS for 10 min and fixed into a Bio-Dot microfiltration apparatus (Bio-Rad). LPS (E. coli strain O127:B8, Sigma-Aldrich) dissolved in PBS (10 μg/100 μl) was added to the wells of the apparatus and the sample was allowed to pass through the membrane under gravity. Serial dilutions of rCsMF6p/HDM (0 to 20 μg) in PBS were applied to each well and incubated for 1 h at room temperature. After incubation, the membrane was removed from the apparatus and blocked with PBST containing 3% skimmed milk for 1 h at room temperature. The membrane was incubated with anti-CsMF6p/HDM diluted 1:1000 in PBST at room temperature for 2 h. After several washes with PBST, the membrane was incubated with 1:1000 diluted HRP-conjugated anti-mouse IgG (Sigma-Aldrich). The immuno-reactive dots were visualized with 4-chloro-1-naphthol (Sigma-Aldrich) and the reaction was terminated by washing the membrane with distilled water.
Cultivation of RAW 264.7 cells and treatment with rCsMF6p/HDM
RAW 264.7 murine macrophage cells (ATCC® TIB-71TM) were purchased from American Type Culture Collection (Manassas, VA, USA) and maintained at 37 °C with 5% CO2 in Dulbecco’s modified Eagle’s medium (DMEM; Welgene Inc., Daegu, Korea) containing 10% heat-inactivated fetal bovine serum (FBS; Gibco, Grand Island, NY, USA), 1% l-glutamine (Welgene Inc.) and 1% penicillin/streptomycin (Gibco). To evaluate the potential cytotoxic effect of rCsMF6p/HDM to RAW 264.7 cells, the cell viability of RAW 264.7 cells after treatment with rCsMF6p/HDM was assessed using the methyl thiazolyl tetrazolium (MTT) method. To remove any LPS contamination in purified rCsMF6p/HDM, the Detoxi-gel endotoxin removing column (Pierce) was used according to the manufacturer’s instruction. No residual amount of endotoxin was detected in the sample when analyzed with Pierce™ LAL Chromogenic Endotoxin Quantitation Kit (Pierce). The LPS-depleted rCsMF6p/HDM was filtered with a syringe filter (0.22 μm; Millipore, Billerica, MA, USA) and used for further experiments. RAW 264.7 cells were continuously deprived of serum by incubation in 1% FBS overnight, followed by incubation in serum-free medium for 12 h. The serum-starved cells were seeded in a 96-well microplate (2 × 104 cells/well) and incubated for 6 h. The cells were treated with different concentrations of rCsMF6p/HDM (0 to 50 μg/ml) and incubated at 37 °C for 24 h. The supernatant was then drained from the plate and 100 μl of MTT (3-(4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazolium Bromide) (2 mg/ml; Sigma-Aldrich) were added in each well. The plate was further incubated for 3 h at 37 °C then followed by the addition of 100 μl of dimethyl sulfoxide (DMSO; Sigma-Aldrich) to each well. The plate was incubated for 30 min at 37 °C and the reaction was read at 595 nm with a Multiskan FC microplate reader (Thermo, Vantaa, Finland). All the experiments were performed in triplicate, and the mean and standard deviation (SD) were calculated. No cellular damage was observed when cells were treated with 20 μg/ml or less of rCsMF6p/HDM.
Cytokine array assay
To analyze the immune response of RAW 264.7 cells induced by rCsMF6p/HDM, the cells were seeded at ~ 70% confluence on a 6-well plate and grown for 24 h under standard culture conditions as described above. Cells were continuously deprived of serum as described above followed by the addition of rCsMF6p/HDM (10 μg/ml) and incubation for 12 h. For negative control, PBS-treated RAW 264.7 cells were used. The expression profiles of cytokines/chemokines in RAW 264.7 cells, which was stimulated by rCsMF6p/HDM, were analyzed compared with the control cells treated with PBS. The culture supernatant was harvested and centrifuged at 1200× rpm for 5 min at 4 °C. The supernatant was subsequently analyzed using the Proteome ProfilerTM Mouse Cytokine Array Panel A (R&D systems, Minneapolis, MN, USA) according to the manufacturer’s instructions.
LPS neutralization assay
To determine the ability of rCsMF6p/HDM to neutralize or modulate LPS-induced inflammatory immune response in RAW 264.7 cells, different concentrations of rCsMF6p/HDM (5 or 10 μg) were incubated with LPS (1 μg) in PBS for 2 h at room temperature. Each mixture was added to serum-starved RAW 264.7 cells cultured in a 6-well plate and incubated for 12 h. The cells were harvested and the expression profiles of IL-6 and TNF-α were analyzed by RT-PCR. Total RNA was extracted from the cells by using RNAiso (Takara) according to the manufacturer’s instructions. RT-PCR was performed using the same method described above. PCR was carried out using specific primer sets for the mouse IL-6, TNF-α and GAPDH genes: IL-6, forward primer (5′-CCG GAG AGG AGA CTT CAC AG-3′) and reverse primer (5′-GGA AAT TGG GGT AGG AAG GA-3′); TNF-α, forward primer (5′-CAT CCT CTC AAA ATT CGA GTG ACA-3′) and reverse primer (5′-TGG GAG TAG ACA AGG TAC AAC CC-3′); and GAPDH, forward primer (5′-ACC AGA GTC CAT GCC ATC AC-3′) and reverse primer (5′-CAC CAC CCT GTT GCT GTA GCC-3′). The amplicon sizes were confirmed as described above and each band was quantified using ImageJ v1.44 software .
Analysis of inflammatory immune responses in RAW 264.7 cells induced by rCsMF6p/HDM
To analyze the inflammatory response in RAW 264.7 cells induced by rCsMF6p/HDM, the effects of MAPK inhibitors on IL-6 and TNF-α production in the cells were analyzed. Inhibitors for p38 (SB203580), c-Jun N-terminal kinase (JNK) (SP600125), extracellular signal-regulated protein kinase (ERK) (U0126), NF-κB (MG132) and AP-1 (SR11302) were used in this study. All inhibitors were purchased from Calbiochem (San Diego, CA, USA). RAW 264.7 cells were seeded in 6-well dishes (2 × 105 cells/well) and cultured to approximately 70% confluence. After changing the media with fresh serum-free media, different concentrations of each inhibitor were added to the cells and incubated for 3 h. Next, rCsMF6p/HDM (10 µg/ml) was added to the cells pretreated with each inhibitor followed by incubation for an additional 12 h. RAW 264.7 cells treated with the same volume of PBS and LPS (1 µg/ml) were used as negative and positive controls, respectively. Cells treated with only rCsMF6p/HDM (10 µg/ml) were also included as a control. The cells were harvested and washed with PBS. The total RNA was isolated using the same protocols described above followed by RT-PCR for IL-6 and TNF-α as described above. The relative levels of IL-6 and TNF-α in each cell supernatant were also quantified by enzyme-linked immunosorbent assay (ELISA) using the mouse Quantikine IL-6 and TNF-α ELISA kits (R&D systems, Minneapolis, MN, USA) according to the manufacturer’s instructions.
Primary sequence similarity and features
Amino acid sequences of CsMF6p/HDM and homologues of Fasciola hepatica (FhMF6p/HDM), Opisthorchis viverrini (OvMF6p/HDM) and Paragonimus westermani (PwMF6p/HDM) were retrieved from the NCBI GenBank database . The signal peptide region was predicted by SignalP v4.0  and removed to develop a three-dimensional (3D) structure model. Multiple sequence alignments were generated using MAFFT v7.407  and displayed by Jalview . Protein physicochemical properties were predicted using the AA-Prop (http://www.biogem.org/tool/aa-prop).
Threading-based modeling and analysis of 3D structure
To obtain the most accurate structures of MF6p/HDM homologues, the top five 3D models of CsMF6p/HDM, FhMF6p/HDM, OvMF6p/HDM and PwMF6p/HDM were predicted using I-TASSER (Iterative Threading ASSEmbly Refinement) . The C-score, ranging from − 5 to 2, corresponds to a high confidence of the model. The structural similarity matrix was generated by measuring the pairwise structural similarities of 20 models based on the Dali Z-scores  using the Dali server . The Dali Z-score higher than 2 indicates fold similarity. Based on the similarity matrix, a structural dendrogram was derived using average linkage clustering . Both the backbone and the side chain of final 3D models were obtained by refinement using GalaxyRefine with the “both mild and aggressive relaxation” option . The quality of the models was evaluated using Ramachandran plot , ERRAT  and ProSA . A structural similarity heatmap image was generated with Plotly (https://plot.ly) and a structural similarity dendrogram was displayed using Evolview v2 . The overall structural superposition was performed using the TM-align . The TM-score value corresponds to the overall fold similarity. The TM-score ranges from 0 to 1, where 1 indicates 100% similarity between the two structures. Structure visualization was performed using PyMOL (http://www.pymol.org) and the hydrophobic surface was applied by running the python script “Color_h” (https://pymolwiki.org/index.php/Color_h).
Molecular docking simulation
The PDB file format was corrected using PDBEditor . LPS (Compound ID 11970143) was retrieved as sdf format from PubChem database  and the format was converted to mol2 format using Open Babel . Molecular docking simulations of MF6p/HDM homologues and LPS were conducted using AutoDock Vina  in DockingServer . A box size of 22 × 40 × 24 Å and box center of 51.6 × 51.5 × 51.3 for CsMF6p/HDM, and a box size of 27 × 34 × 31 Å and box center of 50.2 × 50.1 × 50.3 for FhMF6p/HDM were used for each simulation. Each docking simulation was performed using a setup of 99 different runs (‘ga_run’), 50,000,000 energy evaluations (‘ga_num_evals’) and 540,000 generations (‘ga_num_generation’) in order to terminate when either limit is reached. The remaining options were set as default values. A plot was generated using DataGraph v3.0 (Visual Data Tools Inc., Chapel Hill, NC, USA).
Sequential conservation of trematode MF6p/HDM homologues
CsMF6p/HDM showed the highest sequence identity (82.2%) with OvMF6p/HDM, followed by FhMF6p/HDM (56.6%) and PwMF6p/HDM (50.0%) (Fig. 1). The C-terminal conserved motif, 65LGxKIxxVxxILxxRLTxRxExY87 (x = any), which was reported as potential heme- and LPS-binding site , was conserved in all MF6p/HDM homologues. A putative signal peptide sequence was identified in the N-terminal region of CsMF6p/HDM.
Expression of rCsMF6p/HDM and production of anti-CsMF6p/HDM
A soluble form of rCsMF6p/HDM was expressed in E. coli. The approximate molecular mass of the rCsMF6p/HDM was 9 kDa, consistent with the estimated molecular mass of the deduced amino acid sequence of CsMF6p/HDM with 6 histidine residues (Fig. 2a). A polyclonal antibody against rCsMF6p/HDM (anti-CsMF6p/HDM) was successfully produced in mice immunized with rCsMF6p/HDM and the specificity of the antibody was confirmed by immunoblot (Fig. 2b).
Expression profile of CsMF6p/HDM at the developmental stages
The CsMF6p/HDM mRNA was expressed in both juvenile and adult worms, but not in the metacercariae. The mRNA level increased gradually with the maturation of the worms from juveniles to adults (Fig. 3a). The expression of CsMF6p/HDM protein was also analyzed using immunoblot with the anti-CsMF6p/HDM and C. sinensis worm lysates. CsMF6p/HDM was recognized only in juvenile and adult worm extracts and its expression was increased with maturation, consistent with the mRNA expression pattern in the developmental stages (Fig. 3b).
CsMF6p/HDM binds to LPS
To evaluate LPS-binding ability of rCsMF6p/HDM, a dot blot assay was performed. The rCsMF6p/HDM did not show positive response in the absence of LPS, but showed a positive response in the presence of LPS in a dose-dependent manner, suggesting that rCsMF6p/HDM bound to LPS (Fig. 4).
Overexpression of IL-6 and TNF-α in RAW 264.7 cells stimulated by rCsMF6p/HDM
Expression profiles of diverse cytokines and chemokines in rCsMF6p/HDM-treated and untreated RAW 264.7 cells were comparatively analyzed. The expression of several cytokines and chemokines was increased in RAW 264.7 cells stimulated with rCsMF6p/HDM (Fig. 5a). Compared with negative control (PBS-treated), the levels of several cytokines and chemokines including interferon gamma-induced protein 10 (IP-10), macrophage inflammatory protein 2 (MIP-2), tumor necrosis factor-α (TNF-α), interleukin-1 receptor antagonist (IL-1ra) and interleukin 6 (IL-6) were increased in RAW 264.7 cells stimulated with rCsMF6p/HDM (Fig. 5b). These results suggest that rCsMF6p/HDM induces the synthesis of various pro-inflammatory cytokines and chemokines in RAW 264.7 cells.
rCsMF6p/HDM does not neutralize LPS-mediated pro-inflammatory responses in RAW 264.7 cells
To determine the ability of rCsMF6p/HDM to neutralize LPS-induced pro-inflammatory responses, LPS incubated with different concentrations of rCsMF6p/HDM was treated to RAW 264.7 cells and expression profiles of IL-6 and TNF-α in the cells were analyzed. Compared with LPS control, the expression of IL-6 and TNF-α in RAW 264.7 cells was not affected by rCsMF6p/HDM, indicating that rCsMF6p/HDM did not effectively neutralize the pro-inflammatory effect of LPS in RAW 264.7 cells (Fig. 6). Further, increased expression of IL-6 and TNF-α was observed in RAW 264.7 cells treated with only rCsMF6p/HDM, suggesting that CsMF6p/HDM induces pro-inflammatory responses in RAW 264.7 cells (Fig. 6).
rCsMF6p/HDM induces IL-6 and TNF-α expression via MAPK signaling pathways in RAW 264.7 cells
To investigate the signaling pathways associated with rCsMF6p/HDM-induced IL-6 and TNF-α production, the effects of MAPKs inhibitors on rCsMF6p/HDM-mediated expression of IL-6 and TNF-α in RAW 264.7 cells were analyzed. The cells were pretreated with different concentrations (5 or 10 μM) of each inhibitor, i.e. p38 (SB203580), JNK (SP600125) or ERK (U0126), followed by stimulation with rCsMF6p/HDM. The RT-PCR analysis indicated that the inhibition of p38 and JNK impaired the ability of rCsMF6p/HDM to induce IL-6 expression in RAW 264.7 cells (Fig. 7a). Inhibition of ERK also reduced the expression of IL-6, but less than the inhibitors for p38 and JNK. Treatment with MAPK inhibitors also reduced the expression of TNF-α in RAW 264.7 cells; however, only the high inhibitor concentration (10 μM) attenuated the TNF-α expression in the cells (Fig. 7a). The effects of MAPK inhibitors on the expression of IL-6 and TNF-α were further analyzed by ELISA. The rCsMF6p/HDM-induced IL-6 and TNF-α protein expressions in RAW 264.7 cells were reduced by treatment with the inhibitors, and were consistent with RT-PCR results (Fig. 7b).
rCsMF6p/HDM-induced IL-6 and TNF-α expression is regulated via NF-κB signaling pathway in RAW 264.7 cells
To examine whether the enhanced expression of IL-6 and TNF-α in RAW 264.7 cells stimulated by rCsMF6p/HDM was mediated by NF-κB or AP-1 activation, the effect of NF-κB and AP-1 inhibitors on IL-6 and TNF-α expression was analyzed. The cells were pre-treated with serial concentrations (0.1, 0.5 and 1.0 μM) of inhibitors of NF-κB (MG132) or AP-1 (SR11302) and subsequently stimulated with CsMF6p/HDM. RT-PCR analysis revealed that MG132 abrogated the ability of rCsMF6p/HDM to induce IL-6 and TNF-α expression in RAW 264.7 cells in a dose-dependent manner; however, the expression of IL-6 and TNF-α was not greatly affected by SR11302 (Fig. 8a). Consistent with the RT-PCR analysis, the levels of IL-6 and TNF-α protein in rCsMF6p/HDM-stimulated RAW 264.7 cells were strongly inhibited by MG132, but not by SR11302 (Fig. 8b).
The most appropriate structures of MF6p/HDMs within liver flukes
To elucidate the structural characteristics of trematode MF6p/HDM homologues, the top five 3D models of CsMF6p/HDM, FhMF6p/HDM, OvMF6p/HDM and PwMF6p/HDM were built and comparatively analyzed. A structural similarity dendrogram and a 20 × 20 matrix of pairwise similarities were computed for a set of 20 structures (Additional file 1: Figure S1). Among them, CsMF6p/HDM_5, FhMF6p/HDM_1, OvMF6p/HDM_5 and PwMF6p/HDM_1 were selected as the initial models according to high structural similarity based on pairwise structural comparisons across four species. All of the selected models showed the highest C-score among their individual predicted models: CsMF6p/HDM_5 (− 1.70), FhMF6p/HDM_1 (− 1.03), OvMF6p/HDM_1 (− 1.54) and PwMF6p/HDM_1 (− 1.08). The selected models were further refined and determined as the final models. For example, the final 3D model of CsMF6p/HDM was obtained by sequential refinement (Additional file 2: Figure S2, Additional file 3: Figure S3). Superposition between CsMF6p/HDM and FhMF6p/HDM yielded a TM-score of 0.83, revealing a high structural resemblance of the two molecules (Fig. 9a–c), despite varying hydrophobic surfaces (Fig. 9d, e).
Preferable docking poses of LPS to CsMF6p/HDM and FhMF6p/HDM
To analyze the patterns of preferable LPS binding to CsMF6p/HDM and FhMF6p/HDM, 99 docking simulations and frequency analyses were performed for each binding complex. CsMF6p/HDM and FhMF6p/HDM were predicted to bind LPS with multiple and unique poses with low frequency (1–2%), resulting in a set of 99 binding poses (Additional file 4: Table S1, Additional file 5: Table S2). In the FhMF6p/HDM, a total of 25 poses out of 99 binding poses (mean binding energy of − 3.26E+04 kcal/mol) were predicted. However, no favorable binding affinity was observed for the CsMF6p/HDM. LPS was predicted to bind differently with CsMF6p/HDM and FhMF6p/HDM when all the binding residues of CsMF6p/HDM-LPS complex were normalized to those of FhMF6p/HDM-LPS complex (Additional file 6: Figure S4). An abundance of residues involved in LPS binding were predicted at the N- and C-terminal regions of CsMF6p/HDM, but in the middle and C-terminal regions of FhMF6p/HDM. When the interactions of LPS with CsMF6p/HDM and FhMF6p/HDM were optimized based on their estimated binding affinities, hydrogen bonding interactions were predicted for both proteins (Additional file 7: Figure S5). The optimized binding pose of CsMF6p/HDM with LPS suggested hydrogen bonding interactions between LPS and amino acid residues such as Arg27, Arg31, Lys36, Trp38, Lys68 and Arg79 in the N-terminal region of CsMF6p/HDM (Additional file 7: Figure S5a), but suggested unfavorable binding affinity (9.54E+04 kcal/mol) (Additional file 7: Figure S5c). The optimum binding pose of FhMF6p/HDM with LPS revealed hydrogen bonding between LPS and Arg53 of FhMF6p/HDM but other interactions (polar and hydrophobic bonds) were also predicted in the C-terminal region (Additional file 7: Figure S5b), resulting in highly favorable binding affinity (− 5.69E+04 kcal/mol) (Additional file 7: Figure S5c).
MF6p/HDMs are a broad family of small proteins secreted by helminth parasites . Although the physiological roles of MF6p/HDMs in trematode parasites have not been clearly established, several studies suggested two interesting biological functions of this family of proteins. The molecules act as heme scavengers and transporters to maintain heme homeostasis in trematode parasites, indicating an essential role in parasite survival by detoxifying the toxic free form [9, 14]. Their ability to bind to LPS and prevent LPS-induced inflammation in mice as well as their ability to suppress antigen processing and presentation in macrophages also suggested an immunomodulatory role in suppressing host immune response, which facilitated prolonged parasite survival in the hosts [10,11,12,13, 36, 37].
Here, the biochemical and immunological properties of CsMF6p/HDM were investigated to obtain further insight into the nature of the protein. CsMF6p/HDM is a protein that shares a high level of amino acid sequence similarity with the orthologs from other trematodes. In particular, the C-terminal motif, 65LGxKIxxVxxILxxRLTxRxExY87 (x = any), which has been reported as a potential heme- and LPS-binding site  was well conserved in CsMF6p/HDM. Expression of this protein in various developmental stages of C. sinensis, ranging from juveniles to adults, suggests an essential biological role in the parasitic life-cycle in mammalian hosts. FhMF6p/HDM, the most extensively studied MF6p/HDM, has attracted attention due to its immune modulatory effect against LPS-induced inflammation by preventing the activation of macrophages . Although the C-terminal motif of CsMF6p/HDM showed a high sequence identity (61.8%) with that of FhMF6p/HDM, which was experimentally confirmed to bind with LPS effectively , the rCsMF6p/HDM did not neutralize the pro-inflammatory activity of LPS in RAW 264.7 macrophage cells. Rather, rCsMF6p/HDM induced pro-inflammatory immune responses, characterized by the release of IL-6 and TNF-α (cytokines which play key role in orchestrating the inflammation responses), in RAW 264.7 macrophage cells. The pro-inflammatory responses induced by rCsMF6p/HDM are mediated via the NF-κB-dependent MAPK signaling pathway.
A simple dot blot analysis using rCsMF6p/HDM and LPS revealed that LPS bound to the recombinant protein, which was not consistent with a previous study suggesting that CsMF6p/HDM may not bind to LPS effectively . This disagreement may be attributed to the use of only a small fragment of the C-terminal region of CsMF6p/HDM, and not the full-length protein, in the previous study. The structural integrity of the C-terminal region in the MF6p/HDM has been considered essential for interaction with LPS ; however, this interaction may not be universal across similar peptides derived from MF6p/HDM orthologs in other trematodes. Otherwise, this phenomenon might be attributed to differences in the physicochemical properties of CsMF6p/HDM and FhMF6p/HDM, such as isoelectric point , hydrophobicity  and LPS serotypes (O127:B8 for CsMF6p/HDM; LPS O111:B4 and O55:B5 for FhMF6p/HDM) . To understand the different patterns of CsMF6p/HDM and FhMF6p/HDM in LPS binding and neutralization of LPS-induced inflammatory responses, the docking poses of LPS with CsMF6p/HDM and FhMF6p/HDM were compared. Since the PDB template of MF6p/HDM for 3D homology modeling was unavailable in the present study, a combined method was used to build appropriate 3D models with high quality [39, 40]. A combined approach involving threading-based modeling and all-against-all comparison within the MF6p/HDM subfamily proteins was used to develop the most reliable 3D model of CsMF6p/HDM. Remarkably, although each final 3D model was selected based on the Dali Z-scores via pairwise structural comparison, each 3D model showed the highest C-score value. The 3D model of MF6p/HDM homologues was predicted using an ab initio modeling method, QUARK , previously . Similar structures were also obtained in this study, but the models were not selected due to lower quality and similarity across different species of MF6p/HDM compared with the final structures determined in this study. Based on the final 3D models obtained, the best poses of LPS-CsMF6p/HDM binding were predicted by docking simulation analysis. The N-terminal region of CsMF6p/HDM was most probably involved in binding with LPS despite the unfavorable binding energies. Meanwhile, the C-terminal region of FhMF6p/HDM appears to be a critical region for LPS binding considering that this region showed high affinity binding energies in accordance with a previous study . Moreover, the CsMF6p/HDM and FhMF6p/HDM showed different patterns of binding to LPS. These results collectively suggest that FhMF6p/HDM bind to LPS more preferably compared with CsMF6p/HDM, which was consistent with previous results . Low-affinity binding between CsMF6p/HDM and LPS may result in insufficient capture of LPS by rCsMF6p/HDM. Therefore, the rCsMF6p/HDM failed to effectively neutralize LPS-induced inflammatory response in RAW 264.7 cells.
Interestingly, the rCsMF6p/HDM induced pro-inflammatory responses, which was mediated via NF-κB-dependent MAPK signaling pathway. The pro-inflammatory immune responses in diverse types of mammalian cells induced by C. sinensis excretory and secretory products (ESP) and proteins were reported previously [42,43,44,45]. CsMF6p/HDM is a protein largely identified in the ESP of C. sinensis, which indicates that the protein is actively released from the parasite and continuously interacts with the biliary epithelium in bile ducts. Considering that pro-inflammatory immune responses are induced in RAW 246.7 macrophage cells by rCsMF6p/HDM, this protein may mediate continuous inflammatory response in bile epithelial cells and may be associated with inflammatory events in clonorchiasis.
CsMF6p/HDM is a protein belonging to the MF6p/HDM protein family found in trematode parasites. The expression pattern of CsMF6p/HDM in different developmental stages of C. sinensis suggests its key role in parasite physiology and interactions with host. Although rCsMF6p/HDM binds to LPS, the binding affinity may not be sufficient to maintain a stable complex of rCsMF6p/HDM and LPS. Therefore, rCsMF6p/HDM failed to neutralize LPS-induced inflammatory response in RAW 246.7 cells. Moreover, the rCsMF6p/HDM-induced pro-inflammatory response is characterized by the release of IL-6 and TNF-α in RAW 264.7 macrophage cells. The pro-inflammatory response induced by rCsMF6p/HDM is mediated via NF-κB-dependent MAPK signaling pathway. These results collectively suggest that CsMF6p/HDM mediates C. sinensis-induced inflammation cascades that eventually lead to hepatobiliary diseases.
Availability of data and materials
The data supporting the conclusions of this article are provided within the article and its additional files. The original datasets analyzed in the present study are available from the corresponding author upon request.
polyclonal antibody against CsMF6p/HDM
MF6p/HDM of Clonorchis sinensis
enzyme-linked immunosorbent assay
extracellular signal-regulated kinase
excretory and secretory products
granulocyte-colony stimulating factor
glyceraldehyde 3-phosphate dehydrogenase
MF6p/host defense molecule
interleukin-1 receptor antagonist
interferon gamma-induced protein 10
Iterative Threading ASSEmbly Refinement
c-Jun N-terminal kinase
mitogen-activated protein kinases
macrophage inflammatory protein 2
nuclear factor-kappa B
regulated upon activation normal T cell expressed and secreted
soluble intercellular adhesion molecule-1
tumor necrosis factor-α
Rim HJ. Clonorchiasis: an update. J Helminthol. 2005;79:269–81.
Furst T, Keiser J, Utzinger J. Global burden of human food-borne trematodiasis: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12:210–21.
Wang KX, Zhang RB, Cui YB, Tian Y, Cai R, Li CP. Clinical and epidemiological features of patients with clonorchiasis. World J Gastroenterol. 2004;10:446–8.
Yu SH, Kawanaka M, Li XM, Xu LQ, Lan CG, Rui L. Epidemiological investigation on Clonorchis sinensis in human population in an area of South China. Jpn J Infect Dis. 2003;56:168–71.
Kim TS, Pak JH, Kim JB, Bahk YY. Clonorchis sinensis, an oriental liver fluke, as a human biological agent of cholangiocarcinoma: a brief review. BMB Rep. 2016;49:590–7.
Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F, et al. A review of human carcinogens—part B: biological agents. Lancet Oncol. 2009;10:321–2.
Lee HJ, Lee CS, Kim BS, Joo KH, Lee JS, Kim TS, et al. Purification and characterization of a 7-kDa protein from Clonorchis sinensis adult worms. J Parasitol. 2002;88:499–504.
Zhao QP, Moon SU, Lee HW, Na BK, Cho SY, Kong Y, et al. Evaluation of Clonorchis sinensis recombinant 7-kilodalton antigen for serodiagnosis of clonorchiasis. Clin Diagn Lab Immunol. 2004;11:814–7.
Martínez-Sernández V, Mezo M, González-Warleta M, Perteguer MJ, Gárate T, Romarís F, et al. Delineating distinct heme-scavenging and-binding functions of domains in MF6p/helminth defense molecule (HDM) proteins from parasitic flatworms. J Biol Chem. 2017;292:8667–82.
Robinson MW, Donnelly S, Hutchinson AT, To J, Taylor NL, Norton RS, et al. A family of helminth molecules that modulate Innate cell responses via molecular mimicry of host antimicrobial peptides. PLoS Pathog. 2011;7:e1002042.
Robinson MW, Alvarado R, To J, Hutchinson AT, Dowdell SN, Lund M, et al. A helminth cathelicidin-like protein suppresses antigen processing and presentation in macrophages via inhibition of lysosomal vATPase. FASEB J. 2012;26:4614–27.
Alvarado R, O’Brien B, Tanaka A, Dalton JP, Donnelly S. A parasitic helminth-derived peptide that targets the macrophage lysosome is a novel therapeutic option for autoimmune disease. Immunobiology. 2015;220:262–9.
Martínez-Sernández V, Perteguer MJ, Mezo M, González-Warleta M, Gárate T, Valero MA, et al. Fasciola spp.: mapping of the MF6 epitope and antigenic analysis of the MF6p/HDM family of heme-binding proteins. PLoS One. 2017;12:e0188520.
Martínez-Sernández V, Mezo M, González-Warleta M, Perteguer MJ, Muiño L, Guitián E, et al. The MF6p/FhHDM-1 major antigen secreted by the trematode parasite Fasciola hepatica is a heme-binding protein. J Biol Chem. 2014;289:1441–56.
Toh SQ, Glanfield A, Gobert GN, Jones MK. Heme and blood-feeding parasites: friends or foes? Parasit Vectors. 2010;3:108.
Kang JM, Lee J, Ju HL, Ju JW, Kim JH, Pak JH, et al. Characterization of a gut-associated asparaginyl endopeptidase of Clonorchis sinensis. Exp Parasitol. 2015;153:81–90.
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.
Jenuth JP. The NCBI. Publicly available tools and resources on the Web. Methods Mol Biol. 2000;132:301–12.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview Version 2-a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25:1189–91.
Yang J, Zhang Y. I-TASSER server: new development for protein structure and function predictions. Nucleic Acids Res. 2015;43:W174–81.
Holm L, Kaariainen S, Rosenstrom P, Schenkel A. Searching protein structure databases with DaliLite vol 3. Bioinformatics. 2008;24:2780–1.
Holm L, Laakso LM. Dali server update. Nucleic Acids Res. 2016;44:W351–5.
Heo L, Park H, Seok C. GalaxyRefine: protein structure refinement driven by side-chain repacking. Nucleic Acids Res. 2013;41:W384–8.
Lovell SC, Davis IW, Arendall WB 3rd, de Bakker PI, Word JM, Prisant MG, et al. Structure validation by Calpha geometry: phi, psi and Cbeta deviation. Proteins. 2003;50:437–50.
Colovos C, Yeates TO. Verification of protein structures: patterns of nonbonded atomic interactions. Protein Sci. 1993;2:1511–9.
Wiederstein M, Sippl MJ. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007;35:W407–10.
He Z, Zhang H, Gao S, Lercher MJ, Chen WH, Hu S. Evolview v2: an online visualization and management tool for customized and annotated phylogenetic trees. Nucleic Acids Res. 2016;44:W236–41.
Zhang Y, Skolnick J. TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005;33:2302–9.
Lee J, Kim SH. PDB editor: a user-friendly Java-based protein data bank file editor with a GUI. Acta Crystallogr D Biol Crystallogr. 2009;65:399–402.
Kim S, Thiessen PA, Bolton EE, Chen J, Fu G, Gindulyte A, et al. PubChem substance and compound databases. Nucleic Acids Res. 2016;44:D1202–13.
OʼBoyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR. Open Babel: an open chemical toolbox. J Cheminform. 2011;3:33.
Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31:455–61.
Bikadi Z, Hazai E. Application of the PM6 semi-empirical method to modeling proteins enhances docking accuracy of AutoDock. J Cheminform. 2009;1:15.
Alvarado R, To J, Lund ME, Pinar A, Mansell A, Robinson MW, et al. The immune modulatory peptide FhHDM-1 secreted by the helminth Fasciola hepatica prevents NLRP3 inflammasome activation by inhibiting endolysosomal acidification in macrophages. FASEB J. 2017;31:85–95.
Thivierge K, Cotton S, Schaefer DA, Riggs MW, To J, Lund ME, et al. Cathelicidin-like helminth defence molecules (HDMs): absence of cytotoxic, anti-microbial and anti-protozoan activities imply a specific adaptation to immune modulation. PLoS Negl Trop Dis. 2013;7:e2307.
Martinez-Sernandez V, Orbegozo-Medina RA, Romaris F, Paniagua E, Ubeira FM. Usefulness of ELISA methods for assessing LPS interactions with proteins and peptides. PLoS One. 2016;11:e0156530.
Kim YJ, Yoo WG, Lee MR, Kang JM, Na BK, Cho SH, et al. Molecular and structural characterization of the tegumental 20.6-kDa protein in Clonorchis sinensis as a potential druggable target. Int J Mol Sci. 2017;18:557.
Lu Y, Yoo WG, Dai F, Lee JY, Pak JH, Sohn WM, et al. Characterization of a novel organic solute transporter homologue from Clonorchis sinensis. PLoS Negl Trop Dis. 2018;12:e0006459.
Xu D, Zhang Y. Toward optimal fragment generations for ab initio protein structure assembly. Proteins. 2013;81:229–39.
Nam JH, Moon JH, Kim IK, Lee MR, Hong SJ, Ahn JH, et al. Free radicals enzymatically triggered by Clonorchis sinensis excretory-secretory products cause NF-kappaB-mediated inflammation in human cholangiocarcinoma cells. Int J Parasitol. 2012;42:103–13.
Pak JH, Son WC, Seo SB, Hong SJ, Sohn WM, Na BK, et al. Peroxiredoxin 6 expression is inversely correlated with nuclear factor-kappaB activation during Clonorchis sinensis infestation. Free Radic Biol Med. 2016;99:273–85.
Yan C, Wang YH, Yu Q, Cheng XD, Zhang BB, Li B, et al. Clonorchis sinensis excretory/secretory products promote the secretion of TNF-alpha in the mouse intrahepatic biliary epithelial cells via Toll-like receptor 4. Parasit Vectors. 2015;8:559.
Kim EM, Kwak YS, Yi MH, Kim JY, Sohn WM, Yong TS. Clonorchis sinensis antigens alter hepatic macrophage polarization in vitro and in vivo. PLoS Negl Trop Dis. 2017;11:e0005614.
This study was supported by a grant of the National Research Foundation (NRF) funded by the Korean government (NRF2016R1C1B1009348).
This work was supported by a grant of the National Research Foundation (NRF) funded by the Korean government (NRF2016R1C1B1009348).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
All-against-all comparison of CsMF6p/HDM and homologues.
Additional file 2: Figure S2.
Quality verification of initial 3D model of CsMF6p/HDM.
Additional file 3: Figure S3.
Quality verification of final 3D model of CsMF6p/HDM.
Additional file 4: Table S1.
Docking results of CsMF6p/HDM and LPS (99 simulations).
Additional file 5: Table S2.
Docking results of FhMF6p/HDM and LPS (99 simulations).
Additional file 6: Figure S4.
Difference of binding residues of LPS to between CsMF6p/HDM and FhMF6p/HDM.
Additional file 7: Figure S5.
Interactions of LPS between CsMF6p/HDM and FhMF6p/HDM.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Kang, JM., Yoo, W.G., Lê, H.G. et al. Clonorchis sinensis MF6p/HDM (CsMF6p/HDM) induces pro-inflammatory immune response in RAW 264.7 macrophage cells via NF-κB-dependent MAPK pathways. Parasites Vectors 13, 20 (2020). https://doi.org/10.1186/s13071-020-3882-0
- Clonorchis sinensis
- MF6p/host defense molecule
- Pro-inflammatory immune response