1. Introduction
COVID-19, short for "Coronavirus Disease 2019," is a highly contagious illness caused by the novel coronavirus SARS-CoV-2.
SARS-CoV-2 infection can lead to a range of clinical manifestations, from mild respiratory symptoms to severe respiratory distress and multi-organ failure, these latter frequently associated with an excessive and dysregulated immune response, leading to a state of hyperinflammation. This extreme release of pro-inflammatory cytokines, especially IL-1β, IL-6, IL-8, and TNF-α [
1], in the context of COVID-19, is often associated with a phenomenon called the "cytokine storm”, representing a distinctive sign of the immunopathogenesis of SARS-CoV-2, in addition to being directly related to the severity and mortality of the disease. Even more critical is the fact that some COVID-19 survivors experience the so-called post-acute sequelae of COVID-19 (PASCA), also known as "long COVID" or "long-haul COVID," a range of persistent symptoms and health issues that continue for weeks or months after the acute phase of a COVID-19 infection has resolved. Even individuals who initially experienced mild cases of COVID-19 can develop these long-lasting symptoms [
1]. It is important to note that researchers are still studying the exact mechanisms behind these long-lasting symptoms and how they relate to the initial COVID-19 infection. Some theories suggest that the virus might trigger an inflammatory response that continues even after the virus is no longer present in significant amounts in the body.
In consideration of the absence of any definitive therapeutic treatment for COVID-19, the identification of novel (and potentially selectable as targets) inflammatory pathways that can mitigate the aggressive and uncontrolled inflammatory response in COVID-19 patients, could be crucial to develop effective therapeutic strategies against SARS-CoV-2 infection [
2]. Besides, currently, the exact mechanisms of SARS-CoV-2-induced inflammation are scantly known.
Methylglyoxal (MG) is an extremely reactive metabolite generated as a byproduct of various cellular processes, mainly glucose metabolism [
3]. It is a potent inducer of a process called glycation, occurring when MG reacts with biological macromolecules (proteins, nucleic acids, lipids) in a non-enzymatic manner, leading to the formation of advanced glycation end products (AGEs). AGEs are stable structures that can accumulate in various tissues, a condition known as glycative stress, causing cellular damage and contributing to the development and progression of several diseases [4-6]. MG-H1 (methylglyoxal-derived hydroimidazolone 1), one of the specific products of protein glycation by MG [
3], can activate inflammatory pathways, usually by binding to the receptor RAGE [
4,
7,
8], which is mostly expressed in the lungs, and promote the release of pro-inflammatory molecules [
4,
8].
To counteract pro-inflammatory MG-H1-derived glycative stress, the body has a major defense mechanism, the GSH-dependent glyoxalase 1 (Glo1), that converts MG, the precursor of MG-H1, into the non-toxic D-lactate [
9].
We herein investigated the involvement of MG-H1/RAGE axis in SARS-CoV-2-induced pathogenic inflammation, by exposing human bronchial (BEAS-2B) and alveolar (A549) epithelial cells to the SARS-CoV-2 spike 1 (S1) protein [
1]. We found that S1 triggered inflammation in both cells. Interestingly, we found, specifically in BEAS-2B cells, that do not express ACE2-R, compared to A549 cells, that S1 exerted a pro-inflammatory action through a novel MG-H1/RAGE-based pathway. Moreover, RAGE and Glo1 expression levels in nasopharyngeal swabs from SARS-Cov-2 positive and negative subjects seems to support the results obtained
in vitro. Altogether, our findings reveal a novel mechanism involved in the inflammation triggered by S1, paving the way to the study of MG-H1/RAGE inflammatory axis in SARS-CoV-2 infection as a potential therapeutic target to mitigate COVID-19 associated pathogenic inflammation.
3. Discussion
In the present research, we first observed, in line with the literature [1,14-17], that SARS-CoV-2 S1 spike surface protein is sufficient alone to activate the production of key COVID-19 pro-inflammatory cytokines (IL-1β, IL-6 and IL-8) in alveolar A549 cells, a major COVID-19 cell-type target. Moreover, we observed a comparable pro-inflammatory response also in human bronchial epithelial BEAS-2B cells, which was never investigated before. Altogether, these results are important in confirming the growing body of studies that show that the S1 spike protein, regardless of the whole virus, has a proinflammatory biological activity in lung epithelial cells. Despite epithelial cells, compared to innate immune cells, usually provide feebler inflammatory responses, they release cytokines and/or chemokines able to potently recruit monocytes, lymphocytes and neutrophils in the lungs infected by SARS-CoV-2, thus anyhow making a significant contribution to the “cytokine storm”-related immunopathology underpinning COVID-19 pathogenesis and progression [
14]. Moreover, there is growing evidence supporting that sustained and protracted circulating levels of SARS-CoV-2 S1 spike1 protein are present in PASC [18-21]. Our findings suggest that these plasma levels, being S1 able to induce inflammation itself, may contribute to the sustained inflammatory-associated events, at least at lung level, that characterize post-acute COVID-19 syndrome [
22]. Hence, our data contribute to focus on and strength the proinflammatory role of S1 and its implication in COVID-19 progression, including PACS, whose underlying mechanisms are entirely unknown.
Importantly, in this study we also demonstrated that in human bronchial epithelial BEAS-2B cells, SARS-CoV-2 S1 spike protein induces inflammation through a novel MG-H1/RAGE-based mechanism. MG-H1 is one of the specific advanced glycation end products originated by MG [
3], able to initiate inflammatory pathways, frequently by binding to the receptor RAGE [
4,
7,
8], that is mostly expressed in the lungs [
4,
8]. As previously reported [
23,
24], BEAS-2B cells expresses very low to undetectable levels of ACE2 receptor, to which S1 typically binds in order to allow endocytosis of the viral particle and begin its replicative cycle within the host cell. Hence, our MG-H1/RAGE axis would represent a novel mechanism for S1, and potentially a novel strategy for SARS-Cov-2, to start pathogenic inflammation in cells that do not express the surface ACE2 receptor. In support of this hypothesis, in A549 cells that express ACE2 [
23], MG-H1/RAGE axis is not activated by S1.
Notably, aminoguanidine (AG), by specifically scavenging MG, thus preventing MG-derived MG-H1 formation [
6,
10], reverted RAGE-dependent inflammation. This result, in addition to demonstrate a causative role of MG-H1/RAGE axis in driving inflammation upon S1 challenge, indicates that AG could be promising in developing supportive therapeutics to prevent COVID-19-related pathogenic inflammation. Besides, AG has already proven potentiality in other inflammatory diseases [25-27]. Clearly, further investigation in this ambit is mandatory.
It is now consolidated that RAGE signaling produces pro-inflammatory responses via NF-kB activation [
14,
28] and, indeed, it has been recently described that SARS-CoV-2 spike protein is pro-inflammatory in pulmonary cells using NF-κB pathway activation [14,29-31]. Moreover, in pancreatic islets infected with SARS-CoV-2, elevated MG activates RAGE and NF-kB [
32]. Finally, an upregulation of NF-kB has been observed in COVID-19 patients [
33]. Hence, it is very likely that NF-kB may be involved down-stream RAGE activation also in our model and potentially upon virus infection.
With the aim of unveil how MG-H1 accumulation, leading to MG-H1/RAGE proinflammatory pathway activation, could occur in bronchial BEAS-2B cells, we found that this was ascribed to the down-regulation, at all expression levels, of the enzyme that modulates MG intracellular levels, named Glo1. Again, the mechanism seems to be specific for ACE2-independent BEAS-2B, since ACE2-dependent A549 cells showed increased Glo1 expression (
Figure S3). It has been reported that glycolysis, the main route in the production of MG, is upregulated in SARS-CoV-2-infected cells to provide the substrate required for the virus replication [
34,
35]. Hence, it would be possible that also this metabolic route may contribute to S1-induced MG/MG-H1 accumulation in BEAS-2B cells, even though the fact that in A549 cells where S1 should still stimulate glycolysis, MG-H1 levels does not change upon S1 treatment, would make one conclude that Glo1 is indeed primarily responsible for it.
Moreover, we found that Glo1 down-regulation was dependent on S1-driven Nrf2 desensitization. It is known that the promoter region of Glo1 has a functionally operating antioxidant response element (ARE) [
3] that is bound by Nrf2 to induce Glo1 expression. Here, we found that Nrf2 desensitization causatively led to Glo1 down-regulation, since Nrf2 activation restored Glo1 expression and enzyme activity. In line with previous results [
11,
36,
37], we here confirm Glo1 control by Nrf2. Interestingly, and in agreement with our findings, a very recent
in vivo study has demonstrated that S1 induces Nrf2 down-regulation in rats exposed to this virus subunit [
38]. It is also noteworthy that a very complex interplay exists between Nrf2 and NF-kB pathways [
39]. Nrf2 can decrease the activation of NF-kB, while NF-kB transcription can inhibit Nrf2 activation [
40,
41]. The existence of such a mutual interaction would seem to strength our hypothesis of NF-kB involvement in our pro-inflammatory Nrf2/Glo1/MG-H1/RAGE-(NF-kB) pathway, where Nrf2 desensitization could very likely explain NF-kB activation. Concomitantly, NF-kB activation could inhibit Nrf2 in a circuitry nurturing S1-driven inflammation through Glo1/MG-H1/RAGE axis, which is very intriguing and needs further investigation.
Finally, we found that patients with a moderate/severe disease presented higher levels of RAGE and IL-1β, and lower levels of Glo1 expression, compared with the non-severe group, thus strongly suggesting that the mechanistic
in vitro core results, showing MG-H1/RAGE pathway as a novel SARS-CoV-2 S1-driven pro-inflammatory axis, might occur also during infection. Besides, Alomar
et al. have recently reported significantly higher MG and IL-1β plasma levels in patients who required intensive care unit (ICU) hospitalization compared to uninfected individuals and significantly lower Glo1 levels in patients who died compared to ICU patients that survived [
42]. They also found strong inverse correlations between plasma MG and Glo1 and significant positive correlations between plasma MG and IL-1β. In the same study, authors state, as one limitation of the research, the need of additional methodology to measure MG, Glo1 and IL-1β plasma levels. By detecting MG, Glo1 and IL-1β levels in swabs from the nasopharynx of individuals, by RT-PCR our results appear to support those of Alomar
et al.
4. Materials and Methods
4.1. Reagents
The chemicals used in the present study are analytical grade reagents from various sources. The SARS-CoV-2 S1 spike protein was purchased from Prodotti Gianni S.p.A (Milan, Italy). MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide] and aminoguanidine bicarbonate (AG) were purchased from Merck Spa (Milan, Italy). Nrf2 activator (Nrf2-A) was from EMD Millipore Corporation (Billerica, MA, USA). Nrf2-A was solubilized in dimethyl sulfoxide (DMSO, Merck Spa, Milan, Italy) (final concentration in the essays of DMSO was 0.01%). Controls contained an identical volume of DMSO vehicle. Laemmli buffer was from Invitrogen (Milan, Italy), Roti-Block from Roth (Roth, Germany) and bicinchoninic acid (BCA) kit from Pierce (USA).
4.2. Cell cultures
Human bronchial BEAS-2B and alveolar A549 epithelial cells were from the American Type Culture Collection (ATCC). Both cells were cultured in RPMI medium complemented with antibiotics and antimycotic (Invitrogen, Paisley, UK) at 37°C and 5% CO
2 [
43,
44].
4.3. Cell viability and morphology
Human bronchial BEAS-2B and alveolar A549 cell viability was studied by MTT assay [
6,
45]. Cell morphology was studied by phase-contrast microscopy [
45].
4.4. RNA Isolation, Reverse Transcription, and Real-Time Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR) Analyses
TRIzol Reagent (ThermoFisher Scientific, Milan, Italy) was used to isolate total RNA. RevertAid™ H Minus First Strand cDNA Synthesis Kit (ThermoFisher Scientific, Milan, Italy) was used to produce cDNA, as previously described [
9]. Gene expression of the studied genes versus GAPDH was evaluated by RT-PCR, using MX3000P Real-Time PCR System (Agilent Technology, Milan, Italy). The sequences of the oligonucleotide primers are reported in
Table 1.
A total volume of 20 µL containing 25 ng of cDNA, 1X Brilliant II SYBR® Green QPCR Master Mix, ROX Reference Dye, and 600 nM of specific primers was employed in PCR reactions, thermal cycling conditions being: 1 cycle (95 °C, 5 min); 45 cycles (95 °C, 20 s and 60 °C, 30 s). Comparative analysis of gene expression was performed by the 2−(∆∆CT) method [
6].
4.5. TNF-α, IL-6 and IL-8 detection
TNF-α (cod. KHC3011), IL-6 (cod. BMS213HS) and IL-8 (cod. KHC0081) were measured by ELISA (ThermoFisher Scientific, Milan, Italy).
4.6. Detection of Methylglyoxal (MG)-H1
MG-H1 was measured by ELISA (cat. STA-811, DBA Italia S.r.l., Milan, Italy).
4.7. Cell and nuclear lysis and Western blot
Cell lysis was performed in radioimmunoprecipitation assay (RIPA) lysis buffer [
44]. FractionPREP Cell Fractionation kit (Biovision, Vinci-Biochem, Florence, Italy) was used for nuclear extraction [
44]. Western blot was performed as previously described [
46]. Samples were boiled in Laemmli buffer for 5 min, resolved on SDS-PAGE and blotted onto a nitrocellulose membrane. Roti-Block (room temperature, RT, 1 h) was used to block unspecific binding sites. Membranes incubation with the appropriate primary Abs: mouse anti-Glo1 (D-6) mAb (dilution 1:1000, Santa Cruz, cat. # sc-133144, DBA Italia Srl, Segrate, Italy), mouse anti-RAGE (A-9) mAb (dilution 1:1000, Santa Cruz, cat. # sc-365154, DBA Italia Srl, Segrate, Italy) and mouse anti-GAPDH (6C5) mAb (dilution 1:1000, Santa Cruz, cat. # sc- 32233, DBA Italia Srl, Segrate, Italy), as internal loading control, was performed overnight at 4 °C. Membranes were then incubated with the appropriate HRP-conjugated secondary Ab (RT, 1 h) and ECL was used as revealing system (Amersham Pharmacia, Milan, Italy).
4.8. Nrf2 activation detection
Nrf2 Transcription Factor Assay kit (Colorimetric) (ab207223) (DBA Italia srl, Milan, Italy), was employed to detect Nrf2 activation in nuclear extracts [
44].
4.9. Total protein and enzyme specific activity of Glyoxalase 1 (Glo1) detection
BCA kit (cat. 23225, ThermoFisher Scientific) was used to quantify total protein concentration. Bovine serum albumin was used as a standard. Glo1 enzyme activity was assayed by an established method [
6].
4.10. Patients
Nasopharyngeal swabs (NPS) from individuals referring to the Hospital of Perugia for SARS-CoV-2 molecular testing were collected from December 2022 to May 2023. The study was approved by the Ethical Committee of the Umbria Region (CER Umbria, approval number 3990/21) and of the University of Perugia (15/2023 del 15.03.2023) and carried out in accordance with the Declaration of Helsinki. Informed consent was obtained from the subjects involved in the study. The cohort of SARS-CoV-2 infection cases consisted of 60 patients, stratified into 2 groups, on the basis of clinical symptom severity into non-severe (asymptomatic: n = 20) and moderate/severe (total n = 40; ordinary hospitalization, n = 20 and intensive care, n = 20) groups. Once the specimens were tested for SARS-CoV-2 viral load, the residual sample was used for DNA extraction performed by using STARMag 96 x 4 Universal Cartridge Kit (Seegene, Korea) on automated STARLET (Seegene, Korea), according to the manufacturer's instructions.
4.11. Statistical analysis
Results were expressed as means ± standard deviation (SD) of three independent experiments. One-way analysis of variance with Dunnett’ s correction was employed to determine differences among groups. T-test or χ2 were also applied in the analysis of patients’ data where appropriate. Statistical significance was set at p < 0.05.
Figure 1.
Effect of SARS-CoV-2 S1 spike protein on human bronchial BEAS-2B and alveolar A549 cell viability. Viability of (a,b) BEAS and (c,d) A549 cells unexposed or exposed to S1 was measured by MTT assay. Data report the means of three separate experiments carried out in triplicate and error bars correspond to the standard deviation (SD) of the mean.
Figure 1.
Effect of SARS-CoV-2 S1 spike protein on human bronchial BEAS-2B and alveolar A549 cell viability. Viability of (a,b) BEAS and (c,d) A549 cells unexposed or exposed to S1 was measured by MTT assay. Data report the means of three separate experiments carried out in triplicate and error bars correspond to the standard deviation (SD) of the mean.
Figure 2.
SARS-CoV-2 S1 spike protein induces inflammatory cytokines in human bronchial BEAS-2B and alveolar A549 epithelial cells. BEAS-2B and A549 cells were stimulated with S1 at a concentration of 25 ng/ml. Six and 24 h post-stimulation, the expression of IL-1β was evaluated by real-time RT-PCR, while TNF-α, IL-6 and IL-8 levels, by ELISA. Data represent mean ± SD (n = 3). *p < 0.05, **p < 0.01, ***p < 0.001.
Figure 2.
SARS-CoV-2 S1 spike protein induces inflammatory cytokines in human bronchial BEAS-2B and alveolar A549 epithelial cells. BEAS-2B and A549 cells were stimulated with S1 at a concentration of 25 ng/ml. Six and 24 h post-stimulation, the expression of IL-1β was evaluated by real-time RT-PCR, while TNF-α, IL-6 and IL-8 levels, by ELISA. Data represent mean ± SD (n = 3). *p < 0.05, **p < 0.01, ***p < 0.001.
Figure 3.
SARS-CoV-2 S1 spike protein affects MG-H1 levels and RAGE expression in human bronchial BEAS-2B cells. BEAS-2B and A549 cells were stimulated with S1 at a concentration of 25 ng/ml. 24 h post-stimulation, the levels of MG-H1 (a,c) were evaluated by a specific ELISA kit, while RAGE mRNA and protein expression (b,d), by real-time RT-PCR and western blot (WB), respectively. GAPDH was used as internal loading control for WB normalization. Data represent mean ± SD (n = 3). *p < 0.05, **p < 0.01.
Figure 3.
SARS-CoV-2 S1 spike protein affects MG-H1 levels and RAGE expression in human bronchial BEAS-2B cells. BEAS-2B and A549 cells were stimulated with S1 at a concentration of 25 ng/ml. 24 h post-stimulation, the levels of MG-H1 (a,c) were evaluated by a specific ELISA kit, while RAGE mRNA and protein expression (b,d), by real-time RT-PCR and western blot (WB), respectively. GAPDH was used as internal loading control for WB normalization. Data represent mean ± SD (n = 3). *p < 0.05, **p < 0.01.
Figure 4.
SARS-CoV-2 S1 spike protein induces inflammation in human bronchial BEAS-2B cells through MG-H1/RAGE axis. BEAS-2B cells were pretreated with the specific MG scavenger aminoguanidine (AG) (10 mM for 24 h). 25 ng/ml S1 spike protein was then added and left for further 24 h. The expression of RAGE and IL-1β was evaluated by real-time RT-PCR, while TNF-α, IL-6 and IL-8 levels, by ELISA. Data represent mean ± SD (n = 3). *p < 0.05, **p < 0.01.
Figure 4.
SARS-CoV-2 S1 spike protein induces inflammation in human bronchial BEAS-2B cells through MG-H1/RAGE axis. BEAS-2B cells were pretreated with the specific MG scavenger aminoguanidine (AG) (10 mM for 24 h). 25 ng/ml S1 spike protein was then added and left for further 24 h. The expression of RAGE and IL-1β was evaluated by real-time RT-PCR, while TNF-α, IL-6 and IL-8 levels, by ELISA. Data represent mean ± SD (n = 3). *p < 0.05, **p < 0.01.
Figure 5.
SARS-CoV-2 S1 spike protein controls MG-H1/RAGE proinflammatory pathway through Nrf2-dependent Glyoxalase 1 (Glo1) down-regulation in human bronchial BEAS-2B cells. BEAS-2B cells were exposed to 25 ng/ml S1 for 24 h and (a) Glo1 mRNA, protein expression, specific activity, evaluated by real-time RT-PCR, western blot (WB) and spectrophotometric assay, respectively, and (b) Nrf2 nuclear expression, evaluated by a specific assay, were studied. Pretreatment with 10 µM Nrf2 activator (Nrf2-A) rescued (c) Glo1 expression and activity, reduced (d) MG-H1 intracellular accumulation, measured by a specific ELISA kit, (e) RAGE expression, evaluated by real-time RT-PCR, as well as inflammation, evaluated by IL-1β (f) and IL-6 levels (g), by real-time RT-PCR and ELISA, respectively, compared with S1-challenged cells. GAPDH was used as internal loading control for WB normalization. Data represent mean ± SD (n = 3). *p < 0.05, **p < 0.01, ***p < 0.001.
Figure 5.
SARS-CoV-2 S1 spike protein controls MG-H1/RAGE proinflammatory pathway through Nrf2-dependent Glyoxalase 1 (Glo1) down-regulation in human bronchial BEAS-2B cells. BEAS-2B cells were exposed to 25 ng/ml S1 for 24 h and (a) Glo1 mRNA, protein expression, specific activity, evaluated by real-time RT-PCR, western blot (WB) and spectrophotometric assay, respectively, and (b) Nrf2 nuclear expression, evaluated by a specific assay, were studied. Pretreatment with 10 µM Nrf2 activator (Nrf2-A) rescued (c) Glo1 expression and activity, reduced (d) MG-H1 intracellular accumulation, measured by a specific ELISA kit, (e) RAGE expression, evaluated by real-time RT-PCR, as well as inflammation, evaluated by IL-1β (f) and IL-6 levels (g), by real-time RT-PCR and ELISA, respectively, compared with S1-challenged cells. GAPDH was used as internal loading control for WB normalization. Data represent mean ± SD (n = 3). *p < 0.05, **p < 0.01, ***p < 0.001.
Figure 6.
RAGE, IL-1β and Glo1 expression in nasopharyngeal swabs of SARS-CoV-2 infected patients at different clinical severity of Covid-19. Gene expression of (a) RAGE, (b) IL-1β and (c) Glo1was evaluated by real-time RT-PCR. N: SARS-Cov-2 negative group; NS: SARS-Cov-2 positive with non-severe Covid-19 disease; M/S: SARS-Cov-2 positive with moderate/severe Covid-19 disease. *p < 0.05, **p < 0.01.
Figure 6.
RAGE, IL-1β and Glo1 expression in nasopharyngeal swabs of SARS-CoV-2 infected patients at different clinical severity of Covid-19. Gene expression of (a) RAGE, (b) IL-1β and (c) Glo1was evaluated by real-time RT-PCR. N: SARS-Cov-2 negative group; NS: SARS-Cov-2 positive with non-severe Covid-19 disease; M/S: SARS-Cov-2 positive with moderate/severe Covid-19 disease. *p < 0.05, **p < 0.01.
Table 1.
Sequences of oligonucleotide primers.
Table 1.
Sequences of oligonucleotide primers.
Gene |
Sense primer (5’-3’) |
Antisense primer (5’-3’) |
Glo1 |
AGAAAGCACGGGGTGAAACTG |
TACACCTTCAGTCCCGACTCC |
IL-1β |
GGACCTGGACCTCTGCCCTCTGG |
GCCTGCCTGAAGCCCTTGCTGTAG |
RAGE |
TGAAGGAACAGACCAGGAGACAC |
GCACAGGCTCCCAGACAC |
GAPDH |
CAAGGTCATCCATGACAACTTTG |
GTCCACCACCCTGTTGCTGTAG |