- 1Wannan Medical College, Wuhu, China
- 2Nanjing University of Chinese Medicine, Nanjing, China
Objective: Lianhua Qingwen capsule (LHQW) can attenuate lung injury caused by influenza virus infection. However, it is unclear whether the intestinal microbiota plays a role in LHQW activity in ameliorating viral infectious pneumonia. This study aimed to investigate the role of intestinal microbiota in LHQW activity in ameliorating viral infectious pneumonia and its possible mechanisms.
Research design and methods: A mouse model of influenza A viral pneumonia was established by intranasal administration in BALB/c mice. Detection of influenza virus in the lungs, pathological examination of the lungs and small intestine, and biochemical detection of inflammatory indices were performed. The effects of LHQW on intestinal microbiota were evaluated by 16S rRNA gene sequencing. The key components and targets of LHQW were screened via network pharmacology and verified through molecular docking, molecular dynamics simulation, and free binding energy calculations.
Results: Body weight decreased, inflammatory factor levels were disturbed, and the lung and intestinal mucosal barriers were significantly injured in the infected group. The alpha diversity of the intestinal microbiota decreased, and the abundance of Bacteroidetes, Muribaculaceae_unclassified, and Streptococcus decreased significantly. LHQW treatment reduced the viral load in the lungs, rescued body weight and survival, alleviated lung and intestinal mucosal barrier injury, reversed the reduction in the intestinal microbiota alpha diversity, and significantly increased the abundance of Bacteroidetes and Muribaculaceae. Network pharmacological analysis showed that six active herbal medicinal compounds from LHQW could regulate the intestinal microbiota and inhibit the immune-inflammatory response through the Toll-like receptor (TLR) and nuclear factor-κB (NF-κB) signalling pathways in the lungs.
Conclusion: These results suggest that LHQW is effective for treating influenza A virus infectious pneumonia, and the mechanism is associated with the regulation of the TLR4/NF-κB signalling pathway in the lungs by restoring intestinal microbiota and repairing the intestinal wall.
Introduction
In the past two decades, there have been several respiratory virus pandemics, which seriously threaten human health and negatively impact on the national economy and social life worldwide. Influenza viruses circulating in humans are mainly of two subtypes: type A (HA), H1N1, and H3N2. Human seasonal influenza virus A can generate new strains via antigenic drift. The new strain overcomes existing immunity in humans, leading to a new influenza epidemic (Caton et al., 1982; Treanor, 2004; Boni, 2008). Another significant way of influenza virus evolution leading to pandemics is the antigenic shift, in which the viral genome is rearranged to produce new subtypes with significant antigenic changes (Webster et al., 1982; Zhang et al., 2017). If the new influenza virus spreads efficiently from person to person, it could lead to a global pandemic. Four classes of antiviral drugs are currently approved worldwide for the treatment of influenza: amamantane derivatives, neuraminidase inhibitors, membrane fusion inhibitors, and RNA-dependent RNA polymerase inhibitors (Li et al., 2022). Due to mutations and rearrangements in the influenza virus genome, the antiviral effect of these drugs is also decreasing (Paules and Subbarao, 2017); therefore, it is necessary to explore new therapeutic strategies.
Lianhua Qingwen capsule (LHQW), a traditional Chinese medicine compound, has complex components and various effects (Liang et al., 2020; Huang et al., 2021; Zhou et al., 2021). The LHQW was first used clinically in 2003. Since then, it has been used to treat atypical pneumonia, severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), influenza A H1N1 virus, influenza A H3N2 virus, and influenza H7N9 virus infections (Xing and Liu, 2021). In 2020, the LHQW capsule was recommended for COVID-19 treatment in China (Ni et al., 2020; Wu et al., 2020; Zhuang et al., 2020; Li et al., 2021; Hu et al., 2022). Literature and clinical evidence have indicated that LHQW, which ameliorates symptoms of influenza, regulates the expression of cytokines after viral infection, alleviates lung injury, and has protective effects against COVID-19 (Liang et al., 2020; Huang et al., 2021; Shen and Yin, 2021; Hu et al., 2022). LHQW, combined with routine treatment, can be used as synergistic therapy to significantly improve symptoms and reduce the mortality rate of critically ill patients (Wu et al., 2021; Yang et al., 2021). Animal experiments have also confirmed that LHQW can reduce the viral load in the lungs of mice with viral pneumonia, reduce the expression of inflammatory factors in the lungs, and improve lung injury (Xia et al., 2020; Su et al., 2022). It is well known that viral clearance in the lung depends mainly on lymphocytes and macrophages (Bedi et al., 2022; Cammann et al., 2022; Harris and Borg, 2022; McGee et al., 2022; Verma et al., 2022; Wei et al., 2022; Zhang H. et al., 2022; Zhang M. et al., 2022); however, the present findings do not clarify how LHQW reduces the viral load in the lung and improves lung injury with viral pneumonia.
According to TCM theory, the lung and large intestine interact in terms of physiological, pathophysiological, and immune functions. Both belong to the mucosal immune system, although the intestine and respiratory tract are two separate organs. The latest research has indicated that changes in intestinal microbial composition and function are correlated with the development of lung diseases (“intestine–lung axis”) (Bulanda and Wypych, 2022; Chen et al., 2022; Melo-Gonzalez et al., 2022; Shi et al., 2022; Wang L. et al., 2022; Wang Z. et al., 2022). Moreover, available data indicate that intestinal microbes are strongly correlated with clinical symptoms and inflammatory indices of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (Wang H. et al., 2021; Gutierrez-Castrellon et al., 2022; Sencio et al., 2022; Xavier-Santos et al., 2022). However, the composition and characteristics of the intestinal microbiota altered by LHQW treatment for influenza A virus infection and the relationship between microbiota changes and the antiviral effects of drugs need to be further explored.
This study aimed to investigate the role of intestinal microbiota in LHQW activity in ameliorating viral infectious pneumonia and its possible mechanisms. In this study, 16S rRNA was used to analyze the characteristics of the intestinal microbiota. Computer-aided design technology was used to integrate disease, drug-related genes, and proteins for comprehensive analysis to explore the relationship between microbes and respiratory diseases mediated by LHQW through the “intestine-lung axis” and the potential molecular mechanism of the major active components and targets of LHQW.
Materials and methods
Animal experiment
Reagents and materials
A total of 90 SPF BALB/c mice aged between 6 and 8 weeks (half males and half females) were used in this study. Mice were purchased from Qinglongshan Animal Farm, Jiangning District, Nanjing. The same pathogen-free room was used to accommodate all the mice at 18–25°C and 50–60% humidity. All animal experimental procedures strictly followed the protocol approved by the Ethics Committee of Wannan Medical College (YJS-2020-10-006).
Anti-TLR4 rabbit pAb (A11226, ABclonal Technology Co., Ltd., China), anti-NF-κB p65 rabbit mAb (D14E12, Cell Signaling Technology, United States), anti-phospho-NF-κB p65 (Ser536) anti-rabbit mAb (93H1, Cell Signaling Technology, USA), anti-MyD88 rabbit mAb (D80F5, Cell Signaling Technology, USA), anti-occludin rabbit anti-mouse (ab216327, Abcam, United Kingdom), and anti-ZO1 rabbit anti-mouse (ab216880, Abcam, UK) were used as primary antibodies. The secondary antibody was horseradish enzyme-labelled anti-rabbit IgG [Zsbio ZB-2301]. TRNzol Universal Total RNA Extraction Reagent (DP424, Beijing), Scientific Revertaid First-Strand cDNA Synthesis Kit (10,071,325, Thermo Fisher, USA), and One-Step Prime Script RT–PCR Kit (Perfect Real Time; RR064A, Takara, Japan) were used. Primers were synthesized by Shanghai Shenggong Bioengineering Technology Co., Ltd. Lianhua Qingwen capsules (Z20040063, Yiling, China) and oseltamivir phosphate capsules (H20065415, China Kewei) were used. A nucleic acid protein tester (Biophotometer D30, Germany), high-speed freezing centrifuge (Beckman 64R, USA), and real-time fluorescence quantitative PCR instrument (ABI7500, USA) were used. An enzyme-linked immunosorbent assay (ELISA) test kit for Lipopolysaccharide (LPS) was purchased from Hangzhou Lianke Biotechnology Co., Ltd.
Grouping and administration
The A/Puerto Rico/8/34 (H1N1; PR8) strain was subcultured in chicken embryos. The titer of the amplified virus was 1:440, and the median lethal dose (LD50) was 10–2.1/50 μl. After the mice were lightly anesthetized with ether, 50 μl of 10–2.1 LD50 influenza virus solution was evenly dropped into the nostrils of each mouse to establish the model. The control group(control) received nasal drip of normal saline. Mice were randomly assigned to a control group (18) or an infected group (Buckley and Turner, 2018; Zhang et al., 2021). The mice were fed normally after the nasal drip.
The incidence and death of mice were recorded daily. After nasal drip, the mice were randomly assigned to the virus-infected group (infected), 11 mg d−1 (LHQW-M) Lianhua Qingwen capsule, 22 mg d−1 (LHQW-H) Lianhua Qingwen capsule high-dose group and 395.90 μg d−1 oseltamivir group (OSTW) (Figure 1A). There were 9 female and 9 male mice in each group. The drug group was administered 0.2 ml by gavage at the same time for 7 days, whereas the control and infected groups were given the same amount of distilled water by gavage. The gavage volume of mice is 0.1–0.2 ml/10 g, and what we used is 0.2 ml/20 g. Firstly, opened the Lianhua Qingwen capsule, accurately weighed the powder in the capsule, added it to distilled water, and prepared solutions with different final concentrations. Then, stored the prepared solution in a refrigerator at 4°C, and recovered to room temperature half an hour before each gavage. Food consumption was calculated by weighing the added feed with an electronic balance on the same day and the remaining feed on the next day. The total food intake of the cage mice was the amount of feed added minus the remaining amount (Total food intake of each group of mice).
Figure 1. LHQW improves lung injury and inflammation. (A) Flow diagram; (B) Mouse survival curve: The survival rate of mice in other groups except the infected group (72%) were 100%(n = 90); (C) Weight changes of female mice infected with influenza virus (G) (n = 45); (D) Weight changes of male mice infected with influenza virus (g; n = 45); (E) Daily changes of total food consumption of female mice in each group (g/d; n = 45); (F) Daily changes of total food consumption of male mice in each group (g/d; n = 45); (G) Viral load (copies/mL) in lung tissue of female mice (n = 45); (H) Viral load (copies/mL) in lung tissue of male mice (n = 45); (I) Lung index of female mice (n = 45); (J) Lung index of male mice (n = 45); (K) Lung inhibition rate of female mice (n = 45); (L) Lung inhibition rate in male mice (n = 45); (M) HE staining of mouse lung tissue (200 and 400 times). The mice were randomly assigned to the virus infected group (infected), 11 mg·d−1 (LHQW-M) Lianhua Qingwen capsule, 22 mg·d−1 (LHQW-H) Lianhua Qingwen capsule high-dose group and 395.90 μg·d−1 oseltamivir group (OSTW). Compared with control group, *p < 0.05, **p < 0.01; Compared with infected group, #p < 0.05, ##p < 0.01.
Sample collection and preparation
Mice were anesthetized with ether. Blood was collected from the orbit of each mouse. The mice were sacrificed by cervical dislocation. The abdominal cavity was opened using surgical scissors, and the intestinal canal was dissected. Furthermore, 6 cm jejunum and ileum segments were taken. RNA was extracted from the first 1/2 segment to detect the expression of the inflammatory factor-related mRNAs. The first two segments were fixed in 10% neutral-buffered formalin. The lung tissues of mice were aseptically extracted and evenly divided into three parts. One-third of the total RNA was extracted to detect the expression of inflammatory factor-related mRNA. The other third was homogenized, and the last third was fixed in 10% neutral buffered formalin.
The lung tissue of the mice was aseptically removed, and the connective tissue was trimmed. The blood on the surface was absorbed with filter paper, and the changes in lung organ index were weighed and recorded using an analytical balance. The formulae used in this study were A = B/C*100%, where A represents the lung organ index, B represents the lung weight (mg), and C represents the body weight (g); and D = (E − F)/G*100%, where D represents the inhibition rate of the lung index (%), E represents the average lung index of the infected group, F represents the average lung index of the administration group, and G represents the average lung index of the infected group (Han et al., 2021; Song et al., 2021; Li-Juan et al., 2022; Wei et al., 2022).
Quantitative real-time reverse transcription-PCR was used to detect the expression of related genes
The homogenate supernatants of the lung and small intestine were collected, and the amounts of virus, Interleukin-6 (IL-6), and tumor necrosis factor-α (TNF-α) in the mouse lung were detected using quantitative real-time reverse transcription PCR (qRT-PCR), as was the expression level of Interleukin-1β (IL-1β) and interleukin-10 (IL-10). After homogenization at 4°C and 12,000 RPM, centrifugation was performed for 15 min, after which the supernatant was collected, isopropanol was added at an equal volume to the supernatant, and the mixture was allowed to stand before centrifugation. The supernatant was discarded, and the precipitate remained.
In addition, 1 ml of 75% ethanol was added to each tube and mixed evenly. The supernatant was then centrifuged, discarded, and dried for 5 min. A total of 60 μl of DEPC-treated water was pipetted repeatedly until it was evenly mixed. The sample was dissolved at room temperature, shaken to properly mix, and centrifuged at a low speed for 10 s. After the system was prepared, 20 μl was centrifuged at a low speed for 10 s, vortexed, properly mixed, and then centrifuged again. The reaction parameters were 95°C for 1 min (preheating), 95°C for 15 s, and 60°C for 30 s (40 cycles), followed by the production of a dissolution curve (95°C for 15 s, 60°C for 1 min, and 95°C for 15 s) using β-actin as an internal reference gene. After the reaction, 2−ΔΔCt and the relative gene expression were calculated using the CT method.
Detection of LPS in serum by ELISA
After 30 min of blood sample agglutination and centrifugation, the serum sample was drawn and stored at −80°C. A concentrated standard sample was used for gradient dilution to prepare a standard curve of the serum sample. All the reagents and samples were maintained at room temperature before testing. After soaking the enzyme standard plate, add the standard, add the standard diluent to the blank well, and add 1 × to the sample well detection buffer 90 and 10 μl of sample. Following the kit instructions, double wavelength detection was performed using a microplate reader to determine the optical density (OD) value at the maximum absorption wavelength of 450 nm and a reference wavelength of 570 nm.
Pathological examination
Paraffin sections of the lung and intestine were stained with H&E and alcian blue. After HE staining, the structure of small intestine wall was observed under light microscope. 10 slices of small intestine were taken from each group of mice. Each slice had 10 visual fields, and digital photography was taken. In each photo, the deepest recess depth (subject to the junction of intestinal gland villi to the base of intestinal gland), the thickest mucosal thickness and muscular layer thickness were measured. The morphology and distribution of goblet cells in the epithelium were observed under light microscope, and the number of goblet cell positive cells distributed in each recess was counted. The expression levels of zonula occludens 1 (ZO-1) and occludin were detected by immunohistochemistry. Pathological sections were randomly photographed at magnifications of 100 × and 200 ×. The OD values of the immunohistochemical images were analyzed using ImageJ software to obtain relative expression.
For immunohistochemistry, the fixed small intestine was removed, embedded in paraffin, cut into 5-μm-thick slices, and H2O2 was added and incubated at room temperature for 10 min. After rinsing, 5% bovine serum albumin (BSA) was added and incubated at room temperature for 60 min. Rabbit anti-mouse primary ZO-1 and occludin antibodies were added dropwise without washing before incubation at 4°C overnight. After rinsing, biotinylated sheep anti-rabbit IgG and rabbit anti-goat IgG antibodies were added. The sample was then washed and 3, 3′ diaminobenzidine tetrahydrochloride (DAB) was added at room temperature for 60 min.
For Alcian blue staining, the slices were placed in a drying oven and baked at 66°C for 20–30 min. Three courses of xylene and three courses of ethanol were added successively. Alcian blue staining solution (100 μl) was added dropwise to each slice, which was then dyed in a wet box for 1 h. The staining solution was removed, 100 μl of nuclear solid red staining solution was added to each slice, and sections were washed with water. These samples were passed through three successive passes of ethanol, phenol-xylene, I, and xylene II.
Analysis and characterization of the intestinal microbiota
16S rDNA high-throughput sequencing was performed on the mouse feces collected from each group. The total DNA of bacteria in feces was extracted, and 10 ng of DNA template was used for PCR amplification according to the sequence of the v3–v4 region. The library was constructed using a library building kit, quantified using Qubit and qPCR, and sequenced on a computer. Representative sequences were selected and annotated, and the species were classified using a database.
Western blotting
Fresh small intestinal tissue was collected, and proteins were extracted and quantified. SDS–PAGE, membrane transfer, blocking, and incubation of both primary and secondary antibodies, as well as chemiluminescence and development were carried out. The grayscale of the target protein and internal reference protein was scanned using ImageJ software, and then semiquantitative analysis of the protein content was conducted to obtain NF-κBand p-NF-κB (Phospho NF-κB) levels, as well as the relative expression of Toll-likereceptor4(TLR4) and myeloid differentiation factor 88 (MYD88).
Network pharmacology
Chinese drug compounds were searched in the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) and Traditional Chinese Medicines Integrated Database (TCMID). The PubChem database was used to obtain the structures of the described components, which were imported into the Swiss target prediction database. A target with a prediction score greater than 0 was considered the drug target, and the Batman database was used to obtain the components and targets of gypsum. Oral bioavailability (OB) and drug-likeness (DL) were set at ≥ 30% and ≥ 0.18, respectively, in the TCMSP database to screen for effective components in the LHQW drug group. Finally, the components of Rhodiola in the TCMID database and gypsum in the Batman database were searched (Table 1).
Table 1. Statistics of basic information of traditional Chinese Medicine components - com–onents - targets in drug group.
A total of 253 potentially active components were obtained, and 1,077 drug targets were screened. The GeneCards database was searched using “influenza virus” as the keyword. The targets were selected with an evaluation score greater than 10, and the disease targets were obtained after the removal of weightings.
Cytoscape 3.7.2 software was used to construct the “drug component target disease” network diagram, and a network analyzer was used to analyze the topology of the network diagram. In addition, a PPI cluster analysis diagram was drawn using Cytoscape software. After running the common target in the R language, gene ontology (GO) analysis identified the molecular function, cell composition, biological process, and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway of the top 20 hits.
Molecular docking
The Schrodinger software was used to construct a ligand molecular database for molecular docking. The crystal structure was downloaded from RCSB PDB. The protein structure was imported into Maestro 11.9 platform, and the protein structure was prepared using Schrodinger. The receptor was pretreated, optimized, and minimized (the OPLS3e force field was applied for constraint minimization).
Molecular dynamics simulation
Desmond version 2020 was used for the MD simulation of proteins and compounds. OPLS3e was selected as the molecular force field for the MD simulation, and the TIP3 water model was used to solvate the system. Energy minimization of the entire system was achieved using the OPLS3e force field (all-atomic force field). A Berendsen coupling algorithm was used to create coupling between the temperature and pressure parameters. In the later preparation of the system, 100 ns were run at a time step of 1.2 femtoseconds, and the track was recorded every 10 ps, recording a total of 1,000 frames. The root-mean-square deviation (RMSD) of backbone atoms was calculated, and graphical analysis was performed to illustrate the nature of the interactions between proteins and ligands. The root-mean-square fluctuation (RMSF) of each residue was calculated to determine the major conformational changes in each residue between the initial and dynamic states.
Calculation of MM-GBSA
The basic principle of the molecular mechanics/Poisson Boltzmann (generalized boron) surface area (MM-GBSA) method is to calculate the difference between the binding free energies of two solvated molecules in the binding and free states or to compare the free energies of different solvated conformations of the same molecule.
where includes the bond, bond angle, and dihedral angle energies, is the non-bond van der Waals energy, is the nonbond electrostatic energy, and is the entropy contribution, which can be obtained by normal mode analysis.
Statistical processing
All statistical analyses were performed using the SPSS software (version 18.0). Measurement data are expressed as ± s. Multiple groups were compared using one-way ANOVA. The LSD-t test was used to compare data between the groups.
Results
Lianhua Qingwen capsule reduces lung injury
In this study, the degree of lung injury was evaluated using the mouse lung index, lung inhibition rate, and lung histology (HE staining). Compared with the control group mice, H1N1 infected mice (infected group) showed reduced diet (Figures 1E,F) and weight loss (Figures 1C,D). The infected group had a significantly higher lung index (p < 0.01, Figures 1I,J), widened alveolar septum in the lung tissue, increased inflammatory cell infiltration, and partially collapsed alveolar cavities (Figure 1M). Compared with the infected group, LHQW reduced the viral load in the lungs (p < 0.01, Figures 1G,H) and alleviated the reduction in body weight (p < 0.05, p < 0.01, Figures 1C,D). The LHQW group had a significantly decreased lung index, reduced inflammatory cell infiltration in the pulmonary alveoli and mesenchyme, and a significantly reduced degree of pulmonary interstitial swelling (p < 0.01, Figures 1I–M). The survival rate of mice was significantly higher in the LHQW group (100%) than that in the infected group (72%) (Figures 1B, p < 0.01). Body weight change, total daily food consumption of mice, viral load in mouse lung tissue, lung index, and lung inhibition rate were separately classified and counted for female and male mice. According to the results of statistical calculation, the sex difference is not significant (p > 0.05), so no sex classification statistics was carried out in the subsequent experiments.
Lianhua Qingwen capsule improves the intestinal microbe community
There were differences between the three groups, indicating no errors in the experimental group (Figure 2A). The stacked bar chart of the intestinal microbiota classification at the phylum and genus levels showed that the intestinal microbiota composition was significantly different (Figures 2B,D). Compared with those in the control group, the four indices of alpha diversity of the intestinal microbiota in the infected group decreased to varying degrees, of which the Simpson index decreased significantly (p < 0.05). At the phylum level, the microbiota abundance (Bacteroidetes) decreased significantly, and at the genus level, the microbiota abundance (Muribaculaceae_unclassified, Muribaculum, Odoribacter, Lachnospiraceae_ UCG-006, Prevotellaceae_UCG-001, Anaeroplasma, Absiella, Eubacterium_coprostanoligenes_Group, and Streptococcus) decreased significantly. Compared with the infected group intervention, LHQW reversed the red function in the four indices to varying degrees, and the Shannon index increased significantly (p < 0.05) (Figure 2F). The microbiota abundance (Bacteroidetes, Muribaculaceae_unclassified, Muribaculum, Lachnospiraceae_ UCG-006, and Prevotellaceae_Ucg-001) was significantly higher in the LHQW group and was close to that of the control group (p < 0.05).
Figure 2. LHQW significantly improves the imbalance of intestinal flora induced by H1N1. (A) Beta diversity observed using principal component analysis (PCA) and principal coordinates analysis (PCoA). The same group is circled according to the 95% confidencce interval. The closer the distance between the two samples, the more similar the microbial composition structure between the samples and the smaller the difference. The sp©re analyzed at the phylum and genus levels, and those with significant differences were used to draw a box diagram. Using the genus level as an example, the abscissa represents the species with significant difference with a value of p less than 0.05 in the difference analysis, arranged from left to right according to the abundance from high to low. The ordinate represents the relative abundance. The figure shows the relative abundance of each species in each group. The taxonomic distribution of intestinal bacterial composition at the phylum level, including column stacking diagram (B) and box line diagram (C); The taxonomic distribution of intestinal bacterial composition at the genus level, including column stacking diagram (D) an©ox line diagram (E). (F) Alpha diversity indexes, including Chao1, observed species, Shannon, and Simpson indexes. Chao1 and ACE indexes are mainly used to evaluate flora richness. The larger the value, the higher the flora richnes’. Shannon and Simpson’s indexes are mainly used to evaluate flora diversity. Higher Shannon index and lower Simpson index indicate higher species diversity of the sample.
Lianhua Qingwen capsule reduces intestinal inflammation and improves intestinal mucosal barrier function
HE staining of the small intestine showed that the crypt depth decreased, the muscle layer thinned, and the colon length was shortened in the infected group compared to that in the control group (Figure 3A). Alcian blue staining showed that the number of small intestinal goblet cells and the amount of mucin secreted decreased in the infected group (Figure 3B). Immunohistochemical staining showed reduced expression of ZO-1 (Figure 3C) and occludin (Figure 3D) in the infected group (tan indicates positive expression; Figure 3C). In the LHQW group, the small intestinal crypt depth increased (Figure 3B), the muscle layer thickened (Figure 3C), and the colon length increased compared with that in the infected group (Figure 3A). Alcian blue staining showed that the number of small intestinal goblet cells and the amount of mucin secreted increased in the LHQW group (Figure 3D). LHQW increased ZO-1 and occludin protein expression in the intestinal epithelium. Furthermore, qRT-PCR results showed that IL-6, IL-1β, and TNF-α levels decreased, whereas IL-10 levels increased in the LHQW group compared with those in the infected group (Figure 3E).
Figure 3. LHQW reduces intestinal wall inflammation and improves intestinal mucosal barrier function. (A) Small intestinal HE (100 times, 200 times) and (B) small intestinal recess depth (C) small intestinal muscle layer width according to HE; (A) Colon length; (B,D) the number of goblet cells; Immunohistochemical staining of intestinal mucosal barrier©otein (C) ZO-1 and (E) semi-quantitative analysis of Image J, (D) immunohistochemical staining of Occludin and (F) semi-quantitative analysis of ImageJ. (G) The content of LPS in serum of mic©as detected by ELISA. (E) the expression of inflammatory cytokines in the small intestinal wall was determined by qRT-PCR (IL-6, IL-1β, IL-10, and TNF-α). Compared with control group, *p < 0.05, **p < 0.01; Compared with infected group, #p < 0.05, ##p < 0.01.
Network pharmacology screening results
Network pharmacological screening detected 1,077 drug targets and 421 disease targets. A total of 115 common drug-disease targets were obtained at the intersection of the targets (Figure 4A). The PPI network results were topologically analyzed (Figure 4B), and genes with a degree value greater than the average score were selected as the core targets (Figure 4C).
Figure 4. Network pharmacology screening results. (A) Drug disease comm“n target Wayne diagram and “d”ug component target disease” network diagram. Dark green, light green, pink, and rose red represent drugs, 204 active components in LHQW, 115 common targets, and diseases, respectively. (B) PPI network diagram of protein interaction. (C) Ranking of core targets based on PPI topology analysis (degree ranks the top 30). GO enrichment analysis of LHQW-treated influenza virus: (©Biological process (BP); (E) Cellular component (CC); (F) Molecular function (MF). (G) The top 20 signaling pathways in the KEGG enrichment analysis.
GO results showed that the intersection gene set was enriched in 2333 biological process pathways (Figure 4D), 89 cell components (Figure 4E), and 144 processes related to molecular function (Figure 4F). A total of 172 KEGG pathways were identified using R language. The top 20 terms in the KEGG enrichment bar graph (Figure 4G) were Posey’s sarcoma-associated herpesvirus infection, hepatitis B, lipids and atherosclerosis, AGE-RAGE signalling pathways in diabetes complications, human cytomegalovirus infection, Chagas disease, influenza A virus, toxoplasmosis, measles, C-type lectin receptor signalling pathway, pancreatic cancer, fluid shear stress and atherosclerosis, PD-L1 expression and PD-1 checkpoint pathway in cancer, Toll-like receptor signalling pathway, coronavirus (COVID-19), human T-cell leukemia virus 1 infection, Jerson Prand infection, HIF-1 signalling pathway, TNF signalling pathway, and chemokine signalling pathway. KEGG enrichment analysis showed that infectious diseases mainly caused hair coloring. These results indicate that the Toll-like receptor signalling pathway may be the core mechanism of LHQW-mediated amelioration of influenza A virus infection.
Efficacy analysis of the combination mode and method
AKT1 (PDB ID: 4gv1) target was used to assess the effectiveness of the docking method. The protein had high crystal structure accuracy, no deletion of key residues, clear active sites, and a clear binding mode between small-molecule ligands and active sites (Figure 5C). Protoligand molecules can form strong π-π conjugate bonds with key residues (Ala-230, Met-281, Glu-234, Val-164, and Glu-228). For instance, small ligand molecules form strong hydrophobic bonds with Val-164, thereby stabilizing small molecules. Small molecules can also form strong hydrogen bonds with amino acids (Ala-230, Met-281, Glu-234, and Glu-228) with hydrogen bond distances of 2.0, 2.6, 1.7, and 1.9 Å, respectively. They form short hydrogen bonds with strong binding ability, thus stabilizing small molecules in the active protein sites. In this study, known protoligand pairs were linked to the AKT1 binding site to determine a suitable docking method for screening potential active compounds. The binding conformation overlapped well with that of the ligands in the previous complex (Figure 5D), indicating that the screening method was suitable and effective.
Figure 5. The docking model and analysis of AKT1 with natural ligand. (A) The overall 3D structure of AKT1 complex. The protein backbone was shown as a tube and colored bright blue and pink. The ligand was shown in stick and colored element. (B) A close view of the active site binding with natural ligand. Key residues interacting with ligand were shown as sticks and colored green. (C) The 2D protein-ligand interaction diagram of ligand-AKT1 complex. Protein residues were shown as circles and colored based on their properties: green, hydrophobic residue; purple, polar residue. (D) The re-docking results of active compounds with target.
Virtual screening results
The top 50 compounds were selected based on the energy score obtained from the docking results. 10 compounds were obtained for each gene target by scoring the binding energy and evaluating the key residues of the active sites (Table 2).
Interaction analysis of small protein molecules
The core active components (204) of LHQW were docked with AKT1, STAT3, TLR4, TNF, VEGFA, and IL6 targets. The PyMOL21 software was used to visualize the complexes formed by docking. For each target, the compound was selected with the best score for protein targeting and docking (ligand063, ligand113, ligand103, ligand073, and ligand121). The binding mode between the compounds and proteins was then used to visualize the amino acid residues between the compounds and protein pockets.
Ligand063 has multiple hydrogen bond donors and receptors and thus forms strong hydrogen and hydrophobic bonds with protein active sites (Figure 6A). The hydrogen bond distance was significantly shorter (average of 1.86 Å) than that of the traditional hydrogen bond (3.5 Å). Moreover, the ligand has a strong binding ability and thus anchors small molecules in the protein pockets. Ligand121 contains three benzene rings with strong hydrophobicity and forms strong hydrogen and hydrophobic bonds with the hydrophobic cavity of the protein (Figure 6B). The hydrogen bond distances are 1.6 Å and 1.8 Å. Ligand 113 had a six-membered sugar ring at one end and an oxygen-containing six-membered ring with multiple carboxyl groups at the other end, yielding multiple hydrogen bonds for the receptor (Figure 6C). Ligand103 and ligand073 highly matched the target protein pocket and formed hydrogen, hydrophobic, and other bonds with the protein pocket (Figures 6D-F). These results indicate that ligand063, ligand121, ligand113, and ligand073 are potential active compounds with good binding affinity and matching with AKT1, IL6, STAT3, TLR4, TNF, and VEGFA.
Figure 6. The binding site of target with compound (A-F). The 3D structure of AKT1-ligand063, IL6-ligand121, STAT3-ligand113, TLR4-ligand103, TNF-ligand103, VEGFA-ligand073. (a-f) A close view of the active site binding with AKT1-ligand063, IL6-ligand121, STAT3-ligand113, TLR4-ligand103, TNF-ligand103, VEGFA-ligand073. Key residues interacting with ligands were shown as sticks and colored green. The protein backbone was shown as a tube and colored bright blue.
Molecular dynamics simulation analysis
The stability of reactive proteins and small molecules depends on the RMSD. The larger the RMSD, the more unstable is the protein. The stability of small molecules fluctuated at the beginning and tended to stabilize during movement, reflecting the continuous collision between small molecules and active sites in the protein pocket (Figure 7). This result also showed that the small molecules combined well with the protein pocket and reached a dynamic equilibrium. The average RMSD values of the VEGFA-ligand073, TLR4-ligand103, AKT1-ligand063, and IL6-ligand121 complexes were approximately 2, 5, 10, and 40 ns, respectively. The average RMSD value of TNF-ligand103 was < 2.4 Å and had a slight RMSD fluctuation; thus, the small molecule–protein complex could quickly reach a stable state, reflecting the good stability of the complex. The above results also indicated that the small molecule had a high degree of matching with the protein. The average RMSD of the STAT3-ligand113 complex was < 4.0 Å and showed a slight RMSD fluctuation, possibly due to the unstable chain of the protein. Protein and small molecule complexes quickly reached dynamic equilibrium, with a consistent change trend in the RMSDs of the complexes, indicating that small molecules and proteins can form stable complexes, which change with the conformational state of the proteins in the solvent. Ligand063, ligand121, ligand103, ligand073, and ligand063 formed strong hydrogen and hydrophobic bonds with the active site (Figure 7). In addition, ligand121 exhibited a strong electrostatic interaction with the active site. Ligand 113 formed a strong hydrogen bond with the active site, which induced its binding to the protein. These interactions can improve the stability of small molecules in the protein pockets.
Figure 7. Molecular dynamics simulations. A, B, C, D, E, F represent the RMSD plot during molecular dynamics simulations of AKT1-ligand063, IL6-ligand121, STAT3-ligand113, TLR4-ligand103, TNF-ligand103, VEGFA-ligand073; a, b, c, d, e, f represent the interaction residues of AKT1-ligand063, IL6-ligand121, STAT3-ligand113, TLR4-ligand103, TNF-ligand103, VEGFA-ligand073.
Results of the calculation of free binding energy
The calculation of static molecular docking and molecular mechanics generalized Born surface area (MM-GBSA) free binding energy ensures that the complexes of compounds and targets have sufficient energy for biochemical reactions by providing binding posture and binding free energy. The binding free energy calculated using MM-GBSA supported the molecular docking results. Ligand103 had the highest binding ability to TNF (210.26 ± 22.01) and TLR4 (180.53 ± 19.85; Table 3).
Correlation analysis
Network pharmacological screening and molecular docking results showed that the TNF and Toll receptor signalling pathways targeting TNF and TLR4 might be the main pathways involved in LHQW-mediated amelioration of influenza A virus infection. The 16S rDNA results were combined with the results of molecular biology experiments (qRT–PCR and WB results of lung tissue) to determine the relationship between the amount of virus in the mouse lung and the main target of LHQW in the treatment of influenza A virus and abundance of the intestinal microbiota. The correlation between the abundance of the intestinal microbiota at the phylum and genus levels and the targets screened using computer-aided design technology was also assessed.
At the phylum level, the influenza A virus (Inf A) and TLR4, NF-κB, TNF-α, IL-6, and IL-1β levels were negatively correlated with Bacteroidetes abundance. In contrast, these levels were positively correlated with Verrucomicrobia, Gemmatimonadetes, Rokubacteria, and Chloroflexi. At the genus level, INF-α, TLR4, NF-κB, TNF-α, IL-6, and IL-1β levels were negatively correlated with Prevotellaceae_UCG-001, Muribaculum, and Muribaculaceae_unclassified abundances. In contrast, these levels were positively correlated with Blautia, Klebsiella, Parabacteroides, Roseburia, Bilophila, Eisenbergiella, Citrobacter, Akkermansia, and Clostridiales_unclassified abundance (Figure 8).
Figure 8. Correlation analysis. (A) Toll-like receptor signaling pathway; (B) TLR4, NF-κB, and p-NF-κB levels in mouse lung tissue. MyD88 protein expression level; Correlation network diagram (C) and heat diagram (D) at the phylum level; Correlation network diagram (E) and heat map (F) at the genus level. Compared with control group, *p < 0.05, **p < 0.01; Compared with infected group, #p < 0.05, ##p < 0.01.
Discussion
Oseltamivir (OSTW), a neuraminidase inhibitor commonly used to treat influenza, is an antiviral drug listed in China. OSTW is a classic antiviral drug with good antiviral effects. However, based on the principle of neuraminidase inhibition, the protective effect of this drug on influenza patients has a time window. Taking it 48 h after onset can significantly shorten the course of the disease and inhibit virus replication. Therefore, we selected OSTW as the positive control drug in this study. BALB/c inbred line mice had the same genetic background. The use of Balb/c mice can reduce the influence of genetic background factors among individual mice on the difference in virus infection sensitivity; therefore, it is easier to study and analyze the dynamic characteristics of viral pathogenicity and replication. BALB/c mice were selected for this study according to the standard operating procedures of the Institute for Viral Disease Control and Prevention (2007).
In this study, we found that the inflammatory factors in the lungs of mice with viral pneumonia were disordered, the lungs were damaged, and the TLR4 pathway was significantly activated (Figure 8B). One of the important factors of the influenza virus leading to human death is the abundant expression of pro-inflammatory cytokines. Some studies have shown that lung damage caused by influenza A virus (IVA) is mainly caused by an excessive inflammatory response rather than direct damage to the alveolar epithelium by the virus (Brandes et al., 2013). Excessive immune response can cause a “cytokine storm” that induce acute lung injury (ALI), and further develop acute respiratory distress syndrome (ARDS), leading to respiratory failure and even death. Therefore, in the treatment of influenza virus infection, it is very important to reduce the inflammatory response caused by IVA in the lung (de Jong et al., 2006; Perrone et al., 2008; Bermejo-Martin et al., 2009; Brandes et al., 2013; Li et al., 2014).
During the course of influenza virus infection, three different types of innate immune pattern recognition receptors (PRRs) recognize viral RNA, namely Toll-like receptors (TLRs), retinotide-induced gene I (RIG-I), and NLRP3. After receptor activation, interferon regulators 3/7 (IRF 3/7) and NF-κB are activated to promote the expression of IFNs and pro-inflammatory factors, respectively (Sun et al., 2010). TLRs are the main PRRs for host virus recognition, and during IVA infection, TLR 3/7 is involved in pathogen-associated molecular pattern (PAMPs) identification from IVA by activating myeloid differentiation factor 88 (MYD88), activating downstream transcription factors NF-κB or IRF7, and promoting the expression of pro-inflammatory cytokines such as IL-lβ, IL-18, and type I IFN, respectively (McKinstry et al., 2009).
TNF-α is considered the major pro-inflammatory cytokine capable of causing “cytokine storm” (Cheung et al., 2002), thereby exacerbating the pathogenicity of influenza viruses in humans. (Peper and Van Campen, 1995) found that inhibition of TNF-α reduced lung inflammation but did not affect the clearance of the influenza virus. IL-6 plays a crucial role in resistance to H1N1 influenza infection by promoting neutrophil survival (Dienz et al., 2012). A large body of evidence has shown that the presence of IL-6 may be a marker of persistent inflammation and a direct participant in a comprehensive immune response (Dienz et al., 2012; Spencer et al., 2019; Yousif et al., 2021), IL-6 levels in patients with acute exacerbation were significantly increased compared with those in the stable stage and were positively correlated with viral infection (Hutchinson et al., 2010). IL-1β may be a key inflammatory factor that leads to acute lung injury and ARDS (Fukuyama and Kawaoka, 2011; Niu et al., 2019). It is produced by monocyte-macrophages and endothelial cells and is a pro-inflammatory cytokine involved in host defense responses, such as the immune response to inflammation. It can activate vascular endothelial cells and inflammatory cells such as T lymphocytes and macrophages, release inflammatory mediators, and promote effector cells of the immune system to enter the infection site (Yazdi and Ghoreschi, 2016).
The results of this study showed that the alpha diversity of the intestinal microbiota, the abundance of Bacteroidetes, muribaculaceae_unclassified, and Streptococcus decreased significantly (Figure 2), accompanied by disordered intestinal wall inflammatory factors and injured intestinal mucosal barriers (Figure 3). LHQW treatment increased the survival rate, restored the alpha diversity of intestinal microbiota, recovered the Bacteroidetes and muribaculaceae_unclassified levels, adjusted inflammatory factors of the intestinal walls, and improved the intestinal mucosal barrier (Figure 2).
Intestinal microbiota-derived propionic acid protects against zinc oxide nanoparticle-induced lung injury (Zhang Y. et al., 2022). Baicalin inhibits APEC-induced lung injury by regulating gut microbiota and short-chain fatty acids (SCFA) production (Peng et al., 2021). Mycoplasma gallisepticum baicalin ameliorated-induced inflammatory injury in the chicken lung by regulating the intestinal microbiota and phenylalanine metabolism (Wang J. et al., 2021). The lung and gut microbiota are altered by hyperoxia and contribute to oxygen-induced lung injury in mice (Ashley et al., 2020). Gallacher et al. (2020) found that targeted treatment of the predominant organisms, including those not routinely treated, such as spp., may decrease the development of CLD in preterm-born infants. In summary, these findings suggest that restoration of normal intestinal microbiota can repair lung injury. The lung and intestine have a similar mucosal immune system, which is functionally interconnected. Influenza virus can not only cause lung inflammation but can also be accompanied by gastrointestinal symptoms such as diarrhea and vomiting after infection with influenza H7N9 virus (Gao et al., 2013). In this experiment, we observed numerous inflammatory cell infiltrations in the lungs of mice with pneumonia and significantly increased expression of TNF-α, IL-1β, and IL-6 in the lungs, accompanied by intestinal wall barrier damage, as well as significant changes in the intestinal microbiota (Figures 1C, 2, 3, 8). These phenomena suggest that the pulmonary inflammatory response caused by the virus may lead to cytokines or immune cells entering the intestinal wall through the blood circulation, causing an immune response and resulting in damage to the intestinal wall barrier. Wang et al. (2014) reported that after mice were infected with H1N1 intranasally, the CCR9+ CD4+ T cells in the lungs migrated to the small intestine mediated by CCL25/CCR9 and secreted IFN-γ, changing the composition of the intestinal microbiota and leading to an increase in IL-15 secretion by intestinal epithelial cells, thereby promoting in situ differentiation of CD4+ T cells into Th 17 cells, which eventually triggered intestinal inflammation. Other studies found that in BALB/ C mice inoculated intranasally with influenza H1N1 virus, a significant inflammatory response was observed in the large and small intestine, accompanied by intestinal micromicrobiota disorder (Zhang et al., 2015; Li et al., 2018). Another study reported that after influenza virus infection in an animal model, the abundance of intestinal microbiota, such as Bifidobacterium lactobacillus, was significantly reduced, Enterococcus and Enterobacteria increased significantly, and the bacteria associated with the production of short-chain fatty acids (SCFAs) decreased significantly, which plays an important role in intestinal defense against pathogenic bacteria (Wang et al., 2014; Adak and Khan, 2019).
Bacteroidetes, Muribaculaceae, and Streptococcus are bacteria that produce SCFCs. The intestine has several lines of defense that protect it against the translocation of microorganisms or microbial products into the bloodstream. This multilayer barrier constitutes the largest interface between the external environment and host (Zhang et al., 2021). The gut-vascular barrier, an additional cellular barrier, is situated below the epithelial barrier and is involved in controlling the translocation of microorganisms into the portal vein (Buckley and Turner, 2018). In addition, gut permeability can be altered by gut microbiota dysbiosis, favoring LPS translocation into systemic circulation, with the development of low-grade endotoxemia. LPS is a component of the membrane of gram-negative bacteria present in the gut that can translocate into the systemic circulation, causing non-septic, low-grade endotoxemia (Violi et al., 2022). Gut dysbiosis is a major determinant of low-grade endotoxemia via dysfunction of the intestinal barrier scaffold, which is a prerequisite for LPS translocation into systemic circulation (Johansson et al., 2013; Arbizu et al., 2020). LPS is widely used as an important indicator of intestinal permeability (Li et al., 2021; Liu et al., 2021; Wang Y. et al., 2021; Zhao et al., 2021). Therefore, the detection of LPS in serum by ELISA can be used to evaluate intestinal permeability. In this experiment, the diversity of the intestinal microbiota and the abundance of Bacteroides, Muribaculaceae, and Streptococcus were significantly reduced, resulting in an increase in the abundance of opportunistic bacteria (Figure 2), which then led to intestinal mucosal barrier damage. In this case, enterobacteria metabolites, such as LPS, enter the damaged intestinal wall, aggravate the intestinal wall immune response, and release inflammatory factors, and then amplify the intestinal wall injury, which crosses the injured intestinal wall into the bloodstream and activates the TLR4 signalling pathway in the lungs to aggravate the inflammatory response and injury in the lung, which was also verified in our experiments.
The results of this study showed that LHQW treatment reduced the viral loads and adjusted inflammatory factors in the lungs, alleviated lung injury, and inhibited the activation of the TLR4/NF-κB signalling pathway in viral pneumonia mice (Figures 1C, 8B). Network pharmacological analysis showed that the six active herbal medicine molecules from LHQW could regulate the intestinal microbiota, thereby inhibiting the immune-inflammatory response through the TLR4/NF-κB signalling pathways in the lungs (Table 2). It reveals the relationship between the gut microbiota, intestinal mucosal barrier, and viral pneumonia, and part of the mechanisms of LHQW in the treatment of viral pneumonia. The LHQW capsule, a traditional Chinese medicine compound, has 13 traditional Chinese medicines and many action targets, such as Akt1, P53, Cox-2, and TLR4 (Xia et al., 2020; Su et al., 2022), and has unique advantages in the treatment of influenza virus infection. Formula composition: forsythia suspensa 225 g, honeysuckle 255 g, isatis root 255 g, bitter almond 85 g, menthol 7.5 g, houttuynia cordata 255 g, rhubarb 51 g, patchouli 85 g, Male Fern Rhizome 255 g, Rhodiolae Herba 85 g, Ephedra equisetifolia 85 g, licorice 85 g, gypsum 255 g. Studies have confirmed that the main active ingredients of LHQW are chlorogenic acid, ursolic acid, rutin, luteolin, and quercetin, which exhibit antibacterial activities (Niu et al., 2021; Yang et al., 2021; Zhou et al., 2021).
In this experiment, after treatment with LHQW, the alpha diversity and abundance of probiotics such as Bacteroidetes, Muribaculaceae, and Streptococcus in the intestinal tract of mice with viral pneumonia were significantly increased, and intestinal mucosal barrier injury was significantly alleviated (Figures 2, 3).
As a result, bacterial metabolites (LPS) entering the intestinal wall are reduced and the resulting immune response and inflammatory factors are also decreased. Meanwhile, LPS crossing the intestinal wall to enter the blood stream and eventually reach the lungs was significantly reduced, thus reducing the activation of the LPS/TLR4/NF-κB signalling pathway (Figures 8A,B), thereby decreasing the expression of inflammatory factors in the lungs. In addition, LPS/TLR4 can also be inhibited (Figure 8A), thus reducing lung injury (Figure 9), which is beneficial for the survival of some immune cells and the clearance of the virus in the lungs.
Figure 9. (A) qRT-PCR detection of inflammatory factors (IL-6, IL-1β, IL-10, TNF-α) in the mouse lung tissue. Compared with control group, *p < 0.05, **p < 0.01; Compared with infected group, #p < 0.05, ##p < 0.01. (B) The working model of LHQW capsule against influenza A virus infection. (a) The lung infected with influenza A virus induces an immune response in lung, leading to the release of inflammatory cytokines into blood. (b) The inflammatory cytokines from the blood enter the intestinal wall to induce the immune response, resulting in intestinal wall injury. (c) The LPS from the intestine crosses the damaged intestinal mucosal barrier and enters the blood to activate the TLR4-mediated pathway in lung, causing the lung injury (d).
Conclusion
The present study was designed to determine the effect of intestinal microbiota in the improvement of LHQW activity in viral pneumonia and its possible mechanism. The findings clearly indicate that LHQW is effective for treating influenza A virus infectious pneumonia, and the mechanism is associated with the regulation of the TLR4/NF-κB signalling pathway in the lungs by restoring intestinal microbiota and repairing the intestinal wall. This result provides theoretical support for the effective use of LHQW in the treatment of influenza A virus.
Data availability statement
The data presented in the study are deposited in the NCBI repository, accession number PRJNA889462.
Ethics statement
The animal study was reviewed and approved by all animal experimental procedures strictly followed the protocol approved by the Ethics Committee of Wannan Medical College (YJS-2020-10-006).
Author contributions
PX, ZY, and SD performed the experiments, the data, and wrote the manuscript. ZH and SZ contributed to the study design and overall supervision. All authors have reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was financially supported by the Natural science projects in colleges and universities in Anhui Province (grant number KJ2020ZD56) and the National Natural Science Foundation of China (grant number 81671318). SZ received funding from the Natural science projects in colleges and universities in Anhui Province (grant number KJ2020ZD56). ZH received funding from the National Natural Science Foundation of China (grant number 81671318).
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.
References
Adak, A., and Khan, M. R. (2019). An insight into gut microbiota and its functionalities. Cell. Mol. Life Sci. 76, 473–493. doi: 10.1007/s00018-018-2943-4
Arbizu, S., Chew, B., Mertens-Talcott, S. U., and Noratto, G. (2020). Commercial whey products promote intestinal barrier function with glycomacropeptide enhanced activity in downregulating bacterial endotoxin lipopolysaccharides (LPS)-induced inflammation in vitro. Food Funct. 11, 5842–5852. doi: 10.1039/D0FO00487A
Ashley, S. L., Sjoding, M. W., Popova, A. P., Cui, T. X., Hoostal, M. J., Schmidt, T. M., et al. (2020). Lung and gut microbiota are altered by hyperoxia and contribute to oxygen-induced lung injury in mice. Sci. Transl. Med. 12:eaau9959. doi: 10.1126/scitranslmed.aau9959
Bedi, S., Mudgal, R., Haag, A., and Ono, A. (2022). A Glu-Glu-Tyr sequence in the cytoplasmic tail of the M2 protein renders influenza a virus susceptible to restriction of the Hemagglutinin-M2 association in primary human macrophages. J. Virol. 96:e0071622. doi: 10.1128/jvi.00716-22
Bermejo-Martin, J. F., Ortiz de Lejarazu, R., Pumarola, T., Rello, J., Almansa, R., Ramirez, P., et al. (2009). Th1 and Th17 hypercytokinemia as early host response signature in severe pandemic influenza. Crit. Care 13:R201. doi: 10.1186/cc8208
Boni, M. F. (2008). Vaccination and antigenic drift in influenza. Vaccine 26, C8–C14. doi: 10.1016/j.vaccine.2008.04.011
Brandes, M., Klauschen, F., Kuchen, S., and Germain, R. N. (2013). A systems analysis identifies a feedforward inflammatory circuit leading to lethal influenza infection. Cells 154, 197–212. doi: 10.1016/j.cell.2013.06.013
Buckley, A., and Turner, J. R. (2018). Cell biology of tight junction barrier regulation and mucosal disease. Cold Spring Harb. Perspect. Biol. 10:a029314. doi: 10.1101/cshperspect.a029314
Bulanda, E., and Wypych, T. P. (2022). Bypassing the gut-lung Axis via microbial metabolites: implications for chronic respiratory diseases. Front. Microbiol. 13:857418. doi: 10.3389/fmicb.2022.857418
Cammann, C., Israel, N., Frentzel, S., Jeron, A., Topfstedt, E., Schuler, T., et al. (2022). T cell-specific constitutive active SHP2 enhances T cell memory formation and reduces T cell activation. Front. Immunol. 13:958616. doi: 10.3389/fimmu.2022.958616
Caton, A. J., Brownlee, G. G., Yewdell, J. W., and Gerhard, W. (1982). The antigenic structure of the influenza virus a/PR/8/34 hemagglutinin (H1 subtype). Cells 31, 417–427. doi: 10.1016/0092-8674(82)90135-0
Chen, T. H., Hsu, M. T., Lee, M. Y., and Chou, C. K. (2022). Gastrointestinal involvement in SARS-CoV-2 infection. Viruses 14:1188. doi: 10.3390/v14061188
Cheung, C. Y., Poon, L. L., Lau, A. S., Luk, W., Lau, Y. L., Shortridge, K. F., et al. (2002). Induction of proinflammatory cytokines in human macrophages by influenza a (H5N1) viruses: a mechanism for the unusual severity of human disease? Lancet 360, 1831–1837. doi: 10.1016/S0140-6736(02)11772-7
de Jong, M. D., Simmons, C. P., Thanh, T. T., Hien, V. M., Smith, G. J., Chau, T. N., et al. (2006). Fatal outcome of human influenza a (H5N1) is associated with high viral load and hypercytokinemia. Nat. Med. 12, 1203–1207. doi: 10.1038/nm1477
Dienz, O., Rud, J. G., Eaton, S. M., Lanthier, P. A., Burg, E., Drew, A., et al. (2012). Essential role of IL-6 in protection against H1N1 influenza virus by promoting neutrophil survival in the lung. Mucosal Immunol. 5, 258–266. doi: 10.1038/mi.2012.2
Fukuyama, S., and Kawaoka, Y. (2011). The pathogenesis of influenza virus infections: the contributions of virus and host factors. Curr. Opin. Immunol. 23, 481–486. doi: 10.1016/j.coi.2011.07.016
Gallacher, D., Mitchell, E., Alber, D., Wach, R., Klein, N., Marchesi, J. R., et al. (2020). Dissimilarity of the gut-lung axis and dysbiosis of the lower airways in ventilated preterm infants. Eur. Respir. J. 55:1901909. doi: 10.1183/13993003.01909-2019
Gao, H. N., Lu, H. Z., Cao, B., Du, B., Shang, H., Gan, J. H., et al. (2013). Clinical findings in 111 cases of influenza a (H7N9) virus infection. N. Engl. J. Med. 368, 2277–2285. doi: 10.1056/NEJMoa1305584
Gutierrez-Castrellon, P., Gandara-Marti, T., Abreu, Y. A. A. T., Nieto-Rufino, C. D., Lopez-Orduna, E., Jimenez-Escobar, I., et al. (2022). Probiotic improves symptomatic and viral clearance in Covid19 outpatients: a randomized, quadruple-blinded, placebo-controlled trial. Gut Microbes 14:2018899. doi: 10.1080/19490976.2021.2018899
Han, L., Shi, C., Zeng, X., Cen, L., Mei, X., Fan, J., et al. (2021). A novel bifunctional fusion protein, vunakizumab-IL22, for protection against pulmonary immune injury caused by influenza virus. Front. Immunol. 12:727941. doi: 10.3389/fimmu.2021.727941
Harris, J., and Borg, N. A. (2022). The multifaceted roles of NLRP3-modulating proteins in virus infection. Front. Immunol. 13:987453. doi: 10.3389/fimmu.2022.987453
Hu, K., Guan, W. J., Bi, Y., Zhang, W., Li, L., Zhang, B., et al. (2022). Efficacy and safety of Lianhua Qingwen capsules, a repurposed Chinese herb, in patients with coronavirus disease 2019: a multicenter, prospective, randomized controlled trial. Phytomedicine 85:153242. doi: 10.1016/j.phymed.2021.153800
Huang, K., Zhang, P., Zhang, Z., Youn, J. Y., Wang, C., Zhang, H., et al. (2021). Traditional Chinese medicine (TCM) in the treatment of COVID-19 and other viral infections: efficacies and mechanisms. Pharmacol. Ther. 225:107843. doi: 10.1016/j.pharmthera.2021.107843
Hutchinson, A. F., Black, J., Thompson, M. A., Bozinovski, S., Brand, C. A., Smallwood, D. M., et al. (2010). Identifying viral infections in vaccinated chronic obstructive pulmonary disease (COPD) patients using clinical features and inflammatory markers. Influenza Other Respir. Viruses 4, 33–39. doi: 10.1111/j.1750-2659.2009.00113.x
Institute for Viral Disease Control and Prevention (2007). Chinese Center for Disease Control and Prevention. Standard operating procedures for the National Influenza Center: revised edition (国家流感中心标准操作规程:修订版)[S].2007, May: 189-199. Availalable at: http://www.cnic.org.cn/chn/down/showdown.php?downid=46
Johansson, M. E., Sjovall, H., and Hansson, G. C. (2013). The gastrointestinal mucus system in health and disease. Nat. Rev. Gastroenterol. Hepatol. 10, 352–361. doi: 10.1038/nrgastro.2013.35
Li, D., Feng, Y., Tian, M., Ji, J., Hu, X., and Chen, F. (2021). Gut microbiota-derived inosine from dietary barley leaf supplementation attenuates colitis through PPARgamma signaling activation. Microbiome 9:83. doi: 10.1186/s40168-021-01028-7
Li, H., Liu, X., Chen, F., Zuo, K., Wu, C., Yan, Y., et al. (2018). Avian influenza virus subtype H9N2 affects intestinal microbiota, barrier structure injury, and inflammatory intestinal disease in the chicken ileum, viruses. Viruses 10:270.
Li, Y., Xiao, P., Liu, N., and Zhang, Z. (2022). Efficacy and safety of Chinese medicine Lianhua Qingwen for treating COVID-19: an updated meta-analysis. Front. Pharmacol. 13:888820. doi: 10.3389/fphar.2022.888820
Li, J., Zhang, L., Wu, T., Li, Y., Zhou, X., and Ruan, Z. (2021). Indole-3-propionic acid improved the intestinal barrier by enhancing epithelial barrier and mucus barrier. J. Agric. Food Chem. 69, 1487–1495. doi: 10.1021/acs.jafc.0c05205
Li, Q., Zhou, L., Zhou, M., Chen, Z., Li, F., Wu, H., et al. (2014). Epidemiology of human infections with avian influenza a(H7N9) virus in China. N. Engl. J. Med. 370, 520–532. doi: 10.1056/NEJMoa1304617
Liang, S. B., Zhang, Y. Y., Shen, C., Liang, C. H., Lai, B. Y., Dai, N., et al. (2020). Chinese herbal medicine used with or without conventional Western therapy for COVID-19: an evidence review of clinical studies. Front. Pharmacol. 11:583450. doi: 10.3389/fphar.2020.583450
Li-Juan, L., Kang, S., Zhi-Juan, L., Dan, L., Feng, X., Peng, Y., et al. (2022). Klebsiella pneumoniae infection following H9N2 influenza a virus infection contributes to the development of pneumonia in mice. Vet. Microbiol. 264:109303. doi: 10.1016/j.vetmic.2021.109303
Liu, Y., Cavallaro, P. M., Kim, B. M., Liu, T., Wang, H., Kuhn, F., et al. (2021). A role for intestinal alkaline phosphatase in preventing liver fibrosis. Theranostics 11, 14–26. doi: 10.7150/thno.48468
McGee, M. C., Zhang, T., Magazine, N., Islam, R., Carossino, M., and Huang, W. (2022). PD-1 and ICOS counter-regulate tissue resident regulatory T cell development and IL-10 production during flu. Front. Immunol. 13:984476. doi: 10.3389/fimmu.2022.984476
McKinstry, K. K., Strutt, T. M., Buck, A., Curtis, J. D., Dibble, J. P., Huston, G., et al. (2009). IL-10 deficiency unleashes an influenza-specific Th17 response and enhances survival against high-dose challenge. J. Immunol. 182, 7353–7363. doi: 10.4049/jimmunol.0900657
Melo-Gonzalez, F., Sepulveda-Alfaro, J., Schultz, B. M., Suazo, I. D., Boone, D. L., Kalergis, A. M., et al. (2022). Distal consequences of mucosal infections in intestinal and lung inflammation. Front. Immunol. 13:877533. doi: 10.3389/fimmu.2022.877533
Ni, L., Chen, L., Huang, X., Han, C., Xu, J., Zhang, H., et al. (2020). Combating COVID-19 with integrated traditional Chinese and Western medicine in China. Acta Pharm. Sin. B 10, 1149–1162. doi: 10.1016/j.apsb.2020.06.009
Niu, W. H., Wu, F., Cao, W. Y., Wu, Z. G., Chao, Y. C., and Liang, C. (2021). Network pharmacology for the identification of phytochemicals in traditional Chinese medicine for COVID-19 that may regulate interleukin-6. Biosci. Rep. 41:BSR20202583. doi: 10.1042/BSR20202583
Niu, J., Wu, S., Chen, M., Xu, K., Guo, Q., Lu, A., et al. (2019). Hyperactivation of the NLRP3 inflammasome protects mice against influenza a virus infection via IL-1beta mediated neutrophil recruitment. Cytokine 120, 115–124. doi: 10.1016/j.cyto.2019.04.019
Paules, C., and Subbarao, K. (2017). Influenza. Lancet 390, 697–708. doi: 10.1016/S0140-6736(17)30129-0
Peng, L. Y., Shi, H. T., Tan, Y. R., Shen, S. Y., Yi, P. F., Shen, H. Q., et al. (2021). Baicalin inhibits APEC-induced lung injury by regulating gut microbiota and SCFA production. Food Funct. 12, 12621–12633. doi: 10.1039/d1fo02407h
Peper, R. L., and Van Campen, H. (1995). Tumor necrosis factor as a mediator of inflammation in influenza a viral pneumonia. Microb. Pathog. 19, 175–183. doi: 10.1006/mpat.1995.0056
Perrone, L. A., Plowden, J. K., Garcia-Sastre, A., Katz, J. M., and Tumpey, T. M. (2008). H5N1 and 1918 pandemic influenza virus infection results in early and excessive infiltration of macrophages and neutrophils in the lungs of mice. PLoS Pathog. 4:e1000115. doi: 10.1371/journal.ppat.1000115
Sencio, V., Machelart, A., Robil, C., Benech, N., Hoffmann, E., Galbert, C., et al. (2022). Alteration of the gut microbiota following SARS-CoV-2 infection correlates with disease severity in hamsters. Gut Microbes 14:2018900. doi: 10.1080/19490976.2021.2018900
Shen, X., and Yin, F. (2021). The mechanisms and clinical application of traditional Chinese medicine Lianhua-Qingwen capsule. Biomed. Pharmacother. 142:111998. doi: 10.1016/j.biopha.2021.111998
Shi, C., Zhou, L., Li, H., Shi, X., Zhang, Y., Lu, Y., et al. (2022). Intestinal microbiota metabolizing Houttuynia cordata polysaccharides in H1N1 induced pneumonia mice contributed to Th17/Treg rebalance in gut-lung axis. Int. J. Biol. Macromol. 221, 288–302. doi: 10.1016/j.ijbiomac.2022.09.015
Song, L., Huang, Y., Liu, G., Li, X., Xiao, Y., Liu, C., et al. (2021). A novel immunobiotics Bacteroides dorei ameliorates influenza virus infection in mice. Front. Immunol. 12:828887. doi: 10.3389/fimmu.2021.828887
Spencer, S., Kostel Bal, S., Egner, W., Lango Allen, H., Raza, S. I., Ma, C. A., et al. (2019). Loss of the interleukin-6 receptor causes immunodeficiency, atopy, and abnormal inflammatory responses. J. Exp. Med. 216, 1986–1998. doi: 10.1084/jem.20190344
Su, H., Wu, G., Zhan, L., Xu, F., Qian, H., Li, Y., et al. (2022). Exploration of the mechanism of Lianhua Qingwen in treating influenza virus pneumonia and new coronavirus pneumonia with the concept of “different diseases with the same treatment” based on network pharmacology. Evid. Based Complement. Alternat. Med. 2022:5536266. doi: 10.1155/2022/5536266
Sun, K., Torres, L., and Metzger, D. W. (2010). A detrimental effect of interleukin-10 on protective pulmonary humoral immunity during primary influenza a virus infection. J. Virol. 84, 5007–5014. doi: 10.1128/JVI.02408-09
Treanor, J. (2004). Influenza vaccine--outmaneuvering antigenic shift and drift. N. Engl. J. Med. 350, 218–220. doi: 10.1056/NEJMp038238
Verma, V., Dileepan, M., Huang, Q., Phan, T., Hu, W. S., Ly, H., et al. (2022). Influenza a virus activates cellular Tropomyosin receptor kinase a (TrkA) signaling to promote viral replication and lung inflammation. PLoS Pathog. 18:e1010874. doi: 10.1371/journal.ppat.1010874
Violi, F., Cammisotto, V., Bartimoccia, S., Pignatelli, P., Carnevale, R., and Nocella, C. (2022). Gut-derived low-grade endotoxaemia, atherothrombosis and cardiovascular disease. Nat. Rev. Cardiol. 1–14. doi: 10.1038/s41569-022-00737-2 [Epub ahead of print].
Wang, J., Ishfaq, M., and Li, J. (2021). Baicalin ameliorates Mycoplasma gallisepticum-induced inflammatory injury in the chicken lung through regulating the intestinal microbiota and phenylalanine metabolism. Food Funct. 12, 4092–4104. doi: 10.1039/d1fo00055a
Wang, Z., Li, F., Liu, J., Luo, Y., Guo, H., Yang, Q., et al. (2022). Intestinal microbiota - an unmissable bridge to severe acute pancreatitis-associated acute lung injury. Front. Immunol. 13:913178. doi: 10.3389/fimmu.2022.913178
Wang, J., Li, F., Wei, H., Lian, Z. X., Sun, R., and Tian, Z. (2014). Respiratory influenza virus infection induces intestinal immune injury via microbiota-mediated Th17 cell-dependent inflammation. J. Exp. Med. 211, 2397–2410. doi: 10.1084/jem.20140625
Wang, L., Pelgrim, C. E., Peralta Marzal, L. N., Korver, S., van Ark, I., Leusink-Muis, T., et al. (2022). Changes in intestinal homeostasis and immunity in a cigarette smoke- and LPS-induced murine model for COPD: the lung-gut axis. Am. J. Physiol. Lung Cell. Mol. Physiol. 323, L266–L280. doi: 10.1152/ajplung.00486.2021
Wang, Y., Tian, Y., Zhang, N., Li, X., Wang, X., Wang, W., et al. (2021). Pediococcus pentosaceus PP04 improves high-fat diet-induced liver injury by the modulation of gut inflammation and intestinal microbiota in C57BL/6N mice. Food Funct. 12, 6851–6862. doi: 10.1039/D1FO00857A
Wang, H., Wang, H., Sun, Y., Ren, Z., Zhu, W., Li, A., et al. (2021). Potential associations between microbiome and COVID-19. Front. Med. (Lausanne) 8:785496. doi: 10.3389/fmed.2021.785496
Webster, R. G., Laver, W. G., Air, G. M., and Schild, G. C. (1982). Molecular mechanisms of variation in influenza viruses. Nature 296, 115–121. doi: 10.1038/296115a0
Wei, X., Sun, W., Zhu, P., Ou, G., Zhang, S., Li, Y., et al. (2022). Refined polysaccharide from Dendrobium devonianum resists H1N1 influenza viral infection in mice by activating immunity through the TLR4/MyD88/NF-kappaB pathway. Front. Immunol. 13:999945. doi: 10.3389/fimmu.2022.999945
Wu, L., Chen, Y., Ma, Y., Yang, Z., Yang, N., Deng, W., et al. (2020). Clinical practice guideline on treating influenza in adult patients with Chinese patent medicines. Pharmacol. Res. 160:105101. doi: 10.1016/j.phrs.2020.105101
Wu, J., Wang, Q., Yang, L., Li, Z., and Wang, X. (2021). Potency of Lianhua Qingwen granule combined with paramivir sodium chloride injection in treating influenza and level changes of serum inflammatory factors. Am. J. Transl. Res. 13, 6790–6795.
Xavier-Santos, D., Padilha, M., Fabiano, G. A., Vinderola, G., Gomes Cruz, A., Sivieri, K., et al. (2022). Evidences and perspectives of the use of probiotics, prebiotics, synbiotics, and postbiotics as adjuvants for prevention and treatment of COVID-19: a bibliometric analysis and systematic review. Trends Food Sci. Technol. 120, 174–192. doi: 10.1016/j.tifs.2021.12.033
Xia, Q. D., Xun, Y., Lu, J. L., Lu, Y. C., Yang, Y. Y., Zhou, P., et al. (2020). Network pharmacology and molecular docking analyses on Lianhua Qingwen capsule indicate Akt1 is a potential target to treat and prevent COVID-19. Cell Prolif. 53:e12949. doi: 10.1111/cpr.12949
Xing, D., and Liu, Z. (2021). Effectiveness and safety of traditional Chinese medicine in treating COVID-19: clinical evidence from China. Aging Dis. 12, 1850–1856. doi: 10.14336/AD.2021.0906
Yang, R., Yang, H., Wei, J., Li, W., Yue, F., Song, Y., et al. (2021). Mechanisms underlying the effects of Lianhua Qingwen on sepsis-induced acute lung injury: a network pharmacology approach. Front. Pharmacol. 12:717652. doi: 10.3389/fphar.2021.717652
Yazdi, A. S., and Ghoreschi, K. (2016). The Interleukin-1 family. Adv. Exp. Med. Biol. 941, 21–29. doi: 10.1007/978-94-024-0921-5_2
Yousif, A. S., Ronsard, L., Shah, P., Omatsu, T., Sangesland, M., Bracamonte Moreno, T., et al. (2021). The persistence of interleukin-6 is regulated by a blood buffer system derived from dendritic cells. Immunity 54:e5. doi: 10.1016/j.immuni.2020.12.001
Zhang, H., Alford, T., Liu, S., Zhou, D., and Wang, J. (2022). Influenza virus causes lung immunopathology through down-regulating PPARgamma activity in macrophages. Front. Immunol. 13:958801. doi: 10.3389/fimmu.2022.958801
Zhang, Y., Duan, S., Liu, Y., and Wang, Y. (2021). The combined effect of food additive titanium dioxide and lipopolysaccharide on mouse intestinal barrier function after chronic exposure of titanium dioxide-contained feedstuffs. Part. Fibre Toxicol. 18:8. doi: 10.1186/s12989-021-00399-x
Zhang, M., Li, N., He, Y., Shi, T., and Jie, Z. (2022). Pulmonary resident memory T cells in respiratory virus infection and their inspiration on therapeutic strategies. Front. Immunol. 13:943331. doi: 10.3389/fimmu.2022.943331
Zhang, S., Wei, T., Tianv, H., Cheng, J., Xiao, J., Wang, M., et al. (2015). Small intestinal injury in mice infected with respiratory influenza a virus: evidence for virus induced gastroenteritis. Biotechnol. Lett. 37, 1585–1592. doi: 10.1007/s10529-015-1847-8
Zhang, R., Xu, C., and Duan, Z. (2017). Novel antigenic shift in HA sequences of H1N1 viruses detected by big data analysis. Infect. Genet. Evol. 51, 138–142. doi: 10.1016/j.meegid.2017.03.028
Zhang, Q., Yue, S., Wang, W., Chen, Y., Zhao, C., Song, Y., et al. (2021). Potential role of gut microbiota in traditional Chinese medicine against COVID-19. Am. J. Chin. Med. 49, 785–803. doi: 10.1142/S0192415X21500373
Zhang, Y., Zhang, L., Mao, L., Fan, J., Jiang, X., Li, N., et al. (2022). Intestinal microbiota-derived propionic acid protects against zinc oxide nanoparticles-induced lung injury. Am. J. Respir. Cell Mol. Biol. doi: 10.1165/rcmb.2021-0515OC [Epub ahead of print].
Zhao, L., Xie, Q., Etareri Evivie, S., Liu, D., Dong, J., Ping, L., et al. (2021). Bifidobacterium dentium N8 with potential probiotic characteristics prevents LPS-induced intestinal barrier injury by alleviating the inflammatory response and regulating the tight junction in Caco-2 cell monolayers. Food Funct. 12, 7171–7184. doi: 10.1039/D1FO01164B
Zhou, Y., Niu, M., Zhang, D., Liu, Z., Wu, Q., Chen, J., et al. (2021). Screening for anti-inflammation quality markers of Lianhua Qingwen capsule based on network pharmacology, UPLC, and biological activity. Front. Pharmacol. 12:648439. doi: 10.3389/fphar.2021.648439
Keywords: influenza A virus, intestinal microbiota, network pharmacology, TLR4/NF-κB signalling pathway, Lianhua Qingwen capsule
Citation: Xu P, Yang Z, Du S, Hong Z and Zhong S (2022) Intestinal microbiota analysis and network pharmacology reveal the mechanism by which Lianhua Qingwen capsule improves the immune function of mice infected with influenza A virus. Front. Microbiol. 13:1035941. doi: 10.3389/fmicb.2022.1035941
Edited by:
Claudia Maria Trombetta, University of Siena, ItalyReviewed by:
Iti Saraav, Washington University in St. Louis, United StatesAlexander Larcombe, University of Western Australia, Australia
Copyright © 2022 Xu, Yang, Du, Hong and Zhong. 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: Zongyuan Hong, zyhong@wnmc.edu.cn; Shuzhi Zhong zhongshuzhi2006@163.com
†These authors have contributed equally to this work
‡ORCID: Ping Xu https://orcid.org/0000-0002-0502-3848
Shuzhi Zhong https://orcid.org/0000-0003-0314-9727