Literature DB >> 33644717

Two-pore channels affect EGF receptor signaling by receptor trafficking and expression.

Thomas Müller1, Sonja Grossmann1, Robert Theodor Mallmann1, Carolin Rommel1, Lutz Hein1, Norbert Klugbauer1.   

Abstract

Two-pore channels (TPCs) are key components for regulating Ca2+ current from endosomes and lysosomes to the cytosol. This locally restricted Ca2+ current forms the basis for fusion and fission events between endolysosomal membranes and thereby for intracellular trafficking processes. Here, we study the function of TPC1 and TPC2 for uptake, recycling, and degradation of epidermal growth factor receptor (EGFR) using a set of TPC knockout cells. RNA sequencing analysis revealed multiple changes in the expression levels of EGFR pathway-related genes in TPC1-deficient cells. We propose that a prolonged presence of activated EGFRs in endolysosomal signaling platforms, caused by genetic inactivation of TPCs, does not only affect EGFR signaling pathways but also increases de novo synthesis of EGFR. Increased basal phospho-c-Jun levels contribute to the high EGFR expression in TPC-deficient cells. Our data point to a role of TPCs not only as important regulators for the EGFR transportation network but also for EGFR-signaling and expression.
© 2021 The Author(s).

Entities:  

Keywords:  Biomolecules; Membranes; Molecular Biology

Year:  2021        PMID: 33644717      PMCID: PMC7887427          DOI: 10.1016/j.isci.2021.102099

Source DB:  PubMed          Journal:  iScience        ISSN: 2589-0042


Introduction

Two-pore channels (TPCs) comprise a small family of ion channels with only two representatives in rodents and humans, TPC1 and TPC2. Structurally they are evolutionary intermediates between single-domain TRP and four-domain voltage-gated Na+ or Ca2+ channels (Galione, 2019; Jentsch et al., 2015; Patel, 2015). TPCs are located in membranes of the endolysosomal system and are assumed to play a crucial role for intracellular trafficking processes (Grimm et al., 2017; Marchant and Patel, 2015). As acidic compartments of the endolysosomal system constitute small Ca2+ stores, it can be assumed that TPCs form cytosolic Ca2+ entry pathways, which allow for a transient and locally restricted increase of Ca2+. This elevation of Ca2+ triggers fusion and fission events of endolysosomal membranes and thereby forms the basis for intracellular vesicle trafficking and sorting processes during protein uptake, recycling, and degradation (Grimm et al., 2014; Sakurai et al., 2015). So far, the function of TPCs has mostly been investigated by targeted deletion of single TPC genes in mice. TPC functions have been described in the context of receptor endocytosis, degradation, and recycling; of bacterial protein toxin uptake; and of entry and processing of virus particles (Castonguay et al., 2017; Grimm et al., 2014; Sakurai et al., 2015). Consequently, all these models resulted in distinct phenotypes such as development of fatty liver disease in the case of LDL receptor regulation, an impairment of toxin uptake and severity of intoxication, and endolysosomal trapping of virus particles and prevention of infection. The epidermal growth factor receptor (EGFR) belongs to the ErbB family of growth factor receptors with intrinsic tyrosine kinase activity and is involved in key processes such as cell growth, differentiation, proliferation, and motility (Ceresa and Peterson, 2014). In many tumors EGFR is either upregulated or mutated and numerous recent therapeutic strategies aim at blocking oncogenic EGFR signaling (Shan et al., 2012; Tomas et al., 2014). EGFR signaling is initiated by ligand binding to receptors present at the cell surface, which triggers their dimerization and auto-phosphorylation. Activated receptors are internalized by clathrin-mediated (CME) or clathrin-independent endocytosis (CIE) (reviewed in Bakker et al., 2017). Although ligand binding and initiation of signal transduction occurs at the cell surface, activated EGFR is located for the longest period within endolysosomal membranes where EGFR signaling is still ongoing (Conte and Sigismund, 2016; Sousa et al., 2012; Vieira et al., 1996). The intracellular EGFR transportation network has been well studied over the last decades, and it was shown that ligand concentrations determine the preferential trafficking route of activated receptor (Sigismund et al., 2008). At low EGF concentrations, uptake of EGFRs primarily occurs by CME and EGFRs are transported via early and recycling endosomes back to the plasma membrane (Rappoport and Simon, 2009). At higher EGF concentrations, a saturation effect can be observed and increasing amounts of activated EGFRs are taken up by CIE and are routed via late endosomes to lysosomes for degradation (Bakker et al., 2017). In this context, earlier studies demonstrated a link between EGFR trafficking and TPCs, because deletion of TPC2 causes an accumulation of EGFR in endolysosomal compartments (Grimm et al., 2014; Sakurai et al., 2015). Recent studies indicate that EGFR signaling occurs not only during surface localization (Sousa et al., 2012) but also after internalization of activated EGFR (Conte and Sigismund, 2016) (Wu et al., 2012). For example, transcriptional activity of the ERK1/2 pathway is directly affected by the localization of EGFR in the endolysosomal system (Sousa et al., 2012; Wu et al., 2012). Thus, analysis of EGFR uptake and EGFR trafficking allows one to uncover the potential roles of the endolysosomally localized TPCs within this transportation system. In this study, we took advantage of mouse embryonic fibroblast (MEF) and HeLa cells with single TPC knockouts, and also of a TPC1/2 double knockout to explore its consequences on EGFR signaling. We identify TPCs not only as key regulators for intracellular EGFR trafficking but also for controlling EGFR transcription and surface expression.

Results

Fluorescence microscopy exemplifies an altered EGFR uptake and trafficking in TPC-deficient cells

TPCs play a crucial role for intracellular trafficking processes and affect receptor endocytosis, recycling, and degradation. In this study, we investigated the roles of TPCs for processing and signaling of the EGFR, a well-studied receptor for regulating cell growth not only under physiological but also under pathophysiological conditions. For this purpose, we used MEF cells derived from transgenic mouse lines with a targeted disruption of either TPC1 or TPC2, the only two members of rodent and human TPC family (Arndt et al., 2014; Grimm et al., 2014). Furthermore, a TPC1/2-double knockout MEF cell line was generated in this study by applying the CRISPR-Cas9 system as described in Supplemental information (Figure S1). Initial studies were performed to compare the uptake of EGFR in wild-type and TPC-deficient MEF cells. Cells were incubated for 60 min with 200 ng/mL EGF labeled with Alexa 488. This high concentration favors CIE of the EGF-EGFR complex and initiates trafficking processes that are mainly linked to degradation routes (Sigismund et al., 2008). As a result, TPC-deficient cells internalized higher amounts of labeled EGF than wild-type cells (Figure 1: top panel). In comparison with single knockouts, TPC1/2 double knockouts demonstrated the highest accumulation of EGF. A Rab5 staining highlighting early endosomes did not indicate any differences in the distribution pattern of Rab5-positive compartments in wild-type and knockout cells (Figure 1: middle panel). However, TPC2 and TPC1/2 double knockouts showed numerous examples of co-localization of Rab5- and EGF-positive vesicles, whereas wild-type and TPC1-deficient MEF cells showed only a few (Figure 1: bottom panel and arrowheads in insets in C and D). A quantification of the area covered by EGF/Rab5-positive vesicles indicates significant differences between those groups (Figure S2). The values for wild-type and TPC1-deficient cells were 0.05% ± 0.02% and 0.13% ± 0.05%, respectively. For TPC2 and TPC1/2 double knockouts coverage values were 0.34% ± 0.07% and 0.36% ± 0.11%, respectively. These initial studies indicate that TPCs are involved in different ways in the uptake and intracellular trafficking of the EGF-EGFR complex.
Figure 1

Representative fluorescence microscopic images of MEF cells after 60-min incubation with Alexa 488-labeled EGF

(A–D) (A) Wild-type, (B) TPC1-deficient, (C) TPC2-deficient, and (D) TPC1/2 double knockout MEF cells. Top panel, EGF-Alexa488 fluorescence; middle panel, immunofluorescence staining of Rab5-positive compartments; bottom panel, composition of EGF and Rab5 fluorescence signal together with DAPI staining. Insets in (C and D) indicate examples of EGF- and Rab5-positive vesicles (arrowheads). Scale bar, 10 μm.

Representative fluorescence microscopic images of MEF cells after 60-min incubation with Alexa 488-labeled EGF (A–D) (A) Wild-type, (B) TPC1-deficient, (C) TPC2-deficient, and (D) TPC1/2 double knockout MEF cells. Top panel, EGF-Alexa488 fluorescence; middle panel, immunofluorescence staining of Rab5-positive compartments; bottom panel, composition of EGF and Rab5 fluorescence signal together with DAPI staining. Insets in (C and D) indicate examples of EGF- and Rab5-positive vesicles (arrowheads). Scale bar, 10 μm.

Deletion of TPCs increases EGFR uptake

A fluorescence-activated cell sorting (FACS)-based method was applied to quantify the findings of the initial fluorescence microscopic approach. MEF cells were incubated with high concentrations of EGF-Alexa488 for 10 to 120 min and were subjected to subsequent FACS analysis. The fluorescence intensities correlated with the amount of incorporated EGF and were used to study the kinetics of EGF uptake and receptor trafficking. All Alexa 488 fluorescence histograms showed a rightward shift compared with untreated controls; this effect was strongest for TPC1/2 double knockouts (Figures 2A–2D). Interestingly, TPC2-deficient cells demonstrated the broadest signal distribution. The quantification of the fluorescence signal of all genotypes revealed a continuous increase of fluorescence over the entire experimental period of 120 min and clearly indicated higher values for the knockouts, in particular for the TPC1/2 double knockout (Figure 2E). The single MEF cell knockouts showed a 2- to 3-fold increase, and the double knockout showed a 4- to 6-fold increase over the entire experimental time frame. To generalize and to lay our conclusions on a broader basis we included wild-type and TPC1-deficient HeLa cells into our studies (Castonguay et al., 2017). The FACS analysis was performed in the same way as for MEF cells, and fluorescence histograms again showed a rightward shift when the cells were treated with EGF-Alexa488 (Figures 2F and 2G). The quantification of the fluorescence signals at 10 min, 30 min, and 60 min all resulted in a significant increase in TPC1-deficient HeLa cells when compared with wild-type cells (Figure 2H). In summary, the FACS-based quantification confirmed the observations made in the fluorescence microscopy studies and indicated a faster and appreciably higher uptake of EGF-Alexa488 in TPC-deficient cell lines.
Figure 2

FACS analysis of EGF-Alexa488 uptake in MEF and HeLa cells

(A–D) Representative fluorescence histograms after 60-min incubation with 200 ng/mL EGF-Alexa488 at 37°C for wild-type (A), TPC1-deficient (B), TPC2-deficient (C), and TPC1/2-double knockout (D) MEF cells. Histograms of untreated cells are indicated in light gray and EGF-treated cells in dark gray (wild-type) or in blue colors (knockouts).

(E) Quantification of fluorescence signals for time points t = 10 min, t = 30 min, t = 60 min, and t = 120 min after incubation with EGF-Alexa488. Bar diagrams show mean values and standard errors of the median Alexa 488-fluorescence intensities normalized to an independent wild-type value at t = 180 min. Number of experiments is five (n = 5). Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗∗p ≤ 0.01).

(F and G) Representative fluorescence histograms after 60-min incubation with 200 ng/mL EGF-Alexa488 at 37°C for wild-type (F) and TPC1-deficient (G) HeLa cells. Histograms of untreated cells are indicated in light gray and EGF-treated cells in dark gray (wild-type) or in blue color (TPC1 knockout).

(H) Quantification of fluorescence signals for time points t = 10 min, t = 30 min, and t = 60 min and after incubation with EGF-Alexa488. Bar diagram shows mean values and standard errors of the median Alexa 488-fluorescence intensities normalized to highest values at t = 60 min. The number of experiments is three (n = 3). Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗p ≤ 0.03; ∗∗∗p ≤ 0.01).

(I) Data table showing the calculated fluorescence intensity values and p values.

FACS analysis of EGF-Alexa488 uptake in MEF and HeLa cells (A–D) Representative fluorescence histograms after 60-min incubation with 200 ng/mL EGF-Alexa488 at 37°C for wild-type (A), TPC1-deficient (B), TPC2-deficient (C), and TPC1/2-double knockout (D) MEF cells. Histograms of untreated cells are indicated in light gray and EGF-treated cells in dark gray (wild-type) or in blue colors (knockouts). (E) Quantification of fluorescence signals for time points t = 10 min, t = 30 min, t = 60 min, and t = 120 min after incubation with EGF-Alexa488. Bar diagrams show mean values and standard errors of the median Alexa 488-fluorescence intensities normalized to an independent wild-type value at t = 180 min. Number of experiments is five (n = 5). Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗∗p ≤ 0.01). (F and G) Representative fluorescence histograms after 60-min incubation with 200 ng/mL EGF-Alexa488 at 37°C for wild-type (F) and TPC1-deficient (G) HeLa cells. Histograms of untreated cells are indicated in light gray and EGF-treated cells in dark gray (wild-type) or in blue color (TPC1 knockout). (H) Quantification of fluorescence signals for time points t = 10 min, t = 30 min, and t = 60 min and after incubation with EGF-Alexa488. Bar diagram shows mean values and standard errors of the median Alexa 488-fluorescence intensities normalized to highest values at t = 60 min. The number of experiments is three (n = 3). Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗p ≤ 0.03; ∗∗∗p ≤ 0.01). (I) Data table showing the calculated fluorescence intensity values and p values. To confirm that the above-mentioned effects were indeed caused by deletion of TPCs, we performed rescue experiments. Thus, single TPC knockout MEF cells were transfected either with the corresponding TPC-EGFP-encoding vector or with an EGFP vector as control. After 24 h, cells were stimulated with EGF-TexasRed for 2 h and subjected to FACS analysis. Only cells positive for expression of EGFP constructs were included in the analysis. Both the TPC1 and the TPC2 rescue experiments confirmed that the expression of corresponding TPC-EGFP vector recovered the wild-type phenotype of the MEF cells, i.e., a reduced uptake of EGF-TexasRed (Figure S3).

Lysosomal degradation does not contribute to observed TPC effects

Next, we checked for a possible lysosomal degradation and stability of EGF-Alexa488 conjugate by performing pulse-chase studies. MEF cells of all genotypes were incubated with EGF-Alexa488 for 1 h at 4°C followed by a washing step to remove any unbound fluorescence-labeled EGF. Afterward cells were incubated with non-labeled EGF at 37°C for 30 and 60 min. FACS analysis revealed no significant differences between the genotypes, when comparing fluorescence intensities after 30 and 60 min with initial values (Figure S4). As there was no decline in fluorescence intensities, degradation, inactivation, or outward transfer of the fluorophore can be excluded for at least 60 min. During the entire experiment, labeled EGF either remained within cells or was bound to the cellular surface rendering the occurrence of putative exocytotic processes during previously performed EGF uptake studies unlikely.

TPC-deficient cells express higher amounts of surface-accessible EGFR

Differences in EGFR endocytosis or degradation between wild-type and TPC-deficient MEF cells may also be caused by an altered expression level or membrane distribution of the EGFR. Therefore, EGF-Alexa488 binding experiments and immunoblot-based analyses of EGFR amounts were performed. MEF cells of all genotypes were incubated for 60 min with 200 ng/mL EGF-Alexa488 at 4°C and were evaluated by FACS analysis. As a result, all cells demonstrated a significant binding of labeled EGF (Figures 3A–3D). However, when comparing the different genotypes, we found an increased binding in the order wild-type, TPC1-KO (knockout), TPC2-KO, TPC1/2-doubleKO, highest for the double knockout (Figure 3E). Additionally, cell lysates from all MEF genotypes were analyzed by western blot and quantified using tubulin as loading control.
Figure 3

FACS analysis of EGF-Alexa488 binding to MEF cells

(A–D) Representative fluorescence histograms after 60-min incubation with 200 ng/mL EGF-Alexa488 at 4°C for wild-type (A), TPC1-deficient (B), TPC2-deficient (C), and TPC1/2-double knockout (D) MEF cells. Histograms of untreated cells are indicated in light gray and EGF-treated cells in dark gray (wild-type) or in blue colors (knockouts).

(E) Quantification of fluorescence signals after incubation with EGF-Alexa488 for 1 h. Bar diagram shows mean values and standard errors of the mean (SEM). Alexa 488-fluorescence intensities were normalized to an independent wild-type value. Number of experiments is ten (n = 10). Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗∗p ≤ 0.01). WT: 1.00 ± 0.20. TPC1-KO: 2.88 ± 0.26 (p = 0.007). TPC2-KO: 4.83 ± 0.53 (p < 0.001). TPC1/2-KO: 6.05 ± 0.45 (p < 0.001).

FACS analysis of EGF-Alexa488 binding to MEF cells (A–D) Representative fluorescence histograms after 60-min incubation with 200 ng/mL EGF-Alexa488 at 4°C for wild-type (A), TPC1-deficient (B), TPC2-deficient (C), and TPC1/2-double knockout (D) MEF cells. Histograms of untreated cells are indicated in light gray and EGF-treated cells in dark gray (wild-type) or in blue colors (knockouts). (E) Quantification of fluorescence signals after incubation with EGF-Alexa488 for 1 h. Bar diagram shows mean values and standard errors of the mean (SEM). Alexa 488-fluorescence intensities were normalized to an independent wild-type value. Number of experiments is ten (n = 10). Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗∗p ≤ 0.01). WT: 1.00 ± 0.20. TPC1-KO: 2.88 ± 0.26 (p = 0.007). TPC2-KO: 4.83 ± 0.53 (p < 0.001). TPC1/2-KO: 6.05 ± 0.45 (p < 0.001). The results from these experiments were completely in line with those from the binding studies and indicated the same order of EGFR expression (Figures 4A and 4B). Thus the total amount of EGFR (western blot) as well as surface-accessible EGFR (binding studies) was significantly increased in TPC-deficient MEF cells and was the highest in TPC1/2-doubleKO. Additionally, we investigated EGFR expression in HeLa cells by immunoblot analysis (Figure 4C). The results from these studies indicated roughly a doubling of the EGFR levels in TPC1-deficient HeLa cells and were in line with data from the MEF cells (Figure 4D).
Figure 4

Quantification of EGFR amounts in MEF and HeLa cells by western blot analysis

(A) Representative western blots of wild-type and TPC-KO MEF cell lysates. Top panel, EGFR; bottom panel, tubulin as control.

(B) Quantification of EGFR protein amounts from five independent western blots (n = 5). Bar diagram shows mean values and standard errors of the luminescence signal when compared with wild-type. Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗∗p ≤ 0.01). WT: 100% ± 17%. TPC1-KO: 238% ± 7% (p = 0.005). TPC2-KO: 361% ± 23% (p < 0.001). TPC1/2-KO: 416% ± 32% (p < 0.001).

(C) Representative western blot analysis of wild-type and TPC1-KO HeLa cell lysates. Top panel, EGFR; bottom panel, tubulin as control.

(D) Quantification of EGFR protein amounts from five independent western blots (n = 5). Bar diagram shows mean values and standard errors of the luminescence signals when compared with wild-type. Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗p ≤ 0.03).

Quantification of EGFR amounts in MEF and HeLa cells by western blot analysis (A) Representative western blots of wild-type and TPC-KO MEF cell lysates. Top panel, EGFR; bottom panel, tubulin as control. (B) Quantification of EGFR protein amounts from five independent western blots (n = 5). Bar diagram shows mean values and standard errors of the luminescence signal when compared with wild-type. Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗∗p ≤ 0.01). WT: 100% ± 17%. TPC1-KO: 238% ± 7% (p = 0.005). TPC2-KO: 361% ± 23% (p < 0.001). TPC1/2-KO: 416% ± 32% (p < 0.001). (C) Representative western blot analysis of wild-type and TPC1-KO HeLa cell lysates. Top panel, EGFR; bottom panel, tubulin as control. (D) Quantification of EGFR protein amounts from five independent western blots (n = 5). Bar diagram shows mean values and standard errors of the luminescence signals when compared with wild-type. Datasets were evaluated via ANOVA and Bonferroni post-hoc test (∗∗p ≤ 0.03).

EGFR degradation follows the same kinetics in wild-type and TPC-deficient cells

EGFR levels were investigated by means of western blot analysis in the presence of the protein synthesis inhibitor cycloheximide (10 μg/mL). This approach allows for a quantification of EGFR levels only depending on degradation, but not on de novo synthesis. Serum-starved MEF cells were exposed to 200 ng/mL EGF for up to 120 min, and EGFR levels were determined for each genotype. Control cells were stimulated with EGF in the absence of cycloheximide. As expected, initial EGFR levels varied between genotypes (wild-type [WT]: 100% ± 19%, TPC1-KO: 455% ± 25%, TPC2-KO: 572% ± 21%, TPC1/2-KO: 549% ± 18%) and were normalized to the same starting value to facilitate comparison of degradation kinetics. There were no statistically significant differences in degradation kinetics between wild-type and TPC-deficient MEF cells, indicating that the higher expression levels of EGFR in TPC-deficient cells are not caused by an impaired degradation process (Figure 5).
Figure 5

Degradation of EGFR in the presence of a protein synthesis inhibitor

(A) Representative western blots of wild-type and TPC-deficient MEF cell lysates after different time points of stimulation with 200 ng/mL EGF in the presence of cycloheximide. Control cells were treated with EGF for 120 min in the absence of cycloheximide. Top panel, EGFR; bottom panel, tubulin as control.

(B) Quantification of EGFR protein amounts from four independent western blots (n = 4). Diagram shows decrease in EGFR amount over time in relation to starting levels (0 min). Datasets were evaluated via ANOVA.

Degradation of EGFR in the presence of a protein synthesis inhibitor (A) Representative western blots of wild-type and TPC-deficient MEF cell lysates after different time points of stimulation with 200 ng/mL EGF in the presence of cycloheximide. Control cells were treated with EGF for 120 min in the absence of cycloheximide. Top panel, EGFR; bottom panel, tubulin as control. (B) Quantification of EGFR protein amounts from four independent western blots (n = 4). Diagram shows decrease in EGFR amount over time in relation to starting levels (0 min). Datasets were evaluated via ANOVA.

TPC-deficient cells show higher recovery rates

In the context of regenerative processes, two main routes, recycling of internalized receptor and de novo receptor synthesis, determine the amount of surface-accessible EGFR. To investigate the role of TPCs on these regenerative processes we established protocols that allowed quantifying EGFR surface expression and transcription levels. MEF cells were incubated with a high concentration of 200 ng/mL of non-labeled EGF at 4°C to saturate surface-expressed EGFR. After 1 h, non-labeled EGF was substituted by EGF-Alexa488, and cells were incubated for 10, 30, and 60 min at 37°C. Additionally, to quantify the amount of surface-expressed EGFR at starting time, cells were kept at 4°C for an additional hour with EGF-Alexa488. Quantification of the initial level of surface EGFR was determined by incubation of MEF cells with EGF-Alexa488 alone. Cells of all genotypes were handled in parallel and were analyzed by FACS. The 1-h receptor saturation with non-labeled EGF significantly reduced the subsequent binding of Alexa 488-labeled EGF to MEF cells in all genotypes. In contrast to wild-type cells, which showed a reduction of only 10%, all TPC-deficient cells demonstrated a much stronger drop down of fluorescence intensity (Figure 6A).
Figure 6

Regeneration rates of surface-accessible EGF receptor

(A) Bar diagram shows quantification of EGF-Alexa488 fluorescence signal for different regeneration periods (t = 0 min, t = 10 min, t = 30 min, and t = 60 min). Surface-accessible EGFRs were at first saturated with 200 ng/mL non-labeled EGF and then incubated with EGF-Alexa488. Control cells of each genotype were incubated with EGF-Alexa488 only. The corresponding bars indicate the initial levels of surface-accessible EGFRs. Bar diagram shows mean values and standard errors of the mean (SEM). Alexa 488 fluorescence intensities were normalized to wild-type values. Five independent experiments were performed (n = 5).

(B) Tabular display of measured values. Datasets were compared via ANOVA and Bonferroni post-hoc test for significance versus corresponding control value of each genotype (n.s., p > 0.05).

Regeneration rates of surface-accessible EGF receptor (A) Bar diagram shows quantification of EGF-Alexa488 fluorescence signal for different regeneration periods (t = 0 min, t = 10 min, t = 30 min, and t = 60 min). Surface-accessible EGFRs were at first saturated with 200 ng/mL non-labeled EGF and then incubated with EGF-Alexa488. Control cells of each genotype were incubated with EGF-Alexa488 only. The corresponding bars indicate the initial levels of surface-accessible EGFRs. Bar diagram shows mean values and standard errors of the mean (SEM). Alexa 488 fluorescence intensities were normalized to wild-type values. Five independent experiments were performed (n = 5). (B) Tabular display of measured values. Datasets were compared via ANOVA and Bonferroni post-hoc test for significance versus corresponding control value of each genotype (n.s., p > 0.05). During the following regeneration period from 10 to 60 min, MEF cells of all genotypes showed an increased binding of EGF-Alexa488 (Figures 6A and 6B). Wild-type cells demonstrated only very moderate increases in fluorescence, but reached higher values after 1 h compared with control cells (cells incubated with labeled EGF only) (Figure 6B). All TPC-deficient cells showed a faster regeneration of surface-accessible EGFR. This effect was dependent on the type of TPC knockout and was the highest for TPC1/2-doubleKO (Figure 6B). For wild-type cells, recovery time was approximately 30 min. TPC1-deficient cells did not reach initial levels even after a regeneration window of 60 min. TPC2- and TPC1/2-doubleKO cells only achieved about 60% of recovery after 1 h when compared with starting levels. In view of the absolute amounts of surface-accessible EGFR, all TPC-deficient cells recovered significantly higher EGFR levels than wild-type cells.

Egfr mRNA levels are increased in TPC-deficient cells

The high concentration of 200 ng/mL of EGF favors the non-canonical endosomal and receptor degradation route and points to a de novo receptor synthesis as the main mechanism for receptor recovery. Therefore, we analyzed Egfr transcript levels of all cells by real-time PCR (qPCR). All TPC-deficient cells demonstrated higher transcript levels compared with wild-type cells, in the same order as observed in previous experiments with highest values for TPC1/2-doubleKO showing an about 12-fold increase (Figure S5). Thus, regeneration and qPCR studies clearly indicate that de novo synthesis of EGFR in TPC-deficient cells is the major mechanism responsible for the higher surface expression of EGFR.

RNA sequencing analysis confirms EGFR upregulation and highlights numerous transcriptional changes in EGFR-linked pathways in TPC1-deficient MEF cells

The high expression of EGFR in TPC-deficient MEF cells led us to investigate the consequences of a TPC1 deletion in EGFR-linked pathways. First, RNA sequencing (RNA-seq) coverage tracks for Egfr in wild-type compared with TPC1 knockout cells indicated a more than 3-fold increase of Egfr transcripts in TPC1-KO cells (Figure S6A). These results correspond perfectly to our western blot, qPCR, and FACS data and further substantiates that TPC deficiency causes a strongly increased Egfr expression in MEF cells (Figure S6B). Second, we took advantage of our RNA-seq data and studied the transcriptional changes of major EGFR pathway genes in TPC1-deficient cells. Our analysis focused on pathways involved in cytoskeletal regulation, endocytosis, apoptosis, and protein synthesis and on expression of relevant transcription factors (Table S1). To provide a comprehensive overview, RNA-seq data were transferred into a scheme highlighting the changes of gene expression involved in EGFR pathways (Figure 7). It is evident from this scheme that deletion of TPC1 affects EGFR-dependent pathways in multiple ways.
Figure 7

Differential gene expression of EGFR pathway proteins in TPC1-deficient MEF cells compared with wild-type cells

Color code indicates gradual changes in expression. Upregulated genes are displayed in red and downregulated genes are displayed in blue color. Genes that show no significant differential expression (q > 0.05) in RNA-seq of TPC1-KO versus wild-type MEF cells are shown in white (fold change 1, no difference). Figure was adopted from Brand et al. (2011).

Differential gene expression of EGFR pathway proteins in TPC1-deficient MEF cells compared with wild-type cells Color code indicates gradual changes in expression. Upregulated genes are displayed in red and downregulated genes are displayed in blue color. Genes that show no significant differential expression (q > 0.05) in RNA-seq of TPC1-KO versus wild-type MEF cells are shown in white (fold change 1, no difference). Figure was adopted from Brand et al. (2011).

TPC-deficient cells exhibit prolonged EGFR signaling

Binding of EGF to EGFR results in receptor auto-phosphorylation and initiation of its kinase activity. In accordance with the high surface expression and regeneration of EGFR in MEF cells deficient for TPCs, we investigated possible consequences for EGFR signaling. MEF cells of all genotypes were stimulated with 200 ng/mL EGF for 1, 3, 5, and 10 min and lysed and subjected to SDS-PAGE. Western blots were evaluated for total amounts of ERK1/2, for phosphorylated ERK1/2 (pERK1/2), and for tubulin. Initial ERK1/2 and pERK1/2 levels were determined initially in the absence of EGF. At the beginning of the experiment (t = 0) pERK1/2 was hardly detectable in all probes (Figure 8 and Table S2). EGF stimulation resulted in a fast and transient ERK1/2 phosphorylation from the first to the fifth minute, which was comparable in wild-type and TPC-deficient cells. However, differences emerged after 10 min demonstrating a drop of pERK1/2 compared with the initial levels in wild-type cells, whereas all TPC-deficient cells still showed elevated pERK1/2 levels (Figure 8 and Table S2). As a control, basal ERK1/2 levels were measured for the entire experimental period. Values for ERK1/2 were stable, indicating that availability of non-phosphorylated ERK1/2 was constant for each time point. In summary, the data indicate that deletion of TPCs in MEF cells caused a prolonged activation of EGFR and ERK1/2 signaling.
Figure 8

Quantification of total ERK1/2 and phosphorylated ERK1/2 (pERK1/2) after EGF stimulation

(A) Western blot analysis of pERK1/2 (top panel), total ERK1/2 (middle panel), and tubulin as control (bottom panel) in wild-type, TPC1-, TPC2-, and TPC1/2-deficient MEF cells after incubation with 200 ng/mL EGF for t = 0 min, t = 1 min, t = 3 min, t = 5 min, and t = 10 min.

(B) Time course and quantification of pERK1/2 levels from three independent experiments (n = 3).

(C) Time course and quantification of total ERK1/2 levels from three independent experiments (n = 3). (B and C) Datasets were compared via ANOVA and Bonferroni post-hoc test for significance versus wild-type control (∗p ≤ 0.05; ∗∗∗p ≤ 0.01; n.s. p > 0.05).

Quantification of total ERK1/2 and phosphorylated ERK1/2 (pERK1/2) after EGF stimulation (A) Western blot analysis of pERK1/2 (top panel), total ERK1/2 (middle panel), and tubulin as control (bottom panel) in wild-type, TPC1-, TPC2-, and TPC1/2-deficient MEF cells after incubation with 200 ng/mL EGF for t = 0 min, t = 1 min, t = 3 min, t = 5 min, and t = 10 min. (B) Time course and quantification of pERK1/2 levels from three independent experiments (n = 3). (C) Time course and quantification of total ERK1/2 levels from three independent experiments (n = 3). (B and C) Datasets were compared via ANOVA and Bonferroni post-hoc test for significance versus wild-type control (∗p ≤ 0.05; ∗∗∗p ≤ 0.01; n.s. p > 0.05).

TPC-deficient MEF cells demonstrate higher basal c-Jun phosphorylation levels than wild-type cells

Basal expression of EGFR is mainly regulated by transcription factor SP1, which is typically present in the nucleus in constant amounts. However, several studies indicate that EGFR transcription is further regulated by MAP kinase-activated transcription factor c-Jun (Fang et al., 2014; Mialon et al., 2005; Weston et al., 2004). Therefore, we investigated phosphorylation of c-Jun in wild-type and TPC-deficient cells under basal and EGF-stimulated conditions. Phospho-specific c-Jun antibodies already indicated a strong phosphorylation of c-Jun under basal non-stimulated conditions in TPC1-, TPC2-, and TPC double knockout cells (Figures 9A and 9B). These differences were significant between wild-type and each of the deletion mutants. Stimulation by EGF (200 ng/mL for 5 min) causes higher phospho-c-Jun levels in all MEF cells, with wild-type phospho-c-Jun levels as high as in unstimulated TPC-deficient cells (Figure 9C). This result clearly indicates that TPC-deficient MEF cells already exhibit phospho-c-Jun levels that were achieved in wild-type cells only in the presence of EGF.
Figure 9

Phosphorylation of c-Jun under basal and EGF-stimulated conditions

(A) Representative example of basal phospho-c-Jun levels (top panel), total amount of c-Jun (middle panel), and tubulin loading control (bottom panel) in wild-type, TPC1-, TPC2- and TPC double knockout MEF cells.

(B) Quantification of six independent experiments presented as phospho-c-Jun/c-Jun relation (n = 6).

(C) Bar diagram demonstrating the phospho-c-Jun/c-Jun relation under basal (left) when compared with EGF-stimulated conditions (right) (n = 3).

Phosphorylation of c-Jun under basal and EGF-stimulated conditions (A) Representative example of basal phospho-c-Jun levels (top panel), total amount of c-Jun (middle panel), and tubulin loading control (bottom panel) in wild-type, TPC1-, TPC2- and TPC double knockout MEF cells. (B) Quantification of six independent experiments presented as phospho-c-Jun/c-Jun relation (n = 6). (C) Bar diagram demonstrating the phospho-c-Jun/c-Jun relation under basal (left) when compared with EGF-stimulated conditions (right) (n = 3).

Discussion

The intracellular trafficking network of the EGFR has been subject of numerous studies, and its understanding is of great value for establishing novel and innovative anti-tumor mechanisms (Sigismund et al., 2018). The network forms an ideal model system to study the functions of TPCs for regulation of receptor endocytosis, recycling, and degradation as well as for investigating receptor signaling. In our studies, we used high concentrations of EGF for receptor stimulation to achieve conditions that favor an uptake of EGFRs via CIE leading to lysosomal degradation of most activated receptor proteins (Caldieri et al., 2017). Owing to the preferred localizations of TPC1 in early and TPC2 in late endosomes and lysosomes (Calcraft et al., 2009; Castonguay et al., 2017), our single and double TPC-KO approach allowed us to discriminate between putatively different roles of TPCs in these endolysosomal compartments. The initial fluorescence microscopic studies—using Rab5 as a marker for early endosomes—demonstrated that TPC2- and TPC1/2-deficient MEF cells contained numerous vesicles that were positive for EGF and Rab5. In contrast, wild-type and TPC1-deficient cells showed much less co-localization of EGF and Rab5. These observations indicate that EGFR trafficking is delayed in TPC2-deficient cells. This may be caused by a longer retention time of the receptor in late endolysosomal compartments leading to a partial backlog of EGFR in early endosomes. In previous studies of TPC2-deficient cells, Grimm and colleagues performed co-localization experiments with the lysosomal marker LAMP-1 observing accumulation of EGFR in LAMP1-positive vesicles (Grimm et al., 2014). This also indicates a delay of receptor processing and trafficking and is in line with observations made in our experiments. To gain additional insights into the time-dependent uptake of EGFRs, a FACS-based approach was chosen. Internalization of labeled EGF was quantified over an experimental window of 2 h. For each time point investigated, TPC1 and TPC2 single knockout cell lines demonstrated a 2- to 3-fold higher fluorescence signal than wild-type cells. For the TPC1/2-doubleKO the effect was additive and a 4- to 6-fold increase of EGF internalization was observed. These results are not limited to MEF cells and can be generalized due to our parallel studies with HeLa cells. Rescue experiments confirmed that the phenotype was indeed caused by deletion of TPCs. Re-expression of TPCs in corresponding knockout cells recovered the wild-type phenotype. These observations are in accordance with studies focusing on other model substrates such as LDL to monitor uptake of LDL receptor in TPC2-deficient cells (Grimm et al., 2014). Here, we add the finding that deletion of TPC1 also increases uptake of labeled substrate and that deletion of both channel subtypes causes an additive effect. We hypothesize that TPC1 and TPC2 meet spatial and/or timely diverse functions within the endolysosomal system. Having the specific distribution of the two TPC subtypes in mind (Calcraft et al., 2009; Castonguay et al., 2017) it seems likely that TPC1 is more important for receptor trafficking in early endosomes, whereas TPC2 contributes to transport processes from late endosomes to lysosomes. Stability of the labeled EGF conjugate is of critical importance because premature lysosomal degradation would impair straightforward analysis of our EGFR trafficking studies. Therefore, we performed pulse-chase experiments with EGF-Alexa488 in MEF cells of all genotypes. In none of the cells, degradation or inactivation of the fluorophore occurred within 1 h and was not significantly different between the genotypes. Thus the observed differences between EGF-positive vesicles in wild-type and TPC-deficient cells cannot be attributed to lysosomal degradation of the Alexa 488 conjugate. So far, our studies investigated the consequences of TPC deletion for the uptake of activated EGFR, but did not consider possible differences in the level of surface-expressed EGFR. Therefore, we quantified the amount of surface-expressed EGFRs and total EGFRs in MEF cells by receptor binding experiments and western blot analysis. The results for both approaches were rather comparable with an about 3-fold increase for TPC1-deficient cells, about 4-fold increase for TPC2-deficient cells, and 5- to 6-fold increase for TPC1/2-doubleKO cells. Again, these observations were not limited to MEF cells, but were also found in TPC1-deficient HeLa cells. An altered degradation route might be one reason for the observed differences in Egfr expression levels and pointed us to compare the degradation kinetics of all genotypes. However, in the presence of a protein synthesis inhibitor no differences could be observed between all MEF cell lines. Therefore, the next studies were designed to discriminate between the two major regenerative mechanisms of receptor availability, recycling, and de novo protein biosynthesis. Quantification of surface-accessible EGFR was examined by binding of labeled EGF at low temperatures after saturation of receptors with non-labeled EGF. Regeneration rates were measured for up to 60 min at 37°C. To provide a comprehensive discussion of these studies it is necessary to analyze both the relative changes and the absolute fluorescence values. Wild-type cells had the capacity to reach the initial levels within 30 min, TPC1-deficient cells regenerated to about 80% after 1 h, and TPC2 and TPC1/2-doubleKO cells only achieved about 60% of the starting level. However, when considering the absolute amounts of surface-accessible EGFR levels, regeneration rates of all TPC knockouts were much higher than wild-type rates. Measurement of egfr-transcript levels supported these findings and showed about 4-, 7- and 12-fold increases in TPC1-, TPC2-, and TPC1/2-doubleKO cells. RNA-seq analysis of Egfr transcripts from TPC1-deficient MEF cells also confirmed above-mentioned expression data. Taken together, these results suggest that EGFR de novo synthesis is the main mechanism responsible for the elevated levels of surface-accessible EGFR. The strong increase of EGFR levels in TPC-deficient cells may have distinct consequences on EGFR-dependent signaling pathways. To get a first glimpse of changes in gene expression of EGFR pathway proteins following deletion of TPC1, we performed a detailed RNA-seq analysis when compared with wild-type cells. These data indicated numerous quantitative changes in the gene expression levels of EGFR pathway-related proteins, whereby up- and downregulation was observed. As these results did not favor a single candidate that would explain high EGFR expression, we additionally investigated phosphorylation of selected target proteins. We analyzed ERK1/2 signaling by quantification of phosphorylated ERK1/2 (pERK1/2) in comparison in with total amounts of ERK1/2. Stimulation with EGF resulted in a fast increase of pERK1/2 levels in cells of all genotypes. In TPC-deficient cells ERK1/2 remained in the phosphorylated state for longer than 10 min. At this time, pERK1/2 levels already dropped down to initial non-stimulated levels in wild-type cells. These observations lead to the conclusion that TPC-deficient MEF cells exhibit prolonged EGFR signaling. Total amounts of non-phosphorylated ERK1/2 were unchanged throughout the entire experiment and were not different between the genotypes. Thus, massively increased EGFR expression in TPC-deficient cells does not account for a stronger or faster activation of EGFR, suggesting that the lower levels of EGFR in wild-type cells are already sufficient for maximal activation. A very similar effect was observed by means of a Ned-19 block of TPCs, which significantly increased and extended tyrosine phosphorylation of ERK1/2 (Kilpatrick et al., 2017). Our studies indicate that deletion of TPCs does not only cause dysregulation of endolysosomal trafficking but also affects transcription of the Egfr gene and EGFR signaling. Therefore, we propose a new model linking function of TPCs with EGFR trafficking and signaling (Figure 10). Binding of EGF causes auto-cross-phosphorylation of the EGFR dimer and subsequent receptor activation. As of this time point, EGFR recruits downstream signaling complexes and triggers specific cellular signaling cascades. Noteworthy, EGFR signaling is ongoing as long as the receptor kinase domain is accessible from the cytosolic side (Conte and Sigismund, 2016; Wu et al., 2012). Comparable mechanisms have been reported for numerous other receptors that use endolysosomal compartments as major signaling platforms (Murphy et al., 2009). During maturation from early to late endosomes, EGFR is internalized in multivesicular bodies and receptor signaling is terminated. Ultimately, receptor and ligand degradation occurs in lysosomes. In this context, deletion of TPCs causes a delay of endolysosomal EGFR trafficking and thereby a prolonged EGFR signaling and continuous activation of associated signaling pathways. Exactly that was observed for phosphorylation of ERK1/2: pERK1/2 levels remained unchanged in TPC-deleted MEF cells for longer than 10 min, whereas pERK1/2 levels dropped down in wild-type cells to initial levels. If this mechanism holds true also for other EGFR-activated complexes, it can be hypothesized that further signaling pathways might be affected. Among them c-Jun is of particular interest because it has been shown to be a major factor for Egfr-gene transcription regulation (Fang et al., 2014; Mialon et al., 2005; Weston et al., 2004). We compared phospho-c-Jun levels in wild-type and TPC-deficient MEF cells and identified increased phospho-c-Jun levels in all cells lacking at least one functional TPC gene. Remarkably, this rise was already observable under unstimulated conditions. Thus, we hypothesize that the increased JNK signaling is a major factor that contributes to the high EGFR expression found in TPC-deficient cells. The prolonged EGFR signaling in endolysosomal compartments caused by deletion of TPCs most likely results in increased JNK signaling, which in turn leads to increased Egfr expression, ultimately forming a positive feedback loop. This positive feedback would be an explanation for the strongly increased amounts of surface accessible EGFRs found in TPC-deficient cells.
Figure 10

Model for the function of TPCs for EGFR trafficking and signaling

(A) Dimerization and auto-cross-phosphorylation of EGF receptor after binding of EGF and subsequent initiation of EGFR signaling. Following internalization of EGFR in endocytotic vesicles and fusion of vesicles to early endosomes, EGFR signaling is ongoing as long as intracellular part of EGFR is accessible from cytosol.

(B) Maturation step of early to late endosomes and to multi-vesicular bodies, respectively. Internalization of EGFR into multi-vesicular bodies terminates EGFR signaling.

(C) Late endosomes fuse with lysosomes where EGFR degradation occurs. During maturation from early endosomes to lysosomes the predominant TPC shifts from TPC1 to TPC2.

(D) Activation of ERK1/2 signaling occurs as long as the receptor kinase domain is accessible from cytosol. Deletion of TPCs leads to a prolonged ERK phosphorylation and pERK signaling.

(E) Deletion of TPCs leads to a comparable effect for JNK signaling pathways. Prolonged activation of c-Jun increases transcription of the egfr gene and number of surface-accessible EGFR.

Model for the function of TPCs for EGFR trafficking and signaling (A) Dimerization and auto-cross-phosphorylation of EGF receptor after binding of EGF and subsequent initiation of EGFR signaling. Following internalization of EGFR in endocytotic vesicles and fusion of vesicles to early endosomes, EGFR signaling is ongoing as long as intracellular part of EGFR is accessible from cytosol. (B) Maturation step of early to late endosomes and to multi-vesicular bodies, respectively. Internalization of EGFR into multi-vesicular bodies terminates EGFR signaling. (C) Late endosomes fuse with lysosomes where EGFR degradation occurs. During maturation from early endosomes to lysosomes the predominant TPC shifts from TPC1 to TPC2. (D) Activation of ERK1/2 signaling occurs as long as the receptor kinase domain is accessible from cytosol. Deletion of TPCs leads to a prolonged ERK phosphorylation and pERK signaling. (E) Deletion of TPCs leads to a comparable effect for JNK signaling pathways. Prolonged activation of c-Jun increases transcription of the egfr gene and number of surface-accessible EGFR. In summary, our work sheds light on a novel aspect of TPC function in the endolysosomal system. In addition to the well-known effects on trafficking, TPC deletions also influence gene transcription by ongoing receptor activation in endolysosomal compartments.

Limitations of the study

Our study was performed with MEF and HeLa cells, but not with primary cells. Particularly with regard to tumor cells, the expression of EGFR may vary within wide limits. We cannot rule out the possibility that in some tumors EGFR-related pathways may be affected in different ways.

Resource availability

Lead contact

Further information and requests for resources should be directed to and will be fulfilled by the Lead Contact, Norbert Klugbauer (n.klugbauer@pharmakol.uni-freiburg.de).

Materials availability

Materials generated in this study will be made available upon reasonable request and may require a material transfer agreement.

Data and code availability

Original sequencing data have been deposited in the Short Read Archive at the National Center for Biotechnology Information (NCBI) under the BioProject ID PRJNA694624. RNA-seq dataset generated in this study will be made available upon reasonable request.

Methods

All methods can be found in the accompanying Transparent Methods supplemental file.
  28 in total

1.  Ebola virus. Two-pore channels control Ebola virus host cell entry and are drug targets for disease treatment.

Authors:  Yasuteru Sakurai; Andrey A Kolokoltsov; Cheng-Chang Chen; Michael W Tidwell; William E Bauta; Norbert Klugbauer; Christian Grimm; Christian Wahl-Schott; Martin Biel; Robert A Davey
Journal:  Science       Date:  2015-02-27       Impact factor: 47.728

2.  Suppression of EGFR endocytosis by dynamin depletion reveals that EGFR signaling occurs primarily at the plasma membrane.

Authors:  Leiliane P Sousa; Irit Lax; Hongying Shen; Shawn M Ferguson; Pietro De Camilli; Joseph Schlessinger
Journal:  Proc Natl Acad Sci U S A       Date:  2012-02-27       Impact factor: 11.205

Review 3.  The nuclear epidermal growth factor receptor signaling network and its role in cancer.

Authors:  Toni M Brand; Mari Iida; Chunrong Li; Deric L Wheeler
Journal:  Discov Med       Date:  2011-11       Impact factor: 2.970

4.  Control of EGF receptor signaling by clathrin-mediated endocytosis.

Authors:  A V Vieira; C Lamaze; S L Schmid
Journal:  Science       Date:  1996-12-20       Impact factor: 47.728

5.  Endocytic trafficking of activated EGFR is AP-2 dependent and occurs through preformed clathrin spots.

Authors:  Joshua Z Rappoport; Sanford M Simon
Journal:  J Cell Sci       Date:  2009-04-07       Impact factor: 5.285

6.  Clathrin-mediated internalization is essential for sustained EGFR signaling but dispensable for degradation.

Authors:  Sara Sigismund; Elisabetta Argenzio; Daniela Tosoni; Elena Cavallaro; Simona Polo; Pier Paolo Di Fiore
Journal:  Dev Cell       Date:  2008-08       Impact factor: 12.270

7.  NAADP and the two-pore channel protein 1 participate in the acrosome reaction in mammalian spermatozoa.

Authors:  Lilli Arndt; Jan Castonguay; Elisabeth Arlt; Dorke Meyer; Sami Hassan; Heike Borth; Susanna Zierler; Gunther Wennemuth; Andreas Breit; Martin Biel; Christian Wahl-Schott; Thomas Gudermann; Norbert Klugbauer; Ingrid Boekhoff
Journal:  Mol Biol Cell       Date:  2014-01-22       Impact factor: 4.138

Review 8.  Emerging functions of the EGFR in cancer.

Authors:  Sara Sigismund; Daniele Avanzato; Letizia Lanzetti
Journal:  Mol Oncol       Date:  2017-11-27       Impact factor: 6.603

Review 9.  Two-Pore Channels: Catalyzers of Endolysosomal Transport and Function.

Authors:  Christian Grimm; Cheng-Chang Chen; Christian Wahl-Schott; Martin Biel
Journal:  Front Pharmacol       Date:  2017-02-07       Impact factor: 5.810

Review 10.  EGF receptor trafficking: consequences for signaling and cancer.

Authors:  Alejandra Tomas; Clare E Futter; Emily R Eden
Journal:  Trends Cell Biol       Date:  2013-11-29       Impact factor: 20.808

View more
  2 in total

Review 1.  Regulation of Aging and Longevity by Ion Channels and Transporters.

Authors:  Kartik Venkatachalam
Journal:  Cells       Date:  2022-03-31       Impact factor: 6.600

2.  Two-pore channels: going with the flows.

Authors:  Anthony J Morgan; Lora L Martucci; Lianne C Davis; Antony Galione
Journal:  Biochem Soc Trans       Date:  2022-08-31       Impact factor: 4.919

  2 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.