Next Article in Journal
Human SMILE-Derived Stromal Lenticule Scaffold for Regenerative Therapy: Review and Perspectives
Next Article in Special Issue
Could SGLT2 Inhibitors Improve Exercise Intolerance in Chronic Heart Failure?
Previous Article in Journal
Investigating the Role of GABA in Neural Development and Disease Using Mice Lacking GAD67 or VGAT Genes
Previous Article in Special Issue
Blocking of SGLT2 to Eliminate NADPH-Induced Oxidative Stress in Lenses of Animals with Fructose-Induced Diabetes Mellitus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential In Vitro Effects of SGLT2 Inhibitors on Mitochondrial Oxidative Phosphorylation, Glucose Uptake and Cell Metabolism

1
Institute for Biomedicine and Health Sciences (HEALTH), Joanneum Research Forschungsgesellschaft m.b.H, Neue Stiftingtalstrasse 2, 8010 Graz, Austria
2
Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, AstraZeneca, 431 83 Gothenburg, Sweden
3
Division of Endocrinology and Diabetology, Department of Internal Medicine, Medical University of Graz, Auenbruggerplatz 15, 8036 Graz, Austria
4
Cooperative Centre for Regenerative Medicine (COREMED), Joanneum Research Forschungsgesellschaft m.b.H, Neue Stiftingtalstrasse 2, 8010 Graz, Austria
5
Research Unit for Tissue Regeneration, Repair and Reconstruction, Division of Plastic, Aesthetic and Reconstructive Surgery, Department of Surgery, Medical University of Graz, 8036 Graz, Austria
6
Pharmaceutical Sciences, R&D, AstraZeneca, 415 25 Göteborg, Sweden
7
Late-Stage Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, AstraZeneca, Gaithersburg, MD 20878, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(14), 7966; https://doi.org/10.3390/ijms23147966
Submission received: 21 April 2022 / Revised: 7 July 2022 / Accepted: 9 July 2022 / Published: 19 July 2022

Abstract

:
(1) The cardio-reno-metabolic benefits of the SGLT2 inhibitors canagliflozin (cana), dapagliflozin (dapa), ertugliflozin (ertu), and empagliflozin (empa) have been demonstrated, but it remains unclear whether they exert different off-target effects influencing clinical profiles. (2) We aimed to investigate the effects of SGLT2 inhibitors on mitochondrial function, cellular glucose-uptake (GU), and metabolic pathways in human-umbilical-vein endothelial cells (HUVECs). (3) At 100 µM (supra-pharmacological concentration), cana decreased ECAR by 45% and inhibited GU (IC5o: 14 µM). At 100 µM and 10 µM (pharmacological concentration), cana increased the ADP/ATP ratio, whereas dapa and ertu (3, 10 µM, about 10× the pharmacological concentration) showed no effect. Cana (100 µM) decreased the oxygen consumption rate (OCR) by 60%, while dapa decreased it by 7%, and ertu and empa (all 100 µM) had no significant effect. Cana (100 µM) inhibited GLUT1, but did not significantly affect GLUTs’ expression levels. Cana (100 µM) treatment reduced glycolysis, elevated the amino acids supplying the tricarboxylic-acid cycle, and significantly increased purine/pyrimidine-pathway metabolites, in contrast to dapa (3 µM) and ertu (10 µM). (4) The results confirmed cana´s inhibition of mitochondrial activity and GU at supra-pharmacological and pharmacological concentrations, whereas the dapa, ertu, and empa did not show effects even at supra-pharmacological concentrations. At supra-pharmacological concentrations, cana (but not dapa or ertu) affected multiple cellular pathways and inhibited GLUT1.

1. Introduction

The selective sodium-glucose co-transporter-2 inhibitors (SGLT2) canagliflozin (cana), dapagliflozin (dapa), ertugliflozin (ertu), and empagliflozin (empa) are approved as therapeutic glucose-lowering agents for the treatment of type 2 diabetes patients [1,2,3,4,5,6,7,8]. They show a 250-, 1200-, 2500-, and 2000-fold selectivity for SGLT2 over SGLT1 [9], and were developed based on the naturally occurring SGLT inhibitor, phlorizin. Phlorizin was isolated from apple-tree bark in 1835 [10], and since it was shown to normalize glycemia in 90% of pancreatectomized animals in the 1980s, interest in its potential to treat diabetes has grown [11]. However, phlorizin acts as a nonspecific SGLT inhibitor that is poorly absorbed from the gastrointestinal tract and not sufficiently stable for clinical use; therefore, it has not been considered for further clinical development [12].
Synthetic selective SGLT2 inhibitors, such as cana, empa, and others, have been developed and their clinical development has progressed. Cana, as a highly potent and selective SGLT2 inhibitor, has shown to lower the renal threshold for glucose, increase urinary glucose excretion, improve glycemic control and beta-cell function in rodent models of type 2 diabetes mellitus, and to reduce body-weight gain in rodent models of obesity [13]. Cana, the first SGLT2 inhibitor to be FDA-approved, in 2013 [14], and one of the most commonly used SGLT2 inhibitors, is about 5 to 10 times less selective for SGLT2 over SGLT1 than dapa, ertu, and empa [8]. Thus, cana may inhibit SGLT1, at approved therapeutic doses, in the small intestine [15].
Furthermore, sotagliflozin, a non-selective SGLT inhibitor acting on both SGLT1 and 2 [16], with strong affinity for SGLT2 and mild selectivity for SGLT1, has been developed and, recently, the results of two randomized controlled trials focusing on heart failure and renal outcomes of sotaglifozin-treatment were published [17,18].
SGLT1 is responsible for a major part of the dietary glucose uptake in the intestine [19], and the SGLT1-inhibitory effect gives the added value of delaying glucose absorption in the intestine and preventing post-meal hyperglycemia, [5,6,8,20,21,22,23]. However, although sotagliflozin has been approved in the EU since 2019 for overweight type 1 diabetic patients who do not achieve adequate glycemic control despite optimal insulin therapy [24], it is not yet approved by the FDA due to the increased risk of ketoacidosis [25].
SGLT2 is responsible for the majority of glucose reuptake in the tubular system of the kidneys [19], and selective SGLT2 inhibitors act by blocking the SGLT2-mediated reabsorption of glucose and sodium in the proximal tubule of the kidneys [26]. Apart from the glucose-lowering effect, SGLT2 inhibitors reduce body weight, blood pressure and, more importantly, the risk of cardiovascular events, heart-failure hospitalizations, and the progression of renal diseases [5,6,8,20,21,22,23]. Although the use of SGLT2 inhibitors is generally regarded as safe, cana has been associated with an increased risk of lower-limb amputations in type 2 diabetic patients [23,27,28,29,30]. Hence, discussions are ongoing as to whether the various SGLT2 inhibitors have different off-target effects at the cellular level.
To address potential differences in the clinical profiles between the SGLT2 inhibitors, several in vitro studies have investigated the cellular effects of SGLT2 inhibitors on energy metabolism [31,32,33,34], cellular function, proliferation and migration [35], and cell metabolism [33,36,37]. A common finding of these studies in different cell models is that cana exhibits an inhibitory effect on mitochondrial function [22,23,31,32,33]. Effects of dapa and empa on mitochondrial complex I activity have also been reported, although these are weak and cell-line-dependent [31,32,34].
Endothelial cells make up the lining of the vascular system, and their dysfunction may lead to complications that are observed in diabetes and other metabolic disorders [38,39,40,41]. The autonomous effects of SGLT2 inhibitors on endothelial cells have been investigated for their relevance to vascular function of SGLT2 inhibition [34,35,42]. The data from these studies showed that cana inhibited cell proliferation and glucose uptake (GU), while dapa and empa did not [35,42].
The inhibition of mitochondrial function by cana has also been observed in previous metabolomics studies on several cell lines, which showed alterations in the fatty-acid and purine/pyrimidine metabolism [33,36,37]. Dapa treatment has been shown to increase levels of amino-acid metabolites and to reduce levels of various lipid species [37].
However, currently, no data are available for the effects of different SGLT2 inhibitors on mitochondrial activity, cellular glucose uptake, and metabolic processes in endothelial cells; effects that could potentially influence the clinical profiles of these SGLT2 inhibitors.
Therefore, in this study, we aimed to investigate the effects of cana, dapa, ertu, and empa on mitochondrial oxidative phosphorylation and GU in human-umbilical-vein endothelial cells (HUVECs). To assess the putative impact of SGLT2 inhibitors on the function of the endothelium, we also examined the impact of the SGLT2 inhibitors cana, dapa, and ertu on the metabolomes of these cells.

2. Results

2.1. Effects on the Oxygen-Consumption Rate (OCR)

The HUVEC cells treated with a supra-pharmacological concentration (100 µM) of the SGLT2 inhibitor cana showed inhibition of OCR of nearly 60% (Figure 1A), similar to complex 1 inhibitor rotenone (2 µM). By contrast, dapa (100 µM) reduced the OCR by only 7%, while ertu and empa (both at 100 µM) had no significant effect on the OCR in the HUVECs compared to the DMSO (Figure 1A,B).
Next, we assessed the dose responses of the different SGLT2 inhibitors on the OCR. Cana inhibited the OCR in a concentration-dependent manner with a half-maximal inhibitory concentration (IC50) of 27 µM. The IC50 values for dapa, ertu, and empa were greater than 100 µM (data not shown). The OCR inhibition became apparent at concentrations above 10 µM of cana (Figure 1C), and at 100 µM, the effect on OCR was fully pronounced.

2.2. Effects on the Extracellular Acidification Rate (ECAR)

At supra-pharmacological concentrations (100 µM), cana, dapa, ertu, and empa lowered the ECAR by 45%, 13%, 32%, and 17%, respectively (Figure 2A,B). The IC50 values were 9 µM for cana and 39 µM for ertu. The IC50 values for dapa and empa were greater than 100 µM. At a concentration of 10 µM, the effect of cana on ECAR was less pronounced.

2.3. Effects on the ADP/ATP Ratio

The cana treatment of the HUVECs resulted in an increase in the ADP/ATP ratio at the concentration equal to the Cmax elicited by its pharmacological dose (10 µM), and an even larger statistically significant response at the supra-pharmacological concentration (100 μM). The Dapa treatment, both at the pharmacologically relevant concentration of 0.3 µM and at the supra-pharmacological concentration of 3 μM, failed to cause any significant changes to the ADP/ATP ratio. However, the treatment with ertu at 1 µM (pharmacologically relevant concentration) and 10 μM displayed a weak, concentration-dependent increase in the ADP/ATP ratio (Figure 3). The rotenone increased the ADP/ATP ratio to 25.

2.4. Effects on the GU Using 14C-Deoxyglucose (14C-DOG)

Both cana and ertu inhibited 14C-DOG uptake with IC50 values of 14 µM and 85 µM, respectively (Figure 4A). At 100 µM, dapa and empa showed weak, but significant, increases in 14C-DOG uptake (dapa: p = 0.0411; empa p = 0.0003), with IC50 values > 100 µM (Figure 4A). The GLUT1 inhibitor, BAY-876 (25 µM), inhibited 60% of the 14C-DOG uptake and cana (100 µM) inhibited 97% of the 14C-DOG uptake (Figure 4B).

2.5. Effects of Cana on the GLUT Gene Expression

The treatment of the HUVECs with cana (10 and 100 μM) had no significant effects on the expression of the glucose transporters GLUT 1, 3, 6, 8, and 10, nor detectable cytotoxic effects (data not shown). The rotenone did not cause detectable cytotoxicity.

2.6. Metabolomics

The four main classes of ontology distribution for the MVA_UVA metabolites measured were amino acids, peptides and analogues (16 metabolites), nucleotides and analogues (21 metabolites), carbohydrates and conjugates (ten metabolites), and metabolites related to energy metabolism (seven metabolites). One cana (10 µM) sample was excluded due to a technical problem during the HPLC-MS measurement. Seventeen metabolites were altered significantly (p-value < 0.01) in the cana (100 µM)-treated HUVECs (Figure 5), while five metabolites for the ertu (10 µM)-treated and two metabolites for the dapa (3 µM)-treated cells showed at least low significance levels (p-value < 0.05). The same metabolites and metabolite classes were also changed in the HUVECs treated with a lower concentration of cana (10 µM), although fewer metabolites were affected (three metabolites, p-values ≤ 0.05). All the metabolite specifications are listed in detail in the Supplementary Material.

2.6.1. Effects of Cana on the Levels of Amino Acid Supplying the Tricarboxylic Acid Cycle (TCA) Cycle

Cana at 100 µM significantly increased the levels of most of the amino acids supplying the TCA cycle (Figure 6), including glutamine (p-value: 0.036) and glutamate (p-value: 0.004) (Figure 7). Cana at 10 µM showed similar trends in changes in amino acid levels and increased aspartate levels. Ertu at 10 µM and dapa at 3 µM had no effects on the levels of these amino acids.

2.6.2. Effects of Cana on Glycolysis Activity and Beta-Oxidation

The metabolites used for glycolysis were significantly depleted in the cana (100 µM)-treated HUVECs (p-value < 0.01; Figure 7).
The cana (100 µM)-treated HUVECs showed a trend towards reduced levels of free fatty acids (FFA). Cana at 10 µM showed the same trend, but less pronounced. The acetyl-CoA levels were significantly increased (p-value = 0.005) in the cana (100 µM)-treated HUVECs (Figure 5), whereas the treatment with cana at 10 µM induced a trend, albeit with no significance, towards increased acetyl-CoA levels. The supra-pharmacological concentrations of dapa and ertu (3 µM and 10 µM) had no discernible effect on the metabolites in the glycolysis.

2.6.3. Effects on Redox Equivalents

At 100 µM, the cana significantly increased the NAD+ and decreased the NADH levels. The NAPDH levels were shifted in the direction of the NADP+. At 10 µM, cana depicted a trend towards increased NAD+ and decreased NADH levels, and the trend of the NADPH shifted to the NADP+ levels. The HUVECs treated with cana (100 µM) had a significantly increased ADP/ATP ratio and cana (10 µM) induced the same effect as the trend, but no significance was observed. Furthermore, no significant effect on redox equivalents was induced by the SGLT2 inhibitors ertu and dapa.

2.6.4. Effects on Metabolites Used in Apoptotic Pathways

The levels of the apoptosis-associated metabolites, glutamine, glutamic acid, and n-acetyl-glutamic acid, were significantly increased in the cana (100 µM)-treated cells (p < 0.05). Cana (10 µM) showed a trend of increased glutamic acid and n-acetyl-glutamic acid, but no significance was perceived. No effects on apoptosis-associated metabolites were observed in the dapa- and ertu-treated HUVECs, irrespective of the applied SGLT2 inhibitor concentrations.

2.6.5. Effects on Purine and Pyrimidine Pathway Activity

The cells treated with cana (both 10 and 100 µM) showed a trend towards increased levels of metabolites in the purine/pyrimidine pathways, whereas a trend towards reduced metabolite levels in the purine/pyrimidine pathways was seen for the cells treated with dapa (0.3 and 3 µM) and ertu (1 and 10 µM).

3. Discussion

In the present study, we investigated the differential cell autonomous effects of the selective SGLT2 inhibitors cana, dapa, ertu, and empa on parameters related to mitochondrial activity (OCR, ECAR, ADP/ATP ratio) and on GU in HUVECs. Further, a metabolomics analysis was performed on the cana-, dapa- and ertu-treated HUVECs to assess their effects on the metabolomes of these cells. Our results confirmed the effects of cana on the mitochondrial activity and GU in the endothelial cells, which have been described previously in other cell models, using supra-pharmacological and pharmacologically relevant concentrations. By contrast, dapa, ertu, and empa did not show any of these effects, even when applied at supra-pharmacological concentrations.
In more detail, we reproduced the effects of cana on the OCR in the HUVECs. At 100 µM, the cana inhibited the OCR by 60% in our study, and weak OCR inhibition was also observed after the treatment with cana at 10 µM, which was consistent with the reported results from several other cell lines [31,32,33]. In agreement with the results obtained in the renal proximal tubule epithelial cells [33], dapa and empa did not inhibit the OCR at concentrations of up to 100 μM in the HUVECs in our study.
The reduction in the OCR by cana was likely the result of the electron transport chain (ETC) complex 1 inhibition in our study. Nonetheless, when interpreting the results, one needs to take into account that endothelial cells, such as HUVECs, do not predominantly rely on their mitochondrial energy supply [44,45,46]. In fact, the majority of the energy supply in HUVECs is generated via anaerobic metabolism (e.g., glycolytic pathway), even in the presence of sufficient oxygen [43,45,47]. Thus, lactate production is sustained even during ETC complex 1 inhibition in HUVECs, and proangiogenic signaling is promoted [48]. Cana treatment (10 and 100 µM) seems to additionally diminish the already low-grade usage of oxygen for cellular ATP generation, and to shift the energy production even further in the direction of anaerobic metabolism, i.e., glycolysis. It is unclear whether these effects would translate to the in vivo situation, since they were only significant when applying 100 µM of the SGLT2 inhibitor, while only trends of these effects were observed in the cana (10 µM)-treated HUVECs.
ETC complex 1 inhibition by cana has also been reported in human renal proximal tubule epithelial cells [33]. Shetty et al. have observed that the inhibition of oxidative phosphorylation in liver cells resulted in a rapid stimulation of glucose transport mediated by the GLUT1 transporter, which occurred without a change in the GLUT1 or GLUT1 mRNA content of the cell [49]. In many cell types, GLUT1 is responsible for basal glucose uptake [50], with decreasing intracellular ATP concentrations leading to increased glucose transport via GLUT1 [51]. Rotenone has stimulated GU and glucose consumption in HepG2 and immortalized mouse myoblast cells when mitochondrial oxidative phosphorylation was inhibited [52]. This stimulatory effect led to increased lactate release and thus contributed to the increased acidification of the extracellular space [52]. We found that cana and ertu inhibited the GU in the HUVECs and that the GU decrease in the HUVECs by cana was caused by the inhibition of the GLUT1. The treatment of the HUVECs with cana did not have a significant effect on the expression of the GLUTs, but an effect of cana on the protein level cannot be excluded.
The metabolomics analyses of cana (10 and 100 µM), ertu (1 and 10 µM), and dapa (0.3 and 3 µM)-treated HUVECs showed that cana, in contrast to dapa and ertu, caused a significant decrease in the metabolites (p < 0.01) used in the glycolysis, which was in agreement with the lower GU observed after the cana treatments.
Low glucose levels cause enhanced mitochondrial respiration in endothelial cells, whereas high glucose levels induce reduced mitochondrial respiration, even when enough oxygen is present [42,44,53]. During the cana treatments (both 10 and 100 µM), the mitochondrial oxidative phosphorylation in the HUVECs was impaired, and the reduced glycolysis rate, together with the impaired mitochondrial oxidative phosphorylation, produced a metabolically quiescent state in the cells to conserve the cellular functions.
The decrease in the ECAR in the HUVECs induced by cana at 100 µM was consistent with the increased lactic acid levels observed in the metabolomics analyses. We assume that the lactate was retained in the cells for further use in the regeneration of NAD+ (anaerobic respiration via pyruvic acid).
The levels of the amino acids used for the production of glucose showed a trend towards elevated levels in the cana (10 µM)-treated HUVECs, and were significantly higher in the cana (100 µM)-treated cells, indicating reduced TCA-cycle activity due to its low usage in endothelial cells, as described previously [53], and/or that the metabolites and pathways supplying it were in a state of suspension due to the blockage of the ETC complex 1 and adjoining pathways.
A potential increase in the TCA-cycle activity could also have been caused by the increased activity of the electron transport chain complex II (ETC2), since ETC2 directly links the TCA cycle and the respiratory chain by catalyzing the oxidation of succinate to fumarate [54]. A potential increase in fatty acid oxidation could lead to heightened activity of the ETC2, as has been described by Pfleger et al. [55].
Secker et al. observed the inhibitory effects of cana (50 µM) on glutamate dehydrogenase (GDH) in human renal proximal tubule epithelial cells after treatment [33]. This effect was concurrent with the inhibition of ETC complex I and the replenishment of the TCA-cycle intermediates and intermediates of the amino acid metabolism. The increased levels of glutamine and glutamate also suggest GDH inhibition by cana in the HUVECs. In HUVECs, glutamate is used to produce gamma-aminobutyric acid via glutamic acid decarboxylase, while about one third of the TCA carbons are derived from glutamine [56].
The inhibition of glutamic acid decarboxylase reduces cellular ATP levels and increases pyruvate and fatty-acid oxidation. The blockage of the glutamine metabolism has also been shown to perturb proliferation dynamics and the adoption of a senescence-like phenotypes [57].
In agreement with the results from Hawley et al., who have studied human embryonic kidney cells [31], cana evoked a concentration-dependent increase in the ADP/ATP ratio in the HUVECs in our study. Since glucose is required for mitochondrial ATP production in endothelial cells [58], the inhibition of GU and mitochondrial oxidative phosphorylation would negatively affect the health status of endothelial cells in the long-term.
Trends in a cellular-depleted redox-state caused by lower levels of metabolites, such as NADPH and NADH, in the cana-treated HUVECs were observed. In particular, the NADPH, a crucial cofactor for endothelial nitric oxide synthase [59], was lowered in the cana-treated HUVECs, thus decreasing the functionality and affecting the cellular ability to deal with reactive oxygen species, which in turn would result in endothelial cell dysfunction.
The HUVECs treated with dapa (3 µM) and ertu (10 µM) showed no significant changes regarding glycolysis, the metabolites used in cell respiration and beta-oxidation. However, these cells showed reduced metabolite levels in the purine and pyrimidine pathways, and these effects were more pronounced with increasing SGLT2 inhibitor concentrations. The levels of metabolites from the purine and pyrimidine pathways were significantly lower after the ertu (10 µM)- and dapa (3 µM) treatments, which were similar to the results described in [37,60], than after the cana (10 µM)- and (100 µM)-treatments in the HUVECs. This reflects a significant difference in the mode of action of ertu and dapa compared to cana. A reduction in purine metabolism after the dapa treatment is described by Stack et al. [60], suggesting a reduction in cell death [60] and, potentially, an improvement in endothelial function [37]. The significant inhibitory effect and the observed differences in the mode of action of cana in relation to those of the other three investigated SGLT2 inhibitors in this study could have been caused by cana´s lower selectivity for SGLT2 over SGLT1 compared to dapa, ertu, and empa. This could be one of the reasons why cana inhibited the mitochondrial activity and GU at supra-pharmacological and pharmacological concentrations, whereas dapa, ertu, and empa did not show effects, even at supra-pharmacological concentrations.

4. Materials and Methods

4.1. SGLT2 Inhibitors and GLUT1 Inhibitor

The synthetic selective SGLT2 inhibitors, cana, empa, and ertu [9], were purchased from Ark Pharma (Arlington Heights, IL, USA). The synthetic selective SGLT2 inhibitor, dapa was obtained from AstraZeneca (Gaithersburg, MD, USA). Initial experiments were performed by applying supra-pharmacological concentrations of cana, dapa, ertu, and empa, and the subsequent experiments were performed with concentrations that equaled the plasma peak concentration (Cmax) of the respective SGLT2 inhibitor elicited by its pharmacological doses. The maximal plasma concentration (Cmax) of cana at daily therapeutic dose of 300 mg at steady state is approximately 10 µM (4445 ng/mL) [61,62], Cmax of dapa at daily therapeutic dose of 10 mg at steady state is about 0.3 µM (122 ng/mL) [63], and 1 µM (436 ng/mL) is the Cmax of ertu at daily therapeutic dose of 15 mg at steady state [64]. Cmax of empa at steady state is 0.63 µM (284 ng/mL) at a daily dose of 25 mg [65].
In contrast to sodium-dependent glucose transporters (SGLT), glucose transporters (GLUT) are facilitative glucose transporters [66] that are overexpressed in many tumors [67]. BAY-876 is an orally active and highly selective glucose transporter 1 (GLUT1) inhibitor with an IC50 of 2 nM. BAY-876 is about 130-fold more selective for GLUT1 than GLUT2, GLUT3, and GLUT4, and a potent blocker of glycolytic metabolism and ovarian cancer growth [68].

4.2. Cell Lines and Cell Culture Conditions

HUVECs (Lonza, Switzerland) were plated and expanded in culture medium according to the manufacturer’s instructions (CC-5035 EGM-PLUS BulletKit Medium; CC-5036 EGM-PLUS Basal Media + CC-4542 EGM-PLUS SingleQuots Kit, Lonza, Switzerland).

4.3. Measurement of the OCR and the ECAR

Data were analyzed as change in OCR and ECAR as a function of SGLT2 inhibitor concentration. Effects on cellular respiration and medium acidification was measured with Agilent/Seahorse Extracellular Flux analyzer XF96 at 37 °C. Basal OCR and ECAR were recorded prior to injection of the respective SGLT2 inhibitor to yield final concentrations of 0.2, 1, 5, 25, and 100 µM (n = 4), followed by five repeated measurements of OCR and ECAR (during approximately 35 min). Rotenone (2 µM) was used as reference for the inhibition of the mitochondrial ETC complex I. DMSO was utilized as blank sample. In total, 15,000 and 25,000 HUVEC cells were seeded in 150-microliter cell culture medium per XF96 plate well 18–24 h prior to experiment. XF96 cartridge was equilibrated with calibration buffer overnight according to manufacturer´s instructions. Cell culture medium was changed to XF Base Medium Minimal DMEM (Agilent, Santa Clara, CA, USA) complemented with 5 mM glucose, 1 mM pyruvate, 1 mM glutamine, and 15 mM HEPES 30–40 min prior to first reading. After the last read-out, cells were fixed in 4% formaldehyde, and subsequently washed with 4% formaldehyde prior to adding Hoechst nuclear stain for cell counting with Image Xpress system from Molecular Devices. OCR and ECAR were normalized based on the cell numbers per well.

4.4. Measurement of the ADP/ATP Ratio

The cellular energy status was assessed by measuring the ADP/ATP ratio, following SGLT2 inhibitor treatment for three hours. ADP and ATP levels in the samples were detected using the ADP/ATP Ratio Assay Kit (ab65313, Abcam, Cambridge, UK). Cells were grown as described above. Prior to the measurement, the culture medium was removed from the plate, and Nucleotide Releasing Buffer (50 µL per 103–104 cells) was added. After five minutes of incubation at room temperature while gently shaking, a total of 100 µL of the reaction mix was added and background luminescence was assessed. Fifty µL of cells (103–104) treated with Nucleotide Releasing Buffer were pipetted into a luminometer plate and samples were read after approximately two minutes in a luminometer. A 10 × ADP converting enzyme was diluted 10-fold with Nucleotide Releasing Buffer. ADP levels in cells were read by adding 10 µL of 1 × ADP converting enzyme and reading the output after two minutes. All samples were run as duplicates in four independent experiments. ADP/ATP ratios were calculated based on the manufacturer´s recommendation.

4.5. Measurement of the GU

Effects on the GU were assessed using 14C-deoxyglucose (14C-DOG) as tracer. In total, 40,000 HUVECs per well were seeded in a Cytostar 96w SPA plate and allowed to adhere overnight prior to the experiment. Cytostar plates were pre-coated with collagen I (0.03 µg/uL). The cells were pre-treated for 15 min with 0.2, 1, 5, 25, and 100 µM of the respective SGLT2 inhibitor. Next, they were pre-diluted in assay medium consisting of XF Base MedSium Minimal DMEM (Agilent 102353) and complemented with 5 mM glucose, 1 mM Pyruvate, 1 mM Glutamine, and 15 mM HEPES prior to the addition of 18 µM 14C-DOG tracer (specific activity: 1 µCi/mL). Uptake of 14C-DOG was monitored in real-time using a micro beta scintillation counter (PerkinElmer) with four readings between 10 and 190 min. Cell monolayers were visually inspected for potential loss of adhesion after the last reading. A total of 25 µM BAY-876 (n = 16) was used as reference for effects on GLUT1-mediated uptake against cana (100 µM, n = 8) [68], and DMSO was used as control (n = 32). The 14C-DOG uptake rate was calculated as fmol DOG per minute, and after subtracting the value for 14C-background. Data were fitted to a four-parameter logistic model to calculate IC50.

4.6. RNA Isolation and cDNA Synthesis

Total RNA was extracted from HUVECs treated with cana (10 and 100 µM) after 4 and 24 h, and purified using the Qiagen RNeasy mini kit (Qiagen, Valencia, CA, USA). Cell lysis was performed in an RNases inactivating buffer provided by the manufacturer. Until RNA extraction, the samples were stored at −70 °C. RNA concentrations and 260/280 absorbance ratios were measured spectrophotometrically using Nanodrop ND-1000. Reverse transcription of the purified RNA was performed immediately after extraction. The cDNA was synthesized using a High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, 168 Third Avenue, Waltham, MA, USA). The reaction was set with 800 ng RNA in 20 µL of solution containing RT-random primers, 100 mM deoxyribonucleotide triphosphates (dNTPs) mix, 2.5 U/µL MultiScribe Reverse Transcriptase, and RT buffer. Real-time qPCRs were performed on Applied Biosystems® QuantStudio™ at 50 °C for two minutes and at 95 °C for ten minutes. Then, denaturation cycles (40 cycles of 95 °C for 15 s) were performed. Subsequently, the annealing and extension were performed, during which fluorescence was measured (40 cycles of 60 °C for 60 s). Data expression levels were recorded as quantification cycles, CT. Data were acquired using real-time PCR software (QuantStudio™, Applied Biosystems, USA). The mean CT values of the triplicate reactions were used in the subsequent analysis. Protein Atlas was used for data comparison.

4.7. Metabolomics

Cultivated HUVECs (median cell number 3.12 × 105) were treated with cana (10 µM and 100 µM), ertu (1 µM and 10 µM), and dapa (0.3 µM and 3 µM) and with vehicle solution (1 µM and 10 µM vehicle solution 1 × phosphate-buffered saline 1 × PBS) as controls with six replicates per treatment (n = 6).
Two hours after treatment, the cell culture samples were extracted, dipped into isotonic NaCl solution with 4.5 g/mL glucose and freeze-dried for 9.5 h, as described previously [69,70]. The cells were then reconstituted in 90 µL H2O. Twenty-five µL of each cell sample were used for the measurement and fifteen µL of each cell sample were used for the quality control sample pool (QC). The cell samples were divided into two smaller batches of 24 samples each, and these two batches were measured in one run. Maximum measurement time per batch was 24 h.
Samples, blanks (pure H2O), QCs, and UltimateMix (UM) were measured in a stratified randomized sequence and analyzed with a Dionex Ultimate 3000 HPLC (Thermo Fisher Scientific, USA) equipped with a reversed-phase Atlantis T3 C18 pre- and analytical column (Waters, USA) with an injection volume of 10 µL, as described previously [70,71,72]. Raw data were converted into mzXML (msConvert, ProteoWizard Toolkit v3.0.5) [73], and the detected m/z of known metabolites were identified using PeakScout [70,74,75] with a reference list containing respective mass and retention times.
Principal component analysis (PCA) was performed, centered and scaled to unit variance (R function prcomp). Missing values were imputed by regularized expectation-maximization (R function impute PCA, estim_ncpPCA). For ANOVA analysis, the treatment was taken as the fixed factor (simple ANOVA; R function gls, missing values were not imputed). For specific group comparisons, pairwise post hoc tests (R function lsmeans) were conducted. Technical variability was 15.9% median relative standard deviation (RSD) in the QC samples for 60 multivariate and univariate analysis (MVA_UVA) metabolites and yielded a further 36 univariate-analysis (UVA) metabolites. The MVA_UVA metabolites were classified according to the Human Metabolome Database [43].

4.8. Statistical Analyses

Results for OCR, ECAR, and GU were presented as mean ± SD. Results for ADP/ATP ratios were presented as mean ± SEM.
Data were analyzed for statistical significance using one-way analysis of variance (ANOVA) corrected for multiple comparisons by Tukey statistical hypothesis testing (GraphPad 8.0.1, San Diego, CA, USA). A p-value < 0.05 between groups was considered as statistically significant.
Statistical analysis for metabolomics was performed with R (R Core Team, v3.4.1, R Foundation for Statistical Computing, Vienna, Austria) and TibcoSpotfire (v7.5.0, TIBCO, Palo Alto, CA, USA). All metabolites were controlled for their analytical quality [70] and graded into two classes: (1) suitable for multivariate analysis and univariate analysis (MVA_UVA) and (2) suitable for univariate analysis (UVA). Quality assessment was performed according to [70], and a separate control for outliers due to sequence position, time-point of measurement, replicate number, cell count, date of sampling, measurement batch, and sample extraction events was conducted via PCA.
To correct for cell number differences and reduce technical variability, median normalization was performed. Each metabolite was normalized to the sample median using the following procedure. First, each metabolite´s peak area was scaled between 0 and 1 within one sample. Next, the median of all 0-to-1-scaled metabolites’ peak areas within one sample was calculated, and each peak area was divided by the sample median value.
Finally, normalized peak areas were log-transformed. Kolmogorov–Smirnov and Brown–Forsythe Levene tests were used to discern normality and homoscedasticity of data distribution. The log10-transformed, normalized data were normally distributed according to the Kolmogorov–Smirnov test (93.8% of all metabolites were normally distributed) and homoscedastic according to the Brown–Forsythe-Levene-type test (92.7% of all metabolites were homoscedastic). The p-values corrected for multi-parameter testing via Benjamini–Hochberg are given in the Supplementary Material. All values and additional information are reported in Supplementary Material.

5. Conclusions

Our data confirmed the reported inhibition of mitochondrial activity by cana at supra-pharmacological and, to some extent, at pharmacological concentrations. By contrast, dapa, ertu, and empa did not show any measurable effect on the mitochondrial activity, glucose uptake, and ECAR in the endothelial cells, even at concentrations that were much higher than the Cmax evoked by their pharmacological doses. The supra-pharmacological concentrations of cana, but not dapa, ertu, or empa, affected multiple cellular pathways (glycolysis, beta-oxidation, and cellular energy state) and blocked the GLUT1.
At 100 µM, cana showed a combined reduction in glycolysis and mitochondria respiration. This diminished cellular-energy state could lead to insufficient angiogenesis [52] and could play a role in microvascular complication, through involvement in, for example, impaired wound healing in people with type 2 diabetes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23147966/s1.

Author Contributions

Conceptualization: H.S., C.M., E.Z., T.R.P., H.-C.Y. and X.-R.P.; methodology: E.Z. and T.R.P.; study investigation: E.Z., P.K., H.-C.Y. and S.H.; data curation, visualization and interpretation: E.Z., P.K., B.B., C.M., H.-C.Y., S.H., C.S.E., R.E., S.M., J.O. and X.-R.P.; original draft preparation: B.B., E.Z., H.S., P.K., C.M. and H.-C.Y.; review and editing: B.B., E.Z., P.K., T.R.P., C.M., H.S. and H.-C.Y.; supervision of the study: C.M., T.R.P., H.-C.Y. and X.-R.P.; project administration: E.Z., C.M. and B.B.; funding acquisition, C.M., H.S. and T.R.P.; All authors contributed to the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The study was partly funded by Astra Zeneca.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors.

Acknowledgments

The authors thank Karin Platzer (HEALTH—Joanneum Research Forschungsgesellschaft m.b.H.) for laboratory work.

Conflicts of Interest

This study was partly funded by Astra Zeneca. HS received investigator-initiated study funding from Boehringer Ingelheim, Eli Lilly, MSD, NovoNordisk, and Sanofi-Aventis. HS is on the advisory board of and/or received speaker honoraria from Amgen, Astra Zeneca, Bayer, Boehringer-Ingelheim, Daiichi Sankyo, Eli Lilly, Kapsch, MSD, Novo Nordisk, and Sanofi-Aventis. TRP is an advisory board/panel member of Novo Nordisk A/S, Arecor, Astra Zeneca, Boehringer Ingelheim, Convatec, Eli Lilly, and Sanofi, and received research support from Novo Nordisk A/S and Astra Zeneca. Xiao-Rong Peng, Russell Esterline, and Jan Oscarsson are full-time employees of AstraZeneca. The remaining authors have no conflict of interest to disclose.

References

  1. Hsia, D.S.; Grove, O.; Cefalu, W.T. An update on sodium-glucose co-transporter-2 inhibitors for the treatment of diabetes mellitus. Curr. Opin. Endocrinol. Diabetes Obes. 2017, 24, 73–79. [Google Scholar] [CrossRef] [PubMed]
  2. Nguyen, V.K.; White, J.R. Overview of Ertugliflozin. Clin. Diabetes 2019, 37, 176–178. [Google Scholar] [CrossRef] [PubMed]
  3. Inagaki, N.; Kondo, K.; Yoshinari, T.; Maruyama, N.; Susuta, Y.; Kuki, H. Efficacy and safety of canagliflozin in Japanese patients with type 2 diabetes: A randomized, double-blind, placebo-controlled, 12-week study. Diabetes Obes. Metab. 2013, 15, 1136–1145. [Google Scholar] [CrossRef] [PubMed]
  4. American Diabetes Association. 9. Pharmacologic Approaches to Glycemic Treatment: Standards of Medical Care in Diabetes-2020. Diabetes Care 2020, 43, S98–S110. [Google Scholar] [CrossRef] [Green Version]
  5. Davies, M.J.; D’Alessio, D.A.; Fradkin, J.; Kernan, W.N.; Mathieu, C.; Mingrone, G.; Rossing, P.; Tsapas, A.; Wexler, D.J.; Buse, J.B. Management of Hyperglycemia in Type 2 Diabetes, 2018. A Consensus Report by the American Diabetes Association (ADA) and the European Association for the Study of Diabetes (EASD). Diabetes Care 2018, 41, 2669–2701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Lupsa, B.C.; Inzucchi, S.E. Use of SGLT2 inhibitors in type 2 diabetes: Weighing the risks and benefits. Diabetologia 2018, 61, 2118–2125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Perkovic, V.; Jardine, M.J.; Neal, B.; Bompoint, S.; Heerspink, H.J.L.; Charytan, D.M.; Edwards, R.; Agarwal, R.; Bakris, G.; Bull, S.; et al. Canagliflozin and Renal Outcomes in Type 2 Diabetes and Nephropathy. N. Engl. J. Med. 2019, 380, 2295–2306. [Google Scholar] [CrossRef] [Green Version]
  8. Wiviott, S.D.; Raz, I.; Bonaca, M.P.; Mosenzon, O.; Kato, E.T.; Cahn, A.; Silverman, M.G.; Zelniker, T.A.; Kuder, J.F.; Murphy, S.A.; et al. Dapagliflozin and Cardiovascular Outcomes in Type 2 Diabetes. N. Engl. J. Med. 2019, 380, 347–357. [Google Scholar] [CrossRef]
  9. Cinti, F.; Moffa, S.; Impronta, F.; Cefalo, C.M.A.; Sun, V.A.; Sorice, G.P.; Mezza, T.; Giaccari, A. Spotlight on ertugliflozin and its potential in the treatment of type 2 diabetes: Evidence to date. Drug Des. Dev. Ther. 2017, ume 11, 2905–2919. [Google Scholar] [CrossRef] [Green Version]
  10. Ehrenkranz, J.R.L.; Lewis, N.; Kahn, C.R.; Roth, J. Phlorizin: A review. Diabetes/Metab. Res. Rev. 2005, 21, 31–38. [Google Scholar] [CrossRef]
  11. Rossetti, L.; Smith, D.; Shulman, G.I.; Papachristou, D.; DeFronzo, R.A. Correction of hyperglycemia with phlorizin normalizes tissue sensitivity to insulin in diabetic rats. J. Clin. Investig. 1987, 79, 1510–1515. [Google Scholar] [CrossRef] [PubMed]
  12. Fujimori, Y.; Katsuno, K.; Ojima, K.; Nakashima, I.; Nakano, S.; Ishikawa-Takemura, Y.; Kusama, H.; Isaji, M. Sergliflozin etabonate, a selective SGLT2 inhibitor, improves glycemic control in streptozotocin-induced diabetic rats and Zucker fatty rats. Eur. J. Pharmacol. 2009, 609, 148–154. [Google Scholar] [CrossRef] [PubMed]
  13. Liang, Y.; Arakawa, K.; Ueta, K.; Matsushita, Y.; Kuriyama, C.; Martin, T.; Du, F.; Liu, Y.; Xu, J.; Conway, B.; et al. Effect of Canagliflozin on Renal Threshold for Glucose, Glycemia, and Body Weight in Normal and Diabetic Animal Models. PLoS ONE 2012, 7, e30555. [Google Scholar] [CrossRef] [Green Version]
  14. Kaushal, S.; Singh, H.; Thangaraju, P.; Singh, J. Canagliflozin: A novel SGLT2 inhibitor for type 2 diabetes mellitus. N. Am. J. Med Sci. 2014, 6, 107–113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ohgaki, R.; Wei, L.; Yamada, K.; Hara, T.; Kuriyama, C.; Okuda, S.; Ueta, K.; Shiotani, M.; Nagamori, S.; Kanai, Y. Interaction of the Sodium/Glucose Cotransporter (SGLT) 2 inhibitor Canagliflozin with SGLT1 and SGLT2. J. Pharmacol. Exp. Ther. 2016, 358, 94–102. [Google Scholar] [CrossRef] [Green Version]
  16. Lapuerta, P.; Zambrowicz, B.; Strumph, P.; Sands, A. Development of sotagliflozin, a dual sodium-dependent glucose transporter 1/2 inhibitor. Diabetes Vasc. Dis. Res. 2015, 12, 101–110. [Google Scholar] [CrossRef] [Green Version]
  17. Bhatt, D.L.; Szarek, M.; Steg, P.G.; Cannon, C.P.; Leiter, L.A.; McGuire, D.K.; Lewis, J.B.; Riddle, M.C.; Voors, A.A.; Metra, M.; et al. Sotagliflozin in Patients with Diabetes and Recent Worsening Heart Failure. N. Engl. J. Med. 2021, 384, 117–128. [Google Scholar] [CrossRef]
  18. Bhatt, D.L.; Szarek, M.; Pitt, B.; Cannon, C.P.; Leiter, L.A.; McGuire, D.K.; Lewis, J.B.; Riddle, M.C.; Inzucchi, S.E.; Kosiborod, M.N.; et al. Sotagliflozin in Patients with Diabetes and Chronic Kidney Disease. N. Engl. J. Med. 2021, 384, 129–139. [Google Scholar] [CrossRef]
  19. Rieg, T.; Vallon, V. Development of SGLT1 and SGLT2 inhibitors. Diabetologia 2018, 61, 2079–2086. [Google Scholar] [CrossRef] [Green Version]
  20. Flores, E.; Santos-Gallego, C.G.; Diaz-Mejía, N.; Badimon, J.J. Do the SGLT-2 Inhibitors Offer More than Hypoglycemic Activity? Cardiovasc. Drugs Ther. 2018, 32, 213–222. [Google Scholar] [CrossRef]
  21. Osataphan, S.; Macchi, C.; Singhal, G.; Chimene-Weiss, J.; Sales, V.; Kozuka, C.; Dreyfuss, J.; Pan, H.; Tangcharoenpaisan, Y.; Morningstar, J.; et al. SGLT2 inhibition reprograms systemic metabolism via FGF21-dependent and -independent mechanisms. JCI Insight 2019, 4, e123130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Zinman, B.; E Inzucchi, S.; Lachin, J.M.; Wanner, C.; Ferrari, R.; Fitchett, D.; Bluhmki, E.; Hantel, S.; Kempthorne-Rawson, J.; Newman, J.; et al. Rationale, design, and baseline characteristics of a randomized, placebo-controlled cardiovascular outcome trial of empagliflozin (EMPA-REG OUTCOME™). Cardiovasc. Diabetol. 2014, 13, 102. [Google Scholar] [CrossRef] [Green Version]
  23. Zinman, B.; Wanner, C.; Lachin, J.M.; Fitchett, D.; Bluhmki, E.; Hantel, S.; Mattheus, M.; Devins, T.; Johansen, O.E.; Woerle, H.J.; et al. Empagliflozin, Cardiovascular Outcomes, and Mortality in Type 2 Diabetes. N. Engl. J. Med. 2015, 373, 2117–2128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Europäische Arzneimittelagentur. Zynquista–Zusammenfassung der Merkmale des Arzneimittels. 2021. Available online: https://www.ema.europa.eu/en/documents/product-information/zynquista-epar-product-information_de.pdf. (accessed on 20 April 2022).
  25. U.S. Food and Drug Administration. Proposal To Refuse To Approve a New Drug Application for Sotagliflozin Oral Tablets, 200 Milligrams and 400 Milligrams; Opportunity for a Hearing. Federal Register. 2021, 86, 12471–12473. [Google Scholar]
  26. Ferrannini, E.; Solini, A. SGLT2 inhibition in diabetes mellitus: Rationale and clinical prospects. Nat. Rev. Endocrinol. 2012, 8, 495–502. [Google Scholar] [CrossRef]
  27. Neal, B.; Perkovic, V.; Mahaffey, K.W.; de Zeeuw, D.; Fulcher, G.; Erondu, N.; Shaw, W.; Law, G.; Desai, M.; Matthews, D.R.; et al. Canagliflozin and Cardiovascular and Renal Events in Type 2 Diabetes. N. Engl. J. Med. 2017, 377, 644–657. [Google Scholar] [CrossRef]
  28. Fadini, G.P.; Avogaro, A. SGLT2 inhibitors and amputations in the US FDA Adverse Event Reporting System. Lancet Diabetes Endocrinol. 2017, 5, 680–681. [Google Scholar] [CrossRef] [Green Version]
  29. Birkeland, K.I.; Jørgensen, M.E.; Carstensen, B.; Persson, F.; Gulseth, H.L.; Thuresson, M.; Fenici, P.; Nathanson, D.; Nyström, T.; Eriksson, J.W.; et al. Cardiovascular mortality and morbidity in patients with type 2 diabetes following initiation of sodium-glucose co-transporter-2 inhibitors versus other glucose-lowering drugs (CVD-REAL Nordic): A multinational observational analysis. Lancet Diabetes Endocrinol. 2017, 5, 709–717. [Google Scholar] [CrossRef]
  30. U.S. FDA. FDA Drug Safety Communication: FDA Confirms Increased Risk of Leg and Foot Amputations with the Diabetes Medicine Canagliflozin (Invokana, Invokamet, Invokamet XR). Available online: https://www.fda.gov/media/104870/download (accessed on 20 April 2022).
  31. Hawley, S.A.; Ford, R.J.; Smith, B.K.; Gowans, G.J.; Mancini, S.J.; Pitt, R.D.; Day, E.A.; Salt, I.P.; Steinberg, G.R.; Hardie, D.G. The Na+/Glucose Cotransporter Inhibitor Canagliflozin Activates AMPK by Inhibiting Mitochondrial Function and Increasing Cellular AMP Levels. Diabetes 2016, 65, 2784–2794. [Google Scholar] [CrossRef] [Green Version]
  32. Villani, L.A.; Smith, B.K.; Marcinko, K.; Ford, R.J.; Broadfield, L.A.; Green, A.E.; Houde, V.P.; Muti, P.; Tsakiridis, T.; Steinberg, G.R. The diabetes medication Canagliflozin reduces cancer cell proliferation by inhibiting mitochondrial complex-I supported respiration. Mol. Metab. 2016, 5, 1048–1056. [Google Scholar] [CrossRef]
  33. Secker, P.F.; Beneke, S.; Schlichenmaier, N.; Delp, J.; Gutbier, S.; Leist, M.; Dietrich, D.R. Canagliflozin mediated dual inhibition of mitochondrial glutamate dehydrogenase and complex I: An off-target adverse effect. Cell Death Dis. 2018, 9, 226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Uthman, L.; Homayr, A.; Juni, R.P.; Spin, E.L.; Kerindongo, R.; Boomsma, M.; Hollmann, M.W.; Preckel, B.; Koolwijk, P.; Van Hinsbergh, V.W.M.; et al. Empagliflozin and Dapagliflozin Reduce ROS Generation and Restore NO Bioavailability in Tumor Necrosis Factor α-Stimulated Human Coronary Arterial Endothelial Cells. Cell. Physiol. Biochem. 2019, 53, 865–886. [Google Scholar] [CrossRef] [PubMed]
  35. Behnammanesh, G.; Durante, Z.E.; Peyton, K.J.; Martinez-Lemus, L.A.; Brown, S.M.; Bender, S.B.; Durante, W. Canagliflozin Inhibits Human Endothelial Cell Proliferation and Tube Formation. Front. Pharmacol. 2019, 10, 362. [Google Scholar] [CrossRef] [PubMed]
  36. Nakano, D.; Kawaguchi, T.; Iwamoto, H.; Hayakawa, M.; Koga, H.; Torimura, T. Effects of canagliflozin on growth and metabolic reprograming in hepatocellular carcinoma cells: Multi-omics analysis of metabolomics and absolute quantification proteomics (iMPAQT). PLoS ONE 2020, 15, e0232283. [Google Scholar] [CrossRef]
  37. Mulder, S.; Hammarstedt, A.; Nagaraj, S.B.; Nair, V.; Ju, W.; Hedberg, J.; Greasley, P.J.; Eriksson, J.W.; Oscarsson, J.; Heerspink, H.J.L. A metabolomics-based molecular pathway analysis of how the sodium-glucose co-transporter-2 inhibitor dapagliflozin may slow kidney function decline in patients with diabetes. Diabetes Obes. Metab. 2020, 22, 1157–1166. [Google Scholar] [CrossRef]
  38. Tabit, C.E.; Chung, W.B.; Hamburg, N.M.; Vita, J.A. Endothelial dysfunction in diabetes mellitus: Molecular mechanisms and clinical implications. Rev. Endocr. Metab. Disord. 2010, 11, 61–74. [Google Scholar] [CrossRef] [Green Version]
  39. Paulus, W.J. Endothelial control of vascular and myocardial function in heart failure. Cardiovasc. Drugs Ther. 1994, 8, 437–446. [Google Scholar] [CrossRef]
  40. Vita, J.A.; Hamburg, N. Does endothelial dysfunction contribute to the clinical status of patients with peripheral arterial disease? Can. J. Cardiol. 2010, 26, 45A–50A. [Google Scholar] [CrossRef]
  41. Kopel, T.; Kaufman, J.S.; Hamburg, N.; Sampalis, J.S.; Vita, J.A.; Dember, L.M. Endothelium-Dependent and -Independent Vascular Function in Advanced Chronic Kidney Disease. Clin. J. Am. Soc. Nephrol. 2017, 12, 1588–1594. [Google Scholar] [CrossRef] [Green Version]
  42. Eelen, G.; De Zeeuw, P.; Treps, L.; Harjes, U.; Wong, B.; Carmeliet, P. Endothelial Cell Metabolism. Physiol. Rev. 2018, 98, 3–58. [Google Scholar] [CrossRef]
  43. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vázquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The human metabolome database for 2018. Nucleic Acids Res. 2018, 46, D608–D617. [Google Scholar] [CrossRef] [PubMed]
  44. Koziel, A.; Woyda-Ploszczyca, A.; Kicinska, A.; Jarmuszkiewicz, W. The influence of high glucose on the aerobic metabolism of endothelial EA.hy926 cells. Pflügers Arch.-Eur. J. Physiol. 2012, 464, 657–669. [Google Scholar] [CrossRef] [Green Version]
  45. Davidson, S.M.; Duchen, M.R. Endothelial Mitochondria. Circ. Res. 2007, 100, 1128–1141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. De Bock, K.; Georgiadou, M.; Schoors, S.; Kuchnio, A.; Wong, B.W.; Cantelmo, A.R.; Quaegebeur, A.; Ghesquière, B.; Cauwenberghs, S.; Eelen, G.; et al. Role of PFKFB3-Driven Glycolysis in Vessel Sprouting. Cell 2013, 154, 651–663. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Davidson, S.M. Endothelial mitochondria and heart disease. Cardiovasc. Res. 2010, 88, 58–66. [Google Scholar] [CrossRef] [PubMed]
  48. Ruan, G.-X.; Kazlauskas, A. Lactate Engages Receptor Tyrosine Kinases Axl, Tie2, and Vascular Endothelial Growth Factor Receptor 2 to Activate Phosphoinositide 3-Kinase/Akt and Promote Angiogenesis. J. Biol. Chem. 2013, 288, 21161–21172. [Google Scholar] [CrossRef] [Green Version]
  49. Shetty, M.; Loeb, J.N.; Ismail-Beigi, F. Enhancement of glucose transport in response to inhibition of oxidative metabolism: Pre- and posttranslational mechanisms. Am. J. Physiol. Physiol. 1992, 262, C527–C532. [Google Scholar] [CrossRef]
  50. Zhou, L.; Cryan, E.V.; D’Andrea, M.R.; Belkowski, S.; Conway, B.R.; Demarest, K.T. Human cardiomyocytes express high level of Na+/glucose cotransporter 1 (SGLT1). J. Cell. Biochem. 2003, 90, 339–346. [Google Scholar] [CrossRef]
  51. Shetty, M.; Loeb, J.N.; Vikstrom, K.; Ismail-Beigi, F. Rapid activation of GLUT-1 glucose transporter following inhibition of oxidative phosphorylation in Clone 9 cells. J. Biol. Chem. 1993, 268, 17225–17232. [Google Scholar] [CrossRef]
  52. Hou, W.-L.; Yin, J.; Alimujiang, M.; Yu, X.-Y.; Ai, L.-G.; Bao, Y.-Q.; Liu, F.; Jia, W.-P. Inhibition of mitochondrial complex I improves glucose metabolism independently of AMPK activation. J. Cell. Mol. Med. 2018, 22, 1316–1328. [Google Scholar] [CrossRef] [Green Version]
  53. Krützfeldt, A.; Spahr, R.; Mertens, S.; Siegmund, B.; Piper, H. Metabolism of exogenous substrates by coronary endothelial cells in culture. J. Mol. Cell. Cardiol. 1990, 22, 1393–1404. [Google Scholar] [CrossRef]
  54. Vanova, K.H.; Kraus, M.; Neuzil, J.; Rohlena, J. Mitochondrial complex II and reactive oxygen species in disease and therapy. Redox Rep. 2020, 25, 26–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Pfleger, J.; He, M.; Abdellatif, M. Mitochondrial complex II is a source of the reserve respiratory capacity that is regulated by metabolic sensors and promotes cell survival. Cell Death Dis. 2015, 6, e1835. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Sen, S.; Roy, S.; Bandyopadhyay, G.; Scott, B.; Xiao, D.; Ramadoss, S.; Mahata, S.K.; Chaudhuri, G. γ-Aminobutyric Acid Is Synthesized and Released by the Endothelium. Circ. Res. 2016, 119, 621–634. [Google Scholar] [CrossRef] [PubMed]
  57. Unterluggauer, H.; Mazurek, S.; Lener, B.; Hütter, E.; Eigenbrodt, E.; Zwerschke, W.; Jansen-Dürr, P. Premature senescence of human endothelial cells induced by inhibition of glutaminase. Biogerontology 2008, 9, 247–259. [Google Scholar] [CrossRef] [PubMed]
  58. Yetkin-Arik, B.; Vogels, I.M.C.; Neyazi, N.; Van Duinen, V.; Houtkooper, R.H.; Van Noorden, C.J.F.; Klaassen, I.; Schlingemann, R.O. Endothelial tip cells in vitro are less glycolytic and have a more flexible response to metabolic stress than non-tip cells. Sci. Rep. 2019, 9, 10414. [Google Scholar] [CrossRef] [Green Version]
  59. Leopold, J.A.; Cap, A.; Scribner, A.W.; Stanton, R.C.; Loscalzo, J. Glucose-6-phosphate dehydrogenase deficiency promotes endothelial oxidant stress and decreases endothelial nitric oxide bioavailability. FASEB J. 2001, 15, 1771–1773. [Google Scholar] [CrossRef] [Green Version]
  60. Stack, A.G.; Han, D.; Goldwater, R.; Johansson, S.; Dronamraju, N.; Oscarsson, J.; Johnsson, E.; Parkinson, J.; Erlandsson, F. Dapagliflozin Added to Verinurad Plus Febuxostat Further Reduces Serum Uric Acid in Hyperuricemia: The QUARTZ Study. J. Clin. Endocrinol. Metab. 2021, 106, e2347–e2356. [Google Scholar] [CrossRef]
  61. Devineni, D.; Curtin, C.R.; Polidori, D.; Gutierrez, M.J.; Murphy, J.; Rusch, S.; Rothenberg, P.L. Pharmacokinetics and Pharmacodynamics of Canagliflozin, a Sodium Glucose Co-Transporter 2 Inhibitor, in Subjects With Type 2 Diabetes Mellitus. J. Clin. Pharmacol. 2013, 53, 601–610. [Google Scholar] [CrossRef]
  62. Iijima, H.; Kifuji, T.; Maruyama, N.; Inagaki, N. Pharmacokinetics, Pharmacodynamics, and Safety of Canagliflozin in Japanese Patients with Type 2 Diabetes Mellitus. Adv. Ther. 2015, 32, 768–782. [Google Scholar] [CrossRef] [Green Version]
  63. Komoroski, B.; Vachharajani, N.; Boulton, D.; Kornhauser, D.; Geraldes, M.; Li, L.; Pfister, M. Dapagliflozin, a Novel SGLT2 Inhibitor, Induces Dose-Dependent Glucosuria in Healthy Subjects. Clin. Pharmacol. Ther. 2009, 85, 520–526. [Google Scholar] [CrossRef] [PubMed]
  64. Sahasrabudhe, V.; Terra, S.G.; Hickman, A.; Saur, D.; Shi, H.; O’Gorman, M.; Zhou, Z.; Cutler, D.L. The Effect of Renal Impairment on the Pharmacokinetics and Pharmacodynamics of Ertugliflozin in Subjects With Type 2 Diabetes Mellitus. J. Clin. Pharmacol. 2017, 57, 1432–1443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Heise, T.; Seman, L.J.; Macha, S.; Jones, P.; Marquart, A.; Pinnetti, S.; Woerle, H.J.; Dugi, K. Safety, Tolerability, Pharmacokinetics, and Pharmacodynamics of Multiple Rising Doses of Empagliflozin in Patients with Type 2 Diabetes Mellitus. Diabetes Ther. 2013, 4, 331–345. [Google Scholar] [CrossRef] [Green Version]
  66. Fu, X.; Zhang, G.; Liu, R.; Wei, J.; Zhang-Negrerie, D.; Jian, X.; Gao, Q. Mechanistic Study of Human Glucose Transport Mediated by GLUT1. J. Chem. Inf. Model. 2016, 56, 517–526. [Google Scholar] [CrossRef] [PubMed]
  67. Medina, R.A.; Owen, G.I. Glucose transporters: Expression, regulation and cancer. Biol. Res. 2002, 35, 9–26. [Google Scholar] [CrossRef] [PubMed]
  68. Siebeneicher, H.; Cleve, A.; Rehwinkel, H.; Neuhaus, R.; Heisler, I.; Müller, T.; Bauser, M.; Buchmann, B. Identification and Optimization of the First Highly Selective GLUT1 Inhibitor BAY-876. ChemMedChem 2016, 11, 2261–2271. [Google Scholar] [CrossRef] [Green Version]
  69. Bordag, N.; Rennfahrt, U.; Nachtigall, J.; Maldonado, S.; Reszka, R.; Ramirez-Hernandez, T.; Kamp, H.; Fux, E.; Van Ravenzwaay, B. Fast Sampling of Adherent Cell Cultures for Optimal Metabolomics Results. Metabolomics 2016, 6, 1–4. [Google Scholar] [CrossRef] [Green Version]
  70. Vogel, F.C.; Bordag, N.; Zügner, E.; Trajkovic-Arsic, M.; Chauvistré, H.; Shannan, B.; Váraljai, R.; Horn, S.; Magnes, C.; Siveke, J.T.; et al. Targeting the H3K4 Demethylase KDM5B Reprograms the Metabolome and Phenotype of Melanoma Cells. J. Investig. Dermatol. 2019, 139, 2506–2516.e10. [Google Scholar] [CrossRef]
  71. Buescher, J.M.; Moco, S.; Sauer, U.; Zamboni, N. Ultrahigh Performance Liquid Chromatography−Tandem Mass Spectrometry Method for Fast and Robust Quantification of Anionic and Aromatic Metabolites. Anal. Chem. 2010, 82, 4403–4412. [Google Scholar] [CrossRef]
  72. Trausinger, G.; Gruber, C.; Krahulec, S.; Magnes, C.; Nidetzky, B.; Klimacek, M. Identification of novel metabolic interactions controlling carbon flux from xylose to ethanol in natural and recombinant yeasts. Biotechnol. Biofuels 2015, 8, 157. [Google Scholar] [CrossRef] [Green Version]
  73. Chambers, M.C.; Maclean, B.; Burke, R.; Amodei, D.; Ruderman, D.L.; Neumann, S.; Gatto, L.; Fischer, B.; Pratt, B.; Egertson, J.; et al. A cross-platform toolkit for mass spectrometry and proteomics. Nat. Biotechnol. 2012, 30, 918–920. [Google Scholar] [CrossRef] [PubMed]
  74. Fröhlich, E.E.; Farzi, A.; Mayerhofer, R.; Reichmann, F.; Jačan, A.; Wagner, B.; Zinser, E.; Bordag, N.; Magnes, C.; Fröhlich, E.; et al. Cognitive impairment by antibiotic-induced gut dysbiosis: Analysis of gut microbiota-brain communication. Brain, Behav. Immun. 2016, 56, 140–155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Mueller, K.M.; Hartmann, K.; Kaltenecker, D.; Vettorazzi, S.; Bauer, M.; Mauser, L.; Amann, S.; Jall, S.; Fischer, K.; Esterbauer, H.; et al. Adipocyte glucocorticoid receptor deficiency attenuates aging- and HFD-induced obesity, and impairs the feeding-fasting transition. Diabetes 2017, 66, 272–286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Effects of SGLT2 inhibitors on OCR. Effect of DMSO; DMSO + SGLT2 inhibitor (100 µM); cana, dapa, ertu, and empa on OCR and rotenone (2 µM) in (A) OCR as % average change from baseline and (B) OCR after 40 min (n = 16 (cana and dapa), n = 8 (ertu and empa), and n = 32 (DMSO); **** p < 0.0001, * p = 0.041, ns = not significant, n = replicates of samples). (C) ΔOCR for cana (n = 12).
Figure 1. Effects of SGLT2 inhibitors on OCR. Effect of DMSO; DMSO + SGLT2 inhibitor (100 µM); cana, dapa, ertu, and empa on OCR and rotenone (2 µM) in (A) OCR as % average change from baseline and (B) OCR after 40 min (n = 16 (cana and dapa), n = 8 (ertu and empa), and n = 32 (DMSO); **** p < 0.0001, * p = 0.041, ns = not significant, n = replicates of samples). (C) ΔOCR for cana (n = 12).
Ijms 23 07966 g001
Figure 2. Effects of SGLT2 inhibitors on ECAR. Effect of DMSO and DMSO + SGLT2 inhibitor (100 µM); cana, dapa, ertu and empa on (A) ECAR as % change from baseline (one experiment). (B) ECAR after 40 min (two independent experiments); **** p < 0.0001; ** p = 0.0021; ## p = 0.0026; n = 62 (DMSO), n = 16 (cana), 15 (dapa) and 8 (ertu and empa).
Figure 2. Effects of SGLT2 inhibitors on ECAR. Effect of DMSO and DMSO + SGLT2 inhibitor (100 µM); cana, dapa, ertu and empa on (A) ECAR as % change from baseline (one experiment). (B) ECAR after 40 min (two independent experiments); **** p < 0.0001; ** p = 0.0021; ## p = 0.0026; n = 62 (DMSO), n = 16 (cana), 15 (dapa) and 8 (ertu and empa).
Ijms 23 07966 g002
Figure 3. Effects of SGLT2 inhibitors on the ADP/ATP ratio. Effect of DMSO and DMSO + SGLT2 inhibitor on the ADP/ATP ratio in HUVECs at pharmacologically relevant and supra-pharmacological concentrations of cana (10 and 100 µM), dapa (0.3 and 3 µM), ertu (1 and 10 µM). *** p < 0.004, ns = not significant. Ordinary one-way analysis of variance (ANOVA).
Figure 3. Effects of SGLT2 inhibitors on the ADP/ATP ratio. Effect of DMSO and DMSO + SGLT2 inhibitor on the ADP/ATP ratio in HUVECs at pharmacologically relevant and supra-pharmacological concentrations of cana (10 and 100 µM), dapa (0.3 and 3 µM), ertu (1 and 10 µM). *** p < 0.004, ns = not significant. Ordinary one-way analysis of variance (ANOVA).
Ijms 23 07966 g003
Figure 4. Effect of SGLT2 inhibitors on 14C-DOG uptake (A) and inhibition of GU by cana and BAY-876 (B). (A) 14C-DOG uptake for cana, dapa, ertu, and empa and in controls (dotted lines). ** p = 0.0082 (cana vs. DMSO), **** p < 0.0001 (cana and ertu vs. DMSO), # p = 0.0411 (dapa vs. DMSO), ### p = 0.0003 (dapa vs. DMSO); n = 28 (DMSO), n = 11–12 (cana), n = 11–12 (dapa), n = 7–8 (ertu), and n = 4 (empa). (B) Inhibitory effects on GU in the presence of cana (100 µM) and GLUT1 inhibitor, BAY-876 (25 µM). **** p < 0.0001 (n = 32 (DMSO), n = 8 (cana), and n = 16 (BAY-876)). Effects of cana on GLUT gene expression.
Figure 4. Effect of SGLT2 inhibitors on 14C-DOG uptake (A) and inhibition of GU by cana and BAY-876 (B). (A) 14C-DOG uptake for cana, dapa, ertu, and empa and in controls (dotted lines). ** p = 0.0082 (cana vs. DMSO), **** p < 0.0001 (cana and ertu vs. DMSO), # p = 0.0411 (dapa vs. DMSO), ### p = 0.0003 (dapa vs. DMSO); n = 28 (DMSO), n = 11–12 (cana), n = 11–12 (dapa), n = 7–8 (ertu), and n = 4 (empa). (B) Inhibitory effects on GU in the presence of cana (100 µM) and GLUT1 inhibitor, BAY-876 (25 µM). **** p < 0.0001 (n = 32 (DMSO), n = 8 (cana), and n = 16 (BAY-876)). Effects of cana on GLUT gene expression.
Ijms 23 07966 g004
Figure 5. ANOVA volcano plots for cana (10 µM)-treated (A) and cana (100 µM)-treated HUVECs (B). ANOVA exhibited significant differences between the control groups for 96 metabolites. Levels of metabolites on the left side of the black vertical line decreased; levels of metabolites on the right side of the black vertical line increased. Y-axis: p-values (inverse log scale); x-axis: metabolite ratios (median normalized data; simple treatment). Left black dotted vertical line: ratio of 0.8. Right black dotted vertical line: ratio of 1.2. The red horizontal line denotes a p-value of 0.01, the dashed orange horizontal line denotes a p-value of 0.05, and the gray horizontal line denotes a p-value of 0.1. Metabolites are colored according to the Human Metabolome Database classes [43]. Metabolites supplying the TCA cycle are marked with filled crosses; glycolysis metabolites by filled minus; metabolites used in beta-oxidation with filled circles; metabolites used as energy (co-)factors with filled stars; metabolites used in apoptotic pathways are marked with empty symbols and others with snowflake-like symbols.
Figure 5. ANOVA volcano plots for cana (10 µM)-treated (A) and cana (100 µM)-treated HUVECs (B). ANOVA exhibited significant differences between the control groups for 96 metabolites. Levels of metabolites on the left side of the black vertical line decreased; levels of metabolites on the right side of the black vertical line increased. Y-axis: p-values (inverse log scale); x-axis: metabolite ratios (median normalized data; simple treatment). Left black dotted vertical line: ratio of 0.8. Right black dotted vertical line: ratio of 1.2. The red horizontal line denotes a p-value of 0.01, the dashed orange horizontal line denotes a p-value of 0.05, and the gray horizontal line denotes a p-value of 0.1. Metabolites are colored according to the Human Metabolome Database classes [43]. Metabolites supplying the TCA cycle are marked with filled crosses; glycolysis metabolites by filled minus; metabolites used in beta-oxidation with filled circles; metabolites used as energy (co-)factors with filled stars; metabolites used in apoptotic pathways are marked with empty symbols and others with snowflake-like symbols.
Ijms 23 07966 g005
Figure 6. Scheme of the affected metabolites and pathways in cana-treated HUVECs. Metabolite levels significantly elevated by cana (100 µM) and changed in the same direction by cana (10 µM) are marked with a green up-arrow, metabolite levels lowered by cana (100 µM) are marked with a red down-arrow.
Figure 6. Scheme of the affected metabolites and pathways in cana-treated HUVECs. Metabolite levels significantly elevated by cana (100 µM) and changed in the same direction by cana (10 µM) are marked with a green up-arrow, metabolite levels lowered by cana (100 µM) are marked with a red down-arrow.
Ijms 23 07966 g006
Figure 7. Scheme of affected metabolites and pathways in cana-treated HUVECs. Metabolite levels that were elevated after cana treatment are marked with a green up-arrow, metabolite levels lowered after cana treatment are marked with a red down-arrow. Metabolites that appear more than once in the graph are only marked once with an arrow. Significantly changed metabolite levels in cana (10 µM)- and cana (100 µM)-treated cells are marked with a filled arrow and metabolite levels that are significantly changed only in cana (100 µM)-treated cells are marked with a shaded arrow. The red line at the GLUT 1 transporter indicates transporter inhibition.
Figure 7. Scheme of affected metabolites and pathways in cana-treated HUVECs. Metabolite levels that were elevated after cana treatment are marked with a green up-arrow, metabolite levels lowered after cana treatment are marked with a red down-arrow. Metabolites that appear more than once in the graph are only marked once with an arrow. Significantly changed metabolite levels in cana (10 µM)- and cana (100 µM)-treated cells are marked with a filled arrow and metabolite levels that are significantly changed only in cana (100 µM)-treated cells are marked with a shaded arrow. The red line at the GLUT 1 transporter indicates transporter inhibition.
Ijms 23 07966 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zügner, E.; Yang, H.-C.; Kotzbeck, P.; Boulgaropoulos, B.; Sourij, H.; Hagvall, S.; Elmore, C.S.; Esterline, R.; Moosmang, S.; Oscarsson, J.; et al. Differential In Vitro Effects of SGLT2 Inhibitors on Mitochondrial Oxidative Phosphorylation, Glucose Uptake and Cell Metabolism. Int. J. Mol. Sci. 2022, 23, 7966. https://doi.org/10.3390/ijms23147966

AMA Style

Zügner E, Yang H-C, Kotzbeck P, Boulgaropoulos B, Sourij H, Hagvall S, Elmore CS, Esterline R, Moosmang S, Oscarsson J, et al. Differential In Vitro Effects of SGLT2 Inhibitors on Mitochondrial Oxidative Phosphorylation, Glucose Uptake and Cell Metabolism. International Journal of Molecular Sciences. 2022; 23(14):7966. https://doi.org/10.3390/ijms23147966

Chicago/Turabian Style

Zügner, Elmar, Hsiu-Chiung Yang, Petra Kotzbeck, Beate Boulgaropoulos, Harald Sourij, Sepideh Hagvall, Charles S. Elmore, Russell Esterline, Sven Moosmang, Jan Oscarsson, and et al. 2022. "Differential In Vitro Effects of SGLT2 Inhibitors on Mitochondrial Oxidative Phosphorylation, Glucose Uptake and Cell Metabolism" International Journal of Molecular Sciences 23, no. 14: 7966. https://doi.org/10.3390/ijms23147966

APA Style

Zügner, E., Yang, H. -C., Kotzbeck, P., Boulgaropoulos, B., Sourij, H., Hagvall, S., Elmore, C. S., Esterline, R., Moosmang, S., Oscarsson, J., Pieber, T. R., Peng, X. -R., & Magnes, C. (2022). Differential In Vitro Effects of SGLT2 Inhibitors on Mitochondrial Oxidative Phosphorylation, Glucose Uptake and Cell Metabolism. International Journal of Molecular Sciences, 23(14), 7966. https://doi.org/10.3390/ijms23147966

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop