1. Introduction
As of May 2023, just after the third anniversary of its officially declared outbreak, Coronavirus Disease 2019 (COVID-19) has caused more than 6.9 million deaths worldwide. The causative agent of this disease, the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), has infected over 760 million people globally, and this number is still rising rapidly [
1]. Despite substantial advances in prevention and treatment strategies, this pandemic still poses global health and economic challenges due to the constantly emerging novel variants of the virus and the long-term consequences of the infection, collectively termed post-COVID conditions or long COVID [
2,
3,
4]. COVID-19 was initially viewed as a respiratory disease, but growing evidence supports the critical role of endothelial dysfunction not only in the pulmonary vasculature but in other organs, both in acute cases and during long COVID [
5,
6,
7,
8,
9].
Severe manifestations of COVID-19 are characterized by progressive respiratory failure resulting from diffuse alveolar damage with inflammatory infiltrates and alveolar edema, intra-alveolar fibrin deposition and hemorrhage, endothelialitis, as well as pulmonary and systemic coagulopathy that form obstructive microthrombi in the lung and other organs [
6,
10,
11]. Additional pathological findings in the vasculature of COVID-19 lungs include disruption of intercellular junctions, basal membrane contact loss, neutrophilic capillaritis/endothelialitis, pulmonary thromboembolism, pulmonary infarctions, and venous thrombosis [
12,
13]. The contribution of different factors to COVID-19‒related endotheliopathy is still under debate. A growing body of evidence points to high levels of pro-inflammatory cytokines and chemokines produced in the lung tissue by infected alveolar epithelial cells and alveolar macrophages [
14,
15,
16], platelet activation [
17,
18], as well as to direct contact of endothelial cells (ECs) with SARS-CoV-2 spike protein [
19,
20,
21,
22]. At the same time, recent data indicate that direct viral infection of ECs is less likely to play a major role in these processes [
12,
22,
23,
24,
25].
In the lung, alveolar epithelial cells are surrounded by an extracellular matrix and adjacent pulmonary microvascular ECs, forming the alveolar-capillary endothelial barrier through inter-endothelial junctions [
26,
27]. The permeability properties of this barrier are tightly regulated through interactions between ECs, surrounding tissue, and biologically active molecules in the blood [
26,
28]. Lung microvascular endothelial injury has been linked to the most severe complications of COVID-19, acute respiratory distress syndrome (ARDS), multiorgan failure, and death [
8,
29]. Recent studies have demonstrated that lung endothelial barrier damage and dysfunction - characterized by increased vascular permeability and loss of barrier integrity leading to the leakage of fluid and plasma proteins into the surrounding tissue - play a critical role in the pathophysiology of the disease and contribute to the development of ARDS [
8,
16,
17]. In the current clinical practice, no specific therapeutic strategies aim to restore the endothelial barrier in COVID-19 patients.
Hydrogen sulfide (H
2S) is a gaseous signaling molecule produced in various mammalian cell types, including ECs [
30,
31]. Among a wide array of physiological functions, it plays a fundamental role in vascular homeostasis, modulates inflammatory responses, and reduces vascular leakage [
30,
32,
33,
34,
35]. Several studies have demonstrated that H
2S improves endothelial barrier function in a variety of experimental conditions. H
2S inhalation attenuated pathologically enhanced blood-brain barrier permeability in animal models of cardiac arrest and resuscitation [
36,
37]. H
2S treatment prevented lipopolysaccharide (LPS)-induced hyperpermeability in EC cultures [
38] and protected against LPS inhalation-induced acute lung injury by reducing neutrophil transmigration and inhibiting pro-inflammatory signaling in animal models [
39,
40]. The H
2S donor NaHS attenuated increased endothelial permeability and inflammation in murine lung specimens challenged by particulate matter inhalation [
41]. In contrast, some reports demonstrated that decreased endogenous H
2S production and altered sulfur metabolism reduced vascular permeability [
42,
43]. These currently available studies indicate that H
2S signaling modulates EC barrier function in a context-dependent manner, and its potential disease-specific impact needs to be determined by targeted research [
30]. To our knowledge, no such data exists about H
2S in the context of SARS-CoV-2–mediated lung endothelial barrier disruption.
Recent findings have provided a rationale for considering the therapeutic implications of H
2S donor molecules in COVID-19 therapy [
44,
45]. Our group and others have shown that H
2S significantly attenuates the replication of several respiratory viruses and virus-induced inflammation [
46,
47]. Moreover, it has been speculated that H
2S blocks SARS-CoV-2 entry into host cells by interfering with angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine 2 (TMPRSS2) expression [
45]. In addition, impaired endogenous H
2S availability is linked to cardiovascular [
48] , metabolic [
49], and pulmonary diseases [
50], which are all risk factors for developing severe COVID-19. Potential alterations in the endogenous H
2S plasma level in COVID-19 patients are still under debate, but data suggest that impaired H
2S availability contributes to COVID-19–associated endotheliopathy and a more severe outcome of this disease [
51,
52,
53]. Furthermore, inhalation of the H
2S donor, sodium thiosulfate, elicited protective effects in COVID-19 patients by reducing symptoms and accelerating recovery [
54]. All these results suggest that H
2S may exhibit beneficial effects in the pathomechanism of COVID-19, but the therapeutic potential of slow-releasing H
2S compounds in lung endothelial barrier disruption associated with SARS-CoV-2 infection has never been explored.
In this study, we aimed to assess the effect of the slow-releasing H2S donor GYY4137 on the barrier function of a human lung microvascular EC monolayer in vitro, after challenging the cells with plasma samples from COVID-19 patients or inactivated SARS-CoV-2 virus. To our knowledge, these data are the first of their kind. We also evaluated inflammatory cytokine levels in the patients’ plasma and determined their correlation with disease severity and impact on the endothelial barrier function.
4. Discussion
One of the novel findings of the present study is that treatment with GYY4137, a well-characterized, slow-releasing H
2S donor, ameliorates endothelial barrier disruption caused by plasma samples from COVID-19 patients
in vitro, regardless of disease severity. In our real-time
in vitro assay model, patient plasma altered endothelial barrier permeability in a highly sample-specific manner, causing barrier damage comparable to TNF-α control, a cytokine known to disrupt endothelial junctions [
66,
67,
68]. Surprisingly, plasma-induced barrier disruption did not correlate with disease status based on the patient’s oxygen requirement: Plasma from some patients with mild disease caused as much or even more damage as some critical patients’ plasma, while plasma from some other individuals in each disease severity group did not elicit any increase in endothelial barrier permeability or even triggered a decrease. Additionally, we determined the cytokine/chemokine profile of each plasma sample and found a high correlation between disease severity and the concentration of several biomolecules, most notably IL-8 and IP-10. In fact, their plasma levels mirrored increasing disease severity much closer than any of the three regularly used clinical markers, LDH, CRP, and D-dimer [
71,
72,
73]. On the other hand, none of the cytokines and chemokines we assessed demonstrated correlation between plasma levels and the corresponding
in vitro barrier function assay results. Finally, inactivated SARS-CoV-2 Omicron BA.1 virus particles elicited a very robust endothelial barrier disruption in our assay that was completely reversed by adding GYY4137. Taken together, these data characterize the effects of patient plasma and inactivated virus particles on the lung microvascular endothelial barrier and provide the basis for further efforts to develop novel treatment modalities that specifically target H
2S signaling in COVID-19 patients.
The COVID-19 pandemic has presented an enormous challenge to countries and health systems worldwide, unprecedented in recent history. This challenge is far from being over due to the constantly evolving nature of the causative agent, SARS-CoV-2, and the debilitating long-term effects of the infection [
2,
3]. The two complementary approaches, prevention and treatment, in the global fight against this pandemic have been increasingly successful but require tremendous amounts of money and effort: As of May 2023, more than three years after the first reported cases, the results still leave room for improvement with about 100,000 new cases and 500-1000 deaths recorded daily worldwide [
1]. Vaccine development and distribution, together with other preventive measures, are used as the first line of defense with limited success around the world [
74,
75]. Serving as the second line, the efficacy of therapeutic approaches in reducing morbidity and mortality of the disease has been gradually increasing, but there is still a great need for novel effective and inexpensive drugs to fill in the gaps of currently available treatment modalities, especially outside the highest-income countries [
76,
77]. In light of recent discoveries about its anti-inflammatory, vasculoprotective, and antiviral effects, H
2S has been proposed as a potential defense against COVID-19 [
44,
45]. In this study, we aimed to assess the effect of an H
2S-donor, GYY4137, on the barrier function of lung endothelial cells, after challenging them with plasma samples from COVID-19 patients, including non-survivors or inactivated SARS-CoV-2 virus. In parallel, we also sought to demonstrate any correlation between inflammatory markers in patient plasma and disease severity.
Damaged endothelial barrier and increased microvascular permeability are hallmarks of severe COVID-19 pathology, greatly contributing to disease severity and mortality [
6,
8,
11]. H
2S has been shown to modulate vascular permeability in several reports, and its effects on the endothelial barrier function have been suggested to be potentially disease- and organ-specific [
30,
36,
38]. In this regard, we established a lung microvasculature-specific assay and tested the effects of the H
2S-donor, GYY4137, in the context of COVID-19–associated biological samples (patient’s plasma and inactivated virus particles). We found a significant increase of barrier function (termed CI, NCI, and RNCI as described in the Materials and Methods section) 12 hours after GYY4137 treatment regardless of the existence and nature of pre-treatments. Interestingly, this effect was maintained, diminished or even reversed at later time points depending on the initial treatment. These results are in line with recent reports showing that H
2S inhalation or H
2S-donor treatment reduced pathologically increased vascular permeability in the brains of rats after cardiac arrest [
37] or in the lungs of mice after particulate matter inhalation [
41], respectively. Remarkably, larger initial barrier disruption by COVID-19 patients’ plasma evoked a relatively greater barrier function increase by GYY4137 treatment. Beyond the obvious base-effect (from a lower base, the same nominal increase constitutes a higher percent), we found an unexplained phenomenon that a 25% or larger RNCI decrease caused by plasma treatment was followed by a significantly greater increase evoked by GYY4137 than when the initial barrier damage was smaller. This could be caused by some of the tight or adherens junctions between ECs that are preferentially restored first or faster by GYY4137 and confer different levels of connecting strength between cells [
78,
79], but to fully understand the reason for this bi-phasic effect will require further investigations. Taken together, our data support the notion that an increased level of H
2S in the lung microenvironment, either by increased endogenous production or by pharmaceutical intervention may be beneficial in severe COVID-19.
Another question we addressed was how the cytokine/chemokine profile of patients’ plasma, the endothelial barrier disruption caused by this plasma, and disease severity correlate with each other. As expected, we found higher levels of inflammatory cytokines and chemokines in the plasma of COVID-19 patients with more severe disease, confirming previously reported data [
15,
59,
61]. Moreover, six (TNF-α, MIP-1α, IL-1ra, G-CSF, IL-8, IP-10) of the assessed signaling molecules demonstrated gradually increasing plasma levels consistent with disease severity, much more so than routine clinical markers, LDH, CRP, and D-dimer [
71,
72,
73]. In fact, we found that while these markers (especially LDH and D-dimer) clearly distinguished the critical group from the others, they did not separate the other groups. While there are obvious advantages to recognizing the critical phase in the course of COVID-19 using blood tests designed to detect these molecules, a panel of the six cytokines/chemokines listed above could provide a better resolution to monitor disease progression. Surprisingly, some biomolecules that have been reported to be potential markers for COVID-19 disease severity, including VEGF [
80,
81], MCP-1 [
60], IFN-γ [
82], PDGF [
81] and others, either did not or only partially correlate with disease severity, maybe due to the relatively small number (10-12) of samples in some of the groups. Nevertheless, our data support previous findings that monitoring blood cytokine/chemokine levels, especially for IL-8 [
58,
61] and IP-10 [
60,
83], could be used as additional biomarkers to help identify and manage COVID-19 patients with different disease severity. On the other hand, we found no correlation between plasma cytokine/chemokine levels and endothelial barrier disruption in our
in vitro assay. In light of this finding, it is not surprising that the measured
in vitro barrier disruption does not correlate with disease severity either. This counterintuitive result is, in fact, in line with recent literature [
84] dissecting the factors in COVID-19 patients’ plasma potentially causing endothelial barrier disruption. Kovacs-Kasa et al. verified that the factor(s) in patients’ plasma disrupting microvascular integrity were heat-labile, but no single or set of cytokine(s) could be accounted for enhanced vascular permeability. They also disproved the potential role of ACE2-binding and complement factors C3a and C5a in the phenomenon. Recent studies implicated several molecular mechanisms, including altered expression and function of adhesion and junction proteins (ICAM-1 and 2, VCAM-1, E- and P-selectin, claudins, occludins, VE-cadherin, Connexin-43, and others) [
78,
79,
85], and/or pathologically modified signaling by integrins, TGF-β, complement, the glycocalyx, mitochondria and (most recently) microRNAs [
12,
19,
86,
87], to contribute to lung endothelial barrier disruption, but the causative agents in the plasma initiating these processes are still widely debated. Since the plasma levels of the cytokines/chemokines assessed in this study are several magnitudes lower than necessary to significantly lower CI values in our assay (e.g., for TNF-α, 1-100 pg vs. 1-100 ng respectively), direct effect from these molecules
in vitro could not be expected. However, lung tissue levels of these mediators are estimated to be potentially over 1,000-fold higher than in plasma during severe inflammation [
88,
89], reaching the necessary levels for endothelial activation and barrier disruption
in vivo. It is also worth noting that disease severity groups were solely based on the patients’ oxygen requirements without consideration of any other clinical characteristics. A more complex classification system including several clinical markers and symptoms could yield different results. Therefore, further
in vitro and
in vivo studies will be needed to resolve this debate.
Circulating virus particles and/or viral spike protein in the patient’s blood have also been implicated in inducing increased microvascular permeability. The potential role and significance of the spike protein (or other viral proteins) in endothelial barrier disruption are still highly controversial [
20,
69,
78,
86,
90] and not always assessed when using these assays. For example, endothelial damage has been reported in the lung after using (1) the spike protein that induced degradation of junction proteins [
78] as well as altered integrin and transforming growth factor beta signaling [
19], (2) the nucleoprotein that induced EC activation via Toll-like receptor 2 and mitogen-activated protein kinase signal pathways [
13] and (3) non-infectious pseudovirus expressing the spike protein that compromised mitochondria and impeded endothelial NO synthase activity [
86]. To this end, we tested the effects of inactivated SARS-CoV-2 Omicron BA.1 virus in our
in vitro assay and demonstrated that challenging the HLMVEC monolayer by 5x104 PFU/well infective dose equivalent inactivated virus results in a drop of CI values comparable to the effect of 10 ng/mL TNF-α. We chose to work with a B.1.1.529 variant as this lineage was circulating in the population at the time of experimentation. Most importantly, inactivated virus-induced endothelial barrier disruption was completely reversed by 300 µM GYY4137. We speculate that the barrier-disruptive potential of the virus particle and/or the spike protein could be lineage- and even sublineage-dependent, explaining the inconsistent data about endothelial barrier disruption available from similar studies using proteins derived from Wuhan or WA1/2020 isolates [
69]. For example, it has been well-demonstrated that Omicron linages feature increased ACE2-affinity and immune evasion capabilities due to several mutations, most of which alter the antigenicity of the spike protein and at the same time modify its structure and function as well [
91,
92,
93]. As a result of their unique virological features in comparison to other SARS-CoV-2 strains, Omicron variants exhibit less efficient TMPRSS2 usage, less spike cleavage, lower fusogenicity, and an altered entry mechanism [
94,
95]. Similarly, the use of different endothelial cells could also be the source of inconsistent findings because of the different genetic backgrounds of the original donors [
69]. Further studies addressing virus-host cell interactions with respect to the spectrum of genetic variations of both could be warranted to better assess the clinical relevance of this pathomechanism.
Our study has certain limitations. First, the number of plasma samples per group was somewhat uneven; there were more critical and healthy samples available than samples belonging to the other three disease severity groups. This weakness of the study design may have introduced a bias towards more significant differences between the two larger sample groups than the others, but we do not believe that it fundamentally altered any of our findings. A follow-up study with larger sample numbers could verify our data. Second, we focused on only one aspect of COVID-19–associated endotheliopathies, the alterations of barrier function using ECs only. While a more complex study could put the results in more context, our simplified monoculture-based approach had the advantage of providing clear answers to some of the basic questions: 1. Are there factors in patient plasma capable of altering endothelial barrier function alone? 2. What correlations exist among the cytokine/chemokine profile of plasma samples, the endothelial barrier disruption caused by them and disease severity? 3. Can inactivated virus alone, as a surrogate of using viral proteins, cause endothelial barrier disruption? And most importantly: 4. Does treatment with an H
2S-donor provide beneficial effects against SARS-CoV-2–associated lung microvascular barrier disruption? Third, we only tested one virus variant and primary lung ECs from one donor as proof of concept. As discussed above, these choices have introduced a genetic bias for the virus-host cell interactions, and we believe that our results justify more comprehensive follow-up studies. Finally, we only tested the effects of a one-time treatment with a single H
2S donor molecule, GYY4137. Other H
2S-releasing agents may have more sustained pharmacological effects as recently reviewed by Szabo and Papapetropoulos [
96]. Further studies will be necessary to clarify the effects of repeated treatments using several different H
2S-releasing compounds to verify whether potential clinical trials could be warranted for pharmacological increase/stabilization of the endothelial barrier as a third pillar for the treatment of COVID-19 in addition to immunomodulators and anti-virals.
Author Contributions
Conceptualization, O.E., P. S., J.C.C., A. N. F., and K.M.; methodology, O.E., P. S., E.L., J.C.C., and K.M.; formal analysis, O.E., P.S., and K.M.; investigation, O.E., P.S., G.T., C.L.V., B.J.L., E.L., T.L.J., and K.M.; resources, C.B.L., S.M., J.C.C., A.N.F., and K.M.; data curation, O.E, P.S., G.T., C.L.V., B.J.S., T.L.J., and C.B.L.; writing—original draft preparation, O.E, P.S., and K.M.; writing—review and editing, O.E, P.S., E.L., C.B.L., S.M., J.C.C., A.N.F., and K.M.; funding acquisition, K.M. All authors have read and agreed to the published version of the manuscript.”
Figure 1.
Disease severity positively correlates with plasma cytokine profiles and routine laboratory data. (A-F) Cytokine levels in plasma samples from COVID-19 patients and healthy volunteers were measured using Bio-Plex Pro Human Cytokine 27-plex Assay (Bio-Rad) and plotted grouped by disease severity. Only the statistical differences compared to the healthy control group are highlighted in these panels; (G-I) Clinical laboratory marker measurements (all available longitudinal data) of the same patient cohort were plotted grouped by disease severity. All statistical differences found between patient groups are labeled. All results in this figure are presented as dot plots of individual mean values for each sample. Bars represent group means and standard deviations. The statistical significance was assessed by Kruskal-Wallis non-parametric one-way ANOVA test followed by Dunn's multiple comparisons. *, p < 0.05; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001.
Figure 1.
Disease severity positively correlates with plasma cytokine profiles and routine laboratory data. (A-F) Cytokine levels in plasma samples from COVID-19 patients and healthy volunteers were measured using Bio-Plex Pro Human Cytokine 27-plex Assay (Bio-Rad) and plotted grouped by disease severity. Only the statistical differences compared to the healthy control group are highlighted in these panels; (G-I) Clinical laboratory marker measurements (all available longitudinal data) of the same patient cohort were plotted grouped by disease severity. All statistical differences found between patient groups are labeled. All results in this figure are presented as dot plots of individual mean values for each sample. Bars represent group means and standard deviations. The statistical significance was assessed by Kruskal-Wallis non-parametric one-way ANOVA test followed by Dunn's multiple comparisons. *, p < 0.05; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001.
Figure 2.
Exogenous H2S released by GYY4137 increases endothelial barrier function. HLMVECs were seeded on E-plates and incubated in growth medium for 48 hours to form a confluent monolayer, then starved for 2 hours before treatment. The effects of biomolecules and patient plasma treatment on the barrier function were monitored by electrical impedance measurements (Cell Index) using the xCELLigence Real-Time Cell Analysis system. Higher Cell Index, as well as Normalized and Relative Normalized Cell Index (described in Materials and Methods) values, represent increased barrier function. (A, B) GYY4137 treatment alone raises barrier function in all tested concentrations. Representative ribbon plots (mean of 3-6 replicates of each condition) of different GYY4137 concentrations (A) and histogram (mean ± SD, n=3-9) of normalized data 12 hours after treatment (B) are presented. Statistical differences for each GYY4137 concentration compared to control at 12 hours are labeled. (C) GYY4137 treatment attenuates TNF-α–induced endothelial barrier disruption. Data are shown as mean ± SD of 3-6 measurements 12 hours after treatment. Statistical differences compared to corresponding zero control are labeled. (D) GYY4137 enhances endothelial barrier function at 12 hours after treatment with plasma samples. Ribbon plots show all measurements (as mean of 3-6 replicates of a single sample of each group) for 24 hours after the first treatment. Thirty-minute, 12-hour, and 24-hour time points are marked by blue, dashed vertical lines. The statistical significance was assessed by one-way (B) or two-way (C) ANOVA, followed by Tukey’s multiple comparisons test. ****, p < 0.0001; GYY, GYY4137.
Figure 2.
Exogenous H2S released by GYY4137 increases endothelial barrier function. HLMVECs were seeded on E-plates and incubated in growth medium for 48 hours to form a confluent monolayer, then starved for 2 hours before treatment. The effects of biomolecules and patient plasma treatment on the barrier function were monitored by electrical impedance measurements (Cell Index) using the xCELLigence Real-Time Cell Analysis system. Higher Cell Index, as well as Normalized and Relative Normalized Cell Index (described in Materials and Methods) values, represent increased barrier function. (A, B) GYY4137 treatment alone raises barrier function in all tested concentrations. Representative ribbon plots (mean of 3-6 replicates of each condition) of different GYY4137 concentrations (A) and histogram (mean ± SD, n=3-9) of normalized data 12 hours after treatment (B) are presented. Statistical differences for each GYY4137 concentration compared to control at 12 hours are labeled. (C) GYY4137 treatment attenuates TNF-α–induced endothelial barrier disruption. Data are shown as mean ± SD of 3-6 measurements 12 hours after treatment. Statistical differences compared to corresponding zero control are labeled. (D) GYY4137 enhances endothelial barrier function at 12 hours after treatment with plasma samples. Ribbon plots show all measurements (as mean of 3-6 replicates of a single sample of each group) for 24 hours after the first treatment. Thirty-minute, 12-hour, and 24-hour time points are marked by blue, dashed vertical lines. The statistical significance was assessed by one-way (B) or two-way (C) ANOVA, followed by Tukey’s multiple comparisons test. ****, p < 0.0001; GYY, GYY4137.
Figure 3.
GYY4137 treatment ameliorates endothelial barrier disruption regardless of disease status. (A, B) Barrier function of an HLMVEC monolayer after treatment with COVID-19 or healthy plasma samples and a second treatment using 300 µM GYY4137 was monitored by electrical impedance measurements using the xCELLigence Real-Time Cell Analysis system. Results are expressed as Relative Normalized Cell Index (RNCI, described in Materials and Methods). Data were collected from at least three independent experiments for each sample (n = 3-6). RNCI values of HLMVEC monolayer after 12h incubation with healthy or COVID-19 plasma samples (A) and the effect of GYY4137 on the RNCI of the plasma-treated HLMVEC monolayer at 12 hours (B) were plotted by disease severity. Data are presented as dot plots of mean values from 3-6 repeated measurements of each sample. Bars represent group means and standard deviations. No statistical differences among disease severity groups were found. (C) Cytokine levels in the plasma samples used for barrier function experiments were measured using Bio-Plex Pro Human Cytokine 27-plex Assay (Bio-Rad) and plotted against the corresponding RNCI values at 12 hours in scatterplots. Slope P-values were calculated to determine statistically significant levels of correlation. None were found. (D) Treatment efficacy plot showing RNCI of HLMVEC monolayer after plasma treatment with or without the addition of GYY4137. Dots represent means of 3-6 measurements for each treatment obtained at the 12-hour timepoint. Statistical difference between RNCI values of untreated and GYY-treated samples is shown. (E) The effect of GYY4137 treatment for each plasma sample was plotted grouped by the level of plasma-induced barrier disruption. Bars represent group means + SD of all measurements (3-6 measurements per sample) 12 hours after treatment. Selected statistical differences among groups are labeled. The statistical significance was assessed by Kruskal-Wallis non-parametric one-way ANOVA test followed by Dunn's multiple comparisons (A, B), simple linear regression (C), 2-tailed, non-parametric, Wilcoxon matched-pairs signed rank t-test (D), and one-way ANOVA followed by Tukey’s multiple comparisons test (E). ****, p < 0.0001; ns, p > 0.05; GYY, GYY4137.
Figure 3.
GYY4137 treatment ameliorates endothelial barrier disruption regardless of disease status. (A, B) Barrier function of an HLMVEC monolayer after treatment with COVID-19 or healthy plasma samples and a second treatment using 300 µM GYY4137 was monitored by electrical impedance measurements using the xCELLigence Real-Time Cell Analysis system. Results are expressed as Relative Normalized Cell Index (RNCI, described in Materials and Methods). Data were collected from at least three independent experiments for each sample (n = 3-6). RNCI values of HLMVEC monolayer after 12h incubation with healthy or COVID-19 plasma samples (A) and the effect of GYY4137 on the RNCI of the plasma-treated HLMVEC monolayer at 12 hours (B) were plotted by disease severity. Data are presented as dot plots of mean values from 3-6 repeated measurements of each sample. Bars represent group means and standard deviations. No statistical differences among disease severity groups were found. (C) Cytokine levels in the plasma samples used for barrier function experiments were measured using Bio-Plex Pro Human Cytokine 27-plex Assay (Bio-Rad) and plotted against the corresponding RNCI values at 12 hours in scatterplots. Slope P-values were calculated to determine statistically significant levels of correlation. None were found. (D) Treatment efficacy plot showing RNCI of HLMVEC monolayer after plasma treatment with or without the addition of GYY4137. Dots represent means of 3-6 measurements for each treatment obtained at the 12-hour timepoint. Statistical difference between RNCI values of untreated and GYY-treated samples is shown. (E) The effect of GYY4137 treatment for each plasma sample was plotted grouped by the level of plasma-induced barrier disruption. Bars represent group means + SD of all measurements (3-6 measurements per sample) 12 hours after treatment. Selected statistical differences among groups are labeled. The statistical significance was assessed by Kruskal-Wallis non-parametric one-way ANOVA test followed by Dunn's multiple comparisons (A, B), simple linear regression (C), 2-tailed, non-parametric, Wilcoxon matched-pairs signed rank t-test (D), and one-way ANOVA followed by Tukey’s multiple comparisons test (E). ****, p < 0.0001; ns, p > 0.05; GYY, GYY4137.
Figure 4.
GYY4137 treatment restores endothelial barrier integrity disrupted by inactivated SARS-CoV-2 Omicron BA.1. Barrier function of an HLMVEC monolayer after treatment with inactivated SARS-CoV-2 or 10 ng/mL TNF-α and a second treatment with 300 µM GYY4137 was monitored by electrical impedance measurements using the xCELLigence Real-Time Cell Analysis system and expressed as Normalized Cell Index and Relative Normalized Cell Index (NCI and RNCI, respectively, described in Materials and Methods). Representative plots of three independent experiments are shown. (A) Ribbon plots show all NCI measurements (as mean of 6-12 replicates) for 24 hours after the first treatment. 12-hour time point is marked by dashed vertical line. (B) Box and whiskers plots represent RNCI data 12 hours after treatment (6-12 replicates per treatment). The box extends from the 25th to 75th percentiles of each sample set, the whiskers go down to the smallest value and up to the largest. The line in the middle of the box is plotted at the median. Selected statistical differences among treatments are labeled. The statistical significance was assessed by one-way ANOVA followed by Tukey’s multiple comparisons test using Graph Pad Prism 9. ****, p < 0.0001; ns, p > 0.05; GYY, GYY4137; Omicron, SARS-CoV-2 Omicron B.1.1.529.
Figure 4.
GYY4137 treatment restores endothelial barrier integrity disrupted by inactivated SARS-CoV-2 Omicron BA.1. Barrier function of an HLMVEC monolayer after treatment with inactivated SARS-CoV-2 or 10 ng/mL TNF-α and a second treatment with 300 µM GYY4137 was monitored by electrical impedance measurements using the xCELLigence Real-Time Cell Analysis system and expressed as Normalized Cell Index and Relative Normalized Cell Index (NCI and RNCI, respectively, described in Materials and Methods). Representative plots of three independent experiments are shown. (A) Ribbon plots show all NCI measurements (as mean of 6-12 replicates) for 24 hours after the first treatment. 12-hour time point is marked by dashed vertical line. (B) Box and whiskers plots represent RNCI data 12 hours after treatment (6-12 replicates per treatment). The box extends from the 25th to 75th percentiles of each sample set, the whiskers go down to the smallest value and up to the largest. The line in the middle of the box is plotted at the median. Selected statistical differences among treatments are labeled. The statistical significance was assessed by one-way ANOVA followed by Tukey’s multiple comparisons test using Graph Pad Prism 9. ****, p < 0.0001; ns, p > 0.05; GYY, GYY4137; Omicron, SARS-CoV-2 Omicron B.1.1.529.
Table 1.
Demographic data, clinical characteristics and laboratory findings of COVID-19 patients included in the study. Data sets of clinical measurements at the time of admission (BMI, temperature, oxygen saturation, respiration rate, LDH, CRP, D-dimer) were analyzed for statistical differences between disease severity groups. None were found. The statistical significance was assessed by Kruskal-Wallis non-parametric one-way ANOVA test followed by Dunn's multiple comparisons and one-way ANOVA followed by Tukey’s multiple comparisons test. a highest level of oxygen delivery required; BMI, body mass index; LDH, lactate dehydrogenase; CRP, C-reactive protein.
Table 1.
Demographic data, clinical characteristics and laboratory findings of COVID-19 patients included in the study. Data sets of clinical measurements at the time of admission (BMI, temperature, oxygen saturation, respiration rate, LDH, CRP, D-dimer) were analyzed for statistical differences between disease severity groups. None were found. The statistical significance was assessed by Kruskal-Wallis non-parametric one-way ANOVA test followed by Dunn's multiple comparisons and one-way ANOVA followed by Tukey’s multiple comparisons test. a highest level of oxygen delivery required; BMI, body mass index; LDH, lactate dehydrogenase; CRP, C-reactive protein.
Disease Severity |
Mild (n=7) |
Moderate (n=6) |
Severe (n=4) |
Critical (n=9) |
Room aira (%) |
100.0% |
0.0% |
0.0% |
0.0% |
Nasal cannulaa (%) |
0.0% |
100.0% |
0.0% |
0.0% |
Non-invasive ventilationa (%) |
0.0% |
0.0% |
100.0% |
0.0% |
Invasive ventilationa (%) |
0.0% |
0.0% |
0.0% |
100.0% |
Male (%) |
57.1% |
33.3% |
75.0% |
66.7% |
Age at Admission, median (years) |
55 |
59 |
54 |
55 |
Hispanic ethnicity (%) |
28.6% |
33.3% |
0.0% |
33.3% |
White race (%) |
57.1% |
83.3% |
75.0% |
55.6% |
Black race (%) |
42.9% |
16.7% |
25.0% |
44.4% |
Days admitted, median (days) |
5 |
7 |
21 |
36 |
On dexamethasone (%) |
14.3% |
50.0% |
75.0% |
100.0% |
Taking remdesivir (%) |
57.1% |
83.3% |
100.0% |
100.0% |
No antiviral used (%) |
42.9% |
16.7% |
0.0% |
0.0% |
COVID vaccinated (%) |
0.0% |
0.0% |
0.0% |
0.0% |
Discharged alive (%) |
100.0% |
100.0% |
75.0% |
11.1% |
Transferred to another facility (%) |
0.0% |
0.0% |
25.0% |
33.3% |
Death (%) |
0.0% |
0.0% |
0.0% |
55.6% |
Clinical characteristics at admission (mean ± SD) |
BMI |
34.80 ± 6.41 |
31.34 ± 4.50 |
39.53 ± 15.30 |
37.8 ± 9.24 |
Body weight (kg) |
99.07 ± 19.04 |
85.98 ± 9.45 |
112.71 ± 22.38 |
113.42 ± 32.73 |
Temperature (Degrees, ᵒC) |
37.10 ± 0.43 |
37.46 ± 1.13 |
38.00 ± 0.69 |
37.35 ± 1.00 |
Oxygen saturation (%) |
96.14 ± 2.96 |
95.16 ± 2.22 |
90.75 ± 14.08 |
89.33 ± 4.35 |
Respiration rate (breaths/minute) |
21.85 ± 5.04 |
21.50 ± 3.56 |
27.75 ± 17.63 |
26.77 ± 8.65 |
LDH (U/L) |
750.45 ± 466.14 |
785.25 ± 204.77 |
1167.33 ± 322.76 |
1065.5 ± 389.10 |
CRP (mg/dL) |
8.50 ± 6.14 |
13.67 ± 6.59 |
10.36 ± 5.31 |
14.49 ± 5.77 |
D-dimer (µg/mL) |
1.31 ± 1.48 |
0.57 ± 0.52 |
0.78 ± 0.36 |
1.63 ± 1.77 |