Next Article in Journal
A Case Study of the Morphological and Molecular Variation within a Ciliate Genus: Taxonomic Descriptions of Three Dysteria Species (Ciliophora, Cyrtophoria), with the Establishment of a New Species
Next Article in Special Issue
Involvement of ORAI1/SOCE in Human AML Cell Lines and Primary Cells According to ABCB1 Activity, LSC Compartment and Potential Resistance to Ara-C Exposure
Previous Article in Journal
A Major Intestinal Catabolite of Quercetin Glycosides, 3-Hydroxyphenylacetic Acid, Protects the Hepatocytes from the Acetaldehyde-Induced Cytotoxicity through the Enhancement of the Total Aldehyde Dehydrogenase Activity
Previous Article in Special Issue
Role of Orai3 in the Pathophysiology of Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Developing a Mathematical Model of Intracellular Calcium Dynamics for Evaluating Combined Anticancer Effects of Afatinib and RP4010 in Esophageal Cancer

1
College of Nursing and Health Innovation, The University of Texas at Arlington, Arlington, TX 76019, USA
2
Department of Mathematics, The University of Texas at Arlington, Arlington, TX 76019, USA
3
Institute of Science and Technology Austria, Am Campus 1, 3400 Klosterneuburg, Austria
4
Department of Nutritional Sciences, University of Connecticut, Mansfield, CT 06269, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(3), 1763; https://doi.org/10.3390/ijms23031763
Submission received: 13 October 2021 / Revised: 23 January 2022 / Accepted: 25 January 2022 / Published: 3 February 2022

Abstract

:
Targeting dysregulated Ca2+ signaling in cancer cells is an emerging chemotherapy approach. We previously reported that store-operated Ca2+ entry (SOCE) blockers, such as RP4010, are promising antitumor drugs for esophageal cancer. As a tyrosine kinase inhibitor (TKI), afatinib received FDA approval to be used in targeted therapy for patients with EGFR mutation-positive cancers. While preclinical studies and clinical trials have shown that afatinib has benefits for esophageal cancer patients, it is not known whether a combination of afatinib and RP4010 could achieve better anticancer effects. Since TKI can alter intracellular Ca2+ dynamics through EGFR/phospholipase C-γ pathway, in this study, we evaluated the inhibitory effect of afatinib and RP4010 on intracellular Ca2+ oscillations in KYSE-150, a human esophageal squamous cell carcinoma cell line, using both experimental and mathematical simulations. Our mathematical simulation of Ca2+ oscillations could fit well with experimental data responding to afatinib or RP4010, both separately or in combination. Guided by simulation, we were able to identify a proper ratio of afatinib and RP4010 for combined treatment, and such a combination presented synergistic anticancer-effect evidence by experimental measurement of intracellular Ca2+ and cell proliferation. This intracellular Ca2+ dynamic-based mathematical simulation approach could be useful for a rapid and cost-effective evaluation of combined targeting therapy drugs.

1. Introduction

Esophageal cancer (EC) is the sixth leading cause of cancer mortality worldwide [1]. In 2021, in the United States alone, there is an estimated number of 19,260 new cases of EC and 15,530 deaths from EC (www.seer.cancer.gov/statfacts, accessed on 24 January 2022). Two main types of EC are distinguished by different etiological and pathological characteristics: esophageal squamous cell carcinoma (ESCC) and adenocarcinoma (EAC). While EAC is more prevalent in the USA, ESCC predominates among Asians and male African Americans [2]. Both cancers remain asymptomatic; therefore, patients are usually diagnosed at relatively late stages, with an overall 5-year survival rate of below 20%, according to the SEER Cancer Statistics review [3]. The need for development of accurate and timely treatments for EC, thus, is crucial. In recent years, targeting dysregulated Ca2+ signaling in cancer cells has become an active research area to develop new chemotherapy drugs [4].
Ca2+ signaling plays an important role in cell proliferation, apoptosis, autophagy, migration and cell cycle; thus, its dysregulation is associated with tumor initiation, angiogenesis, progression and metastasis [4]. The intracellular Ca2+ signals have different forms, such as Ca2+ spikes, waves and oscillations [5,6]. They are regulated by both intracellular Ca2+ release and extracellular Ca2+ influx [7]. Store-operated Ca2+ entry (SOCE) is a ubiquitous important extracellular Ca2+ influx, which is mainly mediated by two proteins, i.e., stromal-interacting molecule 1 (STIM1) as endoplasmic reticulum (ER) Ca2+ storage sensor and Orai1 as plasma membrane (PM) Ca2+ channel. During activation of SOCE, the depletion of ER Ca2+ stores triggers translocation of STIM1 to ER-PM junctions, followed by binding and activation of Orai1 channel at PM [8]. Accumulating evidence has shown that STIM1/Orai1-mediated SOCE actively participates in the progression of many cancers [9], such as breast cancer [10], pancreatic adenocarcinoma [11] and prostate cancer [12]. We previously reported that high expression of Orai1 in tumor tissues is associated with poor prognosis in ESCC patients, and SOCE-mediated intracellular Ca2+ oscillations regulate cell proliferation, migration and invasion in ESCC cells [13]. Decreased Orai1-mediated SOCE, either by gene knockdown or pharmacological channel blockers, is able to reduce the frequency of intracellular Ca2+ oscillations in cultured ESCC cells and to inhibit tumor growth in preclinical animal models. Among many SOCE blockers, RP4010 is a recently developed one in clinical trial phase I/IB. Our published study suggested that RP4010 is a promising chemotherapy drug targeting SOCE for ESCC patients [14].
The overexpression and/or mutations of epidermal growth factor receptor (EGFR) family proteins, i.e., EGFR (HER1), ERBB2, ERBB3 and ERBB4, are often found in multiple types of cancer cells and are considered to be a significant prognostic indicator in the clinical intervention of cancer. Upon binding with the EGF or other ligands, EGFR undergoes a conformational change on the extracellular domain, followed by dimerization and trans-phosphorylation of tyrosine kinases in the intracellular domain [15]. Activated EGFR tyrosine kinase triggers multiple signaling pathways. One is through phosphorylation of phospholipase C-γ (PLC-γ)/inositol-1,4,5-trisphosphate (IP3)/IP3 receptor (IP3R) pathway and results in intracellular Ca2+ release from ER [16]. Others include downstream PI3K/AKT and MEK/ERK signaling pathways, which are essential for cancer progression. Accordingly, many tyrosine kinase inhibitors (TKIs) have been developed, and more than 20 have received FDA approval as targeting cancer therapy for head and neck, lung, breast and colon cancers [17]. Afatinib is one of the second-generation TKIs and has an irreversible inhibition on both EGFR and HER2 [18]. As its regulatory approval for use as a treatment for non-small-cell lung cancer and squamous cell carcinoma of the lung, many preclinical and clinical studies support its benefit for patients with other cancers, including recurrent and/or metastatic ESCC [19]. However, a challenge of using afatinib, as well as other TKIs, for chemotherapy is that most patients show TKI resistance after initial response [20]. Consequently, combinational treatments of chemotherapy, immunotherapy or other kinds of targeting therapy have been proposed to improve the treatment effect and patient survival rate. Considering that EGFR stimulates intracellular Ca2+ release and SOCE-mediated Ca2+ influx regulates EGFR downstream AKT and ERK activity, we hypothesize that combined afatinib and RP4010 may achieve an enhanced anticancer effect in ESCC cells. Since the dynamic of intracellular Ca2+ is at the intersection between the two signaling pathways and it can be rapidly simulated by a mathematical model, this study used both mathematical simulation and experimental data of the dynamics of intracellular Ca2+ in cultured ESCC cell lines to evaluate the combined anticancer effects of afatinib and RP4010.

2. Results

2.1. Inhibited Intracellular Ca2+ Oscillations by Afatinib in ESCC Cells

We extracted RNA-Seq data from the TCGA database to compare the expression of EGFR in tumor tissues removed from ESCC patients with normal human esophageal tissues (Figure 1a). The average expression levels, presented as transcripts per million (TMP) in tumor and non-tumor tissues, were 38.4 and 20.4, respectively, suggesting upregulated expression of EGFR in ESCC tumor tissues. Interestingly, more than 16% ESCC patients showed an increased amplification of EGFR, whereas mutation was much less common. In order to reveal the expression of individual EGFR family proteins in ESCC cells, Western blot was employed by using specific antibodies against each EGFR family protein. We compared their expression in four human ESCC cell lines (KYSE-30, KYSE-150, KYSE-70 and KYSE-790) and a non-tumorous esophageal epithelial cell line (Het-1A). All four proteins, i.e., EGFR, ERBB2, ERBB3 and ERBB4, were expressed in ESCC cells, but only EGFR and ERBB2 were significantly upregulated in ESCC cells compared to those in Het-1A cells (Figure 1b,c). As described earlier, afatinib is an effective TKI for EGFR/ERBB2; therefore, we suspected that these ESCC cells should be sensitive to afatinib. The KYSE-150 cell line was selected as a representative ESCC cell model in the following studies.
Our previous study suggested intracellular Ca2+ oscillations as an important Ca2+ code for activation of downstream signaling pathways in ESCC cells [13,14]. Since EGFR regulates Ca2+ release through PLC-γ/IP3/IP3R axis, we next tested whether afatinib could impair intracellular Ca2+ oscillations. KYSE-150 cells were loaded with Fluo-4 AM, a fluorescent Ca2+ indicator, followed by time-lapse live-cell imaging. More than half of the cells underwent intracellular Ca2+ oscillations, as indicated by spontaneous fluctuation in fluorescence (Figure 1c, upper panel, control). However, treatment of 5 μM afatinib for 4 h clearly impaired intracellular Ca2+ oscillations in KYSE-150 cells, as evidenced by almost no spontaneous fluctuation in fluorescence within the 5 min recording window (Figure 1d, lower panel, afatinib).

2.2. Decreased Frequency of Intracellular Ca2+ Oscillations by Treatment of RP4010 or Afatinib

The previous studies from our and other groups have shown that the frequency is the essential parameter of intracellular Ca2+ oscillations to regulate downstream signaling pathways [13]; therefore, we evaluated the impact of RP4010 or afatinib on the frequency of Ca2+ oscillations in KYSE-150 cells. Following the same procedure as described earlier, we performed live-cell Ca2+ imaging in KYSE-150 cells treated with different concentrations of afatinib or RP4010 (Table 1 and Figure 2). The fluorescent intensity changes in each individual cell, reflecting intracellular Ca2+oscillations, were extracted as ΔF/F0 (representative traces in Figure 2b–d). The period between the two peaks of the ΔF/F0 curve was calculated as T (second) in each group, and frequency (Hz) was calculated as f = 1/T. KYSE-150 cells had active intracellular Ca2+ oscillations with a frequency of 0.030 Hz (Figure 2b). Treatment of either RP4010 or afatinib significantly decreased the frequency of intracellular Ca2+ oscillations in KYSE-150 cells (Figure 2c,d). Interestingly, both RP4010 and afatinib affected the periods of Ca2+ oscillations in a dose-dependent manner (Table 1). Treatment of 10 µM RP4010 and 5 µM afatinib extended the period of intracellular Ca2+ oscillations in KYSE-150 cells from 32.8 s to 122.1 and 73.4 s, respectively. Similar results were observed in another ESCC cell line (Supplementary Table S1).

2.3. Mathematical Simulation of Intracellular Ca2+ Oscillations

We simulated the dynamics of [Ca2+]cyto in KYSE-150 cells by using a mathematical model, as described in “Mathematical Modeling”. The parameter values in Table 2 were taken from the model presented by Sneyd et al. [6], with the exceptions of kf (IP3/IP3R-mediated ER Ca2+ release) and α1 (SOCE-mediated Ca2+ influx). By adjusting the key parameters for α1, kf and SERCA pump-mediated ER Ca2+ uptake (Vserca) (Table 2), we were able to simulate the oscillating patterns of [Ca2+]cyto in KYSE-150 cells with the period at 32.8 s (Figure 3a, green trace). Our previous report showed that non-tumorous esophageal epithelial cells contain less active SOCE and ER Ca2+ release [13]. Accordingly, we adjusted the values of α1 and kf by dividing those numbers in cancer cells by 3 and 5, respectively. We were able to obtain extended period of intracellular Ca2+ oscillations at 600.9 s in Het-1A cells (Figure 3a, blue trace); this is consistent with our previous report of quiescent intracellular Ca2+ oscillations in these cells. The Michaelis–Menten kinetics parameters, ci and ki, were determined by fitting to experimental data relating to the period of Ca2+ oscillations to a given dose of one or both drugs.
We then used a global sensitivity analysis, using a multidimensional parameter space to assess the effects of uncertainties of the drug parameters on the period of oscillation. For this purpose, we used a joint Latin hypercube sampling (LHS) method and partial rank correlation coefficient (PRCC) analysis (LHS-PRCC). In this method, we started with the eight important uncertain parameters, namely kf, α1, c1, c2, c3, k1, k2 and k3, which correspond to the dynamics of the SOCE channels and the drug interactions. We then used a Monte Carlo simulation technique to generate M random numbers (using a Weibull distribution) for each of the parameters. Each such set of numbers corresponded to a parameter vector that is used to compute the period of oscillation, and, thus, we obtained M periods of oscillation. We finally computed the PRCC and the p-values (using the Student’s t-statistic) to ascertain the most sensitive parameters that correspond to high PRCC values or p-values < 0.05. The detailed steps of the LHS-PRCC algorithm for sensitivity analysis can be found in [21,22]. The baseline parameter values are listed in Table 3, together with a brief summary of the p-values and the PRCC coefficients.
From Table 3, we can see that, except for kf, c1 and k1, the other five parameters are the sensitive ones that affect the period of oscillation. Out of the sensitive parameters, α1 and c2 are the most sensitive ones. We also note that the corresponding PRCC values for α1 and c2 are −0.831 and 0.812, respectively. This implies that, while an increase in α1 is responsible for the decrease in the period of oscillation, the opposite holds for c1.
Next, we simulated the impact of RP4010 or afatinib on [Ca2+]cyto by introducing variables in the computational model, as described in “Addition of TKI and SOCE Blocker in the Computational Model”. The simulated periods of intracellular Ca2+ oscillations also presented the dose-dependent inhibitory effect of RP410 or afatinib. For example, when 0.5 or 3.0 µM RP4010 was included in the computational simulation, the periods were extended to 46.3 and 83.8 s, respectively (Figure 3b,c). When 0.5 µM afatinib was included, the periods were extended to 47.4 s. The simulated periods for each compound fitted well with the experimental data (Table 1; and Figure 4c, red and blue curves).

2.4. Effects of Combined RP4010 and Afatinib on Intracellular Ca2+ Oscillations

We then simulated [Ca2+]cyto, including both RP4010 and afatinib, into the computational model. Based on the dose-dependent curve for each individual drug, we selected a molar ratio of 2:1 for RP4010 and afatinib for the following studies. As shown in Figure 4a, the computer simulated [Ca2+]cyto displayed oscillations with a period of 92.8 s when 2 µM RP4010 + 1 µM afatinib was included into the calculation. To validate the predicted period of intracellular Ca2+ oscillations, we conducted live-cell intracellular Ca2+ measurements in KYSE-150 cells treated with combined 2 µM RP4010 + 1 µM afatinib (Figure 4b). The calculated period of intracellular Ca2+ oscillations from experimental data was 89.4 s, suggesting that the mathematical model could predict the intracellular Ca2+ oscillations in cells treated with two different drugs. The periods of intracellular Ca2+ oscillations in cells treated with either 2 µM RP4010 or 1 µM afatinib alone were recorded as 56.5 and 52.7 s, respectively (Table 1). The combined treatment of RP4010 and afatinib clearly had a better impact on reducing intracellular Ca2+ oscillations than that caused by one of the single drugs alone. Similarly, we performed simulation and experimental measurement of [Ca2+]cyto treated with combined RP4010/afatinib at other concentrations (all in µM), ranging from 0.125/0.00625 to 5.0/2.5 (Table 4 and Figure 4c). The dose-dependent curve of the period in the combination group (green curve) was above that of both the RP4010 group (blue curve) and afatinib group (red curve), thus indicating that the combined RP4010 and afatinib may have a greater impact on intracellular Ca2+ oscillations than any single agent. A further analysis using the CompuSyn software showed that the combination index (CI) of RP4010 and afatinib is less than 1 (Supplementary Figure S1). Similar results were observed in another ESCC cell line (Supplementary Table S2).

2.5. Synergistical Effects of RP4010 and Afatinib on Cytotoxicity in ESCC Cells

The intracellular Ca2+ oscillations regulate cell proliferation and apoptosis in ESCC cells [13]. We examined the combined effects of RP4010 and afatinib on cancer-cell viability. KYSE-150 cells were treated with 1 µM RP4010, 0.5 µM afatinib or a combination for 48 h. The cells were stained with cell-permeable nuclear dye Hoechst 33342 to visualize the nuclear morphology. While the cells in the vehicle control group showed complete and intact nuclear envelope and normal morphology (Figure 5a, left panel) and 1 µM RP4010 or 0.5 µM afatinib groups showed slightly altered nuclear morphology, such as nuclear shrinkage (middle two panels), combination treatment of 1 µM RP4010/0.5 µM afatinib induced chromatin condensation, DNA fragmentation and apoptotic body formation (right panel, indicated by white arrow heads). All the nuclear characteristic morphological features indicated that these cells underwent apoptosis. Then MTT assay was employed to evaluate the cell viability, and the dose-dependent curves were analyzed in KYSE-150 cells treated with either RP4010 or afatinib individually or in combination with each other (Figure 5b). The concentrations in the RP4010 group (blue curve) were 10, 5, 2.5, 1.25 and 0.625 µM, respectively. The concentrations in afatinib group (red curve) were 5, 2.5, 1.25, 0.625 and 0.3125 µM, respectively. The combination group (green curve) contained RP4010/afatinib (2:1) at concentrations from 5 µM/2.5 µM to 0.625 µM/0.3125 µM (as in Table 2). While RP4010 or afatinib potently reduced cell viability in KYSE-150 cells in a dose-dependent manner (Figure 5b), their combination demonstrated enhanced inhibitory function on cell viability, especially in the low concentrations (Figure 5b, green curve). To further evaluate whether the enhanced function in combined treatment is through simple addition or synergy, we performed a combination index (CI) analysis by using the CompuSyn platform. The CIs at the tested concentrations were less than 1 (CI < 1), indicating that RP4010 and afatinib worked synergistically to inhibit cell proliferation in KYSE-150 cells.

3. Discussion

Targeting SOCE-mediated Ca2+ signaling in cancer cells is an emerging chemotherapy approach, and several SOCE blockers, including RP4010, are currently being evaluated in clinical trials [23]. As a second-generation TKI, afatinib received FDA approval to be used in targeted therapy for patients with EGFR mutation-positive cancers, but not for ESCC. To improve the treatment response and patient survival rate, it is a common practice to combine several treatments or several chemotherapeutic drugs. In this study, we examined whether a combination of afatinib and RP4010 could achieve better anticancer effects in ESCC cells. Using both mathematical simulation and live-cell intracellular Ca2+ measurement, we evaluated the inhibitory effect of afatinib and RP4010 on intracellular Ca2+ oscillations in KYSE-150 cells. This computer-aided mathematical model successfully predicted the frequency of intracellular Ca2+ oscillations in KYSE-150 cells responding to RP4010 and afatinib either individually or in combination. Both experimental and mathematical simulation data showed that a combination of afatinib and RP4010 could synergistically reduce the frequency of intracellular Ca2+ oscillations and helped to predict their effects on cell viability in KYSE-150 cells. This intracellular Ca2+ dynamic-based mathematical simulation approach could be used as a rapid and cost-effective evaluation of combined targeting therapy drugs.
The anticancer function of TKI has been mainly studied in the aspect of inhibiting downstream PI3K/Akt and MEK/ERK pathways; however, its PLC-γ axis has received less attention. During the preparation of this article, Kim et al. reported that gefitinib, a first-generation reversible TKI, inhibits EGF-stimulated intracellular Ca2+ oscillations in non-small-cell lung cancer cells and that restricting extracellular Ca2+ can consequently enhance gefitinib sensitivity [24]. As an irreversible TKI, afatinib potently inhibits signaling from all EGFR family receptor homodimers and heterodimers and the downstream events. Our experimental data showed that afatinib could reduce the frequency of intracellular Ca2+ oscillations in a dose-dependent manner (Figure 1c, Figure 2d and Figure 4c). The working model is that afatinib can regulate the Ca2+ release from ER through the EGFR–PLC–PIP2–IP3–IP3R axis and, thus, reshape the intracellular Ca2+ oscillations (illustration in Figure 6). Another major factor to control the intracellular Ca2+ oscillations is extracellular Ca2+ influx through the SOCE pathway, which is supported by our earlier reports [13,14]. This notion was further confirmed by this study. The data showed that RP4010 could potently inhibit the frequency of intracellular Ca2+ oscillations in ESCC cells in a dose-dependent manner (Figure 4c). It is worthwhile to note that SOCE-mediated Ca2+ influx regulates AKT and ERK1/2, two key molecules in downstream signaling pathways of EGFR, and activation of EGFR tyrosine kinase triggers the phosphorylation of STIM1 at ERK1/2 target sites to active SOCE [25]. Therefore, besides intracellular Ca2+, the crosstalk between EGFR and SOCE signaling pathways may have other interacting points, and these points require further investigation.
The intracellular Ca2+ oscillations are the event downstream of both EGFR and SOCE signaling pathways and control consequent gene transcription through Ca2+-dependent enzymes, including NF-κB, MAPK/ERK, CAMKII, etc. [26] Therefore, the intracellular Ca2+ oscillations could be used as a readout to evaluate any chemotherapeutic drug targeting EGFR or SOCE pathways. Fluorescence-based live-cell imaging can conveniently record the dynamics of intracellular Ca2+ measurement and can be easily up scaled as high-throughput analysis for drug screening. Another advantage using intracellular Ca2+ oscillations as a readout for drug evaluation is that a mathematical-modeling-based toolbox can be established to simulate [Ca2+]cyto. While the experimental data on multiple drugs are difficult to obtain and many of the drugs are expensive to acquire, the computer-aid mathematical simulation could greatly facilitate the mechanistic understanding of cross-talk between two pathways and provide rapid cost-effective test on hundreds of drugs and countless combinations with different ratios and concentrations. In this study, we established a mathematical model to simulate the intracellular Ca2+ oscillations by using a set of parameters obtained from experimental data (Figure 3 and Figure 4). The frequency of intracellular Ca2+ oscillations was the focus to adjust the parameters, using experimental data from a single drug, either RP4010 or afatinib. Then the mathematical model successfully predicted the synergetic action when both drugs were included into consideration, suggesting that this mathematical simulation system could provide reliable evaluation on the anticancer effects of combined TKIs and SOCE blockers.
Both the intracellular Ca2+ oscillations and cell-viability analysis showed that there is a synergetic effect between RP4010 and afatinib (Figure 4 and Figure 5c). This could be explained by the crosstalk between EGFR and SOCE signaling pathways. The aforementioned issue in many failed TKI clinical trials in ESCC is due to short-term adaptive responses and long-term acquired resistant of cancer cells to TKIs. Targeting both EGFR and SOCE simultaneously will provide synergetic anticancer effects to kill the cancer cells and effectively avoid their adaptive responses. Further studies are urgently needed to explore the combined TKIs and SOCE blockers for better chemotherapy, not only for esophageal cancer but also for other cancers, as well.
This study has two limitations. First, this mathematical model was able to simulate the frequency but not the shape and amplitude of Ca2+ oscillations. The model used here was a modified version of one developed by Sneyd et al. [6], with influence from Dupont et al. [27]. The frequency modulation of signal transduction includes NF-κB, MAPK, etc. [24], and our previous studies showed that the main downstream signaling pathway of Ca2+ oscillations in ESCC cells is NF-κB [14]. Therefore, we selected this relatively simple deterministic model to guide the experimental drug combination design. In the future, the model should be improved to include biochemical reactions, stochasticity, and to add spatial effects and other factors to better simulate both frequency and amplitude/shape of Ca2+ oscillations. Second, this study was conducted only in cultured human ESCC cells without in vivo analysis. Future studies on the combined chemotherapy effects of afatinib and RP4010 are warranted to be conducted in ESCC animal models, such as orthotopic xenografted mice or NMBA-induced rats. It will also be interesting to use this mathematical simulation system to evaluate many FDA-approved TKIs other than afatinib and SOCE blockers. Regardless of the limitations, this study demonstrated that a simulation model of the intracellular Ca2+ oscillations has the potential to be developed as an effective and rapid evaluation tool to be used to study whether the effects of the combination of drugs are synergistic, antagonistic or additive and to aid selection of the combinations of drugs producing the desired effects with the minimal dose. It could provide guidance for clinicians to determine optimal combination of TKIs and SOCE blockers.

4. Materials and Methods

4.1. Materials

RP4010 was obtained from Rhizen Pharmaceuticals, S.A (La Chaux-de-Fonds, Switzerland). The compound was dissolved in DMSO to make up a 10 µM stock solution. Afatinib was bought from Selleckchem company (Houston, TX, USA) and dissolved in DMSO to make a stock solution of 10 µM. Human ESCC (KYSE-30, KYSE-150, KYSE-70 and KYSE-790) and normal epithelial cell line (Het-1A) were used in this study [13].

4.2. Cell Culture

ESCC and Het-1A cells were cultured at 37 °C in a 5% CO2 incubator and maintained in a 1:1 mixture of RPMI-1640 medium and Ham’s F12 Medium (Corning, NY, USA) supplemented with 5% fetal bovine serum (VWR, Radnor, PA, USA) and 1% penicillin/streptomycin (Corning, NY, USA).

4.3. Intracellular Ca2+ Oscillations Measurement

KYSE-150 cells were seeded in a black 96-well plate with a clear bottom. After attachment, cells were treated with RP4010 and afatinib, combined or separately, at different concentrations (details are shown in Table 1), for 4 h. For RP4010, the concentration varied from 5 to 0.25 µM. The concentration of afatinib varied from 5 to 0.25 µM. For the combination of RP4010 and afatinib, the ratio of RP4010 to afatinib was 2:1. KYSE-150 cells were loaded with 5 µM Fluo-4 in 96-well imaging plates (BD Biosciences, Franklin Lakes, NJ, USA) at room temperature for 20 min. After washing, cells were kept in culture medium without phenol red. The intensity of fluorescent signals was recorded by using a Hamamatsu digital camera C11440 complemented with DMi8 inverted microscope (Leica, Wetzlar, Germany) with 20x objective (dry lens, NA 0.75). Time-lapse live-cell images were recorded every 5 s, for a total time period of 5 min, and the period between 2 peaks (2 Ca2+ oscillations) was measured in order to calculate the corresponding periods.

4.4. Western Blot

Het-1A, KYSE-30, KYSE-150, KYSE-70 and KYSE-790 cell lines were cultured in a 6-well plate and harvest for Western blot. Cells were lysed with RIPA buffer (150 mM NaCl, 50 mM Tris-Cl, 1 mM EGTA, 1% Triton X-100, 0.1% SDS and 1% sodium deoxycholate, pH 8.0) and then supplemented with proteinase inhibitor cocktail (Sigma-Aldrich, Burlington, MA, USA). After that, protein concentrations were quantified by using a BCA kit (Thermo, Waltham, MA, USA). Primary antibodies used in this study included anti-EGFR (1:1000, Thermo, USA), anti-ERBB2 (1:1000, Thermo, Waltham, MA, USA), anti-ERBB3 (1:1000), anti-ERBB4 (1:1000) and anti-β-actin (1:1000, Proteintech, Rosemont, IL, USA). Secondary antibodies included HRP-labeled goat anti-rabbit IgG (1:5000) and anti-mouse IgG (1:5000, Cell Signaling Technology, Danvers, MA, USA). ECL substrate reagent (Cytiva Amersham, Marlborough, MA, USA) was used to visualize signals on ChemiDoc (BioRad, Hercules, CA, USA).

4.5. Cell Viability Measurement

Cell viability was measured by using MTT assay. KYSE-150 cells were seeded in a 96-well plate at the number of 104 in each well. After attachment, the cells were treated with RP4010 and afatinib, combined or separately, at different concentrations, for 24 h. For RP4010, the concentration varied from 10 to 0.625 µM. The concentration of afatinib varied from 5 to 0.3125 µM. For the combination of RP4010 and afatinib, the ratio of RP4010 to afatinib was 2:1. After 24 h, KYSE-150 cells were incubated with medium containing 10% of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-tetrazolium bromide (MTT, 5 mg/mL), at 37 °C. After 4 h, formazan was dissolved in 150 µL DMSO. Absorbance was measured at 570 nm on SpectraMax i3 (Molecular Devices, San Jose, CA, USA). The survival curve was created by Graphpad Prism 5 (San Diego, CA, USA).

4.6. Mathematical Modeling

The mathematical model used was a modified version of one developed by Sneyd et al., with influence from Dupont et al. [6,27]. The governing dynamical equations for the Ca2+ signaling pathway are as follows:
d c d t = J I P R c , p , h J s e r c a + δ J i n J p m d c e d t = γ J s e r c a J I P R d h d t = h c h 1 τ h c d p d t = ν β p p .
where c denotes the Ca2+ concentration in the cytoplasm, [Ca2+ ]cyto; Ce is the Ca2+ concentration in the ER, [Ca2+ ]ER; p is the concentration of IP3, [IP3]; and h is the rate at which Ca2+ can activate IP3Rs.
The fluxes in the above equations were given as follows:
J I P 3 R = k f P O c e c J s e r c a = V p c 2 K ¯ c e 2 c 2 + k p 2 J i n = α 0 + α 1 K e 4 K e 4 + c e 4 J p m = V p m c 2 K p m 2 + c 2 .
where the J terms represent the Ca2+ fluxes across the cell and ER membranes. JIP3R modeled the flow of Ca2+ from the ER through the IP3R channel and depended upon both [Ca2+]cyto and [Ca2+]ER. As c decreased, the difference between c and Ce increased, which resulted in an increase for JIP3R. The parameter kf was a scaling factor used to control the IP3R density. PO represented the open probability of IP3R and depended on the activation and inactivation rates of IP3R, which were both affected by the binding of Ca2+ and IP3. The SERCA pump moved Ca2+ from the cytoplasm into the ER and changed depending on c and Ce. As Ce decreased with respect to c, the entire term increased. In order to model Ca2+ flow into the cell, Jin combined α0, which accounted for leaks into the cytoplasm through unspecified channels, and α 1 K e 4 K e 4 + c e 4 , which represented SOCE-mediated Ca2+ influx. These channels open and close in response to [Ca2+]ER, Ce. As Ce increased, the denominator in α 1 K e 4 K e 4 + c e 4 increased, causing the entire term to decrease, implying that the cell does not need to continue to fill up ER Ca2+ stores. The opposite applied if Ce decreased, meaning that the cell needs to increase ER Ca2+ stores. A Hill equation was used to model Jpm, the flow of Ca2+ through the plasma membrane pump, with Kpm being the concentration of Ca2+, where half of the binding sites at the pump were occupied, and Vpm is the maximum capacity of the plasma pump. The parameter γ denoted the ratio of the volume of the cytoplasm to the volume of the ER, and the parameter δ was dimensionless, as it was a scale factor relating the fluxes through the plasma membrane and ER membrane. To model IP3 concentration, βp represents the speed it takes p to decay to the steady state ps and ν = βp ps.
The model was further expanded upon with the introduction of the functions α and β:
P O = β β + k β β + α α = A p 1 m ¯ α c h ¯ α c β = B p m ¯ β c h c , t m ¯ α c = m ¯ β c = c 4 K c 4 + c 4 h ¯ α c = h c = K h 4 K h 4 + c 4 1 A p = B p = p 2 K p 2 + p 2 τ h = τ m a x K τ 4 K τ 4 + c 4 .
The function α(c,p) modeled the rate at which [Ca2+ ] inactivating IP3Rs and β(c,p,t) modeled the rate at which [Ca2+ ] activates IP3Rs.

4.7. Addition of TKI and SOCE Blocker in the Computational Model

A new feature in this computational model was the incorporation of the effects of one drug or of the combination of two different drugs acting on different Ca2+ resources. TKI (afatinib) and SOCE blocker (RP4010) were incorporated into the model through D1 and D2, which represent the overall drug effect on each of the channels:
J i n = α 0 + D 2 α 1 K e 4 K e 4 + c e 4 d p d t = D 1 ν β p p ,
As a SOCE channel inhibitor, RP4010 decreased [Ca2+]cyto and was expected to reduce the frequency of the oscillations. As a TKI, afatinib inhibits EGFR tyrosine kinase activity and thus affects the steady state of IP3. It was expected that it negatively affected the ER Ca2+ release through the IP3R channel and, again, reducing [Ca2+]cyto. In order to describe the relationship between the effectiveness of the drug to the dose of drug given to the cell, we used the following equations:
D 1 = 1 f 1 S 1 D 2 = 1 f 2 S 2 f 3 S 1
where the functions fi were described through the Michaelis–Menten kinetics:
f i S = c i   S / k i + S ,   i = 1 , 2 , 3 )
where S1 and S2 represent the doses of afatinib and RP4010, respectively.
All parameters used in the mathematical model were summarized in Table 3.

4.8. Validation of the Mathematical Model

We fitted the parameters of the mathematical model and validated it by using experimental data for different concentrations of drugs individually or in combination with a ratio of RP4010:afatinib = 2:1. The observed experimental data showed that [Ca2+] oscillates in ESCC cells had a period of 32.8 s. We used this value to adjust the model parameters kf, α1 and τmax, accordingly, to fit. In non-tumorous esophageal epithelial cells, which contain much less SOCE and ER Ca2+ release, we divided the values of k f by 5 and α1 by 3 to reflect such differences. The working assumption was that increasing concentrations of TKI or SOCE blockers would gradually decrease the level of intracellular Ca2+ oscillations in tumor cells until to the level in Het-1A cells. Finally, to find the coefficients ci and ki in the Michaelis–Menten equations, we matched the experimental data regarding drug doses and their corresponding periods.
To carry out the numerical simulations, the MATLAB’s ode23s built-in function, based on the explicit Runge–Kutta (2,3) pair of Bogacki and Shampine [28], was used to solve the systems of differential equations. The MATLAB’s lsqcurvefit built-in function was implemented by using the Levenberg–Marquardt algorithm [29] that solves a minimization problem involving a least-squares data-fitting term for fitting the model to the experimental data presented in Table 1.

4.9. Statistical Analysis

In this study, at least 50 cells in each concentration group were included to calculate the period (peak–peak duration) of Ca2+ oscillations. Cancer cells present different patterns of Ca2+ oscillations during cell-cycle progression, and we focused on the G(1)/S phase, a determining phase for cell proliferation that contains majority cells with Ca2+ oscillations. To simplify the model, we trimmed the dataset by removing any data outside 1.5 times the interquartile range. After trimming, the mean values of the period calculated from the rest of the dataset were used for the mathematical simulation, as presented in Table 1 and Table 2.
The combination index (CI) was calculated from the CI equation algorithms by using CompuSyn software (CompuSyn, Inc., Paramus, NJ, USA). Dose–effect Curve and combination index (CI) versus fraction affected (fa) plots were analyzed and exported from CompuSyn [30].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms23031763/s1.

Author Contributions

Conceptualization, Y.C., H.V.K., B.C. and Z.P.; methodology, Y.C., S.C. and Z.P.; software, M.F., E.S., H.V.K. and B.C.; validation, Y.C., M.F., E.S., H.V.K., B.C. and Z.P.; formal analysis, M.F., E.S., H.V.K. and B.C.; investigation, Y.C., M.F., E.S., H.V.K., B.C. and Z.P.; resources, Y.C., M.F., E.S., H.V.K., B.C. and Z.P.; data curation, Y.C., M.F., E.S., H.V.K., B.C. and Z.P.; writing—original draft preparation, Y.C., M.F., E.S., S.R. and Z.P.; writing—review and editing, all authors; visualization, B.C., H.V.K. and Z.P.; supervision, B.C., H.V.K. and Z.P.; project administration, B.C., H.V.K. and Z.P.; funding acquisition, B.C., H.V.K. and Z.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from National Institutes of Health (NIH) (R01 CA185055, S10OD0252300) and The University of Texas System STARs Award (to Z.P.), The University of Texas at Arlington Interdisciplinary Research Program (to B.C., H.V.K. and Z.P.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

This work was partially supported by grants from National Institutes of Health (NIH) (R01 CA185055, S10OD0252300) and The University of Texas System STARs Award (to Z.P.), The University of Texas at Arlington Interdisciplinary Research Program (to B.C., H.V.K. and Z.P.).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Le Bras, G.F.; Farooq, M.H.; Falk, G.W.; Andl, C.D. Esophageal cancer: The latest on chemoprevention and state of the art therapies. Pharmacol. Res. 2016, 113, 236–244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zhang, X.; Pan, Z. Influence of microbiota on immunity and immunotherapy for gastric and esophageal cancers. Gastroenterol. Rep. 2020, 8, 206–214. [Google Scholar] [CrossRef]
  3. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer Statistics, 2017. CA Cancer J. Clin. 2017, 67, 7–30. [Google Scholar] [CrossRef] [Green Version]
  4. Cui, C.; Merritt, R.; Fu, L.; Pan, Z. Targeting calcium signaling in cancer therapy. Acta Pharm. Sin. B 2017, 7, 3–17. [Google Scholar] [CrossRef] [PubMed]
  5. Clapham, D.E. Calcium signaling. Cell 2007, 131, 1047–1058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sneyd, J.; Han, J.M.; Wang, L.; Chen, J.; Yang, X.; Tanimura, A.; Sanderson, M.J.; Kirk, V.; Yule, D.I. On the dynamical structure of calcium oscillations. Proc. Natl. Acad. Sci. USA 2017, 114, 1456–1461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Berridge, M.J.; Lipp, P.; Bootman, M.D. The versatility and universality of calcium signalling. Nat. Rev. Mol. Cell Biol. 2000, 1, 11–21. [Google Scholar] [CrossRef] [PubMed]
  8. Bergmeier, W.; Weidinger, C.; Zee, I.; Feske, S. Emerging roles of store-operated Ca(2)(+) entry through STIM and ORAI proteins in immunity, hemostasis and cancer. Channels 2013, 7, 379–391. [Google Scholar] [CrossRef] [Green Version]
  9. Pan, Z.; Ma, J. Open Sesame: Treasure in store-operated calcium entry pathway for cancer therapy. Sci. China Life Sci. 2015, 58, 48–53. [Google Scholar] [CrossRef] [Green Version]
  10. Yang, S.; Zhang, J.J.; Huang, X.Y. Orai1 and STIM1 are critical for breast tumor cell migration and metastasis. Cancer Cell 2009, 15, 124–134. [Google Scholar] [CrossRef] [Green Version]
  11. Kondratska, K.; Kondratskyi, A.; Yassine, M.; Lemonnier, L.; Lepage, G.; Morabito, A.; Skryma, R.; Prevarskaya, N. Orai1 and STIM1 mediate SOCE and contribute to apoptotic resistance of pancreatic adenocarcinoma. Biochim. Biophys. Acta 2014, 1843, 2263–2269. [Google Scholar] [CrossRef] [Green Version]
  12. Flourakis, M.; Lehen’kyi, V.; Beck, B.; Raphael, M.; Vandenberghe, M.; Abeele, F.V.; Roudbaraki, M.; Lepage, G.; Mauroy, B.; Romanin, C.; et al. Orai1 contributes to the establishment of an apoptosis-resistant phenotype in prostate cancer cells. Cell Death Dis. 2010, 1, e75. [Google Scholar] [CrossRef] [Green Version]
  13. Zhu, H.; Zhang, H.; Jin, F.; Fang, M.; Huang, M.; Yang, C.S.; Chen, T.; Fu, L.; Pan, Z. Elevated Orai1 expression mediates tumor-promoting intracellular Ca2+ oscillations in human esophageal squamous cell carcinoma. Oncotarget 2014, 5, 3455–3471. [Google Scholar] [CrossRef] [Green Version]
  14. Cui, C.; Chang, Y.; Zhang, X.; Choi, S.; Tran, H.; Penmetsa, K.V.; Viswanadha, S.; Fu, L.; Pan, Z. Targeting Orai1-mediated store-operated calcium entry by RP4010 for anti-tumor activity in esophagus squamous cell carcinoma. Cancer Lett. 2018, 432, 169–179. [Google Scholar] [CrossRef]
  15. Bryant, J.A.; Finn, R.S.; Slamon, D.J.; Cloughesy, T.F.; Charles, A.C. EGF activates intracellular and intercellular calcium signaling by distinct pathways in tumor cells. Cancer Biol. Ther. 2004, 3, 1243–1249. [Google Scholar] [CrossRef] [Green Version]
  16. Meisenhelder, J.; Suh, P.G.; Rhee, S.G.; Hunter, T. Phospholipase C-gamma is a substrate for the PDGF and EGF receptor protein-tyrosine kinases in vivo and in vitro. Cell 1989, 57, 1109–1122. [Google Scholar] [CrossRef]
  17. Pottier, C.; Fresnais, M.; Gilon, M.; Jerusalem, G.; Longuespee, R.; Sounni, N.E. Tyrosine Kinase Inhibitors in Cancer: Breakthrough and Challenges of Targeted Therapy. Cancers 2020, 12, 731. [Google Scholar] [CrossRef] [Green Version]
  18. Minkovsky, N.; Berezov, A. BIBW-2992, a dual receptor tyrosine kinase inhibitor for the treatment of solid tumors. Curr. Opin. Investig. Drugs 2008, 9, 1336–1346. [Google Scholar]
  19. Hong, M.H.; Heo, S.G.; Lee, Y.G.; Kim, H.S.; Park, K.U.; Kim, H.G.; Ko, Y.H.; Chung, I.J.; Min, Y.J.; Kim, M.K.; et al. Phase 2 study of afatinib among patients with recurrent and/or metastatic esophageal squamous cell carcinoma. Cancer 2020, 126, 4521–4531. [Google Scholar] [CrossRef]
  20. Ma, P.; Fu, Y.; Chen, M.; Jing, Y.; Wu, J.; Li, K.; Shen, Y.; Gao, J.X.; Wang, M.; Zhao, X.; et al. Adaptive and Acquired Resistance to EGFR Inhibitors Converge on the MAPK Pathway. Theranostics 2016, 6, 1232–1243. [Google Scholar] [CrossRef] [Green Version]
  21. Roy, S.; Pan, Z.; Pal, S. A Fokker-Planck feedback control framework for optimal personalized therapies in colon cancer-induced angiogenesis. arXiv 2021, arXiv:2109.00166. [Google Scholar]
  22. Roy, S.; Pal, S.; Manoj, A.; Kakarla, S.; Villegas, J.; Alajmi, M. A Fokker-Planck framework for parameter estimation and sen-sitivity analysis in colon cancer. AIP Conf. Proc. 2021. to appear. [Google Scholar]
  23. Chang, Y.; Roy, S.; Pan, Z. Store-Operated Calcium Channels as Drug Target in Gastroesophageal Cancers. Front. Pharmacol. 2021, 12, 668730. [Google Scholar] [CrossRef] [PubMed]
  24. Kim, M.S.; Kim, S.H.; Yang, S.H. Restricting extracellular Ca2+ on gefitinib-resistant non-small cell lung cancer cells reverses altered epidermal growth factor-mediated Ca2+ response, which consequently enhances gefitinib sensitivity. PLoS ONE 2020, 15, e0238155. [Google Scholar] [CrossRef] [PubMed]
  25. Casas-Rua, V.; Tomas-Martin, P.; Lopez-Guerrero, A.M.; Alvarez, I.S.; Pozo-Guisado, E.; Martin-Romero, F.J. STIM1 phosphorylation triggered by epidermal growth factor mediates cell migration. Biochim. Biophys. Acta 2015, 1853, 233–243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Smedler, E.; Uhlén, P. Frequency decoding of calcium oscillations. Biochim. Biophys. Acta 2014, 1840, 964–969. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Dupont, G.; Combettes, L.; Bird, G.S.; Putney, J.W. Calcium Oscillations. Cold Spring Harb. Perspect. Biol. 2011, 3, a004226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Bogacki, P.; Shampine, L.F. A 3(2) pair of Runge-Kutta formulas. Appl. Math. Lett. 1989, 2, 321–325. [Google Scholar] [CrossRef] [Green Version]
  29. Marquardt, D. An Algorithm for Least-squares Estimation of Nonlinear Parameters. SIAM J. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  30. Chou, T.C. Drug combination studies and their synergy quantification using the Chou-Talalay method. Cancer Res. 2010, 70, 440–446. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Upregulated expression of EGFR and ERBB2 and inhibited Ca2+ oscillations by afatinib in ESCC cells. (a) RNA-Seq data (TCGA) revealed significant higher expression of EGFR in tumor vs. normal tissues removed from patients. Average expression level (TPM) for tumor and non-tumor is 38.4 and 20.4, respectively. Tumor (n = 182), including ESCC (n = 94), EAC (n = 87) and cystic, mucinous and serous neoplasms (n = 1). Normal (n = 286), including TCGA adjacent normal (n = 13) and GTEx normal esophageal tissue (n = 273); * p-value < 0.0001. (b) Western Blot results of EGFR family members (EGFR, ERBB2, ERBB3 and ERBB4) in ESCC and Het-1A cells. β-actin was used as loading control. (c) Quantification of Western blot result of EGFR family members in ESCC and Het-1A cells; * p-value < 0.05. (d) Time-lapse live-cell fluorescent imaging of intracellular Ca2+ in KYSE-150 cells. Cells were treated with either vehicle control or 5 μM afatinib for 4 h. The intracellular Ca2+ concentrations were presented by a heap map with fluorescent intensity changes (F/F0). Some representative cells with obvious Ca2+ oscillations were indicated by dashed white circle. Scale bar, 20 μm.
Figure 1. Upregulated expression of EGFR and ERBB2 and inhibited Ca2+ oscillations by afatinib in ESCC cells. (a) RNA-Seq data (TCGA) revealed significant higher expression of EGFR in tumor vs. normal tissues removed from patients. Average expression level (TPM) for tumor and non-tumor is 38.4 and 20.4, respectively. Tumor (n = 182), including ESCC (n = 94), EAC (n = 87) and cystic, mucinous and serous neoplasms (n = 1). Normal (n = 286), including TCGA adjacent normal (n = 13) and GTEx normal esophageal tissue (n = 273); * p-value < 0.0001. (b) Western Blot results of EGFR family members (EGFR, ERBB2, ERBB3 and ERBB4) in ESCC and Het-1A cells. β-actin was used as loading control. (c) Quantification of Western blot result of EGFR family members in ESCC and Het-1A cells; * p-value < 0.05. (d) Time-lapse live-cell fluorescent imaging of intracellular Ca2+ in KYSE-150 cells. Cells were treated with either vehicle control or 5 μM afatinib for 4 h. The intracellular Ca2+ concentrations were presented by a heap map with fluorescent intensity changes (F/F0). Some representative cells with obvious Ca2+ oscillations were indicated by dashed white circle. Scale bar, 20 μm.
Ijms 23 01763 g001
Figure 2. (a) Decreased intracellular Ca2+ oscillations in KYSE-150 cells treated with RP4010 or afatinib. (bd) Fluorescent image of KYSE-150 cells loaded with Fluo-4 AM. Scale bar, 50 μm. (b) Representative traces of Ca2+ oscillations in a single KYSE-150 cell treated with vehicle control, (c) RP4010 or (d) afatinib.
Figure 2. (a) Decreased intracellular Ca2+ oscillations in KYSE-150 cells treated with RP4010 or afatinib. (bd) Fluorescent image of KYSE-150 cells loaded with Fluo-4 AM. Scale bar, 50 μm. (b) Representative traces of Ca2+ oscillations in a single KYSE-150 cell treated with vehicle control, (c) RP4010 or (d) afatinib.
Ijms 23 01763 g002
Figure 3. Ca2+ oscillations simulated by mathematical modeling. (a) Simulated Ca2+ oscillations in KYSE-150 (green line) and Het-1A cells (blue line). (b,c) Simulated Ca2+ oscillations in KYSE-150 cells treated with RP4010 at the concentration of 0.5 µM (b) or 3.0 µM (c). (d) Simulated Ca2+ dynamics in KYSE-150 cells treated with afatinib at the concentration of 0.5 µM.
Figure 3. Ca2+ oscillations simulated by mathematical modeling. (a) Simulated Ca2+ oscillations in KYSE-150 (green line) and Het-1A cells (blue line). (b,c) Simulated Ca2+ oscillations in KYSE-150 cells treated with RP4010 at the concentration of 0.5 µM (b) or 3.0 µM (c). (d) Simulated Ca2+ dynamics in KYSE-150 cells treated with afatinib at the concentration of 0.5 µM.
Ijms 23 01763 g003
Figure 4. Effect of combined treatment of RP4010 and afatinib on intracellular Ca2+ oscillations. (a) Mathematical simulation of Ca2+ dynamics in KYSE-150 cells treated with combined RP4010 and afatinib (2:1 ration at 3.0 µM). (b) A representative trace of intracellular Ca2+ concentration measured by time-lapse live-cell imaging in KYSE150 cells treated with combined 2 µM RP4010 and 1 µM afatinib. (c) Dose-dependent curves of simulated periods of intracellular Ca2+ oscillations in KYSE-150 cells treated with either RP4010 (blue), afatinib (red) or combination (green). Individual circle, square or triangle is from experimental data.
Figure 4. Effect of combined treatment of RP4010 and afatinib on intracellular Ca2+ oscillations. (a) Mathematical simulation of Ca2+ dynamics in KYSE-150 cells treated with combined RP4010 and afatinib (2:1 ration at 3.0 µM). (b) A representative trace of intracellular Ca2+ concentration measured by time-lapse live-cell imaging in KYSE150 cells treated with combined 2 µM RP4010 and 1 µM afatinib. (c) Dose-dependent curves of simulated periods of intracellular Ca2+ oscillations in KYSE-150 cells treated with either RP4010 (blue), afatinib (red) or combination (green). Individual circle, square or triangle is from experimental data.
Ijms 23 01763 g004
Figure 5. Synergetic effect of combined RP4010 and afatinib on cell viability in KYSE-150 cells. (a) Nuclei staining by Hoechst 33342 in KYSE-150 cells treated with vehicle control, 1 µM RP4010, 0.5 µM afatinib or 1 µM RP4010 + 0.5 µM afatinib. Apoptotic cells were indicated by nuclear DNA fragmentation and condensation (white arrow heads). Scale bar, 10 μm. (b) Dose-dependent inhibition curves of KYSE-150 cells treated by RP4010 (blue circles), afatinib (red squares) or RP4010 plus afatinib (ratio 2:1, green triangles). Cell inhibition rates were measured by MTT assay. Details in the confined box were amplified as the upper curve. (c) Combination index (CI) plot for combined effect of RP4010 and afatinib. Each blue circle indicated a data point with combined RP4010 and afatinib. Simulation curve is shown in blue line. Average effect values (Fa) are shown on x-axis.
Figure 5. Synergetic effect of combined RP4010 and afatinib on cell viability in KYSE-150 cells. (a) Nuclei staining by Hoechst 33342 in KYSE-150 cells treated with vehicle control, 1 µM RP4010, 0.5 µM afatinib or 1 µM RP4010 + 0.5 µM afatinib. Apoptotic cells were indicated by nuclear DNA fragmentation and condensation (white arrow heads). Scale bar, 10 μm. (b) Dose-dependent inhibition curves of KYSE-150 cells treated by RP4010 (blue circles), afatinib (red squares) or RP4010 plus afatinib (ratio 2:1, green triangles). Cell inhibition rates were measured by MTT assay. Details in the confined box were amplified as the upper curve. (c) Combination index (CI) plot for combined effect of RP4010 and afatinib. Each blue circle indicated a data point with combined RP4010 and afatinib. Simulation curve is shown in blue line. Average effect values (Fa) are shown on x-axis.
Ijms 23 01763 g005
Figure 6. Illustration of intracellular Ca2+ oscillations as the intersection between EGFR and SOCE signaling pathways. [Ca2+]cyto are dependent on EGFR/PLCγ/IP3/IP3R-mediated ER Ca2+ release and Orai1-mediated SOCE. Intracellular Ca2+ oscillations are code, especially the frequency, to control cell proliferation in ESCC cells. Inhibiting either EGFR by Afatinib or Orai1 by RP4010 could reduce the frequency of intracellular Ca2+ oscillations and inhibit cancer-cell proliferations.
Figure 6. Illustration of intracellular Ca2+ oscillations as the intersection between EGFR and SOCE signaling pathways. [Ca2+]cyto are dependent on EGFR/PLCγ/IP3/IP3R-mediated ER Ca2+ release and Orai1-mediated SOCE. Intracellular Ca2+ oscillations are code, especially the frequency, to control cell proliferation in ESCC cells. Inhibiting either EGFR by Afatinib or Orai1 by RP4010 could reduce the frequency of intracellular Ca2+ oscillations and inhibit cancer-cell proliferations.
Ijms 23 01763 g006
Table 1. Dose-dependent drug effects on intracellular Ca2+ oscillations (experimental data).
Table 1. Dose-dependent drug effects on intracellular Ca2+ oscillations (experimental data).
Drug Concentration (µM)Period (s)
Mean ± SEM
RP4010Afatinib
10122 ± 6.6/
5111 ± 6.273.4 ± 4.3
2.560.8 ± 4.061.9 ± 3.5
256.5 ± 2.7/
1.2555.4 ± 3.554.2 ± 2.5
153.5 ± 2.952.7 ± 2
0.62544.1 ± 3.348.2 ± 2.4
0.538.8 ± 1.347.1 ± 1.5
0.3125/41.9 ± 2.6
0.2535.9 ± 1.341.8 ± 1.4
0.12534.4 ± 239.8 ± 1.4
0.0625/38.6 ± 2.2
0.032.8 ± 1.2
Table 2. Parameter values for the mathematical model.
Table 2. Parameter values for the mathematical model.
ParameterDescriptionValueUnits
δUsed to adjust ratio of [Ca2+] across plasma membrane to ER membrane1.5n/a
The concentration of Ca2+ in response to β0.1µM
KcHalf-maximal [Ca2+] for IP3R0.2µM
k fScaling factor that controls [Ca2+] release through IP3R; IP3R density and channel activity3.9s−1
VsercaMaximum capacity of SERCA pump0.9µMs−1
γRatio of cytoplasmic volume to ER volume5.5n/a
VpmMaximum capacity of plasma membrane pump0.11µMs−1
α0Flow of calcium into the cell through an unspecified leak0.0027µMs−1
KPHalf-maximal [IP3] for IP3R0.2µM
τmaxControls rate that β responds to [Ca2+] changes1420s−1
βpRate of decay of p to its steady state0.027s−1
KhThe concentration of Ca2+ activated IP3R0.08µM
Used to adjust the Ca2+ concentration in ER1.9 × 10−5n/a
KeHalf-maximal [Ca2+] ER for SOC channels8µM
KpmHalf-maximal [Ca2+] for plasma membrane pump0.3µM
α1Rate constant for SOC channels0.385µMs−1
Table 3. Summary of sensitivity analysis results.
Table 3. Summary of sensitivity analysis results.
Dose (µM)ParameterBaselinep-ValuePRCC Value
Aft: 0.5
RP: 1.0
kf3.90.6430.034
α10.3851.44 × 10−50−0.831
c10.510.2720.079
c211.74 × 10−460.812
c30.764.6 × 10−16−0.541
k10.340.061−0.1353
k20.64.01 × 10−12−0.4724
k30.88.45 × 10−80.3740
Table 4. Effects on intracellular Ca2+ oscillations by combined RP4010 and afatinib (experimental data).
Table 4. Effects on intracellular Ca2+ oscillations by combined RP4010 and afatinib (experimental data).
Concentration of Combined Drugs (µM)Period (s)
RP4010Afatinib
52.592.3 ± 4.6
2.51.2590.8 ± 4.4
2189.4 ± 6.4
1.50.7576.8 ± 4.6
1.250.62573.1 ± 3.8
10.568.7 ± 3.9
0.6250.312564.2 ± 3.7
0.250.12561.1 ± 3.1
0.1250.062559.6 ± 2.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chang, Y.; Funk, M.; Roy, S.; Stephenson, E.; Choi, S.; Kojouharov, H.V.; Chen, B.; Pan, Z. Developing a Mathematical Model of Intracellular Calcium Dynamics for Evaluating Combined Anticancer Effects of Afatinib and RP4010 in Esophageal Cancer. Int. J. Mol. Sci. 2022, 23, 1763. https://doi.org/10.3390/ijms23031763

AMA Style

Chang Y, Funk M, Roy S, Stephenson E, Choi S, Kojouharov HV, Chen B, Pan Z. Developing a Mathematical Model of Intracellular Calcium Dynamics for Evaluating Combined Anticancer Effects of Afatinib and RP4010 in Esophageal Cancer. International Journal of Molecular Sciences. 2022; 23(3):1763. https://doi.org/10.3390/ijms23031763

Chicago/Turabian Style

Chang, Yan, Marah Funk, Souvik Roy, Elizabeth Stephenson, Sangyong Choi, Hristo V. Kojouharov, Benito Chen, and Zui Pan. 2022. "Developing a Mathematical Model of Intracellular Calcium Dynamics for Evaluating Combined Anticancer Effects of Afatinib and RP4010 in Esophageal Cancer" International Journal of Molecular Sciences 23, no. 3: 1763. https://doi.org/10.3390/ijms23031763

APA Style

Chang, Y., Funk, M., Roy, S., Stephenson, E., Choi, S., Kojouharov, H. V., Chen, B., & Pan, Z. (2022). Developing a Mathematical Model of Intracellular Calcium Dynamics for Evaluating Combined Anticancer Effects of Afatinib and RP4010 in Esophageal Cancer. International Journal of Molecular Sciences, 23(3), 1763. https://doi.org/10.3390/ijms23031763

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

Article Metrics

Back to TopTop