Andrographolide decreased RSV viral load and attenuated inflammation in vitro.
Andrographolide has been found to exert antiviral function over a variety of viruses, but its role on RSV infection has not been reported yet. To clarify whether andrographolide has antiviral effect on RSV, we conducted in vitro experiment using human lung carcinoma epithelial cell line A549. Firstly, we investigated whether andrographolide induced cytotoxicity in A549 cell line and selected appropriate concentrations to conduct in vitro experiment. As shown in Fig. 1A, obvious impact on cell viability was seen in concentration more than 20 µM after 36 h treatment. Therefore, we selected 5µM and 10µM as appropriate concentration to conduct following in vitro experiment. Next, we measured the RSV N gene load both at 24 h and 36 h post-infection. After treatment with andrographolide, RSV viral load was noticeably decreased at 36 h dose-dependently (Fig. 1B). To further elucidate the impact of andrographolide on RSV-induced inflammation in vitro, we measured IL-6 and IL-8 levels in cell supernatants 36 h post-infection. Both cytokine levels were dramatically increased post RSV infection, while andrographolide attenuated them significantly (Fig. 1C and D).
Identification of RSV infection-related genes
To identify RSV regulated cellular genes, RNAseq dataset GSE32139 was downloaded from the Gene Expression Omnibus (GEO) database. As presented in Table 2, there were a total of 1519 DEGs between the mock-infected and the RSV group, and among them, 891 were up-regulated, while 628 were down-regulated. The results are presented as heatmap of all DEGs (Fig. 2A), and volcano plot (Fig. 2B). DEGs represented in red were up-regulated, while those in blue were down-regulated. Detail information about the DEGs of RSV infection was presented as supplementary material (Supplementary Materials Table S1).
Table 2
RSV infection DEGs in GSE32139
Compare
|
All
|
Up
|
Down
|
Threshold
|
RSV vs Mock
|
1519
|
891
|
628
|
DESeq2 pvalue<0.05 |log2FoldChange|>0.5
|
Andrographolide-related target identification via PharmMapper database
PharmMapper Database and PubChem were used as described in the Materials and Methods. The chemical structure of andrographolide (Fig. 3A) and the potential targets of andrographolide along with their gene symbols are presented in Supplementary Table S2.
GO and KEGG pathway enrichment analysis
As illustrated in Fig. 3B, RSV infection-related genes and andrographolide target genes were merged, and the intersection of the two sets was considered as the potential targets of andrographolide for the treatment of RSV infection. There were 25 genes in total, as presented in Table 3. To further understand these targets, GO and KEGG enrichment analysis were performed and shown in Fig. 3C and D. The top five significant functions were regulation of cytokine-mediated signaling pathway, regulation of response to cytokine stimulus, cellular response to interferon-gamma (IFN-γ), regulation of response to IFN-γ, and regulation of IFN-γ-mediated signaling pathway. The KEGG pathways involved were influenza A, pyruvate metabolism, necroptosis, steroid hormone biosynthesis, NOD-like receptor signaling pathway, cytosolic DNA-sensing pathway, chemokine signaling pathway, prolactin signaling pathway, Leishmaniasis, lipid and atherosclerosis, Coronavirus disease COVID-19, PD-L1 expression and PD-1 checkpoint pathway in cancer, Th1 and Th2 cell differentiation, Toll-like receptor signaling pathway, C-type lectin receptor signaling pathway, Th17 cell differentiation, toxoplasmosis, growth hormone synthesis secretion and action, and 2-Oxocarboxylic acid metabolism. Target genes related to these pathways are listed in Table 4.
Table 3
Potential targets of andrographolide for the treatment of RSV infection
Target gene
|
Target protein
|
ADH1C
|
Alcohol dehydrogenase 1C
|
AMY1A
|
Alpha-amylase 1A
|
ANXA5
|
Annexin A5
|
BCAT2
|
Branched-chain-amino-acid aminotransferase, mitochondrial
|
CASP1
|
Caspase-1
|
CCL5
|
C-C motif chemokine 5
|
CFD
|
Complement factor D
|
DUSP6
|
Dual specificity protein phosphatase 6
|
GLO1
|
Lactoylglutathione lyase
|
GMPR
|
GMP reductase 1
|
HNMT
|
Histamine N-methyltransferase
|
JAK2
|
Tyrosine-protein kinase JAK2
|
LTA4H
|
Leukotriene A-4 hydrolase
|
MMP12
|
Macrophage metalloelastase
|
MMP13
|
Collagenase 3
|
NR1H3
|
Oxysterols receptor LXR-alpha
|
SEC14L1
|
SEC14-like protein 1
|
SORD
|
Sorbitol dehydrogenase
|
SPARC
|
SPARC
|
STAT1
|
Signal transducer and activator of transcription 1-alpha/beta
|
STS
|
Steryl-sulfatase
|
SULT1E1
|
Sulfotransferase Family 1E Member 1
|
TNNC1
|
Troponin C, slow skeletal and cardiac muscles
|
TTR
|
Transthyretin
|
TYMS
|
Thymidylate synthase
|
Table 4
Target genes involved in KEGG pathways
ID
|
KEGG Pathway
|
P value
|
Gene Symbol
|
hsa05164
|
Influenza A
|
0.000886596
|
CASP1/CCL5/JAK2/STAT1
|
hsa00620
|
Pyruvate metabolism
|
0.006454215
|
ADH1C/GLO1
|
hsa04217
|
Necroptosis
|
0.007619001
|
CASP1/JAK2/STAT1
|
hsa00140
|
Steroid hormone biosynthesis
|
0.010690573
|
STS/SULT1E1
|
hsa04621
|
NOD-like receptor signaling pathway
|
0.011355888
|
CASP1/CCL5/STAT1
|
hsa04623
|
Cytosolic DNA-sensing pathway
|
0.011373674
|
CASP1/CCL5
|
hsa04062
|
Chemokine signaling pathway
|
0.01274064
|
CCL5/JAK2/STAT1
|
hsa04917
|
Prolactin signaling pathway
|
0.013911961
|
JAK2/STAT1
|
hsa05140
|
Leishmaniasis
|
0.016673315
|
JAK2/STAT1
|
hsa05417
|
Lipid and atherosclerosis
|
0.017247702
|
CASP1/CCL5/JAK2
|
hsa05171
|
Coronavirus disease - COVID-19
|
0.021089305
|
CASP1/CFD/STAT1
|
hsa05235
|
PD-L1 expression and PD-1 checkpoint pathway in cancer
|
0.021902894
|
JAK2/STAT1
|
hsa04658
|
Th1 and Th2 cell differentiation
|
0.023304382
|
JAK2/STAT1
|
hsa04933
|
AGE-RAGE signaling pathway in diabetic complications
|
0.027218185
|
JAK2/STAT1
|
hsa04620
|
Toll-like receptor signaling pathway
|
0.029268752
|
CCL5/STAT1
|
hsa04625
|
C-type lectin receptor signaling pathway
|
0.029268752
|
CASP1/STAT1
|
hsa04659
|
Th17 cell differentiation
|
0.031379942
|
JAK2/STAT1
|
hsa05145
|
Toxoplasmosis
|
0.03355041
|
JAK2/STAT1
|
hsa04935
|
Growth hormone synthesis secretion and action
|
0.037487397
|
JAK2/STAT1
|
hsa01210
|
2-Oxocarboxylic acid metabolism
|
0.048190537
|
BCAT2
|
RT-qPCR verification of potential target genes
To verify change of gene mRNA expression analyzed by network pharmacology, we extracted total RNA at 24 and 36 h post RSV infection on A594 cells. The expression level of target genes (CASP1/CCL5/JAK2/STAT1) and the downstream gene of CASP1 (which also functions as IL-1β-converting enzyme) were measured using RT-qPCR. As shown in Fig. 4, RSV infection dramatically elevated the expression of CASP1, CCL5, JAK2, STAT1 and IL-1β both at 24 and 36 h. Furthermore, andrographolide clearly mitigated the increase of CASP1, CCL5, JAK2, STAT1 and IL-1β mRNA levels at 36 h.
Andrographolide inhibited apoptosis of RSV-infected epithelial cells
There are three major programmed cell death pathways, namely, apoptosis, pyroptosis and necroptosis, which have been documented to be involved in airway epithelial cells infected with viruses[35]. Airway epithelial cell death caused by viral infection has been considered an important defense mechanism that restricts virus replication and spread[18]. Since CASP1 is an important gene involved in pyroptosis signal pathway[41], we conducted flow cytometry analysis using Annexin V-FITC/PI dual staining to figure out whether andrographolide have effect on cell death after RSV infection. A549 cells, intervened or not with andrographolide, were collected after 36h infection with RSV, and was stained with Annexin V-FITC/PI according to the manufacturer’s instruction. We defined Annexin V positive, PI negative quadrant as early apoptotic cells, while dual Annexin V and PI positive quadrant as late apoptotic or necrotic cells[46]. After 36 h of infection, RSV noticeably increased the ratio of apoptosis of A549 cells. Andrographolide obviously inhibited apoptotic ratio of RSV-infected cells, which was dose-dependent responsive. Similar effect was not seen in uninfected cells (Fig. 5).
Andrographolide elevated protein levels of caspase-1, cleaved caspase-1, N terminal of GSDMD and Bcl-2 to suppress apoptosis and promote pyroptosis of RSV infected cells
As the results of RT-qPCR and flow cytometry had indicated, andrographolide may have impact on death pattern of airway epithelial cells infected with RSV. In the process of pyroptosis, formation of inflammasome cleaves pro-caspase-1 into its active form, which then cleaves Gasdermin-D (GSDMD), pro-IL-1β, pro-IL-18, and pro-IL-33 into biologically active forms[15], resulting the production of IL-1β, IL-33, IL-18 and pyroptosis[28, 40]. Therefore, we measured the protein levels of pro-caspase-1, cleaved caspase-1, cleaved IL-1β and N terminal of GSDMD at 36 h post RSV infection to further elucidate the effects of andrographolide on pyroptosis. The results showed that RSV infection decreased the protein levels of all four, but andrographolide elevated them (Fig. 6). Since B-cell lymphoma-2 (Bcl-2) family play a pivotal role in regulating cell apoptosis, we also assessed the expression level of antiapoptotic protein Bcl-2 to identify the effects of andrographolide on apoptosis. RSV suppressed the expression of antiapoptotic protein Bcl-2 at 36 h post infection, while andrographolide increased it (Fig. 6). Based on the above results, we hypothesized that andrographolide elevated protein levels of caspase-1, cleaved caspase-1, N terminal of GSDMD and Bcl-2 to suppress apoptosis and promote pyroptosis of RSV infected cells.