Next Article in Journal
Oral Administration of Fermented Papaya (FPP®) Controls the Growth of a Murine Melanoma through the In Vivo Induction of a Natural Antioxidant Response
Next Article in Special Issue
Exploring Differential Connexin Expression across Melanocytic Tumor Progression Involving the Tumor Microenvironment
Previous Article in Journal
Vonoprazan-Based Third-Line Therapy Has a Higher Eradication Rate against Sitafloxacin-Resistant Helicobacter pylori
Previous Article in Special Issue
Connexins and Integrins in Exosomes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Perspective

A Barter Economy in Tumors: Exchanging Metabolites through Gap Junctions

Department of Physiology, Anatomy & Genetics, University of Oxford, Oxford OX1 3PT, UK
*
Author to whom correspondence should be addressed.
Cancers 2019, 11(1), 117; https://doi.org/10.3390/cancers11010117
Submission received: 30 December 2018 / Revised: 16 January 2019 / Accepted: 18 January 2019 / Published: 20 January 2019

Abstract

:
To produce physiological functions, many tissues require their cells to be connected by gap junctions. Such diffusive coupling is important in establishing a cytoplasmic syncytium through which cells can exchange signals, substrates and metabolites. Often the benefits of connectivity become apparent solely at the multicellular level, leading to the notion that cells work for a common good rather than exclusively in their self-interest. In some tumors, gap junctional connectivity between cancer cells is reduced or absent, but there are notable cases where it persists or re-emerges in late-stage disease. Diffusive coupling will blur certain phenotypic differences between cells, which may seem to go against the establishment of population heterogeneity, a central pillar of cancer that stems from genetic instability. Here, building on our previous measurements of gap junctional coupling between cancer cells, we use a computational model to simulate the role of connexin-assembled channels in exchanging lactate and bicarbonate ions down their diffusion gradients. Based on the results of these simulations, we propose that an overriding benefit of gap junctional connectivity may relate to lactate/bicarbonate exchange, which would support an elevated metabolic rate in hypoxic tumors. In this example of barter, hypoxic cancer cells provide normoxic neighbors with lactate for mitochondrial oxidation; in exchange, bicarbonate ions, which are more plentiful in normoxic cells, are supplied to hypoxic neighbors to neutralize the H+ ions co-produced glycolytically. Both cells benefit, and so does the tumor.

1. Trading within Biological Systems

In economics, barter is the exchange of goods or services without the involvement of money, under the premise that the two sides coincidentally want each other’s offer. As a means of trading, barter is now rare, but is historically significant for its role in enabling the division of labor. In early societies, barter would typically take place between partners who were not very familiar to one another [1], perhaps because a complete and immediate transaction was a way of diffusing concerns about the partner’s trustworthiness at a time when regulated systems of credit were not available [2]. The introduction of money and credit systems transformed trading and allowed societies to accrue considerable wealth.
A biological cell hosts a collection of specialized processes that inherently demonstrate division of labor. Many of these processes require an energetic input, supplied in the form of ATP, which serves as a universally-accepted currency of energy. An analogy can be drawn to an advanced economy, where metabolism generates ATP (‘money’) to power (‘pay’) various specialized biological processes (‘crafts’), so that these collectively support survival and proliferation (‘wealth’). Such a system functions properly because it is in the interest of all constituents that the cell thrives. At the multicellular level, specialization occurs among different tissues which are ‘paid’ by substrates (e.g., glucose) to engage in activities that benefit the whole body. A mechanism of credit (involving signaling and homeostasis) is implemented to ensure that the components of the system do not seek instantaneous returns. Tumors are, however, quite different. Cancer cells have a strong individualistic survival instinct, with little scope for altruistic gestures towards other cells [3,4]; moreover, cancer cells will often parasitically exploit neighboring cells, such as the host stroma [5,6]. Nonetheless, tumors have the potential to benefit from the division of labor because they are not homogenous tissues. Their heterogeneity arises from genetic differences among sub-populations of cancer cells, and also from variation in micro-environmental factors, such as oxygen supply which depends on the distance from the nearest capillary [3]. Reaping the benefits of the division of labor among selfish and suspicious partners may require a primitive exchange arrangement, such as barter. The instantaneous and reciprocal exchange of assets between cancer cells is achievable through inter-cellular conduits called gap junctions, first described in the 1960s [7].

2. Gap Junctions are Conduits for Exchanging Small- to Medium-Sized Molecules

Gap junctions are assembled from the juxtaposition of two hemichannels expressed at the surface of abutting cells [8]. These hemichannels are made of connexin proteins [9] which are coded by genes belonging to a family of at least 24 members [8,10]. The pore that traverses the complete channel offers a low-resistance pathway for exchanging solutes as big as 1 kDa. Essentially all metabolites, second messengers and inorganic ions are within this range of permeating substances, which results in a functional coupling between cells. The overall strength of this coupling depends on the number of gap junctional channels connecting a pair of cells.
Cell-to-cell coupling is essential to the function of many types of tissue (e.g., the spread of electrical excitation in the heart [11]) and, not surprisingly, most cells in the body are junctionally-coupled [12,13]. Many normal epithelia express connexins and manifest cell-to-cell coupling, such as Cx26, 32, 43 in pancreatic ducts and Cx26, 32, 43, 45 in colorectal epithelium [10,14,15]. Notable exceptions to the near-ubiquitous expression of connexins are blood cells, as these are stand-alone units, and fully-developed skeletal muscle, because their contractile units must be electrically-isolated to allow selective recruitment for the purpose of force gradation [16]. Historically, this list of exceptions also included tumors, on the basis that early observations suggested a lack of connectivity [17,18,19], subsequently supported by more recent work [20], as well as evidence for a tumor suppressing effect of connexins [21,22]. Contemporary findings have challenged this dogma [10], and it is now recognized that some cancers do express connexins, at least in some stages of progression [23].
There will be a myriad of roles for connexins in cancer, including some that relate more to their involvement in cell-to-cell adhesion [24,25,26,27], than the provision of a meaningful solute exchange pathway. Clinical studies have not been able to ascertain the direction of correlation between connexin expression and prognosis, as both positive and negative associations have been reported [10]. In colorectal cancer, for example, one study linked high Cx26 expression with poor prognosis and lung metastasis [28], whereas a subsequent study [29] came to the opposite conclusion. On the whole, there appears to be an isoform-specific pattern in the association between expression and prognosis, with Cx26 generally indicating worse prognosis and Cx43 linked to better prognosis [10]. A limitation of the aforementioned studies is that they have not determined the extent to which connexin expression produces meaningful coupling between cancer cells, a principal phenotypic readout of connexin function. Indeed, gap junctional communication is not routinely quantified in most studies of connexins in cancer (with a few exceptions, e.g., [21,25,30,31,32,33,34] and our work [35,36]). With the availability of antibodies and primers, recent studies of connexins in cancer cell lines have tended to quantify their levels in terms of mRNA or protein, which cannot faithfully predict functional coupling. Indeed, mislocalization of connexins to regions other than cell-cell contacts has been described for Cx26 in pancreatic [37] and colorectal cancer [28], putatively. Thus, our appreciation of the physiological role of cell-to-cell diffusive coupling in cancer, particularly in vivo, is still very limited. If, however, a specific type of malignancy manifested strong cell-to-cell diffusive coupling, it would be surprising if solute exchange had no bearing on cancer physiology.
One method of quantifying the degree of coupling between cells is based on fluorescence recovery after photobleaching (FRAP). A fluorescent dye, which can freely pass gap junctions (i.e., of molecular weight lower than 1 kDa), is loaded into a confluent monolayer; an exemplar colorectal cancer (CRC) cell line is shown in Figure 1. Fluorescence in one cell is then bleached (e.g., to 50% of starting level) using a high-power laser, and the signal time course registered in the bleached cell and its direct neighbors. The recovery time course is then fitted with a mono-exponential, and the initial slope (arrow in Figure 1D) provides a rate of fluorescence recovery, dF/dt, which is related to: (i) permeability, P, (ii) the difference in fluorescence between bleached and neighboring cells ΔF (i.e., the diffusion gradient), (iii) the surface area of the bleached cell at the junction with the neighboring cells, s, and (iv) the volume of the bleached cell, v:
dF/dt = P × ΔF × s/v
The s/v ratio can be approximated by quotient ρ of the bleached cell’s perimeter and its area in the xy plane (this assumes constant monolayer height). Thus, permeability, in units of µm/s, is:
P = d F / d t Δ F × ρ
This permeability relates, of course, to the fluorescent marker, which will be bulkier than many endogenous metabolites, signals or ions. Smaller molecules are thus expected to permeate faster by an amount determined by their size, charge and interactions with the pore area [38,39].
At first sight, solute dissipation across a syncytium of cells may appear to be the antithesis of what individual cancer cells need to compete in the process of somatic evolution because diffusion will blur differences in metabolic, signaling and ionic phenotype between connected cells. In late-stage disease, genetic instability and hyper-proliferation result in clonal heterogeneity [3], but since somatic evolution selects phenotype not genotype, the substrate for selection will become less diverse in connected cancer cells [4,40]. The syncytial behaviors of these units limit cell autonomy and constrain cell-to-cell variation, a fundamental determinant of evolutionary change according to game theory [41]. Indeed, mathematical models of somatic evolution routinely consider each cancer cell as a unit [4]. Considering these potential drawbacks, there must be some overriding benefit that intracellular coupling could offer tumors. A plausible beneficiary of coupling is metabolism and pH regulation, two intertwined processes that are fundamental to cancer growth and generate large fluxes of solutes [42,43]. There have been reports for metabolic cooperation between cancer cells, in which substrates are made available to those cells that cannot synthesize them [44]. More recently, we have demonstrated bulk fluxes of lactate [36] and HCO3 ions [35] across gap junctions in spheroids grown form pancreatic ductal adenocarcinoma (PDAC) cell lines.

3. Glycolytic Metabolism Establishes Radial Concentration Gradients of Lactate and Bicarbonate in Solid Tumors

Cancer cells require a steady flow of ATP to service the demands placed by proliferation, invasion and metastasis. Metabolic reprogramming to a glycolytic phenotype is common in many cancer cells (Warburg effect) [42,43,45,46], although a degree of mitochondrial respiration will also take place if oxygen is adequate [43,47]. The end-product of glycolytic metabolism is the ionized form of lactic acid [48], which must be removed from the cell, particularly if there is insufficient mitochondrial capacity to metabolize lactate further in the Krebs cycle. Typically, cancer cell lines produce lactate at a rate of several millimoles per minute per liter of cell volume (mM/min) [49,50,51,52,53], but cancer cell lines can attain rates as high as 20 mM/min [54]. The sole means by which lactate can cross biological membranes is in co-transport with H+ ions aboard H+-monocarboxylate transporter (MCTs) [55,56] or, to a much lesser degree, as lactic acid crossing the lipid bilayer; either way, membrane permeation is electroneutral. Once outside the cell, lactic acid dissociates, and its ionic products diffuse towards blood capillaries to complete the process of venting (Figure 2A).
Membrane permeation and extracellular diffusion can be considered as resistors in series. For cancer cells that are deep at the hypoxic core of solid tumors, the diffusive resistance will be substantial, and no amount of MCT expression could overcome this. Rate-limiting removal of the products of glycolytic metabolism will lead to their retention at the hypoxic core. The extent of this can be simulated mathematically (Figure 2B,C). The equations used in the model are presented in the Appendix A, together with an explanation of the variables and assumptions. For illustrative purposes, we assume the extracellular space occupies 25% of tumor volume, has an intrinsic buffering capacity of 3 mM/min, and its tortuosity reduces ionic diffusivities by a factor of two ([57,58,59,60]). Thus, extracellular lactate, HCO3 and CO2 diffusion coefficients were 650, 900 and 1200 µm2/s, respectively. A similar tortuosity is modelled for the intracellular space; and the apparent H+ diffusion coefficient taken as 200 µm2/s [61]. MCT activity in cells, radius 7 µm and intracellular intrinsic buffering capacity of 15 mM/pH [62], was set to a high level, equivalent to a permeability to lactic acid of 103 µm/s, and lactic acid production fixed at 20 mM/min [54]. CO2/HCO3 buffer is included in the modelling, with freely-permeating CO2 (permeability constant 103 µm/s) [57], and to facilitate its chemical reactions, exofacial carbonic anhydrase (CA) activity increased spontaneous CO2 hydration rate by a factor of 10 (e.g., CAIX/CAXII) [59,60]. The geometry of the modelled tissue is assumed spherical of radius 250 µm to encompass both normoxic regions at the periphery and hypoxic areas at the core. Concentrations at the boundary are set to levels normally found in blood, i.e., CO2 at 1.2 mM, HCO3 at 24 mM, hence extracellular pH (pHe) at 7.4, and negligible lactate. In the model, cells at the tissue boundary were able to fully metabolize lactate (i.e., intracellular lactate is set to nil) and import HCO3 ions by active transport (i.e., intracellular [HCO3] set at 13.5 mM).
Lactic acid produced glycolytically by cancer cells (‘source’) is vented to the boundary of the tissue (‘sink’), down the route of least resistance. Essentially all traffic involves permeation across the membrane of the glycolytic cell, followed by a meandering across the extracellular compartment. Results of the simulation indicate substantial gradients of lactate in both the intra- and extracellular spaces at the steady-state. Retention of lactate is also linked to intra- and extracellular acidity, and hence a radial gradient of [HCO3]. This gradient is notably smaller than the lactate gradient because CO2/HCO3 is not the only acid-neutralizing buffer present inside cells; a sizable pool of buffering is provided by intrinsic buffers (15 mM/pH). The substantial build-up of intra- and extracellular lactate, and the depletion of HCO3 at core regions is consistent with the retention of lactic acid in tumors in vivo [63]. This will likely inhibit glycolysis by the allosteric inhibitory effect of H+ ions of glycolytic enzymes (e.g., phosphofructokinase [64]), and trans-inhibition by lactate. Such a slowing of energy provision would be an Achilles heel to proliferating cancer cells. In response to this challenge, cancer cells are required to regulate their pH, and the canonical mechanism involves active transport [65]. An alternative means of reducing lactic acid build-up, without compromising metabolic rate, is to increase the diffusive capacity for venting lactate and HCO3 ions.

4. Cancer Cells Can Exchange Lactate and Bicarbonate Ions through Gap Junctions

In the modeled scenario shown in Figure 2A, lactate (and HCO3) can only diffuse radially along the extracellular path. In this system, the overall diffusivity for lactate is low, and a large radial gradient is needed to drive an efflux that will balance metabolic production. Transport capacity in the extracellular space is limited because this is the smaller of the two tissue compartments, and after accounting for tortuosity, the diffusion coefficients of small solutes, like lactate and bicarbonate, are of the order of ~103 µm2/s. Lactate and bicarbonate ions are also transported in a chemically converted form (lactic acid and CO2) but this will be a smaller component of flux by >3 and >1 orders of magnitude, respectively. Any means of increasing effective diffusivity would benefit venting, and one way of achieving this is to engage in diffusion across the intracellular syncytium [35,36].
The substantial radial gradient of lactate in the intracellular space provides the energy to drive a dissipating flux, but only if there is sufficient coupling between cells (Figure 3A). On the basis of FRAP experiments, junctionally-coupled CRC and PDAC monolayers can have a permeability constant to calcein of the order of 0.03 µm2/s [36]. Anions such as lactate and bicarbonate are smaller than calcein (7- and 10-fold respectively), thus their permeability can be expected to be two orders of magnitude higher [38,39,66], i.e., in the order of ~2 and ~3 µm/s, respectively. The activity of MCT may produce lactic acid permeability as high as 103 µm2/s, but this refers to the transport of lactic acid, which is >1000-times less abundant than its ionized conjugate. Lactate transport across the syncytial cytoplasmic volume could be described with an effective diffusion coefficient, Deff, given by Equation (3), where D is the diffusion coefficient between permeation barriers of permeability constant P, spaced at intervals d, i.e., approximated by cell diameter:
1 D e f f = 1 D + 1 d × P
Small solutes like lactate may be expected to have effective diffusion coefficients in the cytoplasmic syncytium of ~30 µm2/s. Although this is more than an order of magnitude lower than diffusivity in the extracellular space, the ensuing intracellular flux can be meaningful because of the large volume of the intracellular compartment. Early studies have shown that low pHi uncouples gap junctions [67,68], but more recent work has shown that gap junctions made of Cx43, for example, have a pH optimum that is slightly acidic relative to resting pHi [69,70]; consequently, coupling by Cx43 gap junctions is expected to strengthen as pHi falls modestly, before cells become uncoupled at unphysiologically low pHi.
Adding junctional lactate permeability to the mathematical model (Figure 3A) significantly improves lactic acid venting and reduces lactate build-up at the core (Figure 3B). However, the diffusive dissipation of lactate, a basic substance, will have an effect on pH [36]. At steady-state, lactate-only junctional efflux leaves behind H+ ions that acidify the intracellular space at the core of the tissue (Figure 3D). This generates a substantial radial gradient of [HCO3] (Figure 3C).
To prevent this unwarranted acidification, junctions must also conduct HCO3 ions, which would be expected to have a Deff, in this example, of 40 µm2/s. As shown in the results of a model with junctional lactate and HCO3 permeability (Figure 4A), lactate build-up at the core can be reduced by a fifth (Figure 4B), without significantly acidifying the core (Figure 4C), thanks to a concurrent dissipation of any [HCO3] gradient that develops in response to lactate transport.

5. Concluding Remarks

In this article, we describe a metabolic function for gap junctional coupling between cancer cells in solid tumors. Gradients of hypoxia in solid tumors also produce gradients of other solutes, such as lactate and HCO3. These gradients are the driving forces for diffusion, but transport is only possible if there is adequate conductance. If this were to take place in the intracellular compartment of tumors, it must involve an adequate degree of coupling by gap junctions. In at least some cancers, such coupling is possible if connexin proteins are correctly targeted to points of cell-cell contact. Once a low-resistance pathway is established, it can lead to the diffusive exchange of solutes. Quantitatively the largest fluxes will involve those substances that have large gradients, such as lactate and HCO3 anions. A lactate/HCO3 exchange could benefit both of the coupled cells, as illustrated in Figure 5.
For the more hypoxic cell, gap junctions offer a route for lactate efflux, which reduces the retention of this glycolytic end-product. This was shown experimentally in spheroids of Colo357 cells, a Cx43-coupled PDAC cell line [36]. The recipient (more normoxic) cell can benefit from this lactate delivery by obtaining ATP from the mitochondrial respiration of this substrate [71]. Any excess lactate that is not respired can leave the normoxic cell by MCT, which is a means of removing H+ ions, another benefit to the recipient shown experimentally in PDAC spheroids [36]. Efflux of lactate, a base, from the hypoxic cell could potentially exacerbate its intracellular acidosis, but this is curtailed by HCO3 delivery from the more alkaline normoxic cell, where HCO3 is more plentiful. To sustain HCO3 delivery, the normoxic cell can engage in active uptake, e.g., by Na+-driven HCO3 cotransporters of the SLC4 gene family [35]. From the viewpoint of the hypoxic cell, this active uptake mechanism is an example of remote secondary active transport that we described in PDAC spheroids [35]. Under this arrangement, the energy cost for delivering HCO3 to the hypoxic cell is paid by the normoxic cell. In this example of barter, junctionally coupled cancer cells can benefit from the instantaneous and complete ionic trading.

Funding

This research was funded by the European Research Council, SURVIVE #723997.

Acknowledgments

We acknowledge the many useful discussions with Tobias Dovmark, Alzbeta Hulikova and Sir Walter Bodmer.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

The concentration dynamics of intra- and extracellular solutes were modelled with a diffusion-reaction partial differential equation [57] in a sphere of radius R = 250 µm, which represents a tumor mass that is sufficient for hosting radial gradients:
u t = D · 1 r 2 · x ( r 2 · u x ) + F ( u )
Variable r is the radial distance from the core, u is a vector of concentrations, D is a vector of effective diffusion coefficients, and F is a function describing chemical reactions and membrane permeation.
Vector u describes concentrations in the following order: (1) [H+]i, (2) [HCO3]i, (3) [CO2]i, (4) [Lac]i, (5) [HLac]i, (6) [H+]e, (7) [HCO3]e, (8) [CO2]e, (9) [Lac]e, (10) [HLac]e, where subscripts i and e refer to the intra- and extracellular compartments. The initial concentrations, u0, were 10−7.2 M, 13.5 mM, 1.2 mM, 0 mM, 0 mM, 10−7.4 M, 22 mM, 1.2 mM, 0 mM, 0 mM. The fractional volume of the combined intracellular space vi was 0.75. Cells were assumed to have a radius of 7 µm and thus a surface area/volume radio ρ of 3/7 µm−1.
Intracellular [H+] reactions comprised CA-catalyzed reversible CO2 hydration, reversible lactic acid dissociation, and buffering by intracellular intrinsic buffers:
u 1 t = C A i · ( k h · u 3 k r · u 1 · u 2 ) + k d · u 5 k p · u 1 · u 4 1 + β i / ( ln ( 10 ) · u 1 )
The spontaneous CO2 hydration constant kh was set to 0.18 s−1. The rate constant for the reverse reaction (kr) was thus kh/KCO2, where KCO2 is the equilibrium constant equal to 10−6.15 M. Intracellular CA (CAi) activity was set to 10-fold to provide fast reaction kinetics. Lactate protonation rate constant kp was set to a fast value, 10−10 M−1 s−1, and the reverse rate (kd) equaled kp×Klac, where Klac is the lactic acid dissociation constant, 10−3.86 M. Intracellular buffering capacity βi was set to 15 mM/pH. Intracellular HCO3 reactions comprised CA-catalyzed reversible CO2 hydration:
u 2 t = C A i · ( k h · u 3 k r · u 1 · u 2 )
Intracellular CO2 reactions comprised CA-catalyzed reversible CO2 hydration, and diffusive entry from the extracellular compartment:
u 3 t = C A i · ( k r · u 1 · u 2 k h · u 3 ) +   ρ · P C O 2 · ( u 8 u 3 )
Membrane CO2 permeability PCO2 was set to a value that supported rapid exchange, 103 µm/s. Intracellular lactate reactions comprised reversible lactic acid dissociation:
u 4 t = k d · u 5 k p · u 1 · u 4
Intracellular lactic acid reactions comprised reversible lactic acid dissociation, glycolytic lactic acid production throughout the structure except in the outer rim of thickness δ, and MCT-facilitated entry form the extracellular compartment:
u 5 t = k p · u 1 · u 4 k d · u 5 + J l a c · ( 1 e x p ( R r δ ) ) + ρ · P m c t · ( u 10 u 5 )
The effective membrane lactic acid permeability Pmct was set to a value that supported fast exchange, 103 µm/s, consistent with a high level of MCT activity. The metabolic rate Jlac was set to 20 mM/min. Thickness of the rim δ was set to 10 µm. Extracellular [H+] reactions comprised CA-catalyzed reversible CO2 hydration and reversible lactic acid dissociation, and buffering by extracellular intrinsic buffers:
u 6 t = C A e · ( k h · u 8 k r · u 6 · u 7 ) + k d · u 10 k p · u 6 · u 9 1 + β e / ( ln ( 10 ) · u 6 )
Extracellular CA (CAe) activity was set to 10-fold to provide fast reaction kinetics. Extracellular buffering capacity βe was set to 3 mM/pH. Extracellular HCO3 reactions comprised CA-catalyzed reversible CO2 hydration:
u 7 t = C A e · ( k h · u 8 k r · u 1 · u 2 )
Extracellular CO2 reactions comprised CA-catalyzed reversible CO2 hydration and diffusive entry from the intracellular compartment, with consideration of the intra-to-extracellular volume ratio:
u 8 t = C A e · ( k r · u 6 · u 7 k h · u 8 ) +   v i 1 v i · ρ · P C O 2 · ( u 3 u 8 )
Extracellular lactate reactions comprised reversible lactic acid dissociation
u 9 t = k d · u 10 k p · u 6 · u 9
Extracellular lactic acid reactions comprised reversible lactic acid dissociation, and MCT-facilitated entry form the intracellular compartment, with consideration of the intra-to-extracellular volume ratio:
u 10 t = k p · u 6 · u 10 k d · u 9 + v i 1 v i · ρ · P m c t · ( u 5 u 10 )
Diffusion coefficients for extracellular H+, HCO3, CO2, Lac and HLac were set to half of their values in aqueous media to account for a two-fold tortuosity imposed by the extracellular space, i.e., 6000, 900, 1200, 650 and 650 µm2/s, respectively. The effective intracellular diffusion coefficients, as defined by Equation (3), were set to values that depended on the modelling scenarios. In the absence of junctional coupling, as depicted in Figure 2A, intracellular diffusivity for ions was set to zero, and 1200 µm2/s in the case of membrane-permeable CO2. In the modeling scenario shown in Figure 4A, gap junctions were permeable to all ionic species. In this case, the permeability to HCO3, lactate and lactic acid were in proportion to their diffusivity and set to 3, 2.17 and 2.17 µm/s, giving an effective permeability of 40, 29 and 29 µm2/s. The diffusive properties of H+ ions are unique, as they are set by mobile buffers [72]. Assuming a mobile buffer diffusivity of 200 µm2/s, consistent with being larger molecules than lactate, and a permeability of 0.67 µm/s, their effective diffusivity was 8.9 µm2/s. In the case depicted in Figure 3A, gap junctions were permeable to lactic acid and lactate only, but not H+ or HCO3 ions.
The boundary conditions at r = R for extracellular solutes were set to provide direct coupling with a well-stirred compartment representing blood flow, i.e., concentrations were clamped to levels defined in the initial conditions vector u0. Similarly, the boundary conditions for intracellular solutes was set to values defined in u0 to simulate a scenario of good pHi regulation (and hence constancy of [HCO3]) and zero lactate (i.e., metabolic degradation of any incoming lactate ions). Thus, boundary conditions took the general form:
u i | r = R = u i 0
The equations were solved with the pdepe solver (MATLAB) and presented as radial profiles once the system attained a steady-state.

References

  1. Graeber, D. Toward an Anthropological Theory of Value: The False Coin of our Dreams; Palgrave: Basingstoke, UK, 2001. [Google Scholar]
  2. Humphrey, C. Barter and Economic Disintegration. Man 1985, 20, 48–72. [Google Scholar] [CrossRef]
  3. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed]
  4. Gatenby, R.A.; Gillies, R.J.; Brown, J.S. Evolutionary dynamics of cancer prevention. Nat. Rev. Cancer 2010, 10, 526–527. [Google Scholar] [CrossRef] [PubMed]
  5. Kalluri, R.; Zeisberg, M. Fibroblasts in cancer. Nat. Rev. Cancer 2006, 6, 392–401. [Google Scholar] [CrossRef] [PubMed]
  6. Hanahan, D.; Coussens, L.M. Accessories to the crime: Functions of cells recruited to the tumor microenvironment. Cancer Cell 2012, 21, 309–322. [Google Scholar] [CrossRef] [PubMed]
  7. Loewenstein, W.R.; Socolar, S.J.; Higashino, S.; Kanno, Y.; Davidson, N. Intercellular Communication: Renal, Urinary Bladder, Sensory, and Salivary Gland Cells. Science 1965, 149, 295–298. [Google Scholar] [CrossRef] [PubMed]
  8. Evans, W.H. Cell communication across gap junctions: A historical perspective and current developments. Biochem. Soc. Trans. 2015, 43, 450–459. [Google Scholar] [CrossRef]
  9. Goodenough, D.A. Bulk isolation of mouse hepatocyte gap junctions. Characterization of the principal protein, connexin. J. Cell Biol. 1974, 61, 557–563. [Google Scholar] [CrossRef]
  10. Aasen, T.; Mesnil, M.; Naus, C.C.; Lampe, P.D.; Laird, D.W. Gap junctions and cancer: Communicating for 50 years. Nat. Rev. Cancer 2016, 16, 775–788. [Google Scholar] [CrossRef]
  11. Duffy, H.S.; Fort, A.G.; Spray, D.C. Cardiac connexins: Genes to nexus. Adv. Cardiol. 2006, 42, 1–17. [Google Scholar] [CrossRef]
  12. Kumar, N.M.; Gilula, N.B. The gap junction communication channel. Cell 1996, 84, 381–388. [Google Scholar] [CrossRef]
  13. Aad, G.; Abajyan, T.; Abbott, B.; Abdallah, J.; Abdel Khalek, S.; Abdelalim, A.A.; Abdinov, O.; Aben, R.; Abi, B.; Abolins, M.; et al. Search for magnetic monopoles in sqrt[s]=7 TeV pp collisions with the ATLAS detector. Phys. Rev. Lett. 2012, 109, 261803. [Google Scholar] [CrossRef] [PubMed]
  14. Graham, S.V.; Jiang, J.X.; Mesnil, M. Connexins and Pannexins: Important Players in Tumorigenesis, Metastasis and Potential Therapeutics. Int. J. Mol. Sci. 2018, 19, 1645. [Google Scholar] [CrossRef] [PubMed]
  15. Solan, J.L.; Hingorani, S.R.; Lampe, P.D. Changes in connexin43 expression and localization during pancreatic cancer progression. J. Membr. Biol. 2012, 245, 255–262. [Google Scholar] [CrossRef] [PubMed]
  16. Constantin, B.; Cronier, L. Involvement of gap junctional communication in myogenesis. Int. Rev. Cytol. 2000, 196, 1–65. [Google Scholar] [PubMed]
  17. Loewenstein, W.R.; Kanno, Y. Intercellular communication and the control of tissue growth: Lack of communication between cancer cells. Nature 1966, 209, 1248–1249. [Google Scholar] [CrossRef] [PubMed]
  18. Johnson, R.G.; Sheridan, J.D. Junctions between cancer cells in culture: Ultrastructure and permeability. Science 1971, 174, 717–719. [Google Scholar] [CrossRef]
  19. Borek, C.; Higashino, S.; Loewenstein, W.R. Intercellular communication and tissue growth: IV. Conductance of membrane junctions of normal and cancerous cells in culture. J. Membr. Biol. 1969, 1, 274–293. [Google Scholar] [CrossRef]
  20. Mesnil, M.; Crespin, S.; Avanzo, J.L.; Zaidan-Dagli, M.L. Defective gap junctional intercellular communication in the carcinogenic process. Biochim. Biophys. Acta 2005, 1719, 125–145. [Google Scholar] [CrossRef] [Green Version]
  21. McLachlan, E.; Shao, Q.; Wang, H.L.; Langlois, S.; Laird, D.W. Connexins act as tumor suppressors in three-dimensional mammary cell organoids by regulating differentiation and angiogenesis. Cancer Res. 2006, 66, 9886–9894. [Google Scholar] [CrossRef]
  22. Mesnil, M. Connexins and cancer. Biol. Cell 2002, 94, 493–500. [Google Scholar] [CrossRef] [Green Version]
  23. Jiang, J.X.; Penuela, S. Connexin and pannexin channels in cancer. BMC Cell Biol. 2016, 17 (Suppl. 1), 12. [Google Scholar] [CrossRef]
  24. Lin, J.H.; Takano, T.; Cotrina, M.L.; Arcuino, G.; Kang, J.; Liu, S.; Gao, Q.; Jiang, L.; Li, F.; Lichtenberg-Frate, H.; et al. Connexin 43 enhances the adhesivity and mediates the invasion of malignant glioma cells. J. Neurosci. Off. J. Soc. Neurosci. 2002, 22, 4302–4311. [Google Scholar] [CrossRef]
  25. Pollmann, M.A.; Shao, Q.; Laird, D.W.; Sandig, M. Connexin 43 mediated gap junctional communication enhances breast tumor cell diapedesis in culture. Breast Cancer Res. 2005, 7, R522–R534. [Google Scholar] [CrossRef]
  26. Polusani, S.R.; Kalmykov, E.A.; Chandrasekhar, A.; Zucker, S.N.; Nicholson, B.J. Cell coupling mediated by connexin 26 selectively contributes to reduced adhesivity and increased migration. J. Cell Sci. 2016, 129, 4399–4410. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Naus, C.C.; Laird, D.W. Implications and challenges of connexin connections to cancer. Nat. Rev. Cancer 2010, 10, 435–441. [Google Scholar] [CrossRef]
  28. Ezumi, K.; Yamamoto, H.; Murata, K.; Higashiyama, M.; Damdinsuren, B.; Nakamura, Y.; Kyo, N.; Okami, J.; Ngan, C.Y.; Takemasa, I.; et al. Aberrant expression of connexin 26 is associated with lung metastasis of colorectal cancer. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2008, 14, 677–684. [Google Scholar] [CrossRef]
  29. Nomura, S.; Maeda, K.; Noda, E.; Inoue, T.; Fukunaga, S.; Nagahara, H.; Hirakawa, K. Clinical significance of the expression of connexin26 in colorectal cancer. J. Exp. Clin. Cancer Res. 2010, 29, 79. [Google Scholar] [CrossRef] [Green Version]
  30. Thuringer, D.; Berthenet, K.; Cronier, L.; Solary, E.; Garrido, C. Primary tumor- and metastasis-derived colon cancer cells differently modulate connexin expression and function in human capillary endothelial cells. Oncotarget 2015, 6, 28800–28815. [Google Scholar] [CrossRef] [Green Version]
  31. Bigelow, K.; Nguyen, T.A. Increase of gap junction activities in SW480 human colorectal cancer cells. BMC Cancer 2014, 14, 502. [Google Scholar] [CrossRef]
  32. Valiunas, V.; Polosina, Y.Y.; Miller, H.; Potapova, I.A.; Valiuniene, L.; Doronin, S.; Mathias, R.T.; Robinson, R.B.; Rosen, M.R.; Cohen, I.S.; et al. Connexin-specific cell-to-cell transfer of short interfering RNA by gap junctions. J. Physiol. 2005, 568, 459–468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zhang, Y.W.; Kaneda, M.; Morita, I. The gap junction-independent tumor-suppressing effect of connexin 43. J. Biol. Chem. 2003, 278, 44852–44856. [Google Scholar] [CrossRef]
  34. Al-Ghadban, S.; Kaissi, S.; Homaidan, F.R.; Naim, H.Y.; El-Sabban, M.E. Cross-talk between intestinal epithelial cells and immune cells in inflammatory bowel disease. Sci. Rep. 2016, 6, 29783. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Dovmark, T.H.; Hulikova, A.; Niederer, S.A.; Vaughan-Jones, R.D.; Swietach, P. Normoxic cells remotely regulate the acid-base balance of cells at the hypoxic core of connexin-coupled tumor growths. Faseb. J. Off. Publ. Fed. Am. Soc. Exp. Biol. 2018, 32, 83–96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Dovmark, T.H.; Saccomano, M.; Hulikova, A.; Alves, F.; Swietach, P. Connexin-43 channels are a pathway for discharging lactate from glycolytic pancreatic ductal adenocarcinoma cells. Oncogene 2017, 36, 4538–4550. [Google Scholar] [CrossRef] [Green Version]
  37. Kyo, N.; Yamamoto, H.; Takeda, Y.; Ezumi, K.; Ngan, C.Y.; Terayama, M.; Miyake, M.; Takemasa, I.; Ikeda, M.; Doki, Y.; et al. Overexpression of connexin 26 in carcinoma of the pancreas. Oncol. Rep. 2008, 19, 627–631. [Google Scholar] [CrossRef]
  38. Weber, P.A.; Chang, H.C.; Spaeth, K.E.; Nitsche, J.M.; Nicholson, B.J. The permeability of gap junction channels to probes of different size is dependent on connexin composition and permeant-pore affinities. Biophys. J. 2004, 87, 958–973. [Google Scholar] [CrossRef]
  39. Kanaporis, G.; Brink, P.R.; Valiunas, V. Gap junction permeability: Selectivity for anionic and cationic probes. Am. J. Physiol. Cell Physiol. 2011, 300, C600–C609. [Google Scholar] [CrossRef]
  40. Gillies, R.J.; Verduzco, D.; Gatenby, R.A. Evolutionary dynamics of carcinogenesis and why targeted therapy does not work. Nat. Rev. Cancer 2012, 12, 487–493. [Google Scholar] [CrossRef] [Green Version]
  41. Gatenby, R.A.; Vincent, T.L.; Gillies, R.J. Evolutionary dynamcs in carcinogenesis. Math. Models Methods Appl. Sci. 2005, 15, 1619–1638. [Google Scholar] [CrossRef]
  42. Gatenby, R.A.; Gillies, R.J. Why do cancers have high aerobic glycolysis? Nat. Rev. Cancer 2004, 4, 891–899. [Google Scholar] [CrossRef] [PubMed]
  43. Ward, P.S.; Thompson, C.B. Metabolic reprogramming: A cancer hallmark even warburg did not anticipate. Cancer Cell 2012, 21, 297–308. [Google Scholar] [CrossRef] [PubMed]
  44. Subak-Sharpe, H.; Burk, R.R.; Pitts, J.D. Metabolic co-operation between biochemically marked mammalian cells in tissue culture. J. Cell Sci. 1969, 4, 353–367. [Google Scholar] [PubMed]
  45. Koppenol, W.H.; Bounds, P.L.; Dang, C.V. Otto Warburg’s contributions to current concepts of cancer metabolism. Nat. Rev. Cancer 2011, 11, 325–337. [Google Scholar] [CrossRef]
  46. Vander Heiden, M.G.; Cantley, L.C.; Thompson, C.B. Understanding the Warburg effect: The metabolic requirements of cell proliferation. Science 2009, 324, 1029–1033. [Google Scholar] [CrossRef] [PubMed]
  47. Hu, X.; Chao, M.; Wu, H. Central role of lactate and proton in cancer cell resistance to glucose deprivation and its clinical translation. Signal Transduct. Target. Ther. 2017, 2, 16047. [Google Scholar] [CrossRef]
  48. Hochachka, P.W.; Mommsen, T.P. Protons and anaerobiosis. Science 1983, 219, 1391–1397. [Google Scholar] [CrossRef] [PubMed]
  49. Wu, H.; Ying, M.; Hu, X. Lactic acidosis switches cancer cells from aerobic glycolysis back to dominant oxidative phosphorylation. Oncotarget 2016, 7, 40621–40629. [Google Scholar] [CrossRef]
  50. Yizhak, K.; Le Devedec, S.E.; Rogkoti, V.M.; Baenke, F.; de Boer, V.C.; Frezza, C.; Schulze, A.; van de Water, B.; Ruppin, E. A computational study of the Warburg effect identifies metabolic targets inhibiting cancer migration. Mol. Syst. Biol. 2014, 10, 744. [Google Scholar] [CrossRef]
  51. Bailey, K.M.; Wojtkowiak, J.W.; Cornnell, H.H.; Ribeiro, M.C.; Balagurunathan, Y.; Hashim, A.I.; Gillies, R.J. Mechanisms of buffer therapy resistance. Neoplasia 2014, 16, 354–364.e1-3. [Google Scholar] [CrossRef]
  52. Wu, M.; Neilson, A.; Swift, A.L.; Moran, R.; Tamagnine, J.; Parslow, D.; Armistead, S.; Lemire, K.; Orrell, J.; Teich, J.; et al. Multiparameter metabolic analysis reveals a close link between attenuated mitochondrial bioenergetic function and enhanced glycolysis dependency in human tumor cells. Am. J. Physiol. Cell Physiol. 2007, 292, C125–C136. [Google Scholar] [CrossRef] [PubMed]
  53. Giang, A.H.; Raymond, T.; Brookes, P.; de Mesy Bentley, K.; Schwarz, E.; O’Keefe, R.; Eliseev, R. Mitochondrial dysfunction and permeability transition in osteosarcoma cells showing the Warburg effect. J. Biol. Chem. 2013, 288, 33303–33311. [Google Scholar] [CrossRef] [PubMed]
  54. Cheng, G.; Zielonka, J.; McAllister, D.; Tsai, S.; Dwinell, M.B.; Kalyanaraman, B. Profiling and targeting of cellular bioenergetics: Inhibition of pancreatic cancer cell proliferation. Br. J. Cancer 2014, 111, 85–93. [Google Scholar] [CrossRef] [PubMed]
  55. Halestrap, A.P.; Meredith, D. The SLC16 gene family-from monocarboxylate transporters (MCTs) to aromatic amino acid transporters and beyond. Pflugers Arch. 2004, 447, 619–628. [Google Scholar] [CrossRef] [PubMed]
  56. Ullah, M.S.; Davies, A.J.; Halestrap, A.P. The plasma membrane lactate transporter MCT4, but not MCT1, is up-regulated by hypoxia through a HIF-1alpha-dependent mechanism. J. Biol. Chem. 2006, 281, 9030–9037. [Google Scholar] [CrossRef] [PubMed]
  57. Hulikova, A.; Swietach, P. Rapid CO2 permeation across biological membranes: Implications for CO2 venting from tissue. Faseb. J. Off. Publ. Fed. Am. Soc. Exp. Biol. 2014, 28, 2762–2774. [Google Scholar] [CrossRef] [PubMed]
  58. Swietach, P.; Hulikova, A.; Vaughan-Jones, R.D.; Harris, A.L. New insights into the physiological role of carbonic anhydrase IX in tumour pH regulation. Oncogene 2010, 29, 6509–6521. [Google Scholar] [CrossRef] [Green Version]
  59. Swietach, P.; Patiar, S.; Supuran, C.T.; Harris, A.L.; Vaughan-Jones, R.D. The role of carbonic anhydrase 9 in regulating extracellular and intracellular ph in three-dimensional tumor cell growths. J. Biol. Chem. 2009, 284, 20299–20310. [Google Scholar] [CrossRef]
  60. Swietach, P.; Wigfield, S.; Cobden, P.; Supuran, C.T.; Harris, A.L.; Vaughan-Jones, R.D. Tumor-associated carbonic anhydrase 9 spatially coordinates intracellular pH in three-dimensional multicellular growths. J. Biol. Chem. 2008, 283, 20473–20483. [Google Scholar] [CrossRef]
  61. Hulikova, A.; Black, N.; Hsia, L.T.; Wilding, J.; Bodmer, W.F.; Swietach, P. Stromal uptake and transmission of acid is a pathway for venting cancer cell-generated acid. Proc. Natl. Acad. Sci. USA 2016, 113, E5344–E5353. [Google Scholar] [CrossRef]
  62. Hulikova, A.; Harris, A.L.; Vaughan-Jones, R.D.; Swietach, P. Regulation of intracellular pH in cancer cell lines under normoxia and hypoxia. J. Cell. Physiol. 2013, 228, 743–752. [Google Scholar] [CrossRef] [PubMed]
  63. Griffiths, J.R.; Stevens, A.N.; Iles, R.A.; Gordon, R.E.; Shaw, D. 31P-NMR investigation of solid tumours in the living rat. Biosci. Rep. 1981, 1, 319–325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Webb, B.A.; Forouhar, F.; Szu, F.E.; Seetharaman, J.; Tong, L.; Barber, D.L. Structures of human phosphofructokinase-1 and atomic basis of cancer-associated mutations. Nature 2015, 523, 111–114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Swietach, P.; Vaughan-Jones, R.D.; Harris, A.L. Regulation of tumor pH and the role of carbonic anhydrase 9. Cancer Metastasis Rev. 2007, 26, 299–310. [Google Scholar] [CrossRef] [PubMed]
  66. Swietach, P.; Vaughan-Jones, R.D. Novel method for measuring junctional proton permeation in isolated ventricular myocyte cell pairs. Am. J. Physiol. Heart Circ. Physiol. 2004, 287, H2352–H2363. [Google Scholar] [CrossRef] [Green Version]
  67. Turin, L.; Warner, A. Carbon dioxide reversibly abolishes ionic communication between cells of early amphibian embryo. Nature 1977, 270, 56–57. [Google Scholar] [CrossRef] [PubMed]
  68. Spray, D.C.; Harris, A.L.; Bennett, M.V. Gap junctional conductance is a simple and sensitive function of intracellular pH. Science 1981, 211, 712–715. [Google Scholar] [CrossRef] [PubMed]
  69. Swietach, P.; Rossini, A.; Spitzer, K.W.; Vaughan-Jones, R.D. H+ ion activation and inactivation of the ventricular gap junction: A basis for spatial regulation of intracellular pH. Circ. Res. 2007, 100, 1045–1054. [Google Scholar] [CrossRef]
  70. Garciarena, C.D.; Malik, A.; Swietach, P.; Moreno, A.P.; Vaughan-Jones, R.D. Distinct moieties underlie biphasic H(+) gating of connexin43 channels, producing a pH optimum for intercellular communication. Faseb J. Off. Publ. Fed. Am. Soc. Exp. Biol. 2018, 32, 1969–1981. [Google Scholar] [CrossRef]
  71. Sonveaux, P.; Vegran, F.; Schroeder, T.; Wergin, M.C.; Verrax, J.; Rabbani, Z.N.; De Saedeleer, C.J.; Kennedy, K.M.; Diepart, C.; Jordan, B.F.; et al. Targeting lactate-fueled respiration selectively kills hypoxic tumor cells in mice. J. Clin. Investig. 2008, 118, 3930–3942. [Google Scholar] [CrossRef] [Green Version]
  72. Swietach, P.; Zaniboni, M.; Stewart, A.K.; Rossini, A.; Spitzer, K.W.; Vaughan-Jones, R.D. Modelling intracellular H(+) ion diffusion. Prog. Biophys. Mol. Biol. 2003, 83, 69–100. [Google Scholar] [CrossRef]
Figure 1. (A) Monolayer of SNU1235 colorectal cancer cells loaded with calcein (AM-ester; 7 min/4 µM) and imaged for fluorescence on a Zeiss LSM 700 confocal system (excitation 488 nm/emission > 510 nm). Red outline indicates cell that was selected for bleaching; blue outline highlights its direct neighbors. (B) Bleaching a central cell produces a gradient for calcein diffusion which can be tracked fluorescently. (C) Fluorescence in bleached cell recovers within 3 min, indicative of diffusion through gap junctions. (D) Recovery of fluorescence in the bleached cell is used to calculate junctional calcein permeability (Equation (2)), determined here to be 0.0491 µm/s. Note that fluorescence in neighboring cell decreases as signal from the bleached cell recovers.
Figure 1. (A) Monolayer of SNU1235 colorectal cancer cells loaded with calcein (AM-ester; 7 min/4 µM) and imaged for fluorescence on a Zeiss LSM 700 confocal system (excitation 488 nm/emission > 510 nm). Red outline indicates cell that was selected for bleaching; blue outline highlights its direct neighbors. (B) Bleaching a central cell produces a gradient for calcein diffusion which can be tracked fluorescently. (C) Fluorescence in bleached cell recovers within 3 min, indicative of diffusion through gap junctions. (D) Recovery of fluorescence in the bleached cell is used to calculate junctional calcein permeability (Equation (2)), determined here to be 0.0491 µm/s. Note that fluorescence in neighboring cell decreases as signal from the bleached cell recovers.
Cancers 11 00117 g001
Figure 2. Modeling the radial distribution of lactate, HCO3 and H+ ions in tissues without junctional coupling, i.e., radial diffusion in extracellular spaces only. (A) Schematic representation of the 3-D spherical geometry representing a solid tumor with a hypoxic core. Results of modelling shown as radial concentration gradients at the steady-state for (B) lactate, (C) HCO3 and (D) pH in the extracellular (red) and intracellular spaces (black). Arrows indicate direction of diffusive flux.
Figure 2. Modeling the radial distribution of lactate, HCO3 and H+ ions in tissues without junctional coupling, i.e., radial diffusion in extracellular spaces only. (A) Schematic representation of the 3-D spherical geometry representing a solid tumor with a hypoxic core. Results of modelling shown as radial concentration gradients at the steady-state for (B) lactate, (C) HCO3 and (D) pH in the extracellular (red) and intracellular spaces (black). Arrows indicate direction of diffusive flux.
Cancers 11 00117 g002
Figure 3. Modeling the radial distribution of lactate, HCO3 and H+ ions in tissues with junctional coupling for lactate only. (A) Schematic representation of the 3-D spherical geometry representing a solid tumor with a hypoxic core. Results of modelling shown as radial concentration gradients at the steady-state for (B) lactate, (C) HCO3 and (D) pH in the intracellular space. Arrows indicate direction of diffusive flux. For comparison, dashed lines show results with no junctional coupling, from Figure 2. Note that diffusion of lactate and lactic acid can take place across the intra- and extracellular spaces. Junctional efflux of lactate (the major permeating species) leaves behind H+ ions, and hence a depletion of HCO3; this generates a large radial gradient of [HCO3].
Figure 3. Modeling the radial distribution of lactate, HCO3 and H+ ions in tissues with junctional coupling for lactate only. (A) Schematic representation of the 3-D spherical geometry representing a solid tumor with a hypoxic core. Results of modelling shown as radial concentration gradients at the steady-state for (B) lactate, (C) HCO3 and (D) pH in the intracellular space. Arrows indicate direction of diffusive flux. For comparison, dashed lines show results with no junctional coupling, from Figure 2. Note that diffusion of lactate and lactic acid can take place across the intra- and extracellular spaces. Junctional efflux of lactate (the major permeating species) leaves behind H+ ions, and hence a depletion of HCO3; this generates a large radial gradient of [HCO3].
Cancers 11 00117 g003
Figure 4. Modeling the radial distribution of lactate, HCO3 and H+ ions in tissues with junctional coupling for all ionic species. (A) Schematic representation of the 3-D spherical geometry representing a solid tumor with a hypoxic core. Results of modelling shown as radial concentration gradients at the steady-state for (B) lactate, (C) HCO3 and (D) pH in the intracellular space. Arrows indicate direction of diffusive flux. For comparison, dashed lines show results with no junctional coupling, from Figure 2. In this system, intracellular lactate diffuses radially out of the core, leaving behind H+ ions. This generates an inverse gradient of intracellular [HCO3], which drives a counter-flux of HCO3 i.e., junctional lactate/HCO3 exchange. Overall, lactate venting improves by a fifth.
Figure 4. Modeling the radial distribution of lactate, HCO3 and H+ ions in tissues with junctional coupling for all ionic species. (A) Schematic representation of the 3-D spherical geometry representing a solid tumor with a hypoxic core. Results of modelling shown as radial concentration gradients at the steady-state for (B) lactate, (C) HCO3 and (D) pH in the intracellular space. Arrows indicate direction of diffusive flux. For comparison, dashed lines show results with no junctional coupling, from Figure 2. In this system, intracellular lactate diffuses radially out of the core, leaving behind H+ ions. This generates an inverse gradient of intracellular [HCO3], which drives a counter-flux of HCO3 i.e., junctional lactate/HCO3 exchange. Overall, lactate venting improves by a fifth.
Cancers 11 00117 g004
Figure 5. Schematic diagram of junctional lactate/HCO3 exchange between two coupled cancer cells, one of which is more hypoxic i.e., further away from a capillary. The barter of lactate for HCO3 can benefit both partners.
Figure 5. Schematic diagram of junctional lactate/HCO3 exchange between two coupled cancer cells, one of which is more hypoxic i.e., further away from a capillary. The barter of lactate for HCO3 can benefit both partners.
Cancers 11 00117 g005

Share and Cite

MDPI and ACS Style

Swietach, P.; Monterisi, S. A Barter Economy in Tumors: Exchanging Metabolites through Gap Junctions. Cancers 2019, 11, 117. https://doi.org/10.3390/cancers11010117

AMA Style

Swietach P, Monterisi S. A Barter Economy in Tumors: Exchanging Metabolites through Gap Junctions. Cancers. 2019; 11(1):117. https://doi.org/10.3390/cancers11010117

Chicago/Turabian Style

Swietach, Pawel, and Stefania Monterisi. 2019. "A Barter Economy in Tumors: Exchanging Metabolites through Gap Junctions" Cancers 11, no. 1: 117. https://doi.org/10.3390/cancers11010117

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