Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 23 December 2015
Sec. Membrane Physiology and Membrane Biophysics

Simulation of Cl− Secretion in Epithelial Tissues: New Methodology Estimating Activity of Electro-Neutral Cl− Transporter

  • 1Department of Molecular Cell Physiology, Graduate School of Medical Science, Kyoto Prefectural University of Medicine, Kyoto, Japan
  • 2Department of Health and Sports Sciences, Faculty of Health and Medical Sciences, Kyoto Gakuen University, Kameoka, Japan
  • 3Japan Institute for Food Education and Health, St. Agnes' University, Kyoto, Japan
  • 4Department of Bio-Ionomics, Graduate School of Medical Science, Kyoto Prefectural University of Medicine, Kyoto, Japan

Transcellular Cl secretion is, in general, mediated by two steps; (1) the entry step of Cl into the cytosolic space from the basolateral space across the basolateral membrane by Cl transporters, such as Na+-K+-2Cl cotransporter (NKCC1, an isoform of NKCC), and (2) the releasing step of Cl from the cytosolic space into the luminal (air) space across the apical membrane via Cl channels, such as cystic fibrosis transmembrane conductance regulator (CFTR) Cl channel. Transcellular Cl secretion has been characterized by using various experimental techniques. For example, measurements of short-circuit currents in the Ussing chamber and patch clamp techniques provide us information on transepithelial ion movements via transcellular pathway, transepithelial conductance, activity (open probability) of single channel, and whole cell currents. Although many investigators have tried to clarify roles of Cl channels and transporters located at the apical and basolateral membranes in transcellular Cl secretion, it is still unclear how Cl channels/transporters contribute to transcellular Cl secretion and are regulated by various stimuli such as Ca2+ and cAMP. In the present study, we simulate transcellular Cl secretion using mathematical models combined with electrophysiological measurements, providing information on contribution of Cl channels/transporters to transcellular Cl secretion, activity of electro-neutral ion transporters and how Cl channels/transporters are regulated.

Introduction

Goblet cells located in airway surface epithelia and mucous cells of submucosal glands produce mucins, which are secreted into the airway space (Kim, 2012). Mucus layers formed by this secreted mucins covering airway epithelial surfaces trap pathogens such as bacteria and viruses in the mucus layers, which are removed by mucociliary clearance (Kim, 2012; Komatani-Tamiya et al., 2012). Transcellular Cl secretion in airway epithelial cells followed by paracellular Na+ transport (secretion) from the basolateral space to the apical space increases osmotic pressure in the apical space, resulting in water secretion (Asano et al., 2009; Kim et al., 2014; Marunaka, 2014a; Stanton et al., 2015). Diminution of water secretion driven by NaCl secretion elevates viscosity of mucins, causing dysfunction of mucociliary clearance (Kunzelmann and Schreiber, 2012). Thus, transcellular Cl secretion plays a crucial role in prevention from infection of bacteria and viruses, etc., by producing water secretion into the apical space, covering the apical surface of airway epithelial cells.

Transcellular Cl secretion in epithelial cells is generally mediated by two steps; (1) the entry step of Cl across the basolateral membrane by Cl transporters such as Na+-K+-2Cl cotransporter (NKCC1, an isoform of NKCC), and (2) the releasing step of Cl across the apical membrane via Cl channels such as the cystic fibrosis transmembrane conductance regulator (CFTR) Cl channel (Lee and Foskett, 2010; Li et al., 2012; Stölting et al., 2014; Sun et al., 2014a,b; Stanton et al., 2015). Transcellular Cl secretion across the epithelial tissue has been characterized by using the technique measuring short-circuit currents. Using blockers of Cl channels/transporters combined with the short-circuit current measurement technique, we can determine amounts of transepithelial Cl secretion (Marunaka, 2014a; Sun et al., 2014a,b). However, using this technique, we could not exactly determine how transcellular Cl secretion varies in magnitude and time due to modification of activity of Cl channels and transporters located at the apical and basolateral membranes. In the present study, we tried to clarify how transcellular Cl secretion varies in magnitude and time when the activity of Cl channels and transporters located at the apical and basolateral membranes changes using mathematical simulation with three parameters: (1) the entry step of Cl into the intracellular space from the basolateral space across the basolateral membrane by Cl transporters, (2) the releasing step of Cl from the intracellular space into the apical space across the apical membrane via Cl channels, and (3) the releasing step of Cl from the intracellular space to the basolateral space across the basolateral membrane via Cl channels (recycle, back flow/flux). This simulation method reported in the present study shows us how the activity of Cl channels and transporters located at the apical and basolateral membranes contributes to the transcellular Cl secretion, and the determination of the activity of electro-neutral Cl transporter such as NKCC1. Data have been partly reported in abstract form (Sasamoto et al., 2014).

Methods

Chemicals and Materials

We obtained forskolin, daidzein, apigenin, genistein, protein kinase inhibitor 14–22 amide (PKI), NPPB (5-nitro-2-(3-phenylpropylamino)benzoic acid), nystatin, benzamil, and dimethyl sulfoxide (DMSO) from Sigma-Aldrich (St Louis, MO, USA), and epithelial A6 cells from American Type Culture Collection (ATCC). Forskolin (10 μM), daidzein (100 μM), apigenin (100 μM), genistein (100 μM), PKI (2 μM), NPPB (100 μM), nystatin (50 μM), and benzamil (10 μM) dissolved in DMSO were applied to the solution as the final concentration. The concentrations of forskolin, daidzein, apigenin, genistein, and PKI used in the present study were determined from the observations obtained in previous reports (Niisato et al., 1999), and these compounds were applied into both apical and basolateral solutions.

Cell Culture

A6 cells were derived from the kidney of Xenopus laevis, an amphibian, and cultured (passage 73–84) similar to our previous reports (Niisato and Marunaka, 1999; Sun et al., 2013, 2014a,b). To measure apical and basolateral Cl conductance (GA and GB) and short-circuit current (ISC), we cultured A6 cells for 14 days on 6.5 mm Transwell-Clear permeable supports (0.33 cm2): GA and GB were actually measured in the area of 0.33 cm2. Volume of individual A6 cell was approximately 3.4 × 10−15 m3, and total volume of A6 cells cultured on Transwell-Clear permeable supports (0.33 cm2) was approximately 5.0 × 10−10 m3. The lateral membrane of A6 cells made tight junction expressing claudin-1, the width of which was less than 3 nm, showing the width of the paracellular space was less than 3 nm (Tokuda et al., 2008; Suzuki et al., 2009).

Measurement of Cl Conductance of Apical and Basolateral Membranes (GA and GB)

We transferred monolayers of A6 cells subcultured on tissue culture-treated Transwell filter cups to a modified Ussing chamber (Jim's Instrument, Iowa City, IA, USA) designed to hold the filter cup similar to another type of Ussing chamber (Marques et al., 2013), and continuously measured transepithelial potential difference (PD) by a high-impedance millivoltmeter (VCC-600, Physiologic Instrument, San Diego, CA, USA; Niisato and Marunaka, 1999; Sun et al., 2014b). We applied a pulse of +1 μA constant current every 10 s for 0.5 s to A6 monolayers under open-circuit conditions from the basolateral to the apical space, and calculated the conductance (G) from the ΔPD caused by the 1 μA constant-current pulse using Ohm's law (G = 1 μA/ΔPD mV: ΔPD had a positive value). We applied 100 μM NPPB (a non-selective Cl channel blocker, Niisato and Marunaka, 1999) to the apical or basolateral solution for detection of the NPPB-sensitive conductance of the apical or basolateral membrane used as the Cl conductance. We measured an NPPB-sensitive conductance (Niisato and Marunaka, 1999; Tokuda et al., 2007, 2008, 2009a,b, 2010). To detect the Cl conductance of the apical membrane, we measured the NPPB-sensitive conductance by applying 100 μM NPPB to the apical solution 45 min after addition of 50 μM nystatin to the basolateral membrane that permeabilized the basolateral membrane. To detect the Cl conductance of the basolateral membrane, we measured the NPPB-sensitive conductance by applying 100 μM NPPB to the basolateral solution 45 min after addition of 50 μM nystatin to the apical membrane that permeabilized the apical membrane (Niisato and Marunaka, 1999). We applied forskolin of 10 μM to both apical and basolateral spaces 30 min before addition of nystatin to detect the forskolin action on the membrane conductance. The values of measured PD were within the range of −1 to −30 mV depending on the experimental conditions.

Measurement of Short-Circuit Current (Isc)

We measured Isc in A6 cells according to the method as previously reported (Niisato and Marunaka, 1999; Sun et al., 2014b). To detect transepithelial Cl movements, we applied benzamil of 10 μM into the apical solution to block epithelial Na+ channel (ENaC) contributing to transepithelial Na+ absorption (Niisato and Marunaka, 1999). The value of G·PD (ISC) had a negative value, since PD had a negative value. However, most recently published articles present the Isc showing Cl secretion and also Na+ absorption as a positive current (Ikehara et al., 2014; Marunaka, 2014a,b). Therefore, in the present study, we present the Isc (Cl secretion) as a positive current.

Solutions

The solution used in the present study contained (in mM) 120 NaCl, 3.5 KCl, 1 CaCl2, 1 MgCl2, 5 glucose, 10 HEPES with pH 7.4, since A6 cells used in the present study was derived from the kidney of an amphibian.

Temperature

All experiments shown in the present study were performed at 24–25°C, which are physiological temperatures for amphibian, since A6 cells are amphibian cells.

Data Presentation

Values of ISC and conductance are shown as the mean ± SEM. n means the number of experiments performed in the present study.

Results

Many compounds show various time-dependent patterns in stimulation of transcellular Cl secretion in epithelial cells (Niisato et al., 1999; Hennig et al., 2008; Ao et al., 2013; Luo et al., 2013). Transcellular Cl secretion in epithelial cells is mediated by uptake and release of Cl into and from the intracellular space. To clarify the mechanism on uptake and release of Cl regulated by various types of compounds influencing transcellular Cl secretion, we propose a model of epithelial Cl secretion via the transcellular pathway by comparing this proposed model with experimental data on transcellular Cl secretion measured as ISC in epithelial A6 cells.

Model of Transcellular Cl Secretion in Epithelial Cells

The parameters used in the present study are listed in Table 1. Figure 1 describes a model of transcellular Cl secretion in epithelial tissues. This model contains three Cl moving pathways between the intracellular and extracellular spaces across the apical and basolateral membranes: (1) a Cl releasing pathway from the intracellular space into the apical space, such as Cl channels, across the apical membrane (Pathway A contributing to Cl secretion as a passive Cl moving pathway driven by electrochemical potential of Cl between the intracellular and apical spaces across the apical membrane); (2) a Cl releasing pathway from the intracellular space into the basolateral space, such as Cl channels, across the basolateral membrane (Pathway B not contributing to Cl secretion as a passive Cl moving pathway driven by electrochemical potential of Cl between the intracellular and basolateral spaces across the basolateral membrane); (3) a Cl uptake pathway from the basolateral space into the intracellular space, such as NKCC1, across the basolateral membrane (Pathway C partially, but not all, contributing to Cl secretion as an active Cl moving pathway, such as NKCC1, driven by electrochemical potential of Na+ between the intracellular and basolateral spaces across the basolateral membrane). The transcellular Cl secretion consists of the following pathways: (1) Cl is first taken up into the intracellular space via Pathway C; (2) Cl taken up into the intracellular space by Pathway C is respectively released into the apical and basolateral spaces via Pathways A and B; Cl taken up by Pathway C released into the apical space via Pathway A only contributes to the transcellular Cl secretion.

TABLE 1
www.frontiersin.org

Table 1. Definition of characters.

FIGURE 1
www.frontiersin.org

Figure 1. A model of transcellular Cl secretion in epithelial tissues. F(t) is the intracellular Cl concentration ([Cl-]i) at time = t. JA is the Cl flux through Cl releasing pathways (Cl channels) across the apical membrane, JB is the Cl flux through Cl releasing pathways (Cl channels) across the basolateral membrane, and JC is the Cl flux through Cl transporter uptaking Cl across the basolateral membrane. GA and GB are respectively the apical and basolateral Cl conductance.

We can, in general, estimate the transcellular Cl secretion measuring short-circuit currents in Ussing chamber (Ussing and Zerahn, 1951). Therefore, we consider a model of transcellular Cl secretion measured under a short-circuit (voltage clamp) condition with the apical membrane potential identical to the basolateral membrane potential and all ionic compositions of the apical solution identical to those of basolateral solution.

Based on characteristics of Cl movements described in a model (Figure 1), we defined Cl effluxes as positive values. To clearly show results, we describe how to obtain equations applied in the present study in Appendix: the equation's number shown in text compatible with that in Appendix. We define the intracellular Cl concentration ([Cl-]i) at time = t as f(t) [see Equation (A13) in Appendix], where t is time after application of compounds affecting the Cl movements across the apical and/or basolateral membranes of epithelial cells. We show ISC (transepithelial Cl secretion) at time = t as IA(t) [see Equation (A18) in Appendix]. Table 1 shows meanings of parameters used in the present study.

A case where Cl uptake via a pathway (JC) is constant (JC=JC0: CTAct=CTAct0) and the membrane potential (V) is not changed (V = V0) after application of compounds affecting intracellular environments modifying apical and/or basolateral Cl conductance

These calculations indicate that: (1) an instantaneous Cl secretion (IA(0)) just after application of compounds affecting just GA and/or GB depends on a change in GA but not GB; (2) a Cl secretion at a steady state, IA(∞), depends upon both changes in GA and GB; (3) a transition Cl secretion occurs in response to changes in GA and GB time dependently with a time constant (τ) shown in Equation (A14). We further show some concrete examples of Cl secretion responses to changes in GA and/or GB.

A case under a condition with a relative change in an apical Cl conductance larger than that in the basolateral Cl conductance (GAGA0>GBGB0) without any change in Cl uptake across the basolateral membrane (JC=JC0:CTAct=CTAct0) or membrane potential (V = V0)

A case under a condition with only a change in the apical Cl conductance (GA=20GA0):GA=20GA0, GB=GB0, GB0=30GA0, GAGA0>GBGB0, JC=JC0 (CTAct=CTAct0), V= V0. Figure 2A shows IA (Cl secretion) in A6 cells treated with 10 μM forskolin in the presence of 2 μM protein kinase inhibitor 14–22 amide (PKI). Forskolin induced a transient increase in IA followed by a declining phase, reaching a steady level (Figure 2A). Figure 2B shows simulation of IA shown in Figure 2A: IA increases in response to a 20-fold increase in apical Cl conductance (GA=20 GA0) without any change in basolateral Cl conductance (GB= GB0). The instantaneous Cl secretion (IA (0)) just after a change in GA without any change in GB is 20-fold larger than IA0 [Equation (A18)]. Then, Cl secretion gradually decreases, reaching a steady state; IA (∞) is 12.4-fold larger than IA0 (see Equation A18). If we expect to induce continuous stimulation of Cl secretion, we should apply any compounds providing with larger activation on apical Cl channels than on basolateral Cl channels. τ of Cl secretion transition in response to a change in apical Cl conductance (GB0=30 GA0, GA=20 GA0, GB= GB0) obtained from Equation (A14) is 150CvF(exp(FV0RT)-1)[Cl-]oV0GA0.

FIGURE 2
www.frontiersin.org

Figure 2. The effect of forskolin in the presence of PKI on Cl secretion (IA) and simulation of Cl secretion (IA) and [Cl-]i under a condition with only a change in the apical Cl conductance (GA=20 GA0) without any change in the basolateral Cl conductance (GB=GB0), Cl uptake across the basolateral membrane (JC=JC0) or membrane potential (V = V0) based on GB0=30 GA0: GA=20 GA0, GB=GB0, GB0=30 GA0, GAG0A>GBGB0,JC=JC0,V=V0. (A) The effect of forskolin in the presence of PKI on Cl secretion (IA). Forskolin (10 μM) applied to the apical and basolateral solutions (at time = 0 min in A) 10 min after addition of 2 μM PKI to the apical and basolateral solutions induced an increase in IA followed by a decrease reaching a steady level larger than the basal one. (B) Simulation of IA shown in (A). A stimulant is applied at time = 0 min. IA(0)=20 IA0 [see Equation (A18)]. Then, Cl secretion gradually decreases, reaching a steady state; IA=12.4 IA0 [see Equation (A18)]. (C) Simulation of [Cl-]i in cells secreting Cl secretion (IA) shown in (B). A stimulant is applied at time = 0 min. [Cl-]i gradually decreases.

Simulated Cl secretion under this condition shows a transient increase followed by a decline toward to a steady state (Figure 2B). The simulated Cl secretion pattern shown in Figure 2B was similarly observed in carbachol-stimulated Cl secretion (Hendrick et al., 2014). This means that even if carbachol would activate the Cl channel located at the basolateral membrane, the action of carbachol on the basolateral Cl channel would be smaller than that on the apical Cl channel GAGA0>GBGB0. Combining the simulation and the experimental observation shown in a previous report (Hendrick et al., 2014), we could obtain more information on modulation of the Cl channels/transporters compared with only the experimental observation. Figure 2C shows a simulated change in [Cl-]i of the cells secreting Cl presented in Figure 2B. This simulation indicates that [Cl-]i decreases with a time course (Figure 2C) same as that in the declining phase in Cl secretion following the transient increase (Figure 2B), although [Cl-]i shows no transient change unlike Cl secretion (IA; see Figures 2B,C).

A case where both the apical Cl conductance (GA=2.5GA0) and the basolateral Cl conductance (GB=1.5GB0) increase, but a relative change in an apical Cl conductance is larger than that in the basolateral Cl (GAGA0>GBGB0): GA=2.5 GA0, GB=1.5 GB0, GB0=30GA0, GAGA0>GBGB0, JC= JC0 (CTAct=CTAct0), V = V0.

For example, if a condition, GA=2.5 GA0, GB=1.5 GB0 and GB0=30GA0 (a relative change in GA is larger than that in GB), is considered, the instantaneous Cl secretion (IA (0)) is 2.5-fold larger than IA0 [see Equation (A18)]. Then, Cl secretion gradually decreases, reaching a steady state; IA(∞) is 77.547.5 -fold (about 1.63-fold) larger than IA0 (see Equation A18). τ of Cl secretion transition in response to a change in apical Cl conductance (GA=2.5 GA0,GB=1.5 GB0, and GB0=30 GA0) obtained from Equation (A14) is 147.5CvF(exp(FV0RT)-1)[Cl-]oV0GA0. Hollenhorst et al. have reported an observation on acetylcholine-stimulated Cl secretion (Hollenhorst et al., 2012) similar to this simulation.

In this case, GAGA0>GBGB0 even without any changes in JC (JC=JC0:CTACT=CTACT0), IA(∞) is larger than IA0. This case indicates a time course of transcellular Cl secretion similar to that shown in Figure 2B. Although we could not determine if the basolateral Cl conductance changes, this simulation leads us to conclude that the stimulatory action of acetylcholine on Cl secretion (Hollenhorst et al., 2012) is mediated through activation of the apical Cl channel that is larger than that of the basolateral Cl channel (GAGA0 > GBGB0).

A case under a condition with the same extent increases in both the apical Cl conductance and the basolateral conductance (GAGA0=GBGB0) without any change in Cl uptake across the basolateral membrane (JC=JC0:CTAct=CTAct0) or membrane potential (V = V0) based on GB0=30 GA0; e.g., GA=18 GA0, GB=18 GB0,GB0=30 GA0, GAGA0=GBGB0, JC=JC0(CTAct=CTAct0), V = V0

Figure 3A shows Cl secretion in A6 cells treated with 100 μM daidzein: daidzein induced a transient increase in IA followed by a decrease, reaching a steady state identical to the basal one. Figure 3B shows simulation of IA in daidzein-treated A6 cells: IA increases in response to 18-fold elevation in apical and basolateral Cl conductance (GA=18 GA0, GB=18 GB0; the same extent increases in relative changes in GA and GB). The instantaneous Cl secretion (IA (0)) just after increases in GA and GB is 18-fold larger than IA0 [see Equation (A18)]. Then, Cl secretion gradually decreases, reaching a steady state; IA (∞) is identical to an initial value of Cl secretion [IA0; IA ()=IA0; see Equation (A18)]. τ of Cl secretion transition in response to a change in apical and basolateral Cl conductance (GA=18 GA0, GB= 18 GB0, and GB0=30 GA0) obtained from Equation (A14) is 1558 CvF(exp(FV0RT)-1)[Cl-]oV0GA0.

FIGURE 3
www.frontiersin.org

Figure 3. The effect of daidzein on Cl secretion (IA) and simulation of Cl secretion (IA) and [Cl-]i under a condition with the same extent increases in both the apical Cl and basolateral Cl conductance (GAGA0=GBGB0) without any change in Cl uptake across the basolateral membrane (JC=JC0) or membrane potential (V = V0): e.g., GA=18 GA0, GB=18 GB0, GB0=30 GA0, GAGA0=GBGB0, JC=JC0, V = V0. (A) The effect of daidzein on Cl secretion (IA). Daidzein (100 μM) applied to the apical and basolateral solutions at time = 0 min induced an increase in IA followed by a decrease reaching a steady level identical to the basal one. (B) Simulation of IA shown in (A). A stimulant is applied at time = 0 min. IA (0)=18 IA0 [see Equation (A18)]. Then, Cl secretion gradually decreases, reaching a steady state; IA ()=IA0 [see Equation (A18)]. (C) Simulation of [Cl-]i in cells secreting Cl secretion (IA) shown in (B). A stimulant is applied at time = 0 min. [Cl-]i decreases with time after application of stimulant.

These results indicate an interesting phenomenon that stimulants increasing apical and basolateral Cl conductance to the same extent (GA=18 GA0, GB= 18 GB0) have no influence on a steady-state Cl secretion (IA (∞) = IA (0)) but induces only transient stimulation to Cl secretion (an 18-fold increase). Figure 3C shows a simulated change in [Cl-]i of the cells secreting Cl presented in Figure 3B. The time course of [Cl-]i decrease is same as that at the declining phase of Cl secretion following a transient increase, although [Cl-]i shows no instantaneous change just after a change in GA (see Figures 3B,C). A similar phenomenon is observed in Cl secretion in cells treated with formaldehyde, bile acids, or Pasteurella multocida toxin (Hennig et al., 2008; Ao et al., 2013; Luo et al., 2013), suggesting that formaldehyde, bile acids, or Pasteurella multocida toxin would activate the apical and basolateral Cl channels in the same level (GAGA0=GBGB0).

A case under a condition with a relative change in the apical Cl conductance smaller than that in the basolateral Cl (GAGA0<GBGB0) without any change in Cl uptake across the basolateral membrane (JC=JC0: CTAct=CTAct0) or membrane potential (V = V0): e.g., GA=GA0,GB=2 GB0,GB0=30 GA0,GAGA0<GBGB0,JC=JC0(CTAct=CTAct0),V=V0

Figure 4A shows the time course of Cl secretion in response to two-fold increases in basolateral Cl conductance with no change in apical Cl conductance (GA=GA0, GB= 2 GB0). The instantaneous Cl secretion (IA(0)) just after an increase in GB without any change in GA does not change but remains at an identical level to IA0 (see Equation A18). Then, Cl secretion gradually decreases reaching a steady state; IA(∞) is smaller than an initial value of Cl secretion (IA0) [IA()=3161IA0 (about 0.5-fold of IA0); see Equation (A18)] in response to an increase in GB with no change in GA. τ of Cl secretion transition in response to a change in basolateral Cl conductance (GA=GA0, GB= 2 GB0, and GB0=30 GA0) obtained from Equation (A14) is 161CvF(exp(FV0RT)-1)[Cl-]oV0GA0.

FIGURE 4
www.frontiersin.org

Figure 4. Simulation of Cl secretion (IA) and [Cl-]i under a condition with a relative change in the apical Cl conductance smaller than that in the basolateral Cl (GAGA0<GBGB0) without Cl uptake across the basolateral membrane (JC=JC0) or membrane potential (V = V0): e.g., GA=GA0, GB=2 GB0, GB0=30 GA0, GAGA0<GBGB0, JC=JC0, V = V0. (A) Simulation of IA. A stimulant is applied at time =0 min. IA (0)=IA0 [see Equation (A18)]. Then, Cl secretion gradually decreases reaching a steady state; IA()=3161IA0 [about 0.5-fold of IA0; see Equation (A18)]. (B) Simulation of [Cl-]i in cells secreting Cl secretion (IA) shown in (A). A stimulant is applied at time = 0 min. [Cl-]i gradually decreases.

Figure 4B shows a simulated change in [Cl-]i of the cells secreting Cl presented in Figure 4A. The time course of [Cl-]i decrease (Figure 4B) is same as that of the declining phase of Cl secretion (Figure 4A). These results indicate an interesting phenomenon that stimulants increasing only basolateral Cl conductance have no instantaneous influence on Cl secretion but induce a gradual decrease in Cl secretion (Figure 4A) caused by gradual diminution in [Cl-]i (Figure 4B) due to an increase by an increase in a back flux of Cl to the basolateral space mediated through an increase in the basolateral Cl conductance.

A case where the membrane potential is changed (VV0) by application of compounds affecting the membrane conductance without any effects on Cl uptake via a pathway (JC = JC0:CTAct = CTAct0) or apical or basolateral Cl conductance (GA = GA0 and GB = GB0)

For example, we calculated IA (0) when V changes to −70 mV (V = −70 mV) from −40 mV (V0 = −40 mV) due to activation of K+ channel at the basolateral membrane at 25°C, the instantaneous Cl secretion (IA (0)) is as follows (see Equation A18) using the information on IA0 (0.08 ± 0.02 μA/cm2, n = 9; Table 2). IA (0) = −0.08 -0.07-0.04 1-0.21041-0.0654 +0.12 -0.07 · 3251-0.0654 (0.2104-0.0654)=-0.5418 (μAcm2). As mentioned in the section of Method, it is notable that the value of IA has a negative one, however IA shown in figures is converted to a positive one. This means that only a change in the membrane potential from −40 to −70 mV causes a transient increase in Cl efflux (IA(t)) about 6.8-fold [Equation (A18)] just after the change of the membrane potential followed by a decline toward the value before the change in the membrane potential (IA0) with exponential function with τ of CvF(exp(FVRT)-1) [Cl-]oV(GA0+GB0) (Equation (A14) in Appendix; Figure 5A). This decline of IA(t) with τ of 129CvF(exp(FVRT)-1) [Cl-]oVGA0 [Equation (A14) in Appendix] is due to a decrease in [Cl-]i (Figure 5B) caused by hyperpolarization of the apical and basolateral.

TABLE 2
www.frontiersin.org

Table 2. Cl currents, and apical and basolateral conductance under basal and forskolin (FK)-stimulated conditions.

FIGURE 5
www.frontiersin.org

Figure 5. Simulation of Cl secretion (IA) and [Cl-]i under a condition where the membrane potential is changed (VV0; e.g., V0 = −40 mV and V = −70 mV) by application of compounds affecting the membrane potential without any effects on Cl uptake via a pathway (JC=JC0) or the apical or basolateral Cl conductance (GA=GA0 and GB=GB0; GB0=30 GA0). (A) A change in the membrane potential from −40 to −70 mV at time = 0 min causes a transient increase in Cl secretion (IA) about 6.8-fold [Equation (A18)] just after the change of the membrane potential. Then, Cl secretion (IA) gradually decreases, reaching a steady state; IA (∞) is identical to an initial value of Cl secretion, IA0; IA ()=IA0 [see Equation (A18)]. (B) Simulation of [Cl-]i. [Cl-]i in cells secreting Cl secretion (IA) shown in (A) decreases after a change in the membrane potential from −40 to −70 mV at time = 0 min.

A case where Cl uptake via a pathway (JC) changes (JCJC0: CTActCTAct0) time-dependently after application of compounds affecting intracellular environments modifying Cl secretion without any effects on apical or basolateral Cl conductance (GA=GA0 and GB=GB0) or membrane potential (V = V0)

Figure 6A shows Cl secretion (IA) observed in A6 cells treated with 100 μM apigenin, which induced a gradual increase in IA, reaching a steady level larger than its initial level. Figure 6B shows the time course of Cl secretion under this condition simulated by using Equation (A18) (Appendix) in response to two-fold increases in Cl uptake via a pathway, JC, (JC()= 2 JC0: CTACT()=2CTACT0) after application of a compound affecting intracellular environments modifying Cl secretion with GA=GA0, GB=GB0 and V = V0. IA(0)=IA0 at JC()= 2 JC0 (CTACT()=2CTACT0), GA=GA0, GB=GB0, and V = V0. Then, Cl secretion gradually increases reaching a steady state; IA ()=2 IA0 [Equation (A18) in Appendix]. This simulation (Figure 6B) mimics Cl secretion in A6 cells treated with 100 μM apigenin (Figure 6A). Figure 6C shows a simulated change in [Cl-]i of A6 cells showing Cl secretion shown in Figure 6B. The time course of [Cl-]i increase is same as that of elevating phase of Cl secretion. Similar observations on Cl secretion in cells treated with quercetin or kaemferol (Cermak et al., 1998, 2002; Illek and Fischer, 1998; Asano et al., 2009; Zhang et al., 2011), suggesting that quercetin or kaemferol would induce Cl secretion gradually after their addition by activating NKCC1.

FIGURE 6
www.frontiersin.org

Figure 6. The effect of apigenin on Cl secretion (IA) and simulation of Cl secretion (IA) and [Cl-]i under a condition where Cl uptake via a pathway (JC) changes (JC() = 2JC0: CTACT() = 2CTACT0) after application of compounds affecting intracellular environments modifying Cl secretion without any effects on apical or basolateral Cl conductance (GA=GA0 and GB=GB0) or membrane potential (V = V0). (A) The effect of apigenin on Cl secretion (IA). Apigenin (100 μM) applied to the apical and basolateral solutions at time = 0 min induced a gradual increase in IA reaching a steady level larger than the basal one. (B) Simulation of IA in cells secreting Cl secretion (IA) shown in (A). A stimulant is applied at time = 0 min. IA (0)=IA0 [see Equation (A18)]. Then, Cl secretion (IA) gradually increases reaching a steady state; IA()=2 IA0 [see Equation (A18)]. (C) Simulation of [Cl-]i in cells secreting Cl secretion (IA) shown in (B). A stimulant is applied at time = 0 min. [Cl-]i gradually increases after a two-fold change in JC (JC()=2 JC0: CTACT()=2CTACT0).

A case where Cl uptake via a pathway (JC) is time-dependently changed (JCJC0: CTActCTAct0) after application of compounds that affect intracellular environments modifying Cl secretion associated with changes in apical and basolateral Cl conductance (GAGA0) and (GBGB0) but without any change in the membrane potential (V = V0)

Genistein (100 μM) induced a biphasic increase in IA in A6 cells (Figure 7A): (1) the first phase consisted of a transient increase in IA followed by a decrease, and (2) at the second phase IA gradually increased toward a steady level after reaching the minimum value at the first phase. We simulated this phenomenon of IA observed in A6 cells treated with genistein (Figure 7A; Equation A18). There are two time constants (τ2=1α+βτ1=CvF (exp(FVRT)-1) [Cl-]oV(GA+GB)) in the change in IA(t) as shown in Equations (A11) and (A14). Based on the inequality, τ2=1α+βτ1=CvF (exp(FVRT)-1) [Cl-]oV(GA+GB), the part of IA contributing to the rising phase of IA(t) at the second phase is presented by Equation (A18) with a time constant of τ2, 1α+β [Equation (A11)].

FIGURE 7
www.frontiersin.org

Figure 7. The effect of genistein on Cl secretion (IA) and simulation of Cl secretion (IA) and [Cl-]i under a condition where Cl uptake via a pathway JC changes (JCJC0) time-dependently after application of compounds that affect intracellular environments modifying Cl secretion associated with changes in apical and basolateral Cl conductance (GAGA0 and GBGB0) but without any change in the membrane potential (V = V0). (A) The effect of genistein on Cl secretion (IA). Genistein (100 μM) applied to the apical and basolateral solutions at time = 0 min induced a biphasic increase in IA: (1) the first phase consists of an initial, transient increase in IA followed by a decrease after application of a stimulant at time = 0, and 2) the second phase consists of a gradual increase in IA after IA reaches its minimum value at the first phase toward a steady level larger than the basal level. (B) Simulation of IA in cells secreting Cl secretion (IA) shown in (A). A stimulant is applied at time = 0 min. Cl secretion shows a five-fold increase just after application of a compound modifying Cl uptake via a pathway (JC= 3 JC0), apical Cl conductance (GA=5 GA0) and basolateral Cl conductance (GB=4 GB0; GB0=30 GA0) without any change in membrane potential (V = V0) according to Equation (A18). Then, Cl secretion gradually decreases from its peak. After reaching its minimum value, Cl secretion again increases toward a steady value shown in Equation (A18). (C) Simulated [Cl-]i of cells secreting Cl shown in (B). A stimulant is applied at time = 0 min. [Cl-]i transiently decreases after increases in GA and GB followed by a gradual increase due to an increase in JC.

Figure 7B shows the time course of Cl secretion in response to a compound causing JC()= 3 JC0 (CTACT()=3CTACT0), GA= 5 GA0 and GB= 4 GB0 (GB0= 30 GA0) with no effect on membrane potential (V = V0) according to Equation (A18) in Appendix. IA (0)=5 IA0 (Equation A18 in Appendix). Then, IA gradually decreases from its peak (Figure 7B). After reaching its minimum value, IA again increases toward a steady value. This change in IA shown in Figure 7B mimics Cl secretion presented in Figure 7A. Similar observations have been reported (Niisato et al., 1999; Deachapunya and Poonyachoti, 2013). Figure 7C shows a simulated change in [Cl-]i of the cells secreting Cl presented in Figure 7B. This simulation indicates that [Cl-]i decreases with a time course similar to that in the declining phase in Cl secretion following the initial, transient increase at the first phase due to a [Cl-]i decrease caused by an increase in Cl secretion, and then increases due to elevation of Cl uptake via NKCC1, although [Cl-]i shows no initial, transient increase unlike Cl secretion (IA; see Figures 7B,C). This simulation suggests us that genistein would have stimulatory action on both the apical Cl channel and the basolateral NKCC1.

As shown in Equation (A19), we are able to determine JC()JC0(=CTACT()CTACT0 ) by measuring IAo, IA(), GA0, GB0, GA, and GB. We tried to determine the effect of fosrskolin (10 μM) on JC. After application of forskolin, IA increased reaching an initial peak, then showed a decline phase followed by a gradual increase (Figure 8A). NPPB of 100 μM applied to the apical solution diminished IA (Figure 8A). Statistical results indicate that IAo was 0.08 ± 0.02 μA/cm2 (n = 9) and IA(∞) stimulated by forskolin (10 μM) was 6.40 ± 0.14 μA/cm2 (n = 9; Table 2). To measure the apical Cl conductance (GA), we permeabilized the basolateral membrane using nystatin. Forskolin (10 μM) applied at t = 0 increased IA. Nystatin (50 μM) added to the basolateral membrane at 30 min after application of forskolin induced IA with a negative value (Figure 8B; see Discussion in detail). NPPB (100 μM) added to the apical solution at 45 min after application of nystatin (75 min after forskolin application in Figure 8B) diminished IA. Under this condition, we measured apical Cl conductance (GA) as the NPPB-sensitive conductance (Table 2). Similar to measurement of apical Cl conductance (GA), we also measured the basolateral Cl conductance (GB). As shown in Figure 8C, forskolin increased IA. Nystatin added to the apical membrane induced an increase in IA. NPPB (100 μM) added to the basolateral solution diminished IA. Under this condition, we measured basolateral Cl conductance (GB) as the NPPB-sensitive conductance (Table 2). The nystatin-induced increase in IA was due to elevation of Na+, K+-pump current (see Discussion in detail). NPPB added to the basolateral solution decreased IA by diminishing the Na+, K+-pump activity via blockade of the basolateral Cl channel (conductance; Niisato and Marunaka, 1999). Under this condition, we measured basolateral Cl conductance (GB) as the NPPB-sensitive conductance (Table 2). We also measured GA0 and GB0 without application of forskolin (Table 2) using a similar protocol mentioned above. Applying the values of IAo, IA(), GA0, GB0, GA, and GB shown in Table 2 to Equation (A19), we estimated the forskolin-induced activation of the electro-neutral ion transporter, NKCC1; the activity of NKCC1 was increased to 4.7-fold. Thus, we are able to estimate a relative change in an electro-neutral ion transporter, such as NKCC1, using the electrophysiological techniques combined with mathematical simulation proposed in the present study.

FIGURE 8
www.frontiersin.org

Figure 8. Forskolin action on Cl secretion (IA), apical conductance (GA) and basolateral Cl conductance (GB) under an assumption with no change in membrane potential (V = V0), and simulation of Cl secretion (IA) and [Cl-]i. (A) The effect of foskolin on Cl secretion (IA). Forskolin (10 μM) applied to the apical and basolateral solutions at time = 0 min induced an initial, transient increase in Cl secretion (IA (0)) followed by a decrease. After reaching its minimum level, Cl secretion (IA) again increased toward a steady value. After Cl secretion (IA) reached a steady value, we applied 100 μM NPPB to the apical solution at time = 30 min, which diminished Cl secretion (IA). (B) Apical Cl conductance (GA) of A6 cells treated with forskolin. After application of forskolin (10 μM) at time = 0 min, we added nystatin (50 μM) to the basolateral solution to permeabilize the basolateral membrane at time = 30 min. Then, we applied 100 μM NPPB to the apical solution at time at 75 min. Under this condition, we determined the NPPB-sensitive conductance of apical membrane. (C) Basolateral Cl conductance (GB) of A6 cells treated with forskolin. After application of forskolin (10 μM) at time = 0 min, we added nystatin (50 μM) to the apical solution to permeabilize the apical membrane at time = 30 min. Then, we applied 100 μM NPPB to the basolateral solution at time 75 min. Under this condition, we determined the NPPB-sensitive conductance of basoalteral membrane. (D) Simulated Cl secretion (IA) of A6 cells treated with forskolin. A stimulant is applied at time = 0 min. We obtained IA(0)=50 IA0, IA()=79 IA0, GA=50 GA0 and GB = 1.33 GB0 from this simulation [Equation (A18)]. (E) Simulated [Cl-]i of cells secreting Cl shown in (D). A stimulant is applied at time = 0 min. [Cl-]i transiently decreases after application of forskolin followed by a gradual increase.

Figure 8D shows simulated Cl secretion (IA) of A6 cells treated with forskolin shown in Figure 8A. Figure 8E shows simulated [Cl]i of cells secreting Cl shown in Figure 8D. A transient increase in Cl secretion shown in Figure 8D is due to an increase in apical Cl conductance. A declining phase of Cl secretion following a transient increase in Cl secretion shown in Figure 8D is due to a decrease in [Cl]i (Figure 8E) caused by elevation of Cl secretion. A gradual increase in Cl secretion following the biphasic change in Cl secretion shown in Figure 8D is due to elevation of [Cl]i (Figure 8E) caused by activation of Cl transporter participating in uptake of Cl from the basolateral space.

Discussion

In the present study, we report a mathematical model simulating transcellular Cl secretion combined with electrophysiological measurements. Using this method, we are able to provide information on regulation of the activity of Cl channels and transporters located at the apical and basolateral membranes contributing to the transcellular Cl secretion, and determine the activity of electro-neutral Cl transporter such as NKCC1.

We assumed that activity of Cl transporter, NKCC1, contributing to Cl uptake across the basolateral membrane does not depend on [Cl-]i, although activity of NKCC1 depends on the chemical potential difference between [Cl-]o and [Cl-]i in fact. In the present study, [Cl-]o did not change. This means that we should consider modification of NKCC1 activity by [Cl-]i. However, the NKCC1-mediated ion transport is driven by electrochemical potential of Na+ between the intracellular and basolateral spaces across the basolateral membrane. This driving force of Na+ is much larger than chemical potential of Cl, even that NKCC1 participates in two Cl uptake associated with one Na+ uptake. Therefore, this mathematical model proposed in the present study provides some valuable information, although the mathematical model proposed in the present study has limitation and we should consider a mathematical model including [Cl-]i-dependent regulation of NKCC1 activity like other studies dealing with this point (Weinstein and Krahn, 2010; Weinstein, 2010a,b).

We measured the apical Cl conductance (GA) by pemeabilizing the basolateral membrane applying nystatin to the basolateral solution (Figure 8B). Application of nystatin to the basolateral solution transiently induced IA with a negative value (Figure 8B), which would mean Cl absorption (Cl influx across the apical membrane). This basolateral application of nystatin initially induces depolarization of the basolateral membrane by increasing conductance for monovalent ions including Na+, K+, and Cl, associated with a gradual increase in [Cl-]i. Depolarization of basolateral membrane, in general, induces depolarization of apical membrane in some extent dependent on resistance (conductance) of tight junction (Marunaka, 2014a). This means that nystatin applied to the basolateral solution would initially induce Cl influx (absorption) from the apical solution via depolarization of apical membrane (see an initial, transient change in IA with a negative value shown in Figure 8B). After nystatin initially induces depolarization of basolateral membrane by increasing conductance for monovalent ions including Na+, K+, and Cl, [Cl-]i would gradually increase via increment of Cl influx due to nystatin-induced depolarization of membrane associated with elevation of Cl conductance. This increase in [Cl-]i abolishes Cl movement across the apical membranes (abolishment of Cl movement (IA) observed after the nystatin-induced transient change in IA value shown in Figure 8B). On the contrary, apical application of nystatin had an opposite effect on IA: apical application of nystatin induced an increase in IA with a positive value. As mentioned above, IA increase caused by apical application of nystatin was due to elevation of the Na+, K+-pump current mediated by nystatin-induced increment of the intracellular Na+ concentration via an increase in Na+ influx across the apical membrane (Niisato and Marunaka, 1999).

Since NPPB of 100 μM used in the present might block other types of ion channels in addition to Cl channels (Kucherenko et al., 2013), we should consider a possibility that the Cl conductance measured as the NPPB-sensitive conductance would be overestimated and NPPB would affect epithelial ion transport by blocking some ion channels in addition to blockade of Cl channels. Our previous study using A6 cells (Niisato and Marunaka, 1999) has clearly indicated that the NPPB-sensitive conductance is compatible with the Cl conductance measured by replacement of Cl with an anion such as gluconate much less permeable to Cl channels. Therefore, at least in A6 cells the NPPB-sensitive conductance could be used as the Cl conductance.

A mathematical model of rat ascending Henle limb (Weinstein and Krahn, 2010; Weinstein, 2010a,b) simulates activity of NKCCs depending on concentrations of Na+, K+ and Cl in detail, and Cl flux mediated via NKCCs. However, no information is available on unidirectional epithelial Cl transport (Cl secretion) mediated via NKCCs and Cl channels expressed in polarized epithelial cells from these studies (Weinstein and Krahn, 2010; Weinstein, 2010a,b), although these studies provide activity characteristics of various types of ion transporters and channels dependent on concentrations of Na+, K+, and Cl in renal epithelial cells. Similar to our observation shown in Figure 2A, Hartmann and Verkman indicate that an increase in apical Cl conductance induces a biphasic change in Cl secretion: sudden elevation of apical Cl conductance causes an immediate increase in Cl secretion followed by a slower decrease to a level larger than the baseline at a steady state (Hartmann and Verkman, 1990). They report that the slow decrease in Cl secretion following the immediate increase caused by elevation of apical Cl conductance would be due to a decrease in [Cl]i mediated through the increase in Cl secretion based on elevation of apical Cl conductance (Hartmann and Verkman, 1990). Figure 2C clearly indicates a gradual decrease in [Cl]i, which diminishes a chemical gradient for Cl secretion leading to a decrease in Cl secretion with a slower time course, strongly supporting the report by Hartmann and Verkman (1990). The model proposed by Hartmann and Verkman (1990) contains apical Cl conductance and basolateral NKCC, but not basolateral Cl conductance, while a model proposed in the present study contains basolateral Cl conductance contributing to a back flux of Cl. Simulation of Cl secretion shown in Figure 2 proposes no change in Cl uptake, but a larger increase in apical Cl conductance than that in basolateral Cl conductance (GAGA0>GBGB0). This means that changes in apical and basolateral conductance under a condition of GAGA0>GBGB0 even without any change in Cl uptake increase Cl flux to the apical side associated with a decrease in Cl flux to the basolateral side by redistributing Cl fluxes to apical and basolateral sides dependent on a change in the ratio of GAGA0 and GBGB0. This phenomenon is able to explain the observation that Cl secretion at a steady state after stimulation is larger than the baseline Cl secretion reported by Hartmann and Verkman (1990). Fischer et al. (2007) have reported an observation in airway epithelial cells: application of CdCl2 or DIDS (Cl channel blockers) to the basolateral solution increases Cl secretion stimulated by forskolin. Further, a study (Duta et al., 2006) has indicated that a Cl channel blocker, DIDS, added to the basolateral solution elevates Cl secretion. He et al. (2011) has also reported that a Cl channel inhibitor, CaCCinh-A01, applied to the basolateral solution enhances Cl secretion in colonic epithelial tissues. Thus, these reports (Duta et al., 2006; Fischer et al., 2007; He et al., 2011) indicate that basolateral Cl conductance contributes to a back flux of Cl: a decrease in basolateral Cl conductance with no change in apical Cl conductance (GAGA0>GBGB0) stimulates Cl, supporting the model proposed in the present study. Further, Hartmann and Verkman indicate that a stimulant-induced Cl secretion in cells with smaller cell volume shows a faster change compared to that in cells with larger cell volume (Hartmann and Verkman, 1990). The present study also clearly showed that τ of change in IA (t) is in direct proportion to cell volume [see Equation (A14)], supporting quantitatively the report by Hartmann and Verkman (1990).

Conclusion

The present study provides an idea how transepithelial Cl secretion is modified by changes in activity of Cl channels and transporters and a method for determining changes in the activity of electro-neutral Cl transporters. The ideas and methods demonstrated in the present study provide powerful tools for the clarification of the regulatory mechanism of transepithelial Cl secretion, and very useful for development of new drugs modifying epithelial Cl secretion, although the mathematical model proposed in the present study has limitation to be adaptable to Cl secretion measured in living epithelial cells only under a short-circuit condition.

Funding

This work was supported by Grants-in-Aid from Japan Society of the Promotion of Science (24590283 and 15K08183 to NN, 25893201 and 26713008 to AT, 25670111 and 15K15034 to YM), Salt Science (1235 to YM and NN, 1429 and 1542 to AT), KIT-KPUM-KPU-KPhU Collaborative Research Grant (2013 and 2015) to YM, NN and AT, Society for Research on Umami Taste, Nestlé Nutrition Council, Japan to AT, Kyoto Prefectural Public University Corporation to AT, Kyoto-Funding for Innovation in Health-related R&D Fields to YM and NN, Fuji Foundation for Protein Research to YM and NN, and Cell Research Conference to YM.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Ao, M., Sarathy, J., Domingue, J., Alrefai, W. A., and Rao, M. C. (2013). Chenodeoxycholic acid stimulates Cl secretion via cAMP signaling and increases cystic fibrosis transmembrane conductance regulator phosphorylation in T84 cells. Am. J. Physiol. Cell Physiol. 305, C447–C456. doi: 10.1152/ajpcell.00416.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Asano, J., Niisato, N., Nakajima, K., Miyazaki, H., Yasuda, M., Iwasaki, Y., et al. (2009). Quercetin stimulates Na+/K+/2Cl cotransport via PTK-dependent mechanisms in human airway epithelium. Am. J. Respir. Cell Mol. Biol. 41, 688–695. doi: 10.1165/rcmb.2008-0338OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Cermak, R., Föllmer, U., and Wolffram, S. (1998). Dietary flavonol quercetin induces chloride secretion in rat colon. Am. J. Physiol. 275, G1166–G1172.

PubMed Abstract | Google Scholar

Cermak, R., Kuhn, G., and Wolffram, S. (2002). The flavonol quercetin activates basolateral K+ channels in rat distal colon epithelium. Br. J. Pharmacol. 135, 1183–1190. doi: 10.1038/sj.bjp.0704564

PubMed Abstract | CrossRef Full Text | Google Scholar

Deachapunya, C., and Poonyachoti, S. (2013). Activation of chloride secretion by isoflavone genistein in endometrial epithelial cells. Cell. Physiol. Biochem. 32, 1473–1486. doi: 10.1159/000356584

PubMed Abstract | CrossRef Full Text | Google Scholar

Duta, V., Duta, F., Puttagunta, L., Befus, A. D., and Duszyk, M. (2006). Regulation of basolateral Cl channels in airway epithelial cells: the role of nitric oxide. J. Membr. Biol. 213, 165–174. doi: 10.1007/s00232-006-0062-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fischer, H., Illek, B., Finkbeiner, W. E., and Widdicombe, J. H. (2007). Basolateral Cl channels in primary airway epithelial cultures. Am. J. Physiol. Lung Cell. Mol. Physiol. 292, L1432–L1443. doi: 10.1152/ajplung.00032.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldman, D. E. (1943). Potential, impedance, and rectification in membranes. J. Gen. Physiol. 27, 37–60. doi: 10.1085/jgp.27.1.37

PubMed Abstract | CrossRef Full Text | Google Scholar

Hartmann, T., and Verkman, A. S. (1990). Model of ion transport regulation in chloride-secreting airway epithelial cells. Integrated description of electrical, chemical, and fluorescence measurements. Biophys. J. 58, 391–401. doi: 10.1016/S0006-3495(90)82385-7

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Q., Halm, S. T., Zhang, J., and Halm, D. R. (2011). Activation of the basolateral membrane Cl conductance essential for electrogenic K+ secretion suppresses electrogenic Cl secretion. Exp. Physiol. 96, 305–316. doi: 10.1113/expphysiol.2010.055038

PubMed Abstract | CrossRef Full Text | Google Scholar

Hendrick, S. M., Mroz, M. S., Greene, C. M., Keely, S. J., and Harvey, B. J. (2014). Bile acids stimulate chloride secretion through CFTR and calcium-activated Cl channels in Calu-3 airway epithelial cells. Am. J. Physiol. Lung Cell. Mol. Physiol. 307, L407–L418. doi: 10.1152/ajplung.00352.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hennig, B., Orth, J., Aktories, K., and Diener, M. (2008). Anion secretion evoked by Pasteurella multocida toxin across rat colon. Eur. J. Pharmacol. 583, 156–163. doi: 10.1016/j.ejphar.2008.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodgkin, A. L., and Katz, B. (1949). The effect of sodium ions on the electrical activity of giant axon of the squid. J. Physiol. 108, 37–77. doi: 10.1113/jphysiol.1949.sp004310

PubMed Abstract | CrossRef Full Text | Google Scholar

Hollenhorst, M. I., Lips, K. S., Wolff, M., Wess, J., Gerbig, S., Takats, Z., et al. (2012). Luminal cholinergic signalling in airway lining fluid: a novel mechanism for activating chloride secretion via Ca2+-dependent Cl and K+ channels. Br. J. Pharmacol. 166, 1388–1402. doi: 10.1111/j.1476-5381.2012.01883.x

CrossRef Full Text | Google Scholar

Ikehara, O., Hayashi, H., Waguri, T., Kaji, I., Karaki, S., Kuwahara, A., et al. (2014). Luminal trypsin induces enteric nerve-mediated anion secretion in the mouse cecum. J. Physiol. Sci. 64, 119–128. doi: 10.1007/s12576-013-0302-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Illek, B., and Fischer, H. (1998). Flavonoids stimulate Cl conductance of human airway epithelium in vitro and in vivo. Am. J. Physiol. 275, L902–L910.

PubMed Abstract | Google Scholar

Kim, D., Liao, J., and Hanrahan, J. W. (2014). The buffer capacity of airway epithelial secretions. Front. Physiol. 5:188. doi: 10.3389/fphys.2014.00188

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, K. C. (2012). Role of epithelial mucins during airway infection. Pulm. Pharmacol. Ther. 25, 415–419. doi: 10.1016/j.pupt.2011.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Komatani-Tamiya, N., Daikoku, E., Takemura, Y., Shimamoto, C., Nakano, T., Iwasaki, Y., et al. (2012). Procaterol-stimulated increases in ciliary bend amplitude and ciliary beat frequency in mouse bronchioles. Cell. Physiol. Biochem. 29, 511–522. doi: 10.1159/000338505

PubMed Abstract | CrossRef Full Text | Google Scholar

Kucherenko, Y. V., Wagner-Britz, L., Bernhardt, I., and Lang, F. (2013). Effect of chloride channel inhibitors on cytosolic Ca2+ levels and Ca2+-activated K+ (Gardos) channel activity in human red blood cells. J. Membr. Biol. 246, 315–326. doi: 10.1007/s00232-013-9532-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kunzelmann, K., and Schreiber, R. (2012). Airway epithelial cells–hyperabsorption in CF? Int. J. Biochem. Cell Biol. 44, 1232–1235. doi: 10.1016/j.biocel.2012.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, R. J., and Foskett, J. K. (2010). cAMP-activated Ca2+ signaling is required for CFTR-mediated serous cell fluid secretion in porcine and human airways. J. Clin. Invest. 120, 3137–3148. doi: 10.1172/JCI42992

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Comellas, A. P., Karp, P. H., Ernst, S. E., Moninger, T. O., Gansemer, N. D., et al. (2012). CFTR is required for maximal transepithelial liquid transport in pig alveolar epithelia. Am. J. Physiol. Lung Cell. Mol. Physiol. 303, L152–L160. doi: 10.1152/ajplung.00116.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, Y. L., Guo, H. M., Zhang, Y. L., Chen, P. X., Zhu, Y. X., Huang, J. H., et al. (2013). Cellular mechanism underlying formaldehyde-stimulated Cl secretion in rat airway epithelium. PLoS ONE 8:e54494. doi: 10.1371/journal.pone.0054494

PubMed Abstract | CrossRef Full Text | Google Scholar

Marques, R. D., Praetorius, H. A., and Leipziger, J. (2013). P2Y2 receptor knock-out mice display normal NaCl absorption in medullary thick ascending limb. Front. Physiol. 4:280. doi: 10.3389/fphys.2013.00280

PubMed Abstract | CrossRef Full Text | Google Scholar

Marunaka, Y. (2014a). Characteristics and pharmacological regulation of epithelial Na+ channel (ENaC) and epithelial Na+ transport. J. Pharmacol. Sci. 26, 21–36. doi: 10.1254/jphs.14R01SR

CrossRef Full Text | Google Scholar

Marunaka, Y. (2014b). Importance of expression and function of angiotensin II receptor type 1 in pulmonary epithelial cells. Respir. Physiol. Neurobiol. 196, 39–42. doi: 10.1016/j.resp.2014.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Marunaka, Y., Niisato, N., O'Brodovich, H., and Eaton, D. C. (1999). Regulation of an amiloride-sensitive Na+-permeable channel by a beta2-adrenergic agonist, cytosolic Ca2+ and Cl in fetal rat alveolar epithelium. J. Physiol. 515, 669–683. doi: 10.1111/j.1469-7793.1999.669ab.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyazaki, H., Shiozaki, A., Niisato, N., Ohsawa, R., Itoi, H., Ueda, Y., et al. (2008). Chloride ions control the G1/S cell-cycle checkpoint by regulating the expression of p21 through a p53-independent pathway in human gastric cancer cells. Biochem. Biophys. Res. Commun. 366, 506–512. doi: 10.1016/j.bbrc.2007.11.144

PubMed Abstract | CrossRef Full Text | Google Scholar

Niisato, N., Ito, Y., and Marunaka, Y. (1999). Activation of Cl channel and Na+/K+/2Cl cotransporter in renal epithelial A6 cells by flavonoids: genistein, daidzein, and apigenin. Biochem. Biophys. Res. Commun. 254, 368–371. doi: 10.1006/bbrc.1998.9952

PubMed Abstract | CrossRef Full Text | Google Scholar

Niisato, N., and Marunaka, Y. (1999). Activation of the Na+-K+ pump by hyposmolality through tyrosine kinase-dependent Cl conductance in Xenopus renal epithelial A6 cells. J. Physiol. 518, 417–432. doi: 10.1111/j.1469-7793.1999.0417p.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nilius, B., Sehrer, J., De Smet, P., Van Driessche, W., and Droogmans, G. (1995). Volume regulation in a toad epithelial cell line: role of coactivation of K+ and Cl channels. J. Physiol. 487(Pt 2), 367–378. doi: 10.1113/jphysiol.1995.sp020886

PubMed Abstract | CrossRef Full Text | Google Scholar

Sasamoto, K., Niisato, N., and Marunaka, Y. (2014). Estimation of activity of electro-neutral ion transporter by establishing mathematical model of Cl secretion in epithelial cells based on measurement of transcellular ion transport and membrane ion conductance. J. Physiol. Sci. 64, S227.

Stanton, B. A., Coutermarsh, B., Barnaby, R., and Hogan, D. (2015). Pseudomonas aeruginosa reduces VX-809 stimulated F508del-CFTR chloride secretion by airway epithelial cells. PLoS ONE 10:e0127742. doi: 10.1371/journal.pone.0127742

PubMed Abstract | CrossRef Full Text | Google Scholar

Stölting, G., Fischer, M., and Fahlke, C. (2014). CLC channel function and dysfunction in health and disease. Front. Physiol. 5:378. doi: 10.3389/fphys.2014.00378

PubMed Abstract | CrossRef Full Text

Sun, H., Niisato, N., Inui, T., and Marunaka, Y. (2014a). Insulin is involved in transcriptional regulation of NKCC and the CFTR Cl channel through PI3K activation and ERK inactivation in renal epithelial cells. J. Physiol. Sci. 64, 433–443. doi: 10.1007/s12576-014-0338-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, H., Niisato, N., and Marunaka, Y. (2013). A role of PI3 kinase in insulin-induced enhancement of cAMP-stimulated Cl secretion in renal epithelia A6 cells. J. Physiol. Sci. 63, S197.

Sun, H., Niisato, N., Nishio, K., Hamilton, K. L., and Marunaka, Y. (2014b). Distinct action of flavonoids, myricetin and quercetin, on epithelial Cl secretion: useful tools as regulators of Cl secretion. Biomed. Res. Int. 2014:902735. doi: 10.1155/2014/902735

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, D., Furutani, Y., Inoue, K., Kikukawa, T., Sakai, M., Fujii, M., et al. (2009). Effects of chloride ion binding on the photochemical properties of salinibacter sensory rhodopsin I. J. Mol. Biol. 392, 48–62. doi: 10.1016/j.jmb.2009.06.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Tohda, H., Foskett, J. K., O'Brodovich, H., and Marunaka, Y. (1994). Cl regulation of a Ca2+ activated non selective cation channel in beta agonist treated fetal distal lung epithelium. Am. J. Physiol. 266, C104–C109.

Google Scholar

Tokuda, S., Miyazaki, H., Nakajima, K., Yamada, T., and Marunaka, Y. (2009a). Hydrostatic pressure regulates tight junctions, actin cytoskeleton and transcellular ion transport. Biochem. Biophys. Res. Commun. 390, 1315–1321. doi: 10.1016/j.bbrc.2009.10.144

PubMed Abstract | CrossRef Full Text | Google Scholar

Tokuda, S., Miyazaki, H., Nakajima, K., Yamada, T., and Marunaka, Y. (2010). NaCl flux between apical and basolateral side recruits claudin-1 to tight junction strands and regulate paracellular transport. Biochem. Biophys. Res. Commun. 393, 390–396. doi: 10.1016/j.bbrc.2010.02.002

CrossRef Full Text | Google Scholar

Tokuda, S., Niisato, N., and Marunaka, Y. (2007). Regulation of paracellular ion conductances by NaCl gradients in renal epithelial cells. Biochem. Biophys. Res. Commun. 363, 566–571. doi: 10.1016/j.bbrc.2007.08.187

PubMed Abstract | CrossRef Full Text | Google Scholar

Tokuda, S., Niisato, N., Nagai, T., Taruno, A., Nakajima, K., Miyazaki, H., et al. (2009b). Regulation of paracellular Na+ and Cl conductances by hydrostatic pressure. Cell Biol. Int. 33, 949–956. doi: 10.1016/j.cellbi.2009.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Tokuda, S., Niisato, N., Nakajima, K., and Marunaka, Y. (2008). Regulation of the paracellular Na+ and Cl conductances by the NaCl-generated osmotic gradient in a manner dependent on the direction of osmotic gradients. Biochem. Biophys. Res. Commun. 366, 464–470. doi: 10.1016/j.bbrc.2007.11.167

PubMed Abstract | CrossRef Full Text | Google Scholar

Ussing, H. H., and Zerahn, K. (1951). Active transport of sodium as the source of electric current in the short-circuited isolated frog skin. Acta Physiol. Scand. 23, 110–127. doi: 10.1111/j.1748-1716.1951.tb00800.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinstein, A. M. (2010a). A mathematical model of rat ascending Henle limb. I. Cotransporter function. Am. J. Physiol. Renal. Physiol. 298, F512–F524. doi: 10.1152/ajprenal.00230.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinstein, A. M. (2010b). A mathematical model of rat ascending Henle limb. III. Tubular function. Am. J. Physiol. Renal. Physiol. 298, F543–F556. doi: 10.1152/ajprenal.00232.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinstein, A. M., and Krahn, T. A. (2010). A mathematical model of rat ascending Henle limb. II. Epithelial function. Am. J. Physiol. Renal. Physiol. 298, F525–F542. doi: 10.1152/ajprenal.00231.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, S., Smith, N., Schuster, D., Azbell, C., Sorscher, E. J., Rowe, S. M., et al. (2011). Quercetin increases cystic fibrosis transmembrane conductance regulator-mediated chloride transport and ciliary beat frequency: therapeutic implications for chronic rhinosinusitis. Am. J. Rhinol. Allergy 25, 307–312. doi: 10.2500/ajra.2011.25.3643

PubMed Abstract | CrossRef Full Text | Google Scholar

Appendix

We here show a model of transcellular Cl secretion in epithelial tissues (Figure 1). This model consists of three Cl moving pathways: (1) a releasing pathway of Cl from the intracellular space into the apical one across the apical membrane (Pathway A); (2) a releasing pathway of Cl from the intracellular space into the basolateral one (Pathway B); (3) an uptake pathway of Cl from the basolateral space into the intracellular space, such as NKCC1, across the basolateral membrane (Pathway C). The transcellular Cl secretion is mediated by the following pathways: (1) Cl is taken up into the intracellular space via Pathway C; (2) Cl taken up into the intracellular space by Pathway C is respectively released into the apical and basolateral spaces via Pathways A,B; Cl taken up by Pathway C released into the apical space via Pathway A only contributes to the transcellular Cl secretion. The transcellular Cl secretion is generally measured as short-circuit currents (ISC) in Ussing chamber (Ussing and Zerahn, 1951). Short-circuit currents (ISC) are measured under a voltage clamp condition with the apical membrane potential identical to the basolateral membrane potential (VA = VB) and all ionic compositions of the apical solution identical to those of basolateral solution. A change in [Cl-]i with time (= f(t)) is presented as a differential equation, Equation (A1). Based on Goldman-Hodgkin-Katz equation (Goldman, 1943; Hodgkin and Katz, 1949), Cl fluxes across the apical and basolateral membranes are respectively shown as Equations (A2) and (A3). We applied these equations to the present study on an assumption that the extracellular and intracellular ions and fluids stirred completely, instantly just after their movements across the membrane. Explanation on the parameters used in simulation in the present study is shown in Table 1.

df(t)dt=-SAJA+SBJB+SBJCCv(1)     (A1)
JA=-PAFVRT f(t)-[Cl-]oexp(FVRT) 1-exp(FVRT)     (A2)
JB=-PBFVRT f(t)-[Cl]oexp(FVRT) 1-exp(FVRT)     (A3)

where [Cl-]o is the extracellular Cl concentration, SA and SB are respectively the areas of the apical and basolateral membranes, JA and JB are respectively the passive Cl movements across the unit area of the apical and basolateral membranes driven by electrochemical potential of Cl between the extracellular and intracellular spaces at time = t, JC is the active Cl movement across the unit area of the basolateral membrane via Cl transporters such as NKCC1 driven by electrochemical potential of Na+ between the extracellular and intracellular spaces and is assumed to be independent of [Cl-]i and the membrane potential, PA and PB are respectively permeability of Cl per unit areas of the apical and basolateral membranes, F is Faraday constant, V is the membrane potential (the apical membrane potential = the basolateral membrane potential under the short-circuit condition), R is the gas constant, and T is the absolute temperature. Here, PASA and PBSB are respectively represented as Equations (A4) and (A5) on an assumption that GA and GB shown in the present study (specially measured under nystatin-permeabilized conditions) are subject to [Cl-]o and independent of [Cl-]i based on the relationship between permeability (PA and PB) and conductance (GA and GB) (Goldman, 1943) in order to simplify the derivation of Equation (A13):

PASA[Cl-]o=RTF2 GA    (A4)
PBSB[Cl-]o=RTF2 GB    (A5)

where GA and GB are respectively the apical and basolateral Cl conductance. Combining Equations (A1)–(A5), we obtain Equation (A6).

df(t)dt=1Cv[V(GA+GB)F(exp(FVRT)1)[Cl]of(t)                  +V(GA+GB)F(exp(FVRT)1)exp(FVRT)SBJC(t)]    (A6)

where Cv is cell volume in m3: SB, area of basolateral membrane in m2; V is membrane potential in V (volt); GA and GB are respectively Cl conductance of apical and basolateral membrane in (=1 V Q s ): S is siemens, Q is coulomb and s is second. In the present study, we apply a condition that Cv, V, GA, and GB after application of stimulants are constant without showing time-dependent changes to Equation (A6), even in the cases that stimulants affect V, GA, and GB: when a stimulant affects Cv, V, GA, and/or GB the stimulant-induced changes in Cv V, GA, and/or GB instantaneously reach steady levels keeping constant levels without showing time-dependent changes.

We further tried to establish a mathematical model regarding time-dependent change of JC on assumptions that JC(t) is regulated by biochemical factors but not directly [Cl-]i, although activity of Cl transporters contributing to Cl uptake across the basolateral membrane (see Figure 1) depends on the chemical potential difference between [Cl-]o and [Cl-]i in fact. We considered a case that SBJC(t) (= jCTAct(t)) is determined by the following condition: (1) total amounts (CTT) of Cl transporters (CT) participating in Cl uptake across the basolateral membrane are unchanged, (2) the Cl transporter has active and inactive forms, (3) the inactive form (CTInact) is converted to the active form with a rate of α0 under the basal condition or α after application of stimulant and the active form (CTAct) is converted to the inactive form with a rate of β0 under the basal condition or β after application of stimulant, (4) CTAct(t) and CTInact(t) are respectively the amounts of the active and inactive forms of CT at time = t after application of stimulants and (5) j is the amount of Cl flux mediated by CTAct per second (in mol/s). For example, α0 and α are respectively insertion rates of the Cl transporter into the basolateral membrane from the cytosolic store site under the basal condition and after application of stimulants, and β0 and β are respectively endocytotic rates of the Cl transporter from the basolateral membrane to the cytosolic space under the basal condition and after application of stimulants.

CTInactαβ CTAct    (A7)
dCTAct(t)dt=α CTInact(t)- β CTAct(t)    (A8)
CTT=CTInact(t)+ CTAct(t)    (A9)
SBJC(t)=j CTAct(t)                   =j [αα+βCTT(αα+βCTT                           α0α0+β0CTT)exp(tτ2)]                  =j [CTAct()(CTAct()                                                   CTAct0)exp(tτ2)]    (A10)

where CTACT0=α0α0+β0CTT is the amount of CTAct at a steady state under the basal condition, CTACT()= αα+βCTT is the amount of CTAct at a steady state after application of stimulants and τ2 is a time constant shown in Equation (A11) when CTAct(t) changes with time toward a steady state after application of stimulants.

τ2=1α+β(11)    (A11)

When JC(t) (= CTAct(t)) is constant, α and β are respectively equal to α0 and β0 (α = α0 and β = β0). We obtain Equation (A12) by combining Equations (A6) and (A10).

df(t)dt=1Cv[V(GA+GB)F(exp(FVRT)1)[Cl]of(t)                  +V(GA+GB)F(exp(FVRT)1)exp(FVRT)j CTAct()                  + j(CTAct()CTAct0)exp(tτ2)]    (A12)

Equation (A13) is obtained by resolving Equation (A12).

f(t)=(f(0)[Cl]oexp(FVRT)              +(exp(FVRT)1)[Cl]oV(GA+GB)F j CTAct()              j (CTAct()CTAct0)Cv (1τ11τ2))exp(tτ1)              +​​[Cl]oexp(FVRT)(exp(FVRT)1)[Cl]oV(GA+GB)F j CTAct()              +j (CTAct()CTAct0)Cv (1τ11τ2)exp(tτ2)    (A13)

where τ1 is a time constant shown in Equation (A14).

τ1=CvF(exp(FVRT) -1) [Cl-]oV(GA+GB)    (A14)

We obtain Equation (A15) by combining IA(t) = ZClF SA JA(t) and Equations (A2), (A4), and (A13).

IA(t)=VGA(1exp(FVRT))[Cl]o[(f(0)[Cl]oexp(FVRT)                 +(exp(FVRT)1)[Cl]oV(GA+GB)F j CTAct()                 j (CTAct()CTAct0)Cv (1τ11τ2))exp(tτ1)                (exp(FVRT)1)[Cl]oV(GA+GB)F j CTAct()                +j (CTAct()CTAct0)Cv (1τ11τ2)exp(tτ2)]    (A15)

In the present study, we define IA(t) generated by Cl effluxes (Cl movements from the intracellular space to the extracellular space) as a positive value based on characteristics of Cl movements described in a model (Figure 1), although the IA(t) generated by Cl efflux was originally defined as a negative current. Equation (A1) is written as Equation (A16) because of df(T)dT=0 under the basal condition.

SAJA0+SBJB0+SBJC0=0    (A16)

By combining Equations (A2)–(A5), and (A16), we obtain Equation (A17) on an assumption that [Cl-]i at time = 0 (f(0)) does not change instantaneously just after application of compound affecting transcellular Cl secretion in epithelial cells: f(0) = f0 (f0 is [Cl-]i under the basal condition without application of any compounds affecting transcellular Cl secretion in epithelial cells).

f(0)=F j CTAct0(1-exp(FV0RT))[Cl-]oV0 (GA0+GB0)+[Cl-]oexp(FV0RT)    (A17)

Using a method reported in our previous study (Miyazaki et al., 2008), we measured [Cl-]i (=f0), which was about 40 mM. V0 was measured to be approximately –40 mV by electrophysiological methods reported in previous reports (Tohda et al., 1994; Nilius et al., 1995; Marunaka et al., 1999).

By combining IA0=GA0GA0+GB0F j CTACT0 and Equations (A10), (A15), and (A17), we obtain Equation (A18).

IA(t)​​=IA0[GAGA0(CTAct()CTAct0(GA0+GB0GA+GB)            +(VV0(1exp(FV0RT)1exp(FVRT))CTAct()CTAct0(GA0+GB0GA+GB))            exp(tτ1))+1F(GA0+GB0GA0(VGA(CTAct()CTAct01)(1exp(FVRT))[Cl]o)            (exp(tτ2)exp(tτ1)Cv (1τ11τ2)))]            +[VGA1exp(FVRT)(exp(FV0RT)exp(FVRT))]exp(tτ1)    (A18)

IA(t) (Equation A18) represents the amount of Cl secretion as a current.

A case where Cl uptake via a pathway (JC) is changed (JC(T)JC0) time-dependently after application of compounds that affect intracellular environments modifying Cl secretion associated with changes in apical and basolateral Cl conductance (GAGA0 and GBGB0) but without the membrane potential (V = V0). Equation (A19) is derived from Equation (A18) under a condition with JC(T)  JC0, GA GA0, GB GB0, V = V0, and τ1 ≠ τ2.

JC()JC0=CTAct()CTAct0=IA()IA01+GBGA1+GB0GA0    (A19)

Keywords: Cl channel, Cl transporter, mathematical model, simulation, epithelial Cl secretion

Citation: Sasamoto K, Niisato N, Taruno A and Marunaka Y (2015) Simulation of Cl Secretion in Epithelial Tissues: New Methodology Estimating Activity of Electro-Neutral Cl Transporter. Front. Physiol. 6:370. doi: 10.3389/fphys.2015.00370

Received: 26 June 2015; Accepted: 20 November 2015;
Published: 23 December 2015.

Edited by:

Francisco Javier Alvarez-Leefmans, Wright State University, USA

Reviewed by:

Jorge Fischbarg, Instituto de Investigaciones Cardiológicas CONICET, Argentina
Eric Delpire, Vanderbilt University, USA
Dan R. Halm, Wright State University, USA

Copyright © 2015 Sasamoto, Niisato, Taruno and Marunaka. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yoshinori Marunaka, marunaka@koto.kpu-m.ac.jp

These authors have contributed equally to this work.

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.