# Proteomic Analysis of KCNK3 Loss of Expression Identified Dysregulated Pathways in Pulmonary Vascular Cells

Hélène Le Ribeuz, Florent Dumont, Guillaume>Ruellou, Mélanie Lambert, Thierry Balliau, Marceau Quatredeniers, Barbara Girerd, Sylvia Cohen-Kaminsky, Olaf Mercier, Stéphanie Yen-Nicolaï, et al.

► To cite this version:

Hélène Le Ribeuz, Florent Dumont, Guillaume>Ruellou, Mélanie Lambert, Thierry Balliau, et al.. Proteomic Analysis of KCNK3 Loss of Expression Identified Dysregulated Pathways in Pulmonary Vascular Cells. International Journal of Molecular Sciences, 2020, 21 (19), pp.1-21. 10.3390/ijms21197400. hal-03096431

HAL Id: hal-03096431

https://hal.science/hal-03096431

Submitted on 5 Jan 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

Proteomic Analysis of KCNK3 Loss of Expression
Identified Dysregulated Pathways in Pulmonary
Vascular Cells

Hélène Le Ribeuz <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Florent Dumont <fnref n="4" />, Guillaume Ruellou <fnref n="4" />, Mélanie Lambert <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,
Thierry Balliau <fnref n="5" />, Marceau Quatredeniers <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Barbara Girerd <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Sylvia Cohen-Kaminsky <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,
Olaf Mercier <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Stéphanie Yen-Nicolaÿ <fnref n="4" />, Marc Humbert <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, David Montani <fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID,
Véronique Capuano <fnref n="1" />,<fnref n="2" />,<fnref n="3" /> and Fabrice Antigny <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,*

Received: 11 September 2020; Accepted: 5 October 2020; Published: 7 October 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** The physiopathology of pulmonary arterial hypertension (PAH) is characterized by pulmonary artery smooth muscle cell (PASMC) and endothelial cell (PAEC) dysfunction, contributing to pulmonary arterial obstruction and PAH progression. KCNK3 loss of function mutations are responsible for the first channelopathy identified in PAH. Loss of KCNK3 function/expression is a hallmark of PAH. However, the molecular mechanisms involved in KCNK3 dysfunction are mostly unknown. To identify the pathological molecular mechanisms downstream of KCNK3 in human PASMCs (hPASMCs) and human PAECs (hPAECs), we used a Liquid Chromatography-Tandem Mass Spectrometry-based proteomic approach to identify the molecular pathways regulated by KCNK3. KCNK3 loss of expression was induced in control hPASMCs or hPAECs by specific siRNA targeting KCNK3. We found that the loss of KCNK3 expression in hPAECs and hPASMCs leads to 326 and 222 proteins differentially expressed, respectively. Among them, 53 proteins were common to hPAECs and hPASMCs. The specific proteome remodeling in hPAECs in absence of KCNK3 was mostly related to the activation of glycolysis, the superpathway of methionine degradation, and the mTOR signaling pathways, and to a reduction in EIF2 signaling pathways. In hPASMCs, we found an activation of the PI3K/AKT signaling pathways and a reduction in EIF2 signaling and the Purine Nucleotides De Novo Biosynthesis II and IL-8 signaling pathways. Common to hPAECs and hPASMCs, we found that the loss of KCNK3 expression leads to the activation of the NRF2-mediated oxidative stress response and a reduction in the interferon pathway. In the hPAECs and hPASMCs, we found an increased expression of HO-1 (heme oxygenase-1) and a decreased IFIT3 (interferon-induced proteins with tetratricopeptide repeats 3) (confirmed by Western blotting), allowing us to identify these axes to understand the consequences of KCNK3 dysfunction. Our experiments, based on the loss of KCNK3 expression by a specific siRNA strategy in control hPAECs and hPASMCs, allow us to identify

<footnote n="1">¹ Faculté de Médecine, Université Paris-Saclay, 94270 Le Kremlin-Bicêtre, France;
helene.leribeuz@aol.com (H.L.R.); melanie.lambert91@hotmail.fr (M.L.);
quatredeniers.marceau@gmail.com (M.Q.); barbara.girerd19@gmail.com (B.G.);
sylvia.cohen-kaminsky@universite.paris-saclay.fr (S.C.-K.); olaf.mercier@gmail.com (O.M.);
marc.humbert@bct.aphp.fr (M.H.); davidmontani@gmail.com (D.M.);
veronique.capuano@universite.paris-saclay.fr (V.C.)</footnote>

<footnote n="2">² INSERM UMR_S 999, Hypertension Pulmonaire: Physiopathologie et Innovation Thérapeutique, Hôpital Marie Lannelongue, 92350 Le Plessis-Robinson, France</footnote>

<footnote n="3">³ Assistance Publique - Hôpitaux de Paris (AP-HP), Service de Pneumologie et Soins Intensifs Respiratoires, Centre de Référence de l'Hypertension Pulmonaire, Hôpital Bicêtre, 94270 Le Kremlin-Bicêtre, France</footnote>

<footnote n="4">⁴ UMS Ingénierie et Plateformes au Service de l'Innovation Thérapeutique, Université Paris-Saclay, 92290 Châtenay-Malabry, France; florent.dumont@universite.paris-saclay.fr (F.D.); guillaume.ruellou@universite.paris-saclay.fr (G.R.); stephanie.nicolay@universite-paris-saclay.fr (S.Y.-N.)</footnote>

<footnote n="5">5 PAPPSO-GQE-Le Moulon, INRAE, CNRS, AgroParisTech, Université Paris-Saclay, 91190 Gif-sur-Yvette, France; thierry balliau@inrae.fr</footnote>

<footnote>* Correspondence: fabrice.antigny@universite.paris-saclay.fr; Tel.: +33-140-942-299; Fax: +33-140-942-522</footnote>
differences in the activation of several signaling pathways, indicating the key role played by KCNK3 dysfunction in the development of PAH. Altogether, these results allow us to better understand the consequences of KCNK3 dysfunction and suggest that KCNK3 loss of expression acts in favor of the proliferation and migration of hPASMCs and promotes the metabolic shift and apoptosis resistance of hPAECs.

**Keywords:** PAH; potassium channel; proteomic

# 1. Introduction

*Potassium Channel Subfamily K Member 3 (KCNK3)* was identified in 2013 as a predisposing gene for pulmonary arterial hypertension (PAH) [1]. To date, 12 different mutations in *KCNK3* have been described in PAH patients [2]. PAH is a severe and progressive disease, with an estimated prevalence of 15 to 50 per million of the population [3]. Under the new clinical classification [4], PAH consists of different etiologies leading to precapillary pulmonary hypertension, which is hemodynamically proposed to be defined by a mean pulmonary artery pressure (mPAP) >20 mmHg, a pulmonary artery wedge pressure <15 mmHg, and a pulmonary vascular resistance (PVR) >3 Wood units at rest [4]. *KCNK3* mutations carrier patients have a higher mPAP and are younger than idiopathic PAH (iPAH) patients at diagnosis [2]. All the identified *KCNK3* mutations result in a loss of function in KCNK3 channels [1,2]. KCNK3, also known as Twik-related-acid-sensitive-K⁺ channel (TASK-1) or K2P3.1, displays several electrophysiological characteristics of background K⁺ current, including minimal voltage sensitivity, extracellular pH sensitivity, resistance to classical K⁺ channel inhibitors, and insensitivity to cytoplasmic Ca²⁺ [5]. TASK-1 channels are expressed in pulmonary arterial smooth muscle cells (PASMCs) and endothelial cells (PAECs) [6].

In the absence of KCNK3 mutations, we found a reduced TASK-1 function or expression in lung or isolated pulmonary arteries (PA) from heritable PAH (*BMPR2* carriers patients) and iPAH patients, as well as in isolated PASMCs, demonstrating that reduced KCNK3 function/expression is a hallmark of PAH at the pulmonary vascular level [7]. Otherwise, we recently found that *Kcnk3* deficiency in rats (CRISPR/Cas9 technology) leads to higher susceptibility to monocrotaline- or chronic hypoxia-induced PH, confirming that KCNK3 loss of function predisposes one to the development of PAH. We also found that KCNK3 knockdown enhances proliferation in control human PASMCs (hPASMCs) [8].

In PAH, the PA obstruction is multifactorial, involving at least hPAECs and hPASMCs dysfunctions, leading to excessive PA vasoconstriction, an over-proliferation of PASMCs and PAECs associated with apoptosis resistance, and pro-inflammatory phenotypes. Most of these phenotypes were persistent in vitro in cultured hPASMCs and hPAECs from PAH patients [9].

In pulmonary vascular cells, K⁺ channels play major role in the control of cell contraction, cell proliferation and survival, cell volume, and cell apoptosis [2]. However, the role of KCNK3 dysfunction in the proliferation/apoptosis imbalance of hPASMCs and hPAECs is still unclear, as well as the molecular mechanism link to KCNK3 dysfunction and the aberrant phenotype of hPASMCs or hPAECs in PAH.

To identify the pathological molecular mechanisms induced by KCNK3 loss of expression in hPASMCs and hPAECs, we postulated that a differential proteomic analysis between siControl and siKCNK3 conditions would identify the crucial signaling pathways to better understand the contribution of KCNK3 dysfunction to PAH pathobiology. We used a Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS)-based proteomic approach coupled with bioinformatics analysis to understand the consequences of KCNK3 dysfunction in hPASMCs and hPAECs.
## 2. Material and Methods

### 2.1. Human PASMCs and PAECs Cell Culture

Human lung specimens were obtained at the time of a lobectomy or pneumonectomy for localized lung cancer from control subjects. Pulmonary arteries were excised at a distance from tumor areas. hPAECs and hPASMCs were cultured as previously described [10,11], and were used for the study at passage 4. The patients studied were part of a program approved by the institutional Ethics Committee, and had given their written informed consent (ID RCB: 2018-A01252-53, approved on June 18, 2006). Five patients for the condition of hPAECs and 5 patients for hPASMCs were compared in this experiment, and each of these cohorts was duplicated into siControl vs. siKCNK3 conditions, giving a total of 20 samples. The results and statistics were taken from an average of 5 patients per condition.

### 2.2. Reverse Transcription Quantitative PCR (RT-qPCR)

Total RNA was extracted using TRIzol reagent according to standard procedures. One microgram of total RNA was reverse-transcribed using a QuantiTect Reverse Transcription Kit (Qiagen, Valencia, CA, USA; cat. no. 205311). Gene expression was quantified using qPCR following the standard protocol for ready-to-use TaqMan gene expression assays (Thermo Scientific, Villebon-sur-Yvette, France) on a StepOne Plus Real-Time PCR System (Life Technologies, Carlsbad, CA, USA). The predesigned probe sets used for experiments were human KCNK3 (Hs00605529_m1 Thermo Scientific, Paris, France) and human 18S (Hs99999901_s1 Thermo Scientific, Villebon-sur-Yvette, France).

### 2.3. Western Blot Analyses

Total proteins from hPASMCs and hPAECs or rat lung-tissue samples were prepared as described previously [12]. Anti-HO-1 (Abcam; ab13248) was used at 1/1000, anti-IFIT3 (Clinisciences, Nanterre, France; sc-393396) was used at 1/1000, and anti-β-Actin was used for normalization.

### 2.4. Protein Preparation

A total of 20 µg of proteins from each sample was loaded per land on a 12% SDS-PAGE. A short migration was then performed, allowing the proteins to penetrate the gel, and the gels were stained with Coomassie colloidal blue (EZblue, Sigma-Aldrich, Lyon, France) in order to recover a 2 mm band containing all the proteins.

### 2.5. Protein In-Gel Digestion

A reduction/alkylation step was performed with dithiothreitol (DTT) at 10 mM and iodoacetamide at 55 mM. Pieces of the gels were rehydrated at 4 °C in 12 ng/µL of sequencing-grade modified trypsin (Promega, Charbonnières-Les-Bains, France) solubilized in 25 mM of NH₄HCO₃ for 1 h and then digested at 37 °C overnight. After tryptic digestion, the peptides were extracted by incubating gel pieces in extraction solvent (0.1% AF/50% ACN).

### 2.6. Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS) Analysis

As described in a previous study [13], samples were analyzed on a nano-LC Ultra 2D (eksigent) on-line with a high-resolution mass spectrometer Q-Exactive (Thermo Scientific, France) using a nanoelectrospray interface in positive polarity mode. For each sample, 4 µL of protein digest was loaded onto a Biosphere C18 Precolumn (0.1 × 20 mm, 100 Å, 5 µm; Nanoseparation) at 7.5 µL/min and desalted with 0.1% (*v*/*v*) formic acid and 2% (*v*/*v*) Acetonitrile. After 3 min, the precolumn was connected to a Biosphere C18 Nanocolumn (0.075 × 300 mm, 100 Å, 3 µm; Nanoseparation).

Electrospray ionization was performed at 1.6 kV with an uncoated capillary probe. The buffers were 0.1% formic acid in water (A) and 100% ACN (B). Peptides were separated using a linear gradient
# 3. Results

## 3.1. Validation of KCNK3 Knockdown in hPASMCs and hPAECs

To decipher the changes in the proteome in response to the knockdown of KCNK3, an LC-MS/MS-based proteomics approach was carried in hPASMCs and hPAECs. First, to evaluate the molecular consequences of the KCNK3 loss of expression in hPASMCs and hPAECs, we reduced the KCNK3 expression using siRNA against KCNK3. As presented in Figure 1A, 72 h after cell transfection, the KCNK3 mRNA was decreased by 85% in hPASMCs and by 75% in hPAECs (Figure 1A). Then, we used proteomic approach to identify and measure the relative abundance of proteins differentially expressed in hPASMCs or hPAECs transfected with siControl or with siKCNK3 (Figure 1B).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Experimental procedure for the proteomic approach. (A) siRNA validation by RT-qPCR in hPASMCs (n = 5) and hPAECs (n = 5). \* *p* < 0.05 vs. control (B) Experimental process for the LC-MS/MS approach. (C) Number of proteins differentially expressed in siControl vs. siKCNK3 in hPASMCs and hPAECs. (D) Venn diagram representing the number of proteins identified in hPASMCs only, hPAECs only, and both hPASMCs and hPAECs.</figcaption></figure>

Consistent with the purity of our hPAEC cultures, endothelial specific protein markers were found to be expressed, including Von Willebrand factor, vascular endothelial growth factor receptor 2, platelet endothelial cell adhesion molecule 1, endothelin-converting enzyme 1, angiopoietin-2, and endothelial nitric oxide synthase (not shown).

The purity of hPASMC cultures was confirmed by the absence of specific PAEC markers and the presence of specific smooth muscle markers, including alpha-actinin-1, alpha-actinin-4, myoferlin, actin, alpha cardiac muscle 1, gelsolin, transgelin-2, transgelin, caldesmon, cofilin 2 (muscle) (not shown).

## 3.2. Comparison of Differentially Expressed Proteins Between hPAECs and hPASMCs Treated With siControl or siKCNK3

We identified a total of 2676 proteins differentially expressed in hPAECs transfected with siControl vs. siKCNK3 (Supplementary Table S1), and a total of 1735 proteins differentially expressed in hPASMCs transfected with siControl vs. siKCNK3 (Supplementary Table S2) (Figure 1C). As shown in Figure 1D, 1472 proteins were common in both hPAECs and hPASMCs, 1204 proteins were specific
to hPAECs, and 263 were specific to hPASMCs (Figure 1D). Principal component analysis (PCA) based on the relative protein expression levels showed the differences between each patient in terms of hPASMCs and hPAECs (Figure 2A,B and Figure 2C,D, respectively), highlighting the proteome diversity in each patient, represented in axis 1. PCA also revealed the presence of two protein clusters in the hPASMCs and hPAECs experiments by separating the siControl and siKCNK3 samples under axis 4 and 5. These results demonstrated that KCNK3 knockdown induces proteome modification in hPASMCs and hPAECs.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Principal component analysis (PCA) performed on the relative protein expression level. (A) PCA showing the consequence of KCNK3 knockdown on the hPAEC proteome obtained by spectral counting (SC) (siControl (black) or siKCNK3 (red)). (B) PCA showing the consequence of KCNK3 knockdown on the hPAECs proteome obtained by eXtracted Ion Current (XIC) analysis (siControl (black) or siKCNK3 (red)). (C) PCA showing the consequence of KCNK3 knockdown on the hPASMC proteome obtained by SC analysis (siControl (black) or siKCNK3 (red)). (D) PCA showing the consequence of KCNK3 knockdown on the hPASMCs proteome obtained by XIC analysis and SC analysis (siControl (black) or siKCNK3 (red)). All are represented with axes 1–2 and axes 4–5.</figcaption></figure>

The distribution of all the differentially expressed proteins in the KCNK3-knockdown hPAECs (Figure 3A) and hPASMCs (Figure 3D) were illustrated in a volcano plot. The 10 proteins with the highest fold change in KCNK3-silenced hPAECs were interferon-induced proteins with tetratricopeptide repeats 3 (IFIT3), interferon-induced protein with tetriceopeptide repeats 1(IFIT1), phosphoglycolate phosphate (PGP), actine-like protein 6A (ACL6A), gap junction protein alpha-1 (CXA1), SAM and HD domain containing deoxynucleoside triphosphate triphosphohydrolase 1 (SAMH1), transmembrane emp24 protein transport domain containing 4 (TMED4), MX dynamin-like GTPase 1 (MX1), 2'-5'-oligoadenylate synthase 2(OAS2), and 2'-5'-oligoadenylate synthase 3(OAS3) (Supplementary Table S1). In KCNK3-silenced hPASMCs, the 10 proteins with the highest fold change were G-protein subunit beta 4 (GBB4), DExD-Box helicase 39B (DX39B), proteasome 20S subunit beta 8 (PSB8), SEC24 Homolog A COPII coat complex component (SC24A), heme oxygenase 1 (HMOX1), LDL receptor-related protein 1 (LRP1), heterogeneous nuclear ribonucleoprotein A/B (ROAA), prosaposin (SAP), cyclase-associated actin cytoskeleton regulatory protein 2 (CAP2), microtubule-associated protein 1B (MAP1B) (Supplementary Table S2).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Volcano plot of dysregulated proteins in KCNK3-knockdown hPAECs and hPASMCs, and confirmation by Western blot. (A) Volcano plot representation of dysregulated proteins in hPAECs transfected with siKCNK3 according to their fold change (outside the yellow line interval fold change > 1.5 or < -1.5) and *p*-value (above the yellow line *p*-value < 0.05). (B) Comparison of the Heme Oxygenase-1 (HMOX-1 or HO-1) abundance by LC-MS/MS and Western blot analysis from the same hPAECs lysates. (C) Comparison of the EIF2AK2 abundance by LC-MS/MS and Western blot analysis from the same hPAECs lysates. (D) Comparison of Interferon-induced protein with tetratricopeptide repeats 3 (IFIT3) abundance by LC-MS/MS and Western blot analysis from the same hPAECs lysates. (E) Volcano plot representation of dysregulated proteins in siKCNK3 hPASMCs according to their fold change (outside the yellow line interval fold change > 1.5 or < -1.5) and *p*-value (above the yellow line *p*-value < 0.05). (F) Comparison of the HMOX-1 (HO-1) abundance by LC-MS/MS and Western blot analysis from the same hPASMCs lysates. (G) Comparison of the Eukaryotic Translation Initiation Factor 2 Alpha Kinase 2 (EIF2AK2 or PKR) abundance by LC-MS/MS and Western blot analysis from the same hPASMCs. (H) Comparison of the IFIT3 abundance by LC-MS/MS and Western blot analysis from the same hPASMCs lysates. \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* \*p* < 0.0001 vs. control. Experiments were analyzed using the Mann-Whitney test.</figcaption></figure>
Among these proteins, we selected HMOX1 (or HO-1) and IFIT3 for confirmation by the semi-quantitative Western blot approach to confirm these findings due to the importance of HMOX1 and EIF2AK2 (or PKR) in the pulmonary field [12] and due to the novelty of IFIT3 in the pulmonary hypertension field.

Among these proteins, the HMOX1/HO-1 abundances' behaviors were investigated by Western blot. Consistent with the proteomic analysis (Figure 3B,C and Figure 3E,F), the Western blot analyses showed that the HMOX-1 protein expression was significantly increased in the hPAECs (Figure 3B) and hPASMCs (Figure 3F) transfected with siKCNK3 compared to the siControl conditions, in the same way as the proteomic analysis (Figure 3A,B,E,F, Supplementary Tables S2 and S4). We also analyzed the EIF2AK2 protein expression by Western blot in hPAECs and hPASMCs transfected with siKCNK3, and we found a significant decrease in the EIF2AK2 expression in KCNK3-knockdown hPAECs (Figure 3C) and hPASMCs (Figure 3G) similar to the proteomic results (Figure 3A,C,E,G). Moreover, we analyzed by Western blot the IFIT3 protein expression, and we found a significant decrease in the IFIT3 expression in KCNK3-knockdown hPAECs (Figure 3D) and in hPASMCs (Figure 3H), similar to the proteomic results (Figure 3A,D,E,H). Altogether, these results confirmed our results from the proteomic experiments.

### 3.3. Global Consequences of KCNK3 Knockdown in hPAECs and hPASMCs

As highlighted by the heatmap clustering of differentially expressed proteins of each group of conditions (Figure 4A–D, we found that 157 proteins were significantly upregulated in siKCNK3, while 247 proteins were significantly downregulated in siKCNK3-transfected hPAECs (Figure 4A–D and Supplementary Table S3). A total of 111 proteins were significantly upregulated in siKCNK3, while 154 proteins were significantly downregulated in siKCNK3-transfected hPASMCs (Figure 4C,D and Supplementary Table S3). To identify the canonical pathways, predicted upstream regulators, biological functions, and disease, the significant proteins differentially expressed were analyzed using Ingenuity Pathway Analysis (IPA, Ingenuity Systems, www.ingenuity.com). The IPA analysis revealed differences in 30 canonical pathways (z-score > 1 or −1), 48 upstream regulators, 15 biological functions, and 20 disease and functional networks between KCNK3-knockdown hPAECs. In the hPASMCs, the IPA analysis found significant differences in 14 canonical pathways (z score > 1 or −1), 16 upstream regulators, 9 biological functions, and 13 disease and functional networks. The activation Z-score makes predictions about potential regulators by using information about the direction of gene regulation [16].
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Heatmap of the proteomic results and pathways deregulated in siKCNK3 hPASMCs and hPAECs. (A) Heatmap representation for proteomic results in hPAECs. hPAECs from the same patient between the siControl and siKCNK3 conditions are clustering together. (B) Representation of the 37 major pathways dysregulated in siKCNK3 hPAECs found with IPA (Ingenuity Pathway Analysis) software. (C) Heatmap representation of proteomic results in hPASMCs. hPASMCs from the same patient between the siControl and siKCNK3 conditions are clustering together. (D) Representation of the 37 major pathways dysregulated in hPASMCs transfected with siKCNK3 found with the Ingenuity Pathway Analysis (IPA) software.</figcaption></figure>

### 3.4. Biological Consequences of KCNK3 Knockdown in hPAECs

To identify the most interesting signaling pathways, networks, diseases, and biological functions deregulated by the loss of KCNK3 expression, we separated the hPAECs specific proteome (316 proteins), the hPASMCs specific proteome (222 proteins), and the common proteome of both hPAECs and hPASMCs (53 proteins) (Supplementary Table S3).

For the hPAECs-specific proteome, we selected the five most affected canonical pathways by the combination of the highest -log(*p*-value) in association with the highest z-score. We used the IPA regulation z-score algorithm to identify the canonical pathways, upstream (positive z-score) or downstream (negative z-score) regulators, and biological functions according to our proteomic data. (positive z-score) according to our RNAseq data. In order to enhance the stringency of our analysis, we considered only functions with a z-score < -1.5 or > 1.5. Among them, our results indicate that Eukaryotic Initiation Factor 2 (EIF2) signaling (z score -1.89) was reduced in the KCNK3-knockdown hPAECs, and the mammalian target of rapamycin (mTOR) signaling (z-score 2.23) and superpathway of methionine degradation (z-score 2.82) were activated in the KCNK3-knockdown hPAECs (Supplementary Table S4a). EIF2 signaling integrates diverse stress-related pathways to regulate both the global and specific mRNA translation. mTOR signaling is known to be a master of growth regulator that senses and integrates diverse nutritional and environmental alterations, while
the superpathway of methionine degradation, as indicated by the name, describes the degradation of the amino acid L-methionine in mammals.

In association with deregulated pathways, IPA identified abnormal biological functions such as the splicing of RNA, the processing of mRNA, neuronal cell death, cell survival, the proliferation of EC lines, neoplasia, and the catabolism of protein was found to be inhibited in **KCNK3-knockdown hPAECs** (Figure 5A,B and Supplementary Table S4b), while cell death, the metabolism of membrane lipid derivative, and the cell proliferation of breast cell line were found to be activated in **KCNK3 knockdown-hPAECs** (Figure 5A,B and Supplementary Table S4b).

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Deregulated biological functions in **KCNK3-knockdown hPAECs** and the protein network of interactions. (A) Heatmap representing the up-regulated (orange) or down-regulated (blue) biological functions in hPAECs transfected with siKCNK3. (B) Histogram representing the 10 first deregulated biological functions in hPAECs transfected with siKCNK3, classified according to –log (p-value) (B). (C) Representation of the dysregulated network 1: protein synthesis, molecular transport, hereditary disorder. Red proteins are up-regulated proteins and green proteins are down-regulated proteins.</figcaption></figure>

Moreover, the IPA identified in **KCNK3-knockdown hPAECs** six upstream molecules that should be activated, including lysine demethylase 5 (KDM5), histone desacetylase 2 (HDAC2), transcription factor AP-2 alpha (TFAP2A), and NK2 Homebox 3 (NKX2-3). IPA also identified seven upstream molecules should be inhibited, including interferon alpha, interferon alpha 2 (IFNA2), serine/threonine kinase 11 (STK11), and interferon gamma (IFNG) (Supplementary Table S4c).

The differentially expressed proteins in **KCNK3-knockdown hPAECs** were associated with 19 functional networks, including protein synthesis, molecular transport (network 1) and RNA post-transcriptional modification, and cancer (network 2) (Figure 5C and Supplementary Table S4d). In this network, we found, for example, that ferritin, ferredoxin reductase (FDXR), and ferritin heavy chain 1 (FTH1) were downregulated, while TNF Receptor Superfamily Member 10d (TNFRSF10D) was upregulated (Figure 5C). These results indicate that the cellular iron storage was reduced in the absence of KCNK3, as well as TNFRSF10D signaling, which is related to p53 and ERK signaling.
### 3.5. Biological Consequences of KCNK3 Knockdown in hPASMCs

In hPASMCs, our results indicated that four canonical pathways were reduced in KCNK3-knockdown hPASMCs: EIF2 signaling (z-score = 0.2), Purine Nucleotides De Novo Biosynthesis II (z-Score = -2), and Interleukine-8 (IL-8) signaling (z-score = -2.23), while one canonical pathway was activated: phosphoinositide 3-kinase/serine/threonine kinase Akt (PI3K/AKT) signaling (z-score = 2.45) (Supplementary Table S5a). These results suggest an activation of the PI3K/AKT pathway and an inhibition of EIF2, Purine Nucleotides De Novo Biosynthesis II, and IL-8 pathways in KCNK3-knockdown hPASMCs. Purine Nucleotides De Novo Biosynthesis II is highly conserved among organisms, participating in many aspects of cellular metabolism, including the structure of DNA and RNA, serving as an enzyme cofactor, functioning in cellular signaling, acting as a phosphate group donor, and generating cellular energy. IL-8 is a pro-inflammatory CXC chemokine which is described to regulate gene expression via numerous transcription factors, modulate the cellular proteome by acting at the translational level, and modulate the cytoskeleton cell organization by acting at the posttranslational level. The PI3K/AKT pathway acts as an important regulator of the cell cycle that is directly related to cellular quiescence, proliferation, and cancer, and is described to be over-activated in PAH.

Underlying these alterations, the IPA predicted that four upstream molecules should be activated among them: ten-eleven-translocation 2 (TET2), lysine demethylase 5A (KDM5A), and homeobox A10 (HOXA10). The IPA predicted that seven upstream molecules should be inhibited in KCNK3-knockdown hPASMCs, among them MLX Interacting Protein Like (MLXIPL) and Heat Shock Protein 90 Beta Family Member 1 (HSP90B1) (Supplementary Table S5b). In association with these deregulated pathways, the IPA identified abnormal biological functions, such as the cell spreading and shape change of tumor cell lines, and cell movement, cell attachment, and cell death were found to be significantly activated in KCNK3-knockdown hPASMCs (Figure 6A,B and Supplementary Table S5c).

The differentially expressed proteins in KCNK3-knockdown hPASMCs were associated with 12 functional networks, including RNA post-transcriptional modification, cancer, hematological disease (network 1), RNA damage and repair, protein synthesis, gene expression (network 2) (Figure 6C and Supplementary Table S5d). In this network, we found, for example, that Guanine Nucleotide-Binding protein subunit beta-4 (GNB4) was downregulated, while Centrosomal Protein 170 (CEP170) and (AHNAK) were upregulated (Figure 5C). The reduced expression of GNB4 could have as a consequence a reduced GTPase activity in KCNK3-knockdown hPASMCs, while the over-expression of CEP170 and AHNAK should contribute to the enhancement of microtubule organization, cell structure and migration, and calcium channel regulation in KCNK3-knockdown hPASMCS.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Deregulated biological functions in KCNK3-knockdown hPASMCs and the protein network of interactions. (A) Heatmap representing the up-regulated (orange) or down-regulated (blue) biological functions in hPASMCs transfected with siKCNK3. (B) Histogram representing the 10 first deregulated biological functions in hPASMCs transfected with siKCNK3, classified according to -log (*p*-value). (C) Representation of the dysregulated network 1: RNA post-transcriptional modifications, cancer, and hematological disease. Red proteins are up-regulated proteins and green proteins are down-regulated proteins.</figcaption></figure>

### 3.6. Common Deregulated Pathways to KCNK3-Knockdown hPAECs and hPASMCs

We found that 53 proteins were common to hPAECs and hPASMCs transfected with siKCNK3 (Supplementary Table S6a). Using IPA, we identified two canonical pathways with a z-score value, NRF2 (nuclear factor erythroid 2-related factor 2) and interferon signaling, with a z-score prediction of 1.34 and -2, respectively, in hPASMCs, and with a z-score predicted value of 0.44 and -2, respectively, in hPAECs (Figures 7A and 8A). These results suggest an activation of NRF2 pathway and an inhibition of interferon signaling in KCNK3-knockdown hPAECs and hPASMCs. NRF2 is a transcription factor that regulates several genes (>600) involved in glutathione synthesis, antioxidant proteins, enzyme activity, transporter activity, metabolism, and transcription. The interferon signaling plays a critical role in the immune response, conferring a cellular resistance to viral infection.

Based on our enrichment data, the IPA identified that 10 molecules could explain the decrease in interferon signaling. These inhibited upstream molecules were interferon lambda 1 (IFNL1), interferon regulatory factor 1 (IRF1), interferon beta (IFN Beta), interferon alpha 2 (IFNA2), interferon alpha/beta receptor (IFNAR), interferon regulatory factor 7 (IRF7), interferon gamma IFNG, interferon regulatory factor 3 (IRF3), interferon beta 1(IFNB1), and interferon alpha (Supplementary Table S6b). Moreover, the IPA predicted that only four molecules were activated: bruton tyrosine kinase (BTK), atypical chemokine receptor 1 (ACKR2), mitogen-activated protein kinase 1 (MAPK1), and ETS variant transcription factor 6-RUNX family transcription factor 1 (ETV6-RUNX1) (Supplementary Table S6b). IPA identified that these dysregulated pathways were related to abnormal biological functions, such as the cell death of myeloid cells, and cell cycle progression was found to be significantly inhibited, while DNA damage was found to be significantly activated in both hPAECs and hPASMCs lacking KCNK3 (Figures 7B and 8B and Supplementary Table S6c).
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Common deregulated functions in hPASMCs and hPAECs transfected with siKCNK3. (A) Histogram representing the 2 deregulated canonical pathways in hPAECs and hPASMCs transfected with siKCNK3, using the fold change obtained in hPAECs transfected cells. (B) Heatmap representing the up-regulated (orange) or down-regulated (blue) biological functions in hPAECs and hPASMCs transfected with siKCNK3, using the fold change obtained in hPAECs transfected cells. (C) Representation of dysregulated network 1: immunological disease and organismal injury and abnormalities. Red proteins are up-regulated proteins and green proteins are down-regulated proteins.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Common deregulated functions in hPASMCs and hPAECs transfected with siKCNK3. (A) Histogram representing the 2 deregulated canonical pathways in the hPAECs and hPASMCs transfected with siKCNK3, using the fold change obtained in hPASMCs transfected cells. (B) Heatmap representing the up-regulated (orange) or down-regulated (blue) biological functions in hPAECs and hPASMCs transfected with siKCNK3, using the fold change obtained in hPASMCs transfected cells. (C) Representation of dysregulated network 1: immunological disease and organismal injury and abnormalities. Red proteins are up-regulated proteins and green proteins are down-regulated proteins.</figcaption></figure>
Finally, the common differentially expressed proteins in KCNK3-knockdown hPAECs and hPASMCs were associated with four functional networks, including the networks immunological disease, organismal injury, and abnormalities (network 1) (Figures 7C and 8C, and Supplementary Table S6d). In this network and in accordance with the previous upstream identified molecules, we found, for example, that interferon induced protein 35 (IFI35), IFIT1, IFIT3, and interferon-α induced proteins were downregulated in hPAECs and hPASMCs transfected with siKCNK3, as well as Eukaryotic Translation Initiation Factor 2 Alpha Kinase 2 (EIF2AK2) (Figures 7C and 8C). Common to hPAECs and hPAECs, vesicle trafficking 1 (VTA1) was overexpressed in siKCNK3 conditions (Figures 7C and 8C). The reduced expression of interferon-signaling proteins suggest that, in the absence of KCNK3, hPASMCs and hPAECs should have a drop in their viral infection response. Whereas, the overexpression of VTA1 should contribute to the intensification of lysosomal protein degradation in KCNK3-knockdown hPAECs and hPASMCS. These common deregulated pathways in hPAECs and hPASMCs have highlighted the fact that KCNK3 has an unexpected similar function in hPAECs and hPASMCs, despite the different physiological cell functions, suggesting that these commonly identified pathways are directly related to KCNK3 and independent of cell systems.

## 4. Discussion

KCNK3 reduced function and expression is a hallmark of PAH. KCNK3 dysfunction is described to lead to the plasma membrane depolarization of PASMCs, pulmonary artery vasoconstriction, and increased PASMC proliferation. However, the molecular mechanisms linked to KCNK3 dysfunction in hPASMCs and hPAECs are mostly unknown.

In the present study, beyond the deregulated pattern expected in KCNK3 knockdown in hPASMCs and hPAECs, the proteomic approach has highlighted new relevant deregulated pathways. Indeed, the KCNK3 loss of expression induces the inhibition of interferon signaling, EIF2 signaling (in hPAECs and hPASMCs), and purine nucleotides de novo biosynthesis II and IL-8 signaling (in hPASMCs). In addition, KCNK3 loss of expression induces the activation of the NRF2 pathway (in hPAECs and hPASMCs), glycolysis, apoptosis signaling, the superpathway of methionine degradation and mTOR signaling (in hPAECs), and PI3K/AKT signaling (in hPASMCs). These results help us to better understand the pathophysiological consequences of KCNK3 dysfunction in pulmonary vascular cells and their implications in PAH pathogenesis. In PAH, PAECs dysfunction is characterized by the impairment of endothelial-dependent vasodilatation, by a reduction in anticoagulant properties, by metabolic changes, by the production of reactive oxygen species, and by an unadapted release of pro-inflammatory molecules. These dysfunctions result in PAECs phenotypic changes, including the decreased capacity for vascular tube formation, increased aerobic glycolysis, the loss of several endothelial markers, the acquisition of several mesenchymal markers and a pro-inflammatory phenotype, and increased *in vitro* hPAECs proliferation [9]. Consequently to PAEC dysfunction, hPASMCs are deregulated and characterized by enhanced proliferation, survival, vasoconstriction, and cell movement in PAH [9].

We recently found that KCNK3 knockdown is associated with ERK-1/2 activation, HIF1α stabilization, and the over-phosphorylation of AKT, promoting the proliferation of control hPASMCs [8]. Here, with our analysis we found that *Kcnk3* deficiency is in favor of cell spreading, cell movement, and cell invasion, as well as some canonical pathways linked to cell proliferation (mTOR, PI3K/AKT) and DNA damage. These functions have been well described to promote cancer and tumor development [17]. Based on numerous points of similarity, PAH has emerged as a cancer-like disease, and is characterized by abnormal cell proliferation, aberrant resistance to apoptosis, and the deregulation of cellular metabolism [17], providing new concepts and new therapeutic perspectives for PAH. A large spectrum of growth factors, cytokines, and chemokines has been demonstrated to promote or to initiate PAH and cancer, including platelet-derived growth factor (PDGF) and interleukin 6 (IL-6) [18]. We found that the knockdown of KCNK3 promotes the proliferation of control hPASMCs [8]
In the absence of KCNK3, the control hPAECs proliferation rate was unchanged compared to the siControl conditions (not shown); however, the consequence of the KCNK3 loss of expression in the phenotype of hPAECs is unknown.

## 4.1. HO-1 Relevance in PAH

Common to hPASMCs and hPAECs, we found that *KCNK3* knockdown deregulated the canonical pathway, EIF2 (eukaryotic initiation factor 2) signaling, and eIF4 and p70S6K signaling. EIF2 is crucial for protein synthesis and the initiator binding of tRNA to the ribosome. EIF2 signaling is known to play a critical role in the vascular remodeling and proliferation of PASMCs in chronic hypoxia-induced pulmonary hypertension [19], as well as in pulmonary vein occlusive disease [12]. EIF2 signaling is a part of the integrated stress response (ISR). To restore cellular homeostasis, eukaryotic cells activate a common adaptive pathway in response to stress stimuli. ISR is mediated by several ISR kinases (GCN2, PERK, HRI), and the protein kinase R (PKR)—also named EIF2AK2—is known to phosphorylate the α-subunit of eIF2α, leading to the activation of activating transcription factor 4 (ATF4). ATF4 is known to induce the expression of other transcription factors, including HO-1 [12]. Here, we found that EIF2AK2 is reduced in hPASMCs and hPAECs transfected with siKCNK3. Additionally, in hPASMCs we found that the co-effector of GCN2 (general control nonderepressible 2, another ISR kinase), GCN1, is upregulated in the absence of KCNK3 (Supplementary Table S2). GCN2/GCN1 senses amino acid deficiency and modulates amino acid metabolism as a response to nutrient deprivation [20]. Heme oxygenase-1 (HO-1) is a stress-inducible enzyme HO-1 which maintains cellular homeostasis under stress conditions [21]. HO-1 is well established to play a crucial role in the cell adaptive response to oxidative and cellular stresses, increasing cell proliferation and cell angiogenesis and reducing inflammation and apoptosis [22]. Here, we found that HO-1 was overexpressed in hPASMCs and hPAECs transfected with siKCNK3. HO-1 plays a crucial role in the survival of vascular smooth muscle [23], and could play major role in many human diseases, including pulmonary vein-occlusive disease [12,22]. HO-1 promotes cancer development, since HO-1 overexpression is observed in different types of cancers promoting cell survival [24], while the inhibition of HO-1 reduces cell viability [25]. HO-1 is also described to regulate endothelial-to-mesenchymal transition [26], which is described to contribute to PH pathobiology [27] in relation to *Kcnk3* deficiency in rats [8]. In vascular smooth muscle cells, Liu et al. found that ER stress mediates HO-1 overexpression via the antioxidant response element (ARE), promoting the suppression of apoptosis [28]. NRF2 is the major transcription factor that promotes the ARE pathway [29], and here we found that the NRF2-mediated oxidative stress response is deregulated in *KCNK3*-knockdown hPAECs and hPASMCs, indicating that the loss of KCNK3 is a key signal to activate this pathway and promote cell survival. We also found that *KCNK3* knockdown leads to a reduced EIF2AK2 expression and the overexpression of HO-1, which should lead to increased biliverdin and bilirubin. An elevated serum bilirubin content is found to be a risk factor for death in PAH patients [30]. Moreover, HO-1 activity cleaves the heme to form equimolar ferrous ions, biliverdin-Ixα, and carbon monoxide (CO) [31–33]. The toxic properties of CO are well known in pulmonary medicine; CO replaces O2 from hemoglobin, causing tissue hypoxia [34]. At high CO concentrations, there are other mechanisms of CO-induced toxicity, including apoptosis and lipid peroxidation [35]; however, at low concentrations, CO contributes to several physiological reactions. Indeed, CO activates soluble guanylate cyclase (sGC), inducing the production of cyclic 3':5'-guanosine monophosphate (cGMP) [36].

The sGC/cGMP pathway mediates the effects of CO on smooth muscle cell relaxation [36,37]. CO-induced vascular relaxation could be directly mediated by the activation of calcium-dependent potassium channels (BKCa) [38,39], which could be a compensatory mechanism in the context of KCNK3 dysfunction.
## 4.2. Interferon Signaling Relevance in PAH

The interferon-induced protein with tetratricopeptide repeats (IFIT) gene family has been well studied for their antiviral properties [40]. In most cell types without viral infection, IFIT gene expression is involved in several biological processes, such as virus-induced translation initiation, replication, and double-stranded RNA signaling [41]. The expression of IFIT genes is enhanced by interferon treatment and viral infection [42]. In this study, we found by the LS/MS approach that the IFIT protein family (IFIT1, IFIT3, IFIT5, IFITX) is reduced in hPASMCs and hPAECs transfected with siKCNK3, and by a larger extent in hPAECs. Additionally, we confirmed the reduction in IFIT3 expression by Western blot experiments, indicating that KCNK3 dysfunction leads to IFIT family down-expression. IFIT genes are involved in the response to viral infection as well as KCNK3. In fact, the overexpression of KCNK3 reduces the ability of the HIV-1 accessory protein (Vpu) to release viral particles, while Vpu is described to abolish the TASK-1 current [43]. HIV-1 exploits the host miRNA cellular systems reduce to block the mechanism of innate inhibition. HIV-1 regulates miR-21, which is known to downregulate the KCNK3 expression, allowing a stronger infection efficiency [44]. HIV infection confers an increased risk for a variety of infectious lung diseases, and HIV infection is an established risk factor for PH. The prevalence of PAH is several-fold higher in HIV-infected patients compared with the PH population [45]. Our results suggest that the loss of function of KCNK3 could facilitate PAH in HIV-infected patients.

Due to their antiviral, immunomodulatory, and antitumor properties, interferon type I (IFN1) is naturally overexpressed in the human body in response to pathogens and in tumor cells. We found in a PAH lung increased interferon receptor 1 (IFNAR1) protein levels. We also found that *Ifnar1*<sup>-/-</sup> mice were protected against the development of PH [46], indicating that IFN1 contributes to PAH pathogenesis. The benefits of IFN therapy in combination with chemotherapy and radiotherapy is well documented for several types of cancers. However, it was also reported that IFN therapy has several side effects, triggering PAH [46,47].

## 4.3. KCNK3 Dysfunction and Superpathway of Methionine Degradation

Here, we found an increase in glycolysis and a decrease in the superpathway of methionine degradation in KCNK3-knockdown hPAECs, indicating that the loss of expression of KNCK3 facilitates the glycolysis metabolic shift of hPAECs. In most cancers, enzymes responsible for the Methionine (Met) cycle are altered [48]. Met is also crucial for RNA and protein synthesis, and polyamine [49]. Polyamines increase cell proliferation by reducing cell apoptosis [50].

The effect of Met restriction leads to the inhibition of IGF signaling and mTOR, and a reduction in mitochondrial oxidative stress [51]. In the present study, we found that the superpathway of Met degradation is reduced in hPAECs transfected with siKCNK3, suggesting that the loss of KCNK3 increases the pool of Met by reducing their degradation pathway, which could lead to the activation of the mTOR signaling pathway. Indeed, we also found an activation of the mTOR signaling pathway in KCNK3-knockdown hPAECs. The overactivation of the mTOR pathway is well known to contribute to PAH pathogenesis, and the in vivo inhibition of mTOR by rapamycin reduces the development experimental PH [52]. Moreover, mTOR is a part of the PI3K/AKT/mTOR pathway, which we found to be increased in lungs from *Kcnk3*-mutated rats as well as in human PAH [8,10].

Patients with hypermethioninemia exhibit a lot of symptoms, including neurological problems, sluggishness, muscle weakness, liver problems, and unusual breath [53]. A high level of Met could hypermethylate DNA and then produce undesirable epigenetic effects. DNA methylation contributes to the pathogenesis of pulmonary diseases, including cancer and PAH [54], in particular by promoting phenotypic changes in hPAECs in PAH [55]. Alterations in the DNA methylation of the *SOD2* gene (Superoxide dismutase 2) are known to reduce the expression of SOD2, which is related to pulmonary vascular cell proliferation, apoptosis resistance, and cell metabolic changes [56]. Here, common to hPAECs and hPASMCs, we found that KCNK3 knockdown reduces the expression of SOD2, supporting the role of KCNK3 in promoting PAH via aberrant DNA methylation.
Moreover, in hPAECs (but not in hPASMCs) transfected with siKCNK3, we found a decrease in the expression of gamma-glutamyl hydrolase (GGH). GGH is involved in the turnover of cellular folates, since a reduction in the GGH expression increases the cellular folate content and poly-gamma glutamate. In lung adenocarcinoma, the GGH expression was lower and correlated with higher folate concentrations and increased DNA methylation [57]. In the context of KCNK3 dysfunction and PAH, a downregulation of GGH may have significant implications for the pulmonary vascular cell phenotype.

### 4.4. KCNK3 Dysfunction and Cell Metabolism Changes

In siKCNK3 transfected hPAECs and hPASMCs, we found by the LC/MS approach an activation of glycolysis. PAH pulmonary vascular cells and cancer cells have developed some similar metabolic adaptations [18]. Several studies have shown an increase in glycolysis and a decrease in glucose oxidation in pulmonary vascular cells, similar to the metabolic changes described in proliferating tumor cells [18]. These metabolic changes include the pentose phosphate pathway, glutaminolysis, and fatty acid synthesis and oxidation. In PAH pulmonary vascular cells and in cancer cells, this metabolic shift is mostly driven by HIFs overactivation [18]. Additionally, we recently demonstrated that KCNK3 dysfunction in human and rat PASMCs induces the overexpression of HIF1α [8].

Moreover, we found that KCNK3 knockdown in hPASMCs induced a significant reduction in the glucose-6-phosphate dehydrogenase (G6PD) expression, which is the main producer of NADPH. Clinical reports connect G6PD deficiency with PAH [58]. More, three patients with missense mutations in the G2PD gene were identified in a cohort of PAH patients [59]. Altogether, these reports indicate that G6PD dysfunction could predispose one to the development of PAH. Unexpectedly, an increase in G6PD activity is reported in chronic hypoxia rats [60]. The loss of G6PD was shown to prevent the production of NADPH and a subsequent decrease in NO synthesis, since reduced NADPH is a critical substrate for NO synthesis [61]. These data suggest that the inhibition of G6PD could be deleterious to the PA tone in PH. This is consistent with our recent results, showing that *Kcnk3* deficiency in rats is associated with significant alteration in PA tone, including the impairment of NO-induced relaxation [8].

**Limitations:** Due to the low capability of mass spectrometry to identify hydrophobic channel proteins, which are particularly nonpolar and present in a very limited amount, unfortunately we did not detect a large amount of ion channel families' proteins. However, 33 ion channels or transporters in hPAECs and 28 in hPASMCs (data not shown) were identified without any changes in their expression in the absence of KCNK3.

In the future, the expression of other newly identified signaling pathways will be further investigated by Western blot or PCR analysis in the context of the KCNK3 loss of expression as well as in the more general context of PAH.

In summary, the current study analyses the molecular functional consequences of KCNK3 dysfunction in hPASMCs and hPAECs. Altogether, these findings allow us to better understand the involvement of KCNK3 dysfunction in PAH disease.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/1422-0067/21/19/7400/s1.

**Author Contributions:** Conceptualization: F.A.; Funding acquisition F.A.; Investigation: H.L.R., M.L., T.B., O.M., F.A.; Formal analysis: H.L.R., F.D., G.R., S.Y.-N.; Software: F.D., T.B., G.R. and M.Q.; Supervision: V.C., F.A.; Writing—original draft: F.A.; Writing—review and editing: H.L.R., F.D., B.G., S.C.-K., S.Y.-N., M.H., D.M., V.C., F.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was supported by grants from the French National Institute for Health and Medical Research (INSERM), the Université Paris-Saclay, the Marie Lannelongue Hospital, and the French National Agency for Research (ANR) (grant no. ANR-18-CE14-0023 (KAPAH). H.L.R is supported by ANR-18-CE14-0023. M.L. is supported by the Therapeutic Innovation Doctoral School (ED569). The Region Ile de France has co-funded proteomic instruments (convention n° EX024607).

**Acknowledgments:** The authors thank Yvonne Dutheil from INSERM U999, Hopital Marie Lannelongue, for all her help.
**Conflicts of Interest:** The sponsors had no role in the design, collection, analysis, or interpretation of data, the writing of this article, or the decision to submit it for publication. M.H. and D.M have relationships with drug companies, including Actelion, Bayer, GSK, Novartis, and Pfizer. In addition to being investigators in trials involving these companies, other relationships include consultancy services and memberships to scientific advisory boards. The authors declare no conflict of interest.

## References

1. Ma, L.; Roman-Campos, D.; Austin, E.D.; Eyries, M.; Sampson, K.S.; Soubrier, F.; Germain, M.; Trégouët, D.-A.; Borczuk, A.; Rosenzweig, E.B.; et al. A novel channelopathy in pulmonary arterial hypertension. *New Engl. J. Med.* **2013**, *369*, 351–361. [CrossRef]

2. Lambert, M.; Capuano, V.; Olschewski, A.; Sabourin, J.; Nagaraj, C.; Girerd, B.; Weatherald, J.; Humbert, M.; Antigny, F. Ion channels in pulmonary hypertension: A therapeutic interest? *Int. J. Mol. Sci.* **2018**, *19*, 3162. [CrossRef] [PubMed]

3. Lau, E.M.; Giannoulatou, E.; Celermajer, D.S.; Humbert, M. Epidemiology and treatment of pulmonary arterial hypertension. *Nat. Rev. Cardiol.* **2017**, *14*, 603–614. [CrossRef] [PubMed]

4. Simonneau, G.; Montani, D.; Celermajer, D.S.; Denton, C.P.; Gatzoulis, M.A.; Krowka, M.; Williams, P.G.; Souza, R. Haemodynamic definitions and updated clinical classification of pulmonary hypertension. *Eur. Respir. J.* **2019**, *53*, 1801913. [CrossRef] [PubMed]

5. Olschewski, A.; Veale, E.L.; Nagy, B.M.; Nagaraj, C.; Kwapiszewska, G.; Antigny, F.; Lambert, M.; Humbert, M.; Czirják, G.; Enyedi, P.; et al. TASK-1 (KCNK3) channels in the lung: From cell biology to clinical implications. *Eur. Respir. J.* **2017**, *50*, 1700754. [CrossRef] [PubMed]

6. Enyedi, P.; Czirják, G. Molecular background of leak K + currents: Two-pore domain potassium channels. *Physiol. Rev.* **2010**, *90*, 559–605. [CrossRef] [PubMed]

7. Antigny, F.; Hautefort, A.; Meloche, J.; Belacel-Ouari, M.; Manoury, B.; Rucker-Martin, C.; Péchoux, C.; Potus, F.; Nadeau, V.; Tremblay, E.; et al. Potassium channel subfamily K member 3 (KCNK3) contributes to the development of pulmonary arterial hypertension. *Circulation* **2016**, *133*, 1371–1385. [CrossRef] [PubMed]

8. Lambert, M.; Capuano, V.; Boet, A.; Tesson, L.; Bertero, T.; Nakhleh, M.K.; Remy, S.; Anegon, I.; Péchoux, C.; Hautefort, A.; et al. Characterization of Kcnk3-mutated rat, a novel model of pulmonary hypertension. *Circ. Res.* **2019**, *125*, 678–695. [CrossRef]

9. Humbert, M.; Guignabert, C.; Bonnet, S.; Dorfmüller, P.; Klinger, J.R.; Nicolls, M.R.; Olschewski, A.J.; Pullamsetti, S.S.; Schermuly, R.T.; Stenmark, K.R.; et al. Pathology and pathobiology of pulmonary hypertension: State of the art and research perspectives. *Eur. Respir. J.* **2019**, *53*, 1801887. [CrossRef]

10. Sankhe, S.; Manousakidi, S.; Antigny, F.; Ataam, J.A.; Bentebbal, S.; Ruchon, Y.; Lecerf, F.; Sabourin, J.; Price, L.; Fadel, E.; et al. T-type Ca²⁺ channels elicit pro-proliferative and anti-apoptotic responses through impaired PP2A/Akt1 signaling in PASMCs from patients with pulmonary arterial hypertension. *Biochim. Biophys. Acta (BBA)-Bioenerg.* **2017**, *1864*, 1631–1641. [CrossRef]

11. Perros, F.; Ranchoux, B.; Izikki, M.; Bentebbal, S.; Happé, C.; Antigny, F.; Jourdon, P.; Dorfmüller, P.; Lecerf, F.; Fadel, E.; et al. Nebivolol for improving endothelial dysfunction, pulmonary vascular remodeling, and right heart function in pulmonary hypertension. *J. Am. Coll. Cardiol.* **2015**, *65*, 668–680. [CrossRef] [PubMed]

12. Manaud, G.; Nossent, E.J.; Lambert, M.; Ghigna, M.-R.; Boët, A.; Vinhas, M.-C.; Ranchoux, B.; Dumas, S.J.; Courboulin, A.; Girerd, B.; et al. Comparison of human and experimental pulmonary veno-occlusive disease. *Am. J. Respir. Cell Mol. Biol.* **2020**, *63*, 118–131. [CrossRef]

13. Belouah, I.; Bénard, C.; Denton, A.; Blein-Nicolas, M.; Balliau, T.; Teyssier, E.; Gallusci, P.; Bouchez, O.; Usadel, B.; Zivy, M.; et al. Transcriptomic and proteomic data in developing tomato fruit. *Data Brief.* **2020**, *28*, 105015. [CrossRef] [PubMed]

14. Langella, O.; Valot, B.; Balliau, T.; Blein-Nicolas, M.; Bonhomme, L.; Zivy, M. X!TandemPipeline: A tool to manage sequence redundancy for protein inference and phosphosite identification. *J. Proteome Res.* **2016**, *16*, 494–503. [CrossRef] [PubMed]

15. Valot, B.; Langella, O.; Nano, E.; Zivy, M. MassChroQ: A versatile tool for mass spectrometry quantification. *Proteomics* **2011**, *11*, 3572–3577. [CrossRef] [PubMed]

16. Krämer, A.; Green, J.; Pollard, J.; Tugendreich, S. Causal analysis approaches in ingenuity pathway analysis. *Bioinformatics* **2013**, *30*, 523–530. [CrossRef]
17. Boucherat, O.; Vitry, G.; Trinh, I.; Paulin, R.; Provencher, S.; Bonnet, S. The cancer theory of pulmonary arterial hypertension. *Pulm. Circ.* **2017**, *7*, 285–299. [CrossRef]

18. Humbert, M.; Monti, G.; Brenot, F.; Sitbon, O.; Portier, A.; Grangeot-Keros, L.; Duroux, P.; Galanaud, P.; Simonneau, G.; Emilie, D. Increased interleukin-1 and interleukin-6 serum concentrations in severe primary pulmonary hypertension. *Am. J. Respir. Crit. Care Med.* **1995**, *151*, 1628–1631. [CrossRef]

19. Wang, A.-P.; Li, X.-H.; Yang, Y.-M.; Li, W.-Q.; Zhang, W.; Hu, C.-P.; Zhang, Z.; Li, Y.-J. A critical role of the mTOR/eIF2α pathway in hypoxia-induced pulmonary hypertension. *PLoS ONE* **2015**, *10*, e0130806. [CrossRef]

20. Anda, S.; Zach, R.; Grallert, B. Activation of Gcn2 in response to different stresses. *PLoS ONE* **2017**, *12*, e0182143. [CrossRef]

21. Agarwal, A.; Nick, H.S. Renal response to tissue injury: Lessons from heme oxygenase-1 GeneAblation and expression. *J. Am. Soc. Nephrol.* **2000**, *11*, 965–973. [PubMed]

22. Dulak, J.; Deshane, J.; Jozkowicz, A.; Agarwal, A. Heme oxygenase-1 and carbon monoxide in vascular pathobiology: Focus on angiogenesis. *Circulation* **2008**, *117*, 231–241. [CrossRef] [PubMed]

23. Liu, X.; Wei, J.; Peng, D.H.; Layne, M.D.; Yet, S.-F. Absence of heme oxygenase-1 exacerbates myocardial ischemia/reperfusion injury in diabetic mice. *Diabetes* **2005**, *54*, 778–784. [CrossRef] [PubMed]

24. Nitti, M.; Piras, S.; Marinari, U.M.; Moretta, L.; Pronzato, M.A.; Furfaro, A.L. HO-1 Induction in cancer progression: A matter of cell adaptation. *Antioxidants* **2017**, *6*, 29. [CrossRef] [PubMed]

25. Was, H.; Dulak, J.; Jozkowicz, A. Heme oxygenase-1 in tumor biology and therapy. *Curr. Drug Targets* **2010**, *11*, 1551–1570. [CrossRef] [PubMed]

26. Chen, Y.; Yuan, T.; Zhang, H.; Yan, Y.; Wang, D.; Fang, L.; Lu, Y.; Du, G. Activation of Nrf2 attenuates pulmonary vascular remodeling via inhibiting endothelial-to-mesenchymal transition: An insight from a plant polyphenol. *Int. J. Biol. Sci.* **2017**, *13*, 1067–1081. [CrossRef]

27. Stenmark, K.R.; Frid, M.; Perros, F. Endothelial-to-mesenchymal transition. *Circulation* **2016**, *133*, 1734–1737. [CrossRef]

28. Liu, X.-M.; Peyton, K.J.; Ensenat, D.; Wang, H.; Schafer, A.I.; Alam, J.; Durante, W. Endoplasmic reticulum stress stimulates heme oxygenase-1 gene expression in vascular smooth muscle. *J. Boil. Chem.* **2004**, *280*, 872–877. [CrossRef]

29. Krajka-Kuźniak, V.; Paluszczak, J.; Baer-Dubowska, W. The Nrf2-ARE signaling pathway: An update on its regulation and possible role in cancer prevention and treatment. *Pharmacol. Rep.* **2017**, *69*, 393–402. [CrossRef]

30. Xu, X.-Q.; Lv, Z.-C.; Liu, Q.-Q.; Zhao, Q.-H.; Wu, Y.; Sun, K.; Jiang, X.; Wang, L.; Peng, F.-H.; Jing, Z.-C. Direct bilirubin: A new risk factor of adverse outcome in idiopathic pulmonary arterial hypertension. *Int. J. Cardiol.* **2017**, *228*, 895–899. [CrossRef]

31. Tenhunen, R.; Marver, H.S.; Schmid, R. Microsomal heme oxygenase. Characterization of the enzyme. *J. Boil. Chem.* **1969**, *244*, 6388–6394.

32. Slebos, D.-J.; Ryter, S.W.; Choi, A.M.K. Heme oxygenase-1 and carbon monoxide in pulmonary medicine. *Respir. Res.* **2003**, *4*, 7. [CrossRef] [PubMed]

33. Ryter, S.W. Heme oxygenase-1/carbon monoxide as modulators of autophagy and inflammation. *Arch. Biochem. Biophys.* **2019**, *678*, 108186. [CrossRef]

34. Piantadosi, C.A. Carbon monoxide poisoning. *New Engl. J. Med.* **2002**, *347*, 1054–1055. [CrossRef] [PubMed]

35. Weaver, L.K. Carbon monoxide poisoning. *Crit. Care Clin.* **1999**, *15*, 297–317. [CrossRef]

36. Duckers, H.J.; Boehm, M.; True, A.L.; Yet, S.-F.; San, H.; Park, J.L.; Webb, R.C.; Lee, M.-E.; Nabel, G.J.; Nabel, E.G. Heme oxygenase-1 protects against vascular constriction and proliferation. *Nat. Med.* **2001**, *7*, 693–698. [CrossRef]

37. Morita, T.; Mitsialis, S.A.; Koike, H.; Liu, Y.; Kourembanas, S. Carbon monoxide controls the proliferation of hypoxic vascular smooth muscle cells. *J. Boil. Chem.* **1997**, *272*, 32804–32809. [CrossRef]

38. Zhang, L.; Bouadjel, K.; Manoury, B.; Vandecasteele, G.; Fischmeister, R.; Leblais, V. Cyclic nucleotide signalling compartmentation by PDEs in cultured vascular smooth muscle cells. *Br. J. Pharmacol.* **2019**, *176*, 1780–1792. [CrossRef]

39. Manoury, B.; Idres, S.; Leblais, V.; Fischmeister, R. Ion channels as effectors of cyclic nucleotide pathways: Functional relevance for arterial tone regulation. *Pharmacol. Ther.* **2020**, *209*, 107499. [CrossRef]
40. Fensterl, V.; Sen, G.C. Interferon-induced IFIT proteins: Their role in viral pathogenesis. *J. Virol.* **2014**, *89*, 2462–2468. [CrossRef]

41. Fensterl, V.; Sen, G.C. The ISG56/IFIT1 Gene family. *J. Interf. Cytokine Res.* **2011**, *31*, 71–78. [CrossRef] [PubMed]

42. Pidugu, V.K.; Pidugu, H.B.; Wu, M.-M.; Liu, C.-J.; Lee, T.-C. Emerging functions of human IFIT proteins in cancer. *Front. Mol. Biosci.* **2019**, *6*, 148. [CrossRef] [PubMed]

43. Hsu, K.; Seharaseyon, J.; Dong, P.; Bour, S.; Marbán, E. Mutual functional destruction of HIV-1 Vpu and host TASK-1 channel. *Mol. Cell* **2004**, *14*, 259–267. [CrossRef]

44. Farberov, L.; Herzig, E.; Modai, S.; Isakov, O.; Hizi, A.; Shomron, N. MicroRNA-mediated regulation of p21 and TASK1 cellular restriction factors enhances HIV-1 infection. *J. Cell Sci.* **2015**, *128*, 1607–1616. [CrossRef]

45. Almodovar, S.; Hsue, P.Y.; Morelli, J.; Huang, L.; Flores, S.C. On behalf of the lung HIV study pathogenesis of HIV-associated pulmonary hypertension: Potential role of HIV-1 nef. *Proc. Am. Thorac. Soc.* **2011**, *8*, 308–312. [CrossRef]

46. George, P.M.; Oliver, E.; Dorfmüller, P.; Dubois, O.D.; Reed, D.M.; Kirkby, N.S.; Mohamed, N.A.; Perros, F.; Antigny, F.; Fadel, E.; et al. Evidence for the Involvement of Type I interferon in pulmonary arterial hypertension. *Circ. Res.* **2014**, *114*, 677–688. [CrossRef]

47. Savale, L.; Sattler, C.; Günther, S.; Montani, D.; Chaumais, M.C.; Perrin, S.; Jaïs, X.; Seferian, A.; Jovan, R.; Bulifon, S.; et al. Pulmonary arterial hypertension in patients treated with interferon. *Eur. Respir. J.* **2014**, *44*, 1627–1634. [CrossRef]

48. Cavuoto, P.; Fenech, M. A review of methionine dependency and the role of methionine restriction in cancer growth control and life-span extension. *Cancer Treat. Rev.* **2012**, *38*, 726–736. [CrossRef]

49. The Role of Methionine in Cancer Growth and Control. Available online: https://www.naturalmedicinejournal.com/journal/2015-12/role-methionine-cancer-growth-and-control (accessed on 20 April 2020).

50. Gerner, E.W.; Meyskens, F.L., Jr. Polyamines and cancer: Old molecules, new understanding. *Nat. Rev. Cancer* **2004**, *4*, 781–792. [CrossRef]

51. Santos, J.; Leitão-Correia, F.; Sousa, M.J.; Leão, C. Dietary restriction and nutrient balance in aging. *Oxidative Med. Cell. Longev.* **2015**, *2016*, 1–10. [CrossRef]

52. Houssaïni, A.; Abid, S.; Mouraret, N.; Wan, F.; Rideau, D.; Saker, M.; Marcos, E.; Tissot, C.-M.; Dubois-Randé, J.-L.; Amsellem, V.; et al. Rapamycin reverses pulmonary artery smooth muscle cell proliferation in pulmonary hypertension. *Am. J. Respir. Cell Mol. Biol.* **2013**, *48*, 568–577. [CrossRef] [PubMed]

53. Schweinberger, B.M.; Wyse, A.T. Mechanistic basis of hypermethioninemia. *Amino Acids* **2016**, *48*, 2479–2489. [CrossRef] [PubMed]

54. Gamen, E.; Seeger, W.; Pullamsetti, S.S. The emerging role of epigenetics in pulmonary hypertension. *Eur. Respir. J.* **2016**, *48*, 903–917. [CrossRef] [PubMed]

55. Hautefort, A.; Chesné, J.; Preussner, J.; Pullamsetti, S.S.; Tost, J.; Looso, M.; Antigny, F.; Girerd, B.; Riou, M.; Eddahibi, S.; et al. Pulmonary endothelial cell DNA methylation signature in pulmonary arterial hypertension. *Oncotarget* **2017**, *8*, 52995–53016. [CrossRef] [PubMed]

56. Archer, S.L.; Marsboom, G.; Cipriani, N.; Rehman, J.; Kim, G.; Zhang, H.J.; Toth, P.T.; Svensson, E.C.; Dyck, J.R.; Gomberg-Maitland, M.; et al. Epigenetic attenuation of mitochondrial superoxide dismutase 2 in pulmonary arterial hypertension: A basis for excessive cell proliferation and a new therapeutic target. *Circulation* **2010**, *121*, 2661–2671. [CrossRef] [PubMed]

57. Kim, S.-E.; Hinoue, T.; Kim, M.S.; Sohn, K.-J.; Cho, R.C.; Cole, P.D.; Weisenberger, D.J.; Laird, P.W.; Kim, Y.-I. γ-Glutamyl hydrolase modulation significantly influences global and gene-specific DNA methylation and gene expression in human colon and breast cancer cells. *Genes Nutr.* **2014**, *10*, 444. [CrossRef]

58. Chen, L.Y.; Lai, E.J.; Collins, D.R.; Ostrow, D.N.; Sreenivasan, G.M. A young woman with episodic angioedema, papilledema, and eosinophilia. *Am. J. Hematol.* **2009**, *85*, 124–127. [CrossRef]

59. Kurdyukov, S.; Eccles, C.A.; Desai, A.A.; Gonzalez-Garay, M.L.; Yuan, J.X.-J.; Garcia, J.G.N.; Rafikova, O.; Rafikov, R. New cases of glucose-6-phosphate dehydrogenase deficiency in pulmonary arterial hypertension. *PLoS ONE* **2018**, *13*, e0203493. [CrossRef]
60. Chettimada, S.; Gupte, R.; Rawat, D.; Gebb, S.A.; McMurtry, I.F.; Gupte, S.A. Hypoxia-induced glucose-6-phosphate dehydrogenase overexpression and -activation in pulmonary artery smooth muscle cells: Implication in pulmonary hypertension. *Am. J. Physiol. Cell. Mol. Physiol.* **2014**, *308*, L287–L300. [CrossRef]

61. 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. [CrossRef]

<figure><img src="image_12.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).