Abstract¶
The classic Boron & De Weer (1976) paper provided the first evidence of active regulation of in cells by an energy-dependent acid-base transporter. These authors also developed a quantitative model — comprising passive fluxes of acid-base equivalents across the cell membrane, intracellular reactions, and an active transport mechanism in the cell membrane (modelled as a proton pump) — to help interpret their measurements of intracellular under perturbations of both extracellular and extracellular . This Physiome paper seeks to make that model, and the experimental conditions under which it was developed, available in a reproducible and well-documented form, along with a software implementation that makes the model easy to use and understand. We have also taken the opportunity to update some of the units used in the original paper, and to provide a few parameter values that were missing in the original paper. Finally, we provide an historical background to the Boron & De Weer (1976) proposal for active regulation and a commentary on subsequent work that has enriched our understanding of this most basic aspect of cellular physiology.
1Introduction¶
In 1976 Boron & De Weer published their landmark paper on “Intracellular transients in squid giant axons caused by , , and metabolic inhibitors” Boron & De Weer, 1976. The authors used a squid giant axon preparation and a mathematical model of buffering and the transport of protons, bicarbonate () and to establish the experimental evidence for active regulation of intracellular () by a transporter in the plasma membrane that — at the expense of energy — either moves acid out of the cell, or base into the cell. Today, we refer to such a transporter generically as an acid-extrusion mechanism. For simplicity, Boron & De Weer modelled it as a proton pump, although the result would have been almost indistinguishable had they modelled it as the uptake of or carbonate (). The paper reported on the consequences of adding and then removing extracellular , (where is ammonium), or the metabolic inhibitors, cyanide, azide and dinitrophenol (DNP).
In the first experiment, following exposure of the cell to elevated and , rapidly enters the cell and intracellular equilibrates with the extracellular , and generates intracellular and via the hydration reaction (). The accumulating results in a rapid fall of (see Figure 1A & Figure 1B). To the extent that the membrane is permeable to as well as to , will initially enter the cell passively, along its electrochemical gradient. Soon, however, the accumulation of intracellular reverses the electrochemical gradient and would be expected to lead to the passive efflux of . This loss of cellular would tend to acidify the cell because — to replenish the lost — additional would enter the cell and form even more and (the passive shuttle). Thus, the expectation was that prolonged exposure to would cause to fall rapidly (passive influx of ) and then to drift more slowly in the acidic direction (passive efflux of ). In fact, Boron & De Weer observed an alkaline drift, leading to the postulate of active extrusion of — or an equivalent process[1] — at a rate that exceeds the passive shuttling by the couple (see Figure 1A & Figure 1C).
Figure 1: changes caused by prolonged exposure of a squid giant axon to extracellular in the bulk solution. (A) Original and traces from figure 1 of BDW. Exposing an axon to extracellular causes a rapid fall in followed by a slow and sustained plateau-phase recovery (i.e., rises). Removing extracellular causes to overshoot its initial resting value. Both the plateau-phase recovery (short arrow) and the overshoot (long arrow) are indicative of net acid extrusion during the period of exposure. (B) Cartoon illustrating the processes underlying the initial, rapid acidification phase in (A). The entry of leads to the intracellular production of (and thus to the observed decay) via the reaction . (C) Cartoon illustrating the processes underlying the plateau-phase alkalinisation in (A). After equilibration across the plasma membrane ( nadir in panel (A)), the slow entry of (or, equivalently, the slow exit of ) — which has always been present but was overwhelmed by the influx of — leads to the consumption of (and thus to the observed slow rise) via the reaction . The newly formed then exits the cell. The observed overshoot is the result of the accumulation of during exposure to extracellular . BDW used the mathematical model to postulate the presence of an active acid-extrusion mechanism that would explain both the observed plateau-phase recovery and the overshoot. (A), modified from Boron & De Weer (1976). (B)-(C), modified from Boron (2010).
Following removal of external , intracellular diffuses out, while intracellular combines with to leave the cell as . Thus, the entire intracellular load associated with entry would be removed, returning to its value before the addition of (or to a slightly lower value, to the extent that had passively exited during the exposure). In fact, Boron & De Weer observed that overshoots its resting value by an amount consistent with the net removal of by the active, acid-extrusion mechanism during the exposure (see Figure 1A & Figure 1C).
In the second experiment, following exposure of the cell to extracellular in the form of (ammonium chloride), the intracellular environment rapidly becomes alkaline as enters and combines with to form (equivalent to the hydration of to form and ). If this were the entire story, then would rise monotonically to a relatively alkaline value, and then the subsequent removal of would cause to fall to precisely its initial value. In fact, Boron & De Weer observed that the exposure to causes to rise rapidly and then fall slowly. Moreover, the subsequent removal of causes to undershoot its original value (see Figure 2A). Thus, Boron & De Weer postulated that, during the exposure, passively enters the cell down its electrochemical gradient. Early during the exposure, this influx would oppose the entry and slightly reduce the increase. Later during the exposure, after intracellular rises to match extracellular (), the continued passive influx of would generate intracellular and . The result would be a slow fall of and a rise in intracellular , the latter leading to the passive exit of (the passive shuttle, see Figure 2A & Figure 2B). The experimental data are consistent with the proposed mathematical model.
Figure 2: changes caused by a short and a long exposure of a squid giant axon to extracellular in the bulk solution. (A) Original and traces from figure 2 of BDW. A short exposure of the axon to extracellular causes a rapid rise in , followed by a decay that modestly undershoots (lower short arrow) its initial resting value upon removal of extracellular . A longer exposure of squid giant axons to extracellular causes a rapid rise in , followed by a slow and sustained decay. Removing extracellular causes to undershoot substantially its initial resting value (long arrow). Both the plateau-phase acidification (upper short arrow) and the undershoot (long arrow) are indicative of net acid loading during the period of exposure. (B) Cartoon illustrating the processes underlying the initial alkalinisation phase in (A) for both short and long exposures to extracellular . The initial entry of leads to the intracellular consumption of (and thus to the observed rise) via the reaction . (C) Cartoon illustrating the processes underlying the plateau-phase acidification during the long exposure in (A). After equilibration across the plasma membrane ( peak in panel (A)), the slow entry of — which has always been present but overwhelmed by the influx of — leads to the production of (and thus to the observed slow decay during the plateau phase) via the reaction . The newly formed then exits the cell. The undershoots observed upon removal of extracellular , during both short and long exposures, are the result of the accumulation of during exposure to extracellular . BDW used the mathematical model to postulate the above sequence of events, including both the plateau-phase acidification and the undershoot. (A), modified from Boron & De Weer (1976). (B)-(C), modified from Boron (2010).
Finally, exposure of the cells — in turn — to cyanide, DNP and azide resulted in intracellular acidosis, consistent with the accumulation of acid metabolites.
In the present paper, we re-formulate the models from Boron & De Weer (1976), henceforth referred to as ‘BDW’, and specify the simulation using the Physiome modelling standards CellML Cuellar et al., 2003 and SED-ML Bergmann et al., 2017 in order to ensure that the model reproduces the graphs in the original paper and that the model is fully curated.[2] Note that this effort requires the specification of some parameters used in BDW’s simulations, but not described in the BDW paper. The curated and annotated model is made available in a form that users can run with OpenCOR[3] to understand the model and to explore the effect of changes in parameter values.
2pH Buffering by Weak Acids and Weak Bases¶
We begin by reviewing a few rudimentary concepts of buffering by weak acids and bases Roos & Boron, 1981Bevensee & Boron, 2013Boron & Boulpaep, 2016 to provide the background for understanding the derivation and implementation of the BDW model.
Buffers.¶
According to Brönsted’s definition Brönsted, 1923, an acid is any substance that can donate a . Conversely, a base is any substance that can accept a . A buffer is any substance that can reversibly consume or produce , thereby minimising changes in .
The dissociation of the uncharged weak acid () to the anionic weak base () is described by the equilibrium reaction:
which is governed by the equilibrium constant[4]
An example is the carbonic acid () dissociation reaction,
The total weak acid concentration, , is the sum of and . Note that is one of the two main unknowns in the BDW model for weak acids.
The dissociation of the cationic weak acid () to the uncharged weak base () is described by the equilibrium reaction,
where the equilibrium constant is
An example is the dissociation reaction,
The total weak base concentration, , is the sum of and . Note that is one of the two main unknowns in the BDW model for weak bases.
The buffer pair.¶
The formation of and from by hydration is given by the equilibrium reaction
where the equilibrium constant is
Taking logarithms of both sides of (5), and recognising from Henry’s law that
where is the solubility coefficient for and is the partial pressure of , we obtain the familiar Henderson-Hasselbalch equation
Here, [5] and .
In terms of the nomenclature above, one might regard as the weak acid [6], and as its conjugate base .
Buffering power ().¶
By definition, is the amount of strong base (e.g., ), or the negative of the amount of strong acid (e.g., ), that one must add to of solution to raise by one unit:
The units of are . For additional details, refer to Roos & Boron (1981)Boron & Boulpaep (2016). Note that BDW defined as a negative number, as did Koppel and Spiro in their original definition of buffering Koppel, 1914Roos & Boron, 1980, rather than as a now-conventional positive quantity, as did Van Slyke in his later work Van Slyke, 1922. BDW’s definition, which they consistently applied, has no effect on the outcome of their simulations. In the present paper, we will follow the definition of Van Slyke — defining as a positive number — and make appropriate sign changes to the derived equations.
3The Boron & De Weer Model for the Permeation by an Uncharged Weak Acid and its Conjugate, Anionic Weak Base¶
The BDW mathematical model consists of two time-dependent ordinary differential equations (ODEs), one describing the time-course of the concentration of total intracellular buffer () and the other the time-course of the intracellular free concentration (which is related to ). BDW derived these two equations for the general cases in which any buffer pair , or any buffer pair , can move passively across the plasma membrane of a prototype cell. Then, they applied these two general equations to their specific experimental conditions, namely exposure of a cell (a squid giant axon) to equilibrated extracellular or to equilibrated extracellular .
Here, following BDW’s approach, we begin by deriving the equations for . In the next section, we apply the same general formalism to .
Derivation for weak acids.¶
Imagine that a cell is exposed to a solution containing equilibrated and that both and initially move into the cell — because of the chemical gradient in the case of , and because of the electrochemical gradient in the case of .
An integrated form of Fick’s first law of diffusion describes the net passive influx[7] of ()
where is flux () and () is the membrane permeability to the uncharged weak acid . Note that this is a passive diffusion equation because is uncharged.
The constant field equation — also known as the Goldman, Hodgkin, Katz (GHK) Goldman, 1943Hodgkin & Huxley, 1952 equation — describes the net passive influx of ():
where () is the membrane permeability to the charged conjugate base , is the membrane potential (intracellular relative to extracellular potential), and is a shorthand for . Note that and have units of .
Although and can interconvert in the cytosol, BDW assumed that the intracellular concentration of total weak acid only can change due to the transmembrane fluxes of and (see Figure 3). Thus, the time rate of change of is
where () is the area-to-volume ratio for the cell, and converts the transmembrane flux per unit area (in units of ) to a time rate of change per unit cell volume ( or ). (11) is the first of two ODEs of the BDW model for the buffer pair .
![Cartoon illustrating the main assumptions in the BDW model of permeating uncharged weak acid \mathrm{HA} and its conjugate anionic weak base \mathrm{A^-}. The BDW model consists of two time-dependent ODEs. The first one describes the time-course of the intracellular concentration of total weak acid \mathrm{[TA]_i}, and the second one describes the time-course of \mathrm{[H^+]_i}. BDW assumed that \mathrm{[TA]_i} changes in time because of the transmembrane fluxes of \mathrm{HA} (J_\mathrm{HA}) — modelled according to Fick’s first law of diffusion — and \mathrm{A^-} (J_\mathrm{A^-}) — modelled according to the Goldman, Hodgkin, Katz (GHK) equation. According to BDW, the time rate of change of \mathrm{[H^+]_i} depends on the net rate \mathrm{d}Q/\mathrm{d}t at which acids are added into the cytosol. BDW assumed that \mathrm{d}Q/\mathrm{d}t depends on (i) the release of \mathrm{H^+} by some fraction x of the entering \mathrm{HA} (i.e., xJ_\mathrm{HA}), (ii) the consumption of \mathrm{H^+} by some fraction y of the entering \mathrm{A^-} (i.e., yJ_\mathrm{A^-}), and (iii) the additional rate of intracellular \mathrm{H^+} consumption via metabolism or active acid extrusion (J_\mathrm{H^+}).](https://pub.curvenote.com/9aa0a502-2696-45fc-93bc-78a2d8af4b7d/public/Figure 3-73c27f0c2fe4c79d60d871eacb50375d.png)
Figure 3:Cartoon illustrating the main assumptions in the BDW model of permeating uncharged weak acid and its conjugate anionic weak base . The BDW model consists of two time-dependent ODEs. The first one describes the time-course of the intracellular concentration of total weak acid , and the second one describes the time-course of . BDW assumed that changes in time because of the transmembrane fluxes of () — modelled according to Fick’s first law of diffusion — and () — modelled according to the Goldman, Hodgkin, Katz (GHK) equation. According to BDW, the time rate of change of depends on the net rate at which acids are added into the cytosol. BDW assumed that depends on (i) the release of by some fraction of the entering (i.e., ), (ii) the consumption of by some fraction of the entering (i.e., ), and (iii) the additional rate of intracellular consumption via metabolism or active acid extrusion ().
Later, Bevensee and Boron defined the time rate of change per unit volume (e.g., ) as a ‘pseudoflux’ , with the area-to-volume ratio folded into the value of Bevensee & Boron, 2013. Physiologists sometimes prefer to present experimental data in terms of pseudoflux because most mammalian cells often have complex geometries that make it difficult to estimate surface area.
In deriving the second ODE of their model, BDW started by noting that the time rate of change of free protons, , depends on the rate at which acids are added into the cytosol per unit volume and per unit time — denoted () where is the total intracellular acid content. Like , both and are pseudofluxes.
In their simple system of a squid giant axon exposed to (i.e., ) and (i.e., ), BDW assumed that only three general processes affect : (i) the release of by some fraction () of the entering (i.e., ), (ii) the consumption of by some fraction of the entering (i.e., ), and (iii) the “additional” rate of intracellular consumption or active extrusion of (; see Figure 3) above the fixed background rate of extrusion necessary to balance the fixed background rate of acid loading (i.e., addition of or equivalent acid, or removal of or equivalent base) in the absence of . Thus,
A critical insight by BDW is that during each infinitesimal increment in time during which a bolus of enters the cell, the entering redistributes itself between (and ) vs , according to the pre-existing ratios and . Thus, the fraction of entering that remains is
This is also the fraction of entering that combines with and becomes . Combining the above expression with (2),
which BDW defined as . Conversely, the fraction of that remains is
This is also the fraction of entering that dissociates to form and . Combining the above expression with (2),
In summary, (12) becomes:
BDW modelled () as the additional proton-extrusion rate above the fixed background rate
where () is the proton-pumping rate constant, is the additional flux of above the background flux of , which occurs at the resting of (i.e., resting of ). Note that has units of (), consistent with the membrane permeability terms and in (9) and (10).[8]
The BDW authors used the definition of buffering power, in its infinitesimal form, to derive the relation between and , as shown in the following steps.
Our first goal is to obtain an expression for in terms of . According to the chain rule:
By definition (see (8)), (), or equivalently
Our next goal is to obtain an expression for in terms of . According to the chain rule:
By definition, , so that:
Combining (22) and (23), we have
or equivalently,
Substituting (21) into (25), we obtain
Finally, substituting (17) into (26),
which is the second equation of the BDW model.
Substituting for (from (14)), , in (11) and (27), we obtain the two ODEs of the BDW model in terms of and :
where (from (9)), and (from (10)) are given by:
and is given by (18).
The numerical solution of the above two equations yields the time courses of and , which in turn yield the time-courses of and via:
where
Simulation for experiments.¶
BDW employed (28) and (29) to simulate the experiments in which they exposed a squid giant axon to a solution containing equilibrated . Their simulation protocol was a step change in (a) extracellular from to ( or, with , ) and (b) extracellular from to (the value that has in a solution containing at of )[9]. The step change is applied for () at constant .
Table 1 and Table 2 report the parameter values used by BDW. Table 1 provides parameter values that are common to both the and the experiments. Table 2 provides parameter values exclusive to the experiments only.
Table 1:Parameter values used in both simulations of squid-axon experiments and experiments.
| Symbol | Name | BDW Value | Unit | New Value | Unit |
|---|---|---|---|---|---|
| temperature | () | º (º ) | |||
| gas constant | |||||
| Faraday constant | |||||
| area/volume ratio | [10] | ||||
| extracellular |
Table 2:Parameter values for simulations of squid-axon experiments.
| Symbol | Name | BDW Value | Unit | New Value | Unit |
|---|---|---|---|---|---|
| buffering power | |||||
| solubility constant for | [11] | ||||
| partial pressure of | |||||
| partial pressure of | |||||
| extracellular | |||||
| extracellular | |||||
| membrane permeability | |||||
| membrane permeability | |||||
| acid dissociation constant | |||||
| acid dissociation constant | |||||
| membrane potential | [12] | ||||
| pump rate constant | [13] | ||||
| intracellular | |||||
| basal | [14] |
In the present work, the differential (28) and (29) — when coded in CellML and solved with OpenCOR — produce the plots in Figure 4. The simulation file Boron-CO2.sedml contains the computational setting for running the model. Open the .sedml file in OpenCOR and click Run Simulation. The initial conditions are and . Note that Figure 4 illustrates the time courses not only of — as presented by BDW — but also of quantities (e.g., various solute concentrations and fluxes) not displayed in the original paper; these values are useful for understanding the processes that contribute to the transient. Moreover, our curated and annotated version of the BDW model also allows one to alter the parameter values from those originally chosen by BDW, thereby extending the ability of the user to investigate the predictive power of the computational model.
![Solution of the BDW model during and following a 2700~\mathrm{s} period of externally applied \mathrm{CO_2}. In these simulations \mathrm{pH_o}=7.70 and \mathrm{[HCO_3^-]_o} is determined from the equilibrium with \mathrm{[H^+]_o} and \mathrm{CO_2} (footnote 9). Note that, during the plateau phase, \mathrm{[HCO_3^-]_i} continues to rise as \mathrm{pH_i} rises at a constant \mathrm{[CO_2]_i} (the proton pumping rate k is set to 300 \mathrm{s^{-1}}, thus k/\rho= 0.0375 \mathrm{m\cdot s^{-1}}). Note also that, after the removal of \mathrm{CO_2}/\mathrm{HCO_3^-}, \mathrm{pH_i} rises to a higher value (\sim 8.15) than its starting value (\sim 7.4), indicating the net extrusion of acid from the cell during the \mathrm{CO_2}/\mathrm{HCO_3^-} exposure.](https://pub.curvenote.com/9aa0a502-2696-45fc-93bc-78a2d8af4b7d/public/Figure 4-b527c01fcde0bf138feda7fc2762a7fc.png)
Figure 4:Solution of the BDW model during and following a period of externally applied . In these simulations and is determined from the equilibrium with and (footnote 9). Note that, during the plateau phase, continues to rise as rises at a constant (the proton pumping rate is set to , thus ). Note also that, after the removal of , rises to a higher value () than its starting value (), indicating the net extrusion of acid from the cell during the exposure.
4The Boron & De Weer Model for the Permeation by an Uncharged Weak Base and its Conjugate, Cationic Weak Acid¶
Following an approach analogous to the one outlined above for weak acids, BDW derived two time-dependent ODEs. The first describes the time-course of the concentration of total intracellular buffer (), and the other the time-course of the intracellular free , for any buffer pair .
Derivation for weak bases.¶
Imagine that a cell is exposed to a solution containing equilibrated , and that both and initially move into the cell — because of the chemical gradient in the case of , and because of the electrochemical gradient in the case of .
![Cartoon illustrating the main assumptions in the BDW model of permeating uncharged weak base \mathrm{B} and its conjugate anionic weak acid \mathrm{BH^+}. The BDW model consists of two time-dependent ODEs. The first one describes the time-course of the intracellular concentration of total weak base \mathrm{[TB]_i}, and the second one describes the time-course of \mathrm{[H^+]_i}. BDW assumed that \mathrm{[TB]_i} changes in time because of the transmembrane fluxes of \mathrm{HA} (J_\mathrm{B}) — modelled according to Fick’s first law of diffusion — and \mathrm{BH^+} (J_\mathrm{BH^+}) — modelled according to the Goldman, Hodgkin, Katz (GHK) equation. According to BDW, the time rate of change of \mathrm{[H^+]_i} depends on the net rate \mathrm{d}Q/\mathrm{d}t at which acids are added into the cytosol. BDW assumed that \mathrm{d}Q/\mathrm{d}t depends on (i) the release of \mathrm{H^+} by some fraction x of the entering \mathrm{BH^+} (i.e., xJ_\mathrm{BH^+}), (ii) the consumption of \mathrm{H^+} by some fraction y of the entering \mathrm{B} (i.e., yJ_\mathrm{B}), and (iii) the additional rate of intracellular \mathrm{H^+} consumption via metabolism or active acid extrusion (J_\mathrm{H^+}).](https://pub.curvenote.com/9aa0a502-2696-45fc-93bc-78a2d8af4b7d/public/Figure 5-60f01f6992b5873c6bc8e57ceafb9194.png)
Figure 5:Cartoon illustrating the main assumptions in the BDW model of permeating uncharged weak base and its conjugate anionic weak acid . The BDW model consists of two time-dependent ODEs. The first one describes the time-course of the intracellular concentration of total weak base , and the second one describes the time-course of . BDW assumed that changes in time because of the transmembrane fluxes of () — modelled according to Fick’s first law of diffusion — and () — modelled according to the Goldman, Hodgkin, Katz (GHK) equation. According to BDW, the time rate of change of depends on the net rate at which acids are added into the cytosol. BDW assumed that depends on (i) the release of by some fraction of the entering (i.e., ), (ii) the consumption of by some fraction of the entering (i.e., ), and (iii) the additional rate of intracellular consumption via metabolism or active acid extrusion ().
Assuming, as in Figure 5, that only can change due to the transmembrane fluxes of () and (), the time rate of change of — analogous to (11) above — is
where () is again the area-to-volume ratio for the cell. The equation
is an integrated form of Fick’s first law of diffusion that describes the net passive flux of B, and
describes the net passive influx of according to the GHK equation. In the previous two equations, () is the membrane permeability to the uncharged weak base , () is the membrane permeability to the charged conjugate weak acid , and is a shorthand for . (33) is the first of two ODEs of the BDW model for the buffer pair .
The second equation of the BDW model for a weak base — analogous to (27) above — is
where is the same as in (18) and
and
Substituting for , , , in (33) and (36), we obtain the two ODEs of the BDW model in terms of and
where
and is given by (18).
Numerically integrating the above two equations yields the time courses of and , from which we can compute the time-courses of and from
where
Simulation for experiments.¶
BDW employed (39) and (40) to simulate the experiments in which they exposed a squid giant axon to equilibrated . Their simulation protocol was a step change in extracellular from to (that is, a step change in from to , and in from to ) applied for () at constant .[15]
Table 1 and Table 3 report the parameter values used by BDW. Note that in the simulations, is always zero, that is, does not affect these processes.
Table 3:Parameter values for simulations of squid-axon experiments.
| Symbol | Name | BDW Value | Unit | New Value | Unit |
|---|---|---|---|---|---|
| buffering power | |||||
| extracellular total ammonia | [16] | ||||
| extracellular | |||||
| extracellular | |||||
| membrane permeability | |||||
| membrane permeability | [17] | ||||
| acid dissociation constant | |||||
| acid dissociation constant | |||||
| membrane potential | [18] | ||||
| pump rate constant | |||||
| intracellular | [19] |
The differential (39) and (40) — when coded in CellML and solved with OpenCOR — produce the plots in Figure 6. The simulation file Boron-NH3.sedml contains the computational setting for running the model. Open the .sedml file in OpenCOR and click Run Simulation. The initial conditions are and .
![Solution of the BDW model during and following a 1500~\mathrm{s} period of externally applied \mathrm{NH_4Cl}. In these simulations \mathrm{[NH_4^+]_o}=9 \mathrm{mM}. The intracellular fluid becomes alkaline as \mathrm{NH_3} enters (note the J_\mathrm{NH_3} time course) and hydrates to form \mathrm{NH_4^+} and \mathrm{OH^-}. Additional passive \mathrm{NH_4^+} entry (note J_\mathrm{NH_4^+} time course) down its electrochemical gradient opposes the effect of the \mathrm{NH_3} entry and slightly reduces the \mathrm{pH_i} increase. Upon removal of \mathrm{NH_4Cl}, \mathrm{[NH_3]_i} and \mathrm{[NH_4^+]_i} decay towards their original values, but \mathrm{pH_i} drops well below its original value of 7.40.](https://pub.curvenote.com/9aa0a502-2696-45fc-93bc-78a2d8af4b7d/public/Figure 6-1047425d7036d5b8aa4aa54608386890.png)
Figure 6:Solution of the BDW model during and following a period of externally applied . In these simulations . The intracellular fluid becomes alkaline as enters (note the time course) and hydrates to form and . Additional passive entry (note time course) down its electrochemical gradient opposes the effect of the entry and slightly reduces the increase. Upon removal of , and decay towards their original values, but drops well below its original value of .
5Discussion¶
The publication of this retrospective paper provides an opportunity to clarify some concepts in the original paper that have benefitted from subsequent experimental and theoretical advances. We also provide some additional parameters missing from BDW and correct some minor errors. Most importantly, the curated model is now freely available on the Physiome website[20] in standardised form (CellML) that can be run in the open source software OpenCOR. A follow-up paper will be written that recasts the equations in bond graph form to facilitate their incorporation into more complex models where regulation is coupled with other cellular processes.
5.1Historical Context¶
The 1976 Boron & De Weer paper introduced the first models to simulate the time course of . By developing a predictive mathematical model based on first principles, BDW provided a quantitative basis for interpreting their new data on the time-dependent response of to step changes in extracellular (and pairs in general) and (and pairs in general). The models also provided a clear, quantitative basis for interpreting BDW’s new data on how cells regulate their , which BDW modelled as a -dependent -extrusion mechanism. Below, we will introduce a broader concept termed “acid extrusion” Boron, 1977. The first of the two BDW models elucidates how the transmembrane fluxes of a neutral weak acid and its anionic conjugate weak base affects , with the acid-extrusion becoming increasingly important as falls. The second model simulates how the transmembrane fluxes of a neutral base and its cationic conjugate weak acid affects .
Of course, BDW were not the first to undertake a quantitative assessment of how acids or bases affect, or are affected by, the of a solution. Below, we divide the earlier work into two major categories, (a) analyses of how neutral weak acids (and their anion conjugate weak bases) or neutral weak bases (and their cationic conjugate weak acids) affect in simple systems, and (b) analyses of how the distribution of (or ) across a barrier, such as a cell membrane, are affected by and .
5.2Development of the Concept of Buffering in Simple and Complex Systems¶
Buffering power.¶
In 1914, Koppel (1914) introduced the first modern definition of the chemical buffering (i.e., “magnitude of moderation”, or ) of by a weak-acid/weak-base conjugate pair, and — based on first principles — derived an expression for buffering power. Because these authors defined in terms of the amount of strong acid that one must add to a solution to produce a change, is a negative number. Roos & Boron (1980) translated the Koppel-Spiro paper from its original German, and provided a historical context.
Initially unaware of the work of Koppel and Spiro, Van Slyke (1922) independently defined buffering power — to which he assigned the Greek letter . Because he defined in terms of the amount of strong base that one must add to a solution to produce a change (see (8)), is a positive number. Although modest differences exist between the efforts of Koppel and Spiro on the one hand and Van Slyke on the other, they are quite similar. Nevertheless, it is Van Slyke’s definition of that has become the modern convention throughout chemistry and physiological chemistry.
Koppel & Spiro and Van Slyke quantitatively described how — in a one-compartment solution — weak acids, weak bases, ampholytes, and weak-acid/base mixtures can buffer added strong acid or strong base. In their analysis, the system both begins and ends in an equilibrium state. Of course, in their pioneering work, these authors had no reason to contemplate time courses or barriers separating more than one compartment.
In their work, BDW defined as a negative number, as Koppel and Spiro defined their .
The “Davenport” diagram.¶
This nomogram Boron & Boulpaep, 2016 is a powerful tool for graphically computing the effects of respiratory acid-base disorders (caused by changes in in a system open to ) and metabolic acid-base disorders (caused by the addition or removal of or a strong acid or base). The underlying assumption for the Davenport diagram is that the system is in equilibrium. The first component of a Davenport diagram (see Figure 7) is a plot of vs for one or more values of — these are the isopleths that describe the equilibrium among , , and . At any on any isopleth, the slope is the open-system buffering power (). The second component is a linear plot, on the same axes, of the concentration of all protonated forms of all non- buffers vs . At any , the slope is the buffering power of all non- buffers (), and the line is termed the non- buffer line. Its linearity implies that is insensitive to changes in . The intersection of an isopleth with the non- buffer line describes the current state of the system, when both buffer and non- buffers are simultaneously in equilibrium. Davenport developed a series of rules for using this paradigm to interpret acid-base disorders, and these rules are well founded in physical chemistry.
![A Davenport diagram (from figure 5A of BDW). This nomogram consists of two kinds of plots. The first kind of plot is represented by the four \mathrm{CO_2} isopleths that slope upwards from left to right. Each isolpleth represents all possible combinations of \mathrm{[HCO_3^-]_i} and \mathrm{pH_i} for a given \mathrm{[CO_2]_i} (described here as the \% of the air phase that is \mathrm{CO_2}). The second kind of plot is represented by the two lines that slope downwards from left to right. The slope of these parallel lines describes the buffering power of non-\mathrm{CO_2}/\mathrm{HCO_3^-} buffers. BDW assumed that the experiment starts at point A, at \mathrm{pH_i} = 7.32 and 0.1\% \mathrm{CO_2}. The addition of 5\% \mathrm{CO_2} causes the \mathrm{pH_i} at equilibrium to fall to the point represented by point B. The extrusion of acid during the \mathrm{CO_2}/\mathrm{HCO_3^-} exposure causes the system to move along the 5\% \mathrm{CO_2} isopleth from point B to point C. Finally, upon removal of \mathrm{CO_2}/\mathrm{HCO_3^-}, the system returns to the original \mathrm{CO_2}/\mathrm{HCO_3^-} isopleth, but now at point D. The difference between points D and A represents the \mathrm{pH_i} overshoot. In the BDW paper, \beta — the slope of the lines in this figure — appeared to be 25 \mathrm{mM}/\mathrm{pH} unit. In fact, the value of \beta determined in the \mathrm{NH_3}/\mathrm{NH_4^+} experiments (also shown as a Davenport-like diagram in figure 5B of BDW) was 9 \mathrm{mM}/\mathrm{pH} unit. The reason for this discrepancy was probably that BDW delivered the \mathrm{CO_2}/\mathrm{HCO_3^-} solution to the axon using a peristaltic pump and Silastic tubing, which they later realised has a high \mathrm{CO_2} permeability. Thus, the \mathrm{[CO_2]} reaching the axon was <5\%, accounting for the artificially inflated value for \beta. In their follow-up paper , BDW delivered the \mathrm{CO_2}/\mathrm{HCO_3^-} solutions from glass syringes and through Saran tubing, which has an extremely low \mathrm{CO_2} permeability.](https://pub.curvenote.com/9aa0a502-2696-45fc-93bc-78a2d8af4b7d/public/Figure 7-959f86e0fa6716eed8e8f81852ae85ee.png)
Figure 7:A Davenport diagram (from figure 5A of BDW). This nomogram consists of two kinds of plots. The first kind of plot is represented by the four isopleths that slope upwards from left to right. Each isolpleth represents all possible combinations of and for a given (described here as the of the air phase that is ). The second kind of plot is represented by the two lines that slope downwards from left to right. The slope of these parallel lines describes the buffering power of non- buffers. BDW assumed that the experiment starts at point A, at and . The addition of causes the at equilibrium to fall to the point represented by point B. The extrusion of acid during the exposure causes the system to move along the isopleth from point B to point C. Finally, upon removal of , the system returns to the original isopleth, but now at point D. The difference between points D and A represents the overshoot. In the BDW paper, — the slope of the lines in this figure — appeared to be unit. In fact, the value of determined in the experiments (also shown as a Davenport-like diagram in figure 5B of BDW) was unit. The reason for this discrepancy was probably that BDW delivered the solution to the axon using a peristaltic pump and Silastic tubing, which they later realised has a high permeability. Thus, the reaching the axon was , accounting for the artificially inflated value for . In their follow-up paper Boron & De Weer, 1976, BDW delivered the solutions from glass syringes and through Saran tubing, which has an extremely low permeability.
The Davenport diagram traces its origins to the analyses of blood by several eminent investigators about a century ago. It was Henderson (1921) — as far as we are able to ascertain — who in figure 4 of his paper was the first to plot vs for two (rather than ) isopleths, and for two different values of (i.e., those produced by and ).
The power of the Davenport approach is that, knowing the initial conditions and the dependence of , one can compute with fair accuracy (using the nomogram) or great accuracy (using a computer to solve the equations numerically) the result of virtually any acid-base disorder in the pathophysiological range for a system containing and a mixture of non- buffers. In their figure 5A (reproduced here as Figure 7), BDW used a Davenport approach to describe the initial steady state of a squid giant axon (point A at , ), the initial effect of an exposure to increased (point B, an example of intracellular respiratory acidosis), the effect of the plateau-phase recovery (point C, an example of intracellular compensatory metabolic alkalosis), and finally the effect of removing extracellular (point D, an example of metabolic alkalosis) to account for the overshoot. If one does not know , the Davenport diagram allows one to compute it from the initial and final . The numerical integration of the BDW equations — when and the acid extrusion rate are both zero — should in principle yield, at infinite time, results consistent with the Davenport diagram.
In their paper (their figure 5B), BDW introduced a novel Davenport-like diagram for the buffer system, with on the ordinate (replacing ), isopleths (replacing isopleths), and a line describing non- buffering power (replacing the line describing non- buffering power). Like the classical Davenport diagram, this one (or others like it, constructed for other buffer pairs) can be a useful tool for interpreting — in the steady state — problems in acid-base chemistry.
In the Davenport analysis, the initial and final conditions represent equilibria. Davenport-like diagrams provides no information about the time course of between the initial and final states. Nor can Davenport-like diagrams deal with time course, barriers (e.g., cell membranes), permeabilities to substances other than the neutral molecule (e.g., , ), or active transport. Of course, in their pioneering work, Henderson, Davenport, and other authors contributing to this nomogram had no reason to contemplate these future complexities.
5.3Pre-BDW Analyses of Transmembrane Distributions of Weak Acids and Bases¶
As summarised by Roos & Boron (1981), about a century ago, several authors — who assumed that equilibrates across the plasma membrane but that is impermeant — used the sum to compute the steady-state of several cell types. Then, beginning in 1940, a series of authors introduced three successively more sophisticated mathematical analyses for the steady-state transmembrane distribution of a neutral weak acid and its anionic conjugate weak base (i.e., ), and three more for the distribution of a neutral weak base and its cationic, conjugate weak acid (i.e., ). We will now present these analyses in order of increasing complexity, and according to their sequence in time (see Figure 8). They all have in common the assumption that the system is either in equilibrium or at least in a steady state supported by the input of energy.
Figure 8:Timeline of acid-base/ models prior to BDW. Note that the time-dependent BDW model of weak acid/conjugate weak base () collapses to Roos (1965)Roos (1975) steady-state model for . The Roos (1965)Roos (1975) model, in turn, reduces to the Milne et al. (1958) model for when the membrane potential () approaches zero. Finally, the Milne et al. (1958) model reduces to the Jacobs (1940) model — where only one uncharged species is permeant — as the permeability of () approaches zero. Similarly, the time-dependent BDW model of weak base/conjugate weak acid () collapses to Boron & Roos (1976) steady-state model for . The Boron & Roos (1976) model in turn reduces to the Orloff & Berliner (1956) model for when approaches zero. Finally, the Orloff & Berliner (1956) model reduces to the Jacobs-like model — where only one uncharged species is permeant — as () approaches zero.
Jacobs (, neutral weak acid).¶
Jacobs (1940) presented a general mathematical model that describes the equilibrium transmembrane distribution of total weak acid (), assuming that only the neutral species (), but not the conjugate weak base () can cross the membrane:
In an accompanying document, we show the derivation of the above equation[21], the linear form of which is:
We also present the derivation of this linear version in an accompanying document[22]. A major assumption in the derivations of both (44) and (45) is that does not merely move but fully equilibrates across the cell membrane. That is, the system is in equilibrium.
The notion that only the neutral member of the buffer species (i.e., ) can traverse the membrane is known as nonionic diffusion. (44) and (45) tell us that, as rises (i.e., falls), rises steeply because rises exponentially with — a concept known as “trapping” (of the ).
Jacobs-like equation (, neutral weak base).¶
Although Jacobs did not present the analogous equation for total weak base (), others have derived it, including Roos & Boron (1981):
In an accompanying document, we show the derivation of (46),[23] the linear form of which is:
We also present a derivation of (47).[24] A major assumption in the derivations of (46) and (47) — analogous to the situation for (44) and (45) — is that fully equilibrates across the cell membrane. That is, the system is in equilibrium.
The notion that only the neutral member of the buffer species (i.e., ) can traverse the membrane is another example of nonionic diffusion. (46) and (47) tell us that, as falls (i.e., rises), rises steeply because rises exponentially with the decrease in — another example of “trapping” (of the ). Such trapping is especially important in renal physiology, where acidic fluid in renal tubules can trap the cationic form of a buffer pair (e.g., ).
(44) and (46) provide the theoretical foundation for using ratios for permeant weak acids (e.g., , above, and the later 5,5’-dimethyl-2,4-oxazolidinedione [DMO] technique) or ratios for permeant weak bases (e.g., for methylamine; see Boron & Roos (1976)) for computing steady-state . For example, equation tells us that, as rises, will rise nearly exponentially; this occurs because rises in a precisely exponential fashion.
Orloff & Berliner (, ).¶
Orloff & Berliner (1956) extended the Jacobs-like model to a neutral weak base and its cationic conjugate weak acid by permitting not just but also to permeate a barrier separating compartments 1 and 2. They avoided the complication of electrodiffusion by assuming a transmembrane voltage of zero (i.e., ). Because they recognised that the flux of across the barrier would cause to drift in opposite directions in the two compartments, they assumed that independent, energy-requiring processes would stabilise in the two compartments and establish a steady state described by
In an accompanying document[25], we show the derivation of (48). If we put (48) into the same form as (47), which expresses the buffer concentrations in terms of and then — for cells — we have
As approaches zero, (49) reduces to (47) — the Jacobs-like equation for a weak base. The accompanying document[26] also shows the derivation of (49), as well as the mathematics that shows the limit of this expression as approaches zero.
Note that the flux of — which leads to a flux of in the opposite direction — tends to push the system off equilibrium. As noted above, the Orloff-Berliner equation requires that the system be in a steady-state, which can be achieved, as they recognised, only by an input of energy to maintain all relevant concentrations constant over time.
Milne et al (, ).¶
Milne and colleagues (1958) developed a steady-state expression similar to that of Orloff & Berliner, but for the transmembrane distribution of a neutral weak acid and its anionic conjugate weak:
An accompanying document[27] shows the derivation of (50), which we can also put in the form of (49), which expresses buffer concentrations in terms of and :
As approaches zero, (51) reduces to (45) — the Jacobs equation for a weak acid. The accompanying document[28] also shows the derivation of (51), as well as the mathematics that shows the limit of this expression as approaches zero.
Roos (, non-zero ).¶
In 1965 and 1975, Roos extended the Irvine model by allowing to assume non-zero values Roos, 1965Roos, 1975, and derived the following equation,
where has the same meaning as in the derivation of the BDW equations: . An accompanying document[29] shows the derivation of (52). This document also shows that, at the limits of certain parameters, (52) simplifies to the expected equation:
- As , the Roos equation simplifies to the equation of Milne et al (which assumes ).
- As , the Roos equation simplifies to the Jacobs equation (which assumes ).
- As , the Roos equation simplifies to the Nernst equation (which assumes permeability to only ).
The Roos equation was important historically because it allowed one to assess possible errors in values computed from the transmembrane distribution of a neutral weak acid (e.g., , DMO) and its monovalent anion conjugate weak base. These errors could in principle arise from membrane permeability to (as already considered by Milne et al), as influenced by .
Boron & Roos (, non-zero ).¶
Finally, Boron & Roos (1976) derived an equation similar to the Roos equation, but for the distribution of a neutral weak base and its monovalent cationic conjugate weak acid:
where has the same meaning as in the derivation of the BDW equations: . An accompanying document[30] shows the derivation of (53). This document also shows that, at the limits of certain parameters, (53) simplifies to the expected equation:
- As , the Boron-Roos equation simplifies to the equation of Orloff & Berliner (which assumes ).
- As , the Boron-Roos equation simplifies to the Jacobs-like equation for a neutral weak base (which assumes ).
- As , the Boron-Roos equation simplifies to the Nernst equation (which assumes permeability to only ).
The Boron-Roos equation was important historically because it allowed one to assess possible errors in values computed from the transmembrane distribution of a neutral weak base (e.g., methylamine) and its monovalent anion conjugate weak base. These errors could in principle arise from membrane permeability to , as influenced by . In their paper, Boron and Roos used the transmembrane distribution of methylamine/methylammonium to monitor a downward drift in caused by the passive influx of methylammonium. This was the first use of a chemical-distribution technique to follow changes over time.
We have already noted for the Orloff-Berliner equation that the derivation requires that the system be in a steady-state, which can be achieved only by an input of energy to maintain all relevant concentrations constant over time. The same is true for the equations of Milne et al, Roos, and Boron & Roos.
5.4Comparison of the Pre-BDW Analyses with the BDW Equations¶
The BDW equations build on the earlier work on buffering and transmembrane distributions of weak acids and bases, presented in the previous two sections. Of course, the work of Koppel and Spiro, and Van Slyke, as well as their predecessors who developed the physico-chemical principles of acid-base chemistry, is at the heart of the BDW model.
An important aspect of the Davenport diagram is its predictive power. For example, given as well as the initial and , the Davenport approach can predict the effect of an increase in on the final equilibrium conditions. However, the Davenport approach makes no statement about mechanism or time course. The BDW approach has all the predictive power of Davenport, but also addresses mechanism and time course.
The six pre-BDW approaches for assessing transmembrane distributions of and all start with the weak acid/base present and the system in an equilibrium or at least in a steady state. Somehow, the system — the cell and its surrounding fluid — went from a condition with no weak acid/base present to a condition with the weak acid/base present at equilibrium/steady state. The older models make no attempt to describe how and how fast the system achieved the new state, and — unlike the Davenport approach — have no predictive value for relating initial and final conditions. It is worth noting that the investigators who developed these six approaches were interested mainly in using tracer levels of weak acids/bases to compute . The intention was that tracer levels would have minimal effects on the state of the system — hence, the minimal interest in the prediction. Note that, at infinite time (and with no acid extrusion), the BDW equations reduce to those six transmembrane-distribution analyses presented in the previous section.
To some extent, the BDW models represent a merger of the Davenport and the six pre-BDW approaches for assessing transmembrane distributions of and . Like the Davenport approach, the BDW approach is predictive. However, unlike Davenport’s approach, the BDW models provide insight into mechanism and time course, and are applicable even when the system is far from equilibrium/steady state. Like the six transmembrane-distribution models, the BDW models provide insight into how and affect , , , and . Unlike the six transmembrane-distribution models, the BDW models are predictive and provide insight into mechanism and time course.
It is worth noting that BDW derived their equations under the influence of Albert Roos, who had derived the transmembrane-distribution model for Roos, 1975 and who inspired Boron’s derivation of the model Boron & Roos, 1976.
5.5Post-BDW Development of Concepts¶
Fundamental law of regulation.¶
Recall that one of the intermediate steps of the derivation of the BDW equations for the system was (21), which we reproduce here:
Another intermediate step was the description of in (17), which we reproduce here:
Combining the above two equations yields a primitive form of the fundamental law of regulation:
In the later literature, Boron and colleagues in effect dissected the term — the net rate at which appear in the cytosol () — into two concepts that are more general than those considered by BDW in (56): the intracellular acid-loading rate () and the intracellular acid-extrusion rate ().
In the narrow definition of (56) the only term is . In the post-BDW literature, Boron and colleagues defined to comprise every process that adds the equivalent of to or removes the equivalent of from the cytosol, including channels (mediating passive influx), a variety of transporters (e.g., - exchangers mediating efflux), and the metabolic production of acid. Boron et al. (1979) measured and introduced the term “rate of acid introduction”. Roos & Boron (1981) later replaced this term when they coined “acid-loading rate”.
In the narrow definition of (56), the terms are and . Note that BDW defined as the -extrusion rate above baseline. In the post-BDW literature, Boron and colleagues defined to comprise every process that removes the equivalent of from or adds the equivalent of to the cytosol, including a variety of transporters (e.g., pumps, - exchangers) that mediate efflux and others (e.g., -driven transporters, cotransporters) that mediate base efflux. It appears that Boron in 1977 provided the first definition of acid extrusion.
Recasting (56) in terms of and ,
The modern version of the fundamental law of regulation is:
Here, with the explicit inclusion of , and have the units ().
(58) tells us that is stable when , rises when , and falls with . BDW did not derive the above equation in their 1976 paper. The first statement of the concept of (58) was in sentence form in the review by Roos & Boron (1981), who also provided a graphical example in their figure 12. By 1989, Boron presented a version of (58) that lacks the term . Thus, he implicitly defined and in terms of (). This 1989 paper also included examples of how (58) could help one interpret the time course of in many experimental settings. For example, a bolus introduction of acid into a cell — an “acute” intracellular acid load, coined by Boron by the time of this 1989 review — rapidly lowers , but also increases and decreases . Because now exceeds , recovers to its initial value. By 1992, Boron coined the term “fundamental law of regulation” to describe (58).
Note that an acute or bolus intracellular acid load is to be distinguished from a chronic intracellular load (the rate of which is ). The amount of an acute acid load can instantly be quantitated in acid equivalents. The amount of a chronic acid load can only be quantitated in acid equivalents by integrating over time.
Later, Boron (2004) provided a more detailed description of regulation, as well as a tongue-in-cheek analogy — which he had been using for years in lectures — between regulation and the temperature regulation of a house, complete with multiple furnaces (acid-extruders), multiple air conditioners (acid loaders), heat capacity (buffering power), a thermostat ( sensitivity of the transporters), and weather radar (extracellular sensors for , , and ).
Flux vs pseudoflux.¶
One can define the acid-loading and acid-extrusion rates as strict fluxes, with units of . Such of use of or — for example (the acid-extruding flux mediated by a cotransporter, ) or (the acid-loading flux mediated by a - exchanger) — is most appropriate in cases in which the cell has a simple geometry (e.g., a squid axon or Xenopus oocyte). However, for cells with complex geometries — where surface-to-volume ratios are difficult to define — physiologists often present experimental data in the units “”. To avoid confusion between the two systems of measurement, Bevensee & Boron (1995) introduce the term “pseudoflux” and the symbol (rather than ). Thus, in their paper, the authors referred to and (rather than and ), so that the fundamental law of regulation becomes
Effect of isodirectional fluxes of or fluxes on near a membrane.¶
The term in (14) and (37), and the term in (16) and (38) inspired an insight that yields equations analogous to the familiar Henderson-Hasselbalch equation. The starting point was the following question: if both and cross the membrane in the same direction, how will their isodirectional fluxes affect on the cis side (the side from which they depart) and the trans side (the side to which they go)?
In the discussion of their paper, Musa-Aziz et al. (2009) showed that if the flux ratio is the same as the concentration ratio , the flux will have no effect on because and are appearing (or disappearing) at the inner surface of the membrane in a proportion equal to their respective, pre-existing concentrations. Similarly, if the flux ratio is the same as the concentration ratio , the flux will have no effect on . Stated somewhat differently, as in equation 4 of Musa-Aziz et al. (2009),
and
Here, is the at which the isodirectional fluxes and will have no effect on . Similarly, is the at which the isodirectional fluxes and will have no effect on . If = , then their null values are the same (i.e., at this , the isodirectional fluxes of and will have no effect on either ). One can write equations similar to (60) and (61), but for and (or and ), which we do here for the first time:
and
The previous four equations can be valuable for interpreting, for example, the effects of isodirectional fluxes of and . The equations cannot predict the speed or extent of the change, only the direction. For instance, if we expose a cell to a solution and falls, does mean that ? The intuitive answer would be, “yes”. The actual answer is, “not necessarily”. Imagine that — in a -free environment — we have a -free cell at a of . Although no is present, we assume that the of the equilibrium would be if were present. We now suddenly add to the extracellular fluid — the precise ratio is of no consequence here because we will focus on the intracellular fluid. and now begin to enter the cell by any route. If the ratio , then will not change from its original value of because and are entering the cell precisely in the correct ratio for a of . If (in this case), will fall because the alkalinising effects of the influx are less than the acidifying effects of the influx. This approach cannot tell us how fast or how far will fall, only that, initially at least, it must fall. The same analysis can be applied — simultaneously — to the extracellular fluid.
5.6Post-BDW Models of Acid-Base Fluxes/Chemistry¶
Following the development of the BDW model, the field of regulation has seen several modelling efforts. Here, we summarise several, with emphasis on those models that, in our opinion and to the best of our knowledge, have contributed to advance the field, either by extending the BDW model or by introducing new modelling paradigms.
Keifer and Roos.¶
In 1981, Keifer and Roos refined the BDW model by modifying the BDW assumption that addition of an infinitesimal amount of weak acid (or , , or ) during an infinitesimal increment in time does not alter the pre-existing equilibrium ratio or Keifer & Roos, 1981. These authors were interested in the transmembrane fluxes of the neutral weak acid DMO and its conjugate weak base. Keifer & Roos assumed that — at the end of the infinitesimal time increment during which and entered/left the cell — the cytosolic , and re-equilibrated. They used this variation on the BDW approach to estimate and , finding that the plasma membrane permeability to is higher than to . In the appendix of their paper, Keifer & Roos derived their refined version of the BDW equations for :
where
The first of these two equations is identical to that presented by BDW, whereas the second includes the new terms ().
In the appendix of their paper, Keifer and Roos noted that it is possible to derive a comparable pair of equations for . In their 1981 review, Roos & Boron reported the two equations of the refined BDW model for (see their (25) and (26)):
where
Note that as approaches zero, or as approaches infinity, the new terms () approach zero, and thus the refined BDW model collapses to the original BDW model. Keifer & Roos did not consider net efflux (i.e., active acid extrusion).
In order to test whether the Keifer-Roos refinement really improves the predictions of the original BDW model, we implemented the BDW model with and without the Keifer-Roos refinement and employed the two models to simulate how changes when (a) only enters the cell (i.e., permeability is zero) and (b) no proton pumping is present (i.e., is zero). With assumptions ‘a’ and ‘b’, the BDW model will take the system to equilibrium, where we can compare the predicted final with that produced by the Davenport diagram — our gold standard for the value of at equilibrium. Using the parameter values of Table 1 and Table 2, we find that both the refined BDW model and the Davenport diagram predict final values of , whereas the original BDW model predicts a slightly lower of — lower because the original BDW does not incorporate as much self-buffering (as Keifer & Roos termed it).
Extending BDW to an epithelial cell.¶
As part of their study of the permeability of gastric-gland cells, Waisbren et al. (1994) extended the BDW model in three ways. First, rather than integrating two equations (i.e., and ), they integrated three equations (i.e., , , and ). Second, after each step of the integration, they re-equilibrated , , and in the cytosol — a maneuver equivalent to the Keifer-Roos extension of the BDW model. Third, they modelled two separate extracellular fluids (each an infinite reservoir), the equivalent of a luminal solution facing one half of the cell that represented the apical membrane, and a basolateral solution facing the other half of the cell that represented the basolateral membrane. The geometry of the cell was still cylindrical, like the squid axon. In their review, Boron et al. (1994) provide additional information about this modelling, which showed that — to account for the physiological data in the paper by Waisbren et al. (1994) — the (membrane area) ( permeability) product must be -fold greater for the basolateral than the apical side of the epithelial cell. Thus, this modelling was a critical part of the main conclusion of the paper by Waisbren et al. (1994), which identified the apical membranes of the gastric chief and parietal cells as the first known membranes with negligible permeability.
Model of regulation by the Vaughan-Jones group.¶
Following BDW’s approach for modelling the net rate of acid addition into the cytosol, Leem et al. (1999) developed the first comprehensive mathematical model of regulation in cardiac myocytes. They simulated the experimentally observed recovery from an acid load (obtained with the ammonium prepulse technique, introduced by BDW) or a base load (obtained with the acetate-prepulse technique). In their model, Leem and coworkers incorporated acid-base fluxes mediated by four sarcolemma transporters, two acid extruders (i.e., - exchangers and cotransporters) and two acid loaders (i.e., - exchangers and hypothetical - exchangers) as well as a time-dependent intracellular buffering by the system. Their simulations nicely reproduced the changes that they observed in rather complex physiological experiments.
Model of diffusion by the Vaughan-Jones group.¶
Richard Vaughan-Jones and his colleagues have extensively investigated the spatial variability of intracellular diffusion in a series of experimental studies (e.g., confocal measurements of ) complemented by mathematical modelling Vaughan-Jones et al., 2002Zaniboni et al., 2003Swietach et al., 2003. These authors developed two-dimensional diffusion models of intracellular diffusion from a constant source to estimate the apparent diffusion constant of in cardiac ventricular myocytes. They found that the apparent mobility is low in cardiac cells because of the presence of mobile and immobile buffers. Later, Swietach and coworkers developed two-dimensional reaction-diffusion models, which they solved numerically in the longitudinal direction only Swietach et al., 2003Swietach et al., 2005. This approach allowed them to investigate further the spatial variability of intracellular diffusion via a proposed shuttling among intracellular mobile buffers. More recently, this group developed a reaction-diffusion model of regulation in tumor spheroids to study the role of carbonic anhydrase (CA) IX in facilitating excretion from tumor cells, thus favoring tumor-cell survival and proliferation Swietach et al., 2008Swietach et al., 2009.
Models of the Gros group.¶
Gerolf Gros and his colleagues have employed mathematical modelling of acid-base physiology to estimate the apparent membrane permeability of the red blood cell (RBC) membrane to () and (), as well as extra- and intracellular CA activity, from mass-spectrometric data and the O-exchange technique Wunder et al., 1997. Their compartmental model comprises a system of six ordinary differential equations that account for the O-exchange reactions among , and in both the intracellular and extracellular compartments, and for transfer of and across the intracellular and extracellular compartments.
In 2005, Endeward & Gros used this modelling framework to estimate and of the apical membrane of the colonic epithelium in guinea pigs. They found that the apical membrane has a quite low , possibly because of the absence of membrane proteins (e.g., aquaporin 1, AQP1) that act as conduits (i.e., gas channels) for the movement of .
In 2009, these same authors extended their previous compartmental model of O-exchange to a one-dimensional reaction-diffusion model of an RBC surrounded by an extracellular unconvected layer that can exchange solutes with a well-stirred bulk solution Endeward & Gros, 2009. They employed the model to assess the influence of intracellular and extracellular unconvected layers in the estimation of RBC membranes. The combination of physiological experiments and modelling allowed them to estimate the extra- and intracellular unstirred layers, and confirm their earlier conclusions that the of RBC membranes is low in the absence of gas channels.
Models of the CWRU group.¶
With the goal of assessing the role of channels in producing decreases in and transient alkaline transients in extracellular-surface () as enters a Xenopus oocyte Endeward et al., 2006Musa-Aziz et al., 2009Musa-Aziz et al., 2010, the group at Case Western Reserve University developed a three-dimensional reaction-diffusion model of influx into a spherical cell Somersalo et al., 2012. The model accounts for a multitude of buffer reactions, as well as solute diffusion within the unconvected intracellular fluid (ICF) and the extracellular unconvected fluid (EUF) that surrounds the cell. Because the electrophysiologists established experimental conditions in which only can cross the plasma membrane, the modellers allowed only to diffuse between the ICF and EUF. However, all solutes can diffuse within the ICF as well as between the EUF and the bulk extracellular fluid (bECF) that surrounds the EUF. Somersalo and coworkers employed the model to investigate a variety of theoretical conditions. For example, they investigated how changes in the width of the EUF or in the value of affect and transients. The model suggests that, in oocytes, the background (i.e., in the absence of channels) must be low in order for channels to have the observed effect on the maximal rate of intracellular acidification, , or on the maximal height of the alkaline transient, .
In 2014, Occhipinti and coworkers extended the above theoretical model and employed it to investigate the role that CA II and CA IV play in enhancing transmembrane fluxes Musa-Aziz et al., 2014Musa-Aziz et al., 2014Occhipinti et al., 2014. The model was able to mimic the and experiments under a variety of experimental conditions. Moreover, it provided a novel observation that the effects on cytosolic CA II and extracellular CA IV on and are supra-additive. More recently, the CWRU group further extended the oocyte model to include the special microenvironment that the electrode creates when pushed against the oocyte membrane to record the transient Calvetti et al., 2020. This model, solved using finite-element method, predicts that the special microenvironment between the blunt tip of the electrode and the oocyte membrane greatly amplifies the alkaline as enters the cell.
6The Future¶
The advice, “It’s tough to make predictions, especially about the future” is attributed to Yogi Berra, amateur philosopher and legend of American baseball. His advice certainly applies here. We imagine that the grand goal of the acid-base segment of world-wide biomedical science would be to model whole-body acid-base regulation on a cell-by-cell basis for a wide range of physiological and pathophysiological challenge. Of course, we all expect modelling to provide new insights that we can test in physiological experiments. However, the models should ultimately provide medical professionals with powerful insights into pathophysiology. These goals represent an enormous challenge in multi-scale modelling — multiscale in both time and space. That lofty goal is probably many decades away, even with appropriate grant support and continued development of computer hardware and software.
The process could begin at the single-cell level. Even here, the challenge is daunting because, as pointed out by Occhipinti & Boron (2015), the movements of acid-base equivalents across the plasma membrane create complex interdependencies among buffer reactions, diffusion processes and transporter mechanisms. Intuition alone cannot explain or predict the consequences of these numerous simultaneous processes on . Advanced mathematical models — including details on processes or more complex geometries — performed in conjunction with physiological experiments can provide a useful tool to make predictions and provide mechanistic explanations on observed changes. The community has begun taking the first steps with cells of simple geometry, like oocytes (i.e., spheres). Even so, current models only simulate passive diffusion of substances like and lactic acid. Current reaction-diffusion models do not yet incorporate electrodiffusion, let alone electrodiffusion through particular channels. The next step might be to incorporate specific acid-base transport systems (e.g., cotransporters, - exchangers) with characteristic kinetic properties. Ultimately, one would like to include all traffic (i.e., including non-acid-base traffic) across the cell membrane, and the regulation of this traffic. Model validation would require extensive physiological studies on these simple cell systems to verify that simulations of ion conductances (and ultimately channels) and transporters are reasonable. Moreover, even these relatively simple first steps would require efficient computational methods for solving the governing equations.
From the level of the single cell, the field must move in two opposite directions. In the reductionist direction, we must move from whole cells to nanodomains (i.e., the mesoscopic level) and ultimately to single molecules. At first, we guessed at even the number of different kinds of such proteins. Now we know the identities of specific proteins and their amino-acid sequences, and sometimes even their structures. Even so, cells are often capable of making many protein variants from one gene, and regulating these in response to myriad influences. Even the protein structures are merely an early step in understanding transporter mechanism.
In the opposite, integrative direction, the community must extend the models to more complex geometries, like the cells of a simple epithelium, spindle-shaped cells (e.g., muscle), cylinders (e.g., axons, dendritic processes), and ultimately cells of very complex geometry (e.g., neurons, astrocytes). In real life, even a simple epithelium is far more complex than a cuboid. The apical membrane (nearest the tight junctions) of a cell of the proximal tubule (PT) has microvilli. The basolateral membrane (nearest the blood) has infoldings. These complexities create nanoenvironments on both the cytosolic and extracellular sides of the membrane. The tight junctions that separate the PT cells allow the paracellular diffusion of water and certain solutes. The lateral interspaces between adjacent cells also creates nanoenvironments. All of these nanoenvironments are likely to be extremely important for the physiology of the epithelium. Further complexities include the myriad cellular organelles. These not only affect diffusion of water and solutes through the cytosol (by creating tortuosity), but also can engage in their own collection of transport processes that can affect and contribute to buffering. Several groups have developed either steady-state or time-dependent compartmental models of acid-base transport in PTs and, more generally, solute reabsorption in various segments of the nephron Thomas & Dagher, 1994Krahn & Weinstein, 1996Thomas et al., 2006Weinstein et al., 2007Weinstein & Sontag, 2009Layton & Edwards, 2014Edwards & Layton, 2017. Bransen and coworkers have used a finite-element approach to solve the partial differential equations that describe their reaction-diffusion model of - exchange in the microvilli of the PT Brasen et al., 2014.
Beyond cells, we must create whole tissues (e.g., a cylindrical PT) from many cells, and create whole organs (e.g., a kidney) from these tissues. Models of a whole kidney, for example, would have to consider the interstitial fluid and capillaries, the complex exchanges of substances among them, and the changes in composition that occur as fluids flow along the lumens of the tubule and capillaries. Comparable modelling must extend to every tissue and organ, and consider the complex intercommunications among them via the circulatory and neuro-endocrine systems as well as metabolism. In the case of acid-base physiology, one must model the specific roles that certain organs — the brain, the pulmonary system, and the kidneys — play in whole-body regulation.
All along the way, community members must cooperate because no one group of investigators can possibly accomplish the ultimate goal single-handedly. The models must be open source, modular and sharable — and the community must share. Finally, in order for modelling approaches to be consistent, the community must establish standard approaches for gathering physiological data, and agree on standard values for physiological parameters. Here, an organisation like the International Union of Physiological Sciences (IUPS)[31] can play a critical role.
8Conflict of Interests¶
None declared.
Acknowledgments¶
Supported by National Institute of Health (NIH) grants K01-DK107787 (to RO), U01GM111251 and DK113197 (to WFB), Office of Naval Research (ONR) grant N00014-15-1-2060 (to WFB), and a Multidisciplinary University Research Initiative (MURI) grant N00014-16-1-2535 from the DoD (to WFB). WFB gratefully acknowledges the support of the Myers/Scarpa endowed chair. SS acknowledges the financial support provided by the Aotearoa Foundation.
Of course, other energy-requiring processes — as yet undiscovered at the time — could also have accounted for the increase.
Note that denotes a concentration in units of or .
The definition of as a scale based on powers of 10 was introduced by Sørensen in an attempt to simplify the notation of and to avoid having to resort to decimals for tiny amounts of Sørensen, 1909.
Note that although is often regarded as an acid, the true weak acid is , the product of the reaction of with .
BDW used to denote a transmembrane flux. In the present paper, we use the more widely used notation .
In the BDW paper, the final equation before the Appendix had a typographical error, omitting the in the following equation: , where is a flux (today represented as ). Indeed, WFB had included in his extant Fortran code.
BDW arrived at the value by rearranging the equilibrium relation for outside the cell: , when (or equivalently, ), and .
Note that BDW did not report the value of , but rather the value of the fiber (i.e., axon) diameter, which is equal to and corresponds to a of or .
Note that BDW did not report the value of . This value is inferred from a of (reported in the legend of BDW’s figure 6) and a of (in the original Fortran code). Although BDW reported a value for of (taken from Harned & Davis (1943), referring to at º ), they used this value only for the Davenport diagram in their figure 5A.
BDW did not report the value of , but in the Fortran code, used , which matched the measured mean value.
In figure 6A of the BDW paper, the proton pumping rate constant () had values of , , , , or .
BDW did not report this value, but in their code used .
BDW arrived at the value by rearranging the equilibrium relation outside the cell: , when , , and (or equivalently, ). The extracellular concentration can be obtained as
In their original Fortran code that generated the plots in their figure 6B, BDW used (the value used in some of their early experiments), and not as (the value used in their later experiments) indicated in the figure and legend for figure 6B in their original paper.
In figure 6B of the BDW paper, had values of , , , and .
BDW did not report the value of , but in their code used , which matched the measured mean value.
In their code BDW used the value of and not , which they used when constructing the Davenport diagram in their figure 5B.
See equation 24 in Derivation of the Jacobs Neutral Weak Acid Equation (Log) in Supplementary Files
See equation 25 in Derivation of the Jacobs Neutral Weak Acid Equation (Linear) in Supplementary Files
See equation 33 in Derivation of a Jacobs-like Equation (Log) for a Neutral Weak Base in Supplementary Files
See equation 21 in Derivation of a Jacobs-like Equation (Linear) for a Neutral Weak Base in Supplementary Files
See equation 26 in Derivation of the Orloff-Berliner Equation in Supplementary Files
See equation 55 in Derivation of the Orloff-Berliner Equation in Supplementary Files
See equation 26 in Derivation of the Equation of Milne et al for a Neutral Weak Acid in Supplementary Files
See equation 62 in Derivation of the Equation of Milne et al for a Neutral Weak Acid in Supplementary Files
See equation 45 in Derivation of the Roos Equation in Supplementary Files
See equation 44 in Derivation of the Boron-Roos Equation in Supplementary Files
- Boron, W. F., & De Weer, P. (1976). Intracellular pH transients in squid giant axons caused by CO₂, NH₃, and metabolic inhibitors. Journal of General Physiology, 67(1), 91–112.
- Boron, W. F. (2010). Sharpey-Schafer lecture: gas channels. Experimental Physiology, 95(12), 1107–1130.
- Cuellar, A. A., Lloyd, C. M., Nielsen, P. F., Bullivant, D. P., Nickerson, D. P., & Hunter, P. J. (2003). An Overview of CellML 1.1, a Biological Model Description Language. Simulation, 79(12), 740–747.
- Bergmann, F. T., Nickerson, D., Waltemath, D., & Scharm, M. (2017). SED-ML web tools: generate, modify and export standard-compliant simulation studies. Bioinformatics (Oxford, England), 33(8), 1253–1254.
- Roos, A., & Boron, W. F. (1981). Intracellular pH. Physiological Reviews, 61(2), 296–434.
