spacer gif spacer gif spacer gif spacer gif spacer gif
 QUICK SEARCH:   [advanced]


spacer gif
     Home     Help     Feedback     Subscriptions     Archive     Search     Table of Contents    

This Article
Right arrow Summary Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Burton, R. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Burton, R. F.
The Journal of Experimental Biology 205, 3587-3600 (2002)
Copyright © 2002 The Company of Biologists Limited


Perspective

Temperature and acid—base balance in ectothermic vertebrates: the imidazole alphastat hypotheses and beyond

Richard F. Burton

Institute of Biomedical and Life Sciences, Thomson Building, University of Glasgow, Glasgow G12 8QQ, UK

e-mail: R.F.Burton{at}bio.gla.ac.uk

Accepted 22 August 2002

Summary

The `imidazole alphastat hypothesis' states that intracellular and extracellular pH, partly via buffering by imidazole groups, change with temperature in a way that keeps imidazole and protein ionization constant, thus maintaining cell function and minimizing shifts of base equivalents and total CO2, while adjustment of PCO2 involves imidazole-based receptors. `The hypothesis', which is actually several hypotheses, has been variously perceived and judged, but its underlying conceptual framework remains largely valid, and is reformulated using differential equations requiring less information input than their integral equivalents. Their usefulness is illustrated with published data on temperature responses in fish cells and whole tetrapods. Mathematical modelling allows general principles to be explored with less immediate concern for uncertainties in experimental data and other information. In tetrapods, it suggests that warming is followed by a loss of base equivalents from the body, and that this loss is due to metabolic adjustments that are not part of pH homeostasis. Uncertainties include intracellular buffer values, local variations in PCO2 within the body, the possible role of buffering by bone mineral, and the temperature dependence of pK values for CO2/HCO3- and imidazole groups. The equations utilize a single, notional, temperature-dependent pK value for all non-bicarbonate buffers in a given body compartment. This approximates to the `passive component' of pH adjustment to temperature change as measured by the homogenate technique. Also discussed are the diversity of cell responses within individual animals, relevant aspects of the control of ventilation, metabolism and transmembrane transport, and the basis of optimum pH—temperature relationships.

Key words: acid—base balance, alphastat, imidazole alphastat, temperature, pH, carbon dioxide tension, vertebrate, fish, amphibian, reptiles, cell pH

Introduction

The `imidazole alphastat hypothesis' of Reeves (1972Go) has inspired many studies on the effects of temperature on acid—base balance. Some writers accept it without question, while others (e.g. Heisler, 1986cGo) are dismissive. However, it is not one hypothesis, but several, and some of these may apply to varying degrees in different species and tissues. It should be judged accordingly. The purpose of this paper is not to review the extensive literature, but more to analyze and comment on the constituent hypotheses. To help in this, a simple mathematical model is presented. Its use is illustrated, sometimes with new light shed on old data. Discussion of the model's limitations raises other issues that are equally relevant to some published experimental studies. This paper is not concerned with effects of temperature beyond critical limits, where steady-state function is disturbed.

What Reeves (1972Go) describes as a "general conceptual framework" was put forward mainly as a basis for computing relationships between pH, CO2 tension (PCO2) and total CO2 content (CCO2). It can be regarded as having the following components, of which those marked with asterisks are highlighted by Reeves himself as conclusions. Components (1) and (2) are not so much proposals of `the hypothesis' as fundamental facts underlying it. (12) is a natural addition (Reeves, 1985Go).

  1. Equilibrium constants (pK) of chemical reactions are generally temperature-dependent, including those for the protonation of imidazole groups (pKIm), i.e. for the equilibrium: Im+H+=ImH+. The fractional dissociation, [Im]/([Im]+[ImH+]), is the `alpha' in `alphastat'.
  2. Intracellular and extracellular non-bicarbonate buffering is dominated by imidazole groups, notably those of protein histidines.
  3. * In a solution buffered by comparable amounts of imidazole and CO2/HCO3- in closed system, imidazole buffering dominates, so that {Delta}pH/{Delta}T is close to {Delta}pKIm/{Delta}T, where T is temperature.
  4. * The temperature dependence of pH in extracellular and intracellular fluids in vivo is similar to the temperature dependence of pKIm (i.e. {Delta}pH/{Delta}T is close to {Delta}pKIm/{Delta}T). Equality implies constant imidazole alpha, and thus the preservation of protein net charge.
  5. * Preservation of protein net charge assures "optimal enzyme activities, protein conformation and structural stability, as well as stable ion and water distribution across cell membranes based on Donnan principles". (Unlike protein net charge, Donnan equilibria are not now seen as very relevant here. This is because the distribution of most ionic species across cell membranes is influenced by special transport mechanisms.)
  6. * "If ventilation is regulated to maintain a constant alpha for one compartment, say blood, alpha imidazole will be maintained for other compartments, regardless of the quantitative makeup of the total intracellular buffer value."
  7. * "Alphastat control ensures that changes in total carbon dioxide stores are small with body temperature changes, and thus the transients of loading and unloading such stores, and the associated disturbances of intracellular acid—base state are minimized." The stores consist of HCO3-, CO32- and dissolved CO2 (but see Discussion) and the transmembrane movements of HCO3- and CO2 are best treated here as distinct processes (see items 9-10).
  8. Within any body compartment, regulation of PCO2 to produce this constant imidazole ionization as temperature varies should result in nearly constant [HCO3-].
  9. Therefore adaptation to varying temperature should not require transfer across cell membranes of acid—base relevant substances other than CO2 (i.e. of HCO3-, H+, OH-, or weak acids and bases). With constant imidazole ionization in all body compartments, no such exchanges would be required with the environment either.
  10. The generation and breakdown of HCO3- involves one-for-one loss and production of CO2. Preservation of imidazole ionization therefore minimizes transient changes in the body's content of dissolved CO2.
  11. * Key protein imidazole groups might function as respiratory receptors governing ventilation.
  12. Key protein imidazole groups might also function as receptors involved in regulating transmembrane transport of HCO3- and H+, and hence both cell pH and the acid—base function of kidneys and gills.

Underlying the whole scheme is an attractive idea that is implied rather than stated: to the extent that protein properties are all similarly affected by pH and temperature through the dominant role of imidazole groups, each imidazole group facilitates pH homeostasis for all others.

Reeves (1972Go) did not propose that all conditions apply exactly. Thus item (7) refers to changes that are "small" or "minimized"; in bullfrogs he found no dependence of CO2 content on temperature in blood and cardiac muscle, but an increase with warming in liver, and possibly (0.05<P<0.1) in striated muscle. He acknowledged that inorganic phosphate, and to a small extent N-terminal {alpha}-amino groups, contribute to intracellular buffering (item 2). He did not refer to the small imidazole compounds, e.g. carnosine, that abound in some cells (Crush, 1970Go; Burton, 1983Go).

Items (4) or (9) are most often singled out as `the hypothesis'. Reeves (1976aGo) put the focus more clearly on preservation of protein net charge (Z), as had Stadie et al. (1925Go) and Austin et al. (1927Go) in relation to serum and blood. Many authors have since done the same. Indeed, Cameron (1989Go) suggested that a `Z-stat' model might be more appropriate than an imidazole alphastat model. But is Z exactly the key variable? Where buffering is the issue, i.e. in relation to changes in [HCO3-], it is the ionization and net charge of all the non-bicarbonate buffers that matters. The properties of enzymes may depend largely on the states of specific imidazole, and other, groups, but enzyme properties are also much affected by the crowding effect of other proteins present, and that may itself be affected by net charge (Garner and Burg, 1994Go; Elcock and McCammon, 2001Go).

The temperature dependence of pK values
Individually, imidazole groups on proteins have diverse values of {Delta}pKIm/{Delta}T, ranging from -0.01 to -0.02°C-1 in myoglobin (Bhattacharya and Lecomte, 1997Go). The value for one histidine in ribonuclease even changes, in vitro, from -0.01 to -0.05°C-1 above 32°C (Roberts et al., 1969Go). Therefore values of {Delta}pKIm/{Delta}T can hardly match those of {Delta}pH/{Delta}T in every instance (item 4). Despite this variation, one may define a single quantity, {Delta}pKprot/{Delta}T, that characterizes the temperature dependence of protein net charge. This is not a simple mean of individual values, but, for proteins with net charge Z and combined buffer value ßprot, is such as to satisfy the equation:

(1)
For individual buffer groups, and small temperature changes, the value of {Delta}pK/{Delta}T is calculable from the standard enthalpy change, {Delta}H°, as {Delta}H°/(2.303RT2), where R is the gas constant (1.99 cal °C-1 mol-1; 1 cal=4.186J) and T is the absolute temperature.

As discussed by Heisler (1986cGo), various values of {Delta}pKprot/{Delta}T have been assumed in the interpretation of experimental results, notably -0.018°C-1 as calculated for histidine at 15°C (Edsall and Wyman, 1958Go), and -0.021°C-1, which lies between values for imidazole and 4-methylimidazole (i.e. -0.020 and -0.022°C-1, respectively), from the same source and again at 15°C. Measurements on actual proteins are preferable, but few are available, not all of these representative. For plasma proteins near pH7, Reeves (1976bGo) found {Delta}H° to be 6940 cal mol-1, corresponding to a value of {Delta}pKprot/{Delta}T of -0.018°C-1 at 20°C. For oxyhaemoglobin near pH7, he found respective values of 7300 cal mol-1 and -0.019°C-1. For oxyhaemoglobin at pH 6.8-8.0, Wyman had found {Delta}H° to be near 6000-7000 cal mol-1 (Edsall and Wyman, 1958Go). For proteins extracted from white skeletal muscle of the eelpout Zoarces viviparus, van Dijk et al. (1997Go) obtained a less negative value of -0.013°C-1.

Some vertebrate tissues are also buffered by substantial amounts of low molecular mass imidazole compounds, such as L-histidine itself and the dipeptides carnosine, anserine and balenine (=ophidine) (Crush, 1970Go; Burton, 1983Go; Abe, 2000Go). In fish white muscle, containing amounts between 0 and 148 mmol kg-1 wet mass, these compounds may be more concentrated than in red muscle, be responsible for much of the buffer capacity, and show a strong correlation with lactic dehydrogenase activity (Abe, 2000Go). As calculated from figure 2 of Abe and Okuma (1991Go), the value of {Delta}pKIm/{Delta}T for carnosine is -0.018°C-1, while values for anserine and balenine are both -0.013°C-1. From the same figure, the value for histidine is about -0.020°C-1 (compared with -0.018°C-1 as given above). Hitzig et al. (1994Go) obtained values of -0.0166°C-1 for imidazole (compared with -0.02°C-1 as given above) and -0.0154°C-1 for carnosine. For several of these compounds there is thus significant disagreement.

Although buffering is dominated by imidazole groups, and concentrations of inorganic phosphate are often very low (e.g. Marjanovic et al., 1998Go), the latter may sometimes contribute significantly to intracellular buffering (but note that concentrations are sometimes artefactually raised in tissue samples by hydrolysis of organic phosphates; Heisler and Neumann, 1980Go). Values of {Delta}pK/{Delta}T for inorganic phosphate are small, i.e. -0.006°C-1 for 0-5°C and -0.001 for 35-40°C-1 (Seo et al., 1983Go), or -0.003°C-1 overall, as also found by Alberty (1972Go). Adenosine triphosphate (ATP), because it is present mainly as MgATP, contributes little to buffering. Phosphocreatine, with its pK of 4.6 (Edsall and Wyman, 1958Go), is not a significant buffer at cell pH. Although taurine has a high pK (8.8 at 37°C), some cells contain enough to make this a significant buffer, and {Delta}pK/{Delta}T is probably strongly negative (Bevans and Harris, 1999Go).

Temperature effects are modelled below in terms of a single, notional, non-bicarbonate buffer with a temperature-dependent pK denoted pK*. It stands for all the different non-bicarbonate buffers and buffer groups that might be present. With these denoted by subscripts A, B etc., {Delta}pK*/{Delta}T is in principle calculable as follows (Burton, 1973Go):

(2)

In general, data now available give only approximate values in this way, but estimates have been obtained more directly from the effects of temperature on tissue homogenates. Thus, data of Heisler and Neumann (1980Go) for the dogfish, Scyliorhinus stellaris, imply mean values for {Delta}pK*/{Delta}T of -0.017°C-1 in white skeletal muscle, -0.012°C-1 in red skeletal muscle and -0.016°C-1 in heart muscle. (I have derived these by recombining separate estimates for phosphate-like and imidazole-like buffers, themselves originally calculated using assumed values for {Delta}pK/{Delta}T of -0.002°C-1 and -0.021°C-1, respectively.) The concentrations of inorganic phosphate in the homogenates used could have been raised by hydrolysis of ATP and phosphocreatine, leading to underestimates of {Delta}pK*/{Delta}T. Heisler and Neumann argued against that and it is also possible that some inorganic phosphate was precipitated by the added calcium (approximately 3.3 mmol l-1).

In a simple solution of non-bicarbonate buffers, free of HCO3-, {Delta}pH/{Delta}T={Delta}pK*/{Delta}T. Therefore {Delta}pK*/{Delta}T may be estimated by measuring {Delta}pH/{Delta}T. A similar method has been applied to tissue homogenates, with the value of {Delta}pH/{Delta}T described as the `passive component' of the pH adjustment to temperature change (Pörtner, 1990Go; Pörtner et al., 1990Go; Pörtner and Sartoris, 1999Go). In fact the original CO2 and HCO3- are retained in the closed system, but this has little effect on {Delta}pH/{Delta}T (Burton, 1973Go; van Dijk et al., 1997Go), so that this should still provide a good estimate of {Delta}pK*/{Delta}T (see item 3 of the hypothesis). An important aspect of this method is that continuing metabolism and the hydrolysis of organic phosphates are prevented by rapid dissection and freeze-clamping of the tissue, followed by homogenization in a solution containing KF and nitrilotriacetic acid that precipitate and chelate Mg2+ and Ca2+. The method is also used for measuring non-bicarbonate buffer values and in vivo pH.

For skeletal muscle of the black racer snake, Coluber constrictor, the value for {Delta}pK*/{Delta}T thus found averages -0.013°C-1 (Stinner et al., 1998Go). Corresponding values in both the cane toad Bufo marinus and the bullfrog Rana catesbeiana average -0.011°C-1 (Stinner and Hartzler, 2000Go), and a similar value applies in tissues of the eelpout Pachycara brachycephalum (Pörtner and Sartoris, 1999Go). For white skeletal muscle of Z. viviparus the value of {Delta}pK*/{Delta}T is about -0.006°C-1 (van Dijk et al., 1997Go). It is not known what component makes this value so much less negative than that of {Delta}pKprot in the same muscle (-0.013°C-1). All these estimates of {Delta}pK*/{Delta}T could have been made slightly less negative by the conversion of MgATP, a poor buffer at intracellular pH, to ATP, which is more like inorganic phosphate in its buffer properties (Alberty, 1972Go), but the concentration in Z. viviparus averages only 3.6 mmol kg-1 fresh mass (van Dijk et al., 1997Go).

There is evidently variation amongst tissues and species, but, for modelling of hypothetical, generalized cells, it seems that representative values of {Delta}pK*/{Delta}T, at least for muscle, may be taken as about -0.011 to -0.018°C-1. Muscle containing much carnosine may be expected to yield values in the more negative part of this range. Values need to be more accurately known for calculations on real experimental data.

Other background information
It is a general finding that arterial blood pH falls with increasing temperature. According to the compilation of Heisler (1986cGo), values of {Delta}pH/{Delta}T in plasma are mostly -0.005 to -0.018°C-1 in reptiles, -0.011 to -0.018°C-1 in amphibians and -0.007 to -0.014°C-1 in fish, with respective mean values of approximately -0.011, -0.015 and -0.012°C-1. Separate studies on a given species have sometimes yielded very disparate results. Intracellular values of {Delta}pH/{Delta}T, again as compiled by Heisler (1986cGo), have been found to vary between 0 and -0.029°C-1, with a mean of approximately -0.012°C-1. Again there are marked differences amongst tissues and species, and values may vary with temperature range (Boutilier et al., 1987Go; Stinner and Hartzler, 2000Go). Heisler (1986cGo) discusses sources of error. The compilation of Ultsch and Jackson (1996Go) includes more recent data on both plasma and cells.

Within cells, there is thus some similarity between values of {Delta}pH/{Delta}T and overall {Delta}pKIm/{Delta}T, but exact comparisons are not yet generally possible. Despite this, it is evident that overall imidazole ionization, dependent on the difference, must increase with cooling in some cells and decrease in others. In the case of white skeletal muscle of Z. viviparus, values of both {Delta}pKprot/{Delta}T for dialyzed homogenates and intracellular {Delta}pH/{Delta}T have been measured, averaging -0.013°C-1 (as already noted) and -0.016°C-1, respectively, a difference that is not statistically significant. Thus the net charge on these proteins, and the dissociation state of their histidine residues, hardly varies with temperature (van Dijk et al., 1997Go). The fractional dissociation of carnosine imidazole has been measured directly, by proton NMR, in tail muscle of intact unanaesthetized newts (Notophthalmus viridescens); it was independent of temperature between 10 and 30°C, thus showing alphastat regulation for carnosine (Hitzig et al., 1994Go).

Concentrations of HCO3- in arterial plasma are typically 14-40 mmol l-1 in reptiles, with the higher values mostly occurring in chelonians, 10-30 mmol l-1 in amphibians, and 3-15 mmol l-1 in water-breathing fish (Toews and Boutilier, 1986Go; Heisler, 1986bGo; Ultsch and Jackson, 1996Go). The generally lower concentrations in water breathers relate to the fact that PCO2 in these is typically lower than in air breathers (Rahn, 1967Go). There is a general tendency for plasma [HCO3-] to fall with warming, especially in fish, but as discussed more fully below, the trend is usually small or absent in tetrapods (Heisler, 1986bGo,cGo; Ultsch and Jackson, 1996Go). It is reversed in adults, but not juveniles, of S. stellaris (Heisler et al., 1980Go).

It is harder to generalize about intracellular HCO3-, especially since absolute concentrations and temperature effects both differ amongst tissues. For modelling, however, reasonable representative values may often suffice. Figures 3-7 of Ultsch and Jackson (1996Go) suggest that the cytoplasmic pH of muscle and liver is commonly about 0.3-0.6 unit lower than arterial plasma pH. If PCO2 is assumed to be not much higher in the cells than in arterial plasma (but see below), then the Henderson—Hasselbalch equation suggests that [HCO3-] in cells should often be about a quarter or half that in plasma.

Intracellular non-bicarbonate buffer values vary with tissue and species, with measurements in skeletal and cardiac muscle of ectothermic vertebrates being typically 25-110 mequiv l-1 cell water pH unit-1 (Castellini and Somero, 1981Go; Heisler, 1986aGo; Milligan and Wood, 1986Go; Abe, 2000Go). Many such values must have been significantly raised artefactually by the release of inorganic phosphate from phosphocreatine and ATP (Pörtner, 1989Go, 1990Go). With the homogenate technique (Pörtner et al., 1990Go), this hydrolysis is avoided (see above) and in gastrocnemius muscle of B. marinus the non-bicarbonate buffer value has been found in that way to be 29.8 mequiv l-1 cell water pH unit-1 (Pörtner, 1990Go). Measured similarly, values in white muscle of Z. viviparus average 31 mequiv l-1 cell water pH unit-1, being slightly higher at 12°C than at 0°C (van Dijk et al., 1997Go). A small positive error results from the conversion of MgATP to ATP (see above) but, in Z. viviparus, the ATP concentration of 3.6 mmol kg-1 fresh mass (van Dijk et al., 1997Go) implies a maximum contribution to the buffer value of 2.1 mequiv l-1 fresh mass pH unit-1 (i.e. 3.6x0.575; Burton, 1973Go).

Buffer values in separated plasma, about 6-8 mequiv l-1 pH unit-1 in mammals, are usually lower in plasma of ectothermic vertebrates and lower still in their interstitial fluid (Cameron and Kormanik, 1982Go; Heisler, 1986aGo; Tufts and Perry, 1998Go).

Other general findings are discussed below. These include the tendency for PCO2 to rise with increasing temperature and, in many tetrapod species, for whole-body CO2 content to fall.

Modelling the effects of temperature
In thinking generally about the effects of temperature on acid—base balance, as opposed to fully analyzing experimental data, it is helpful to start with a single-compartment model described mainly in terms of differential equations expressing changes per degree rise in temperature. This reduces the number of variables and allows easy exploration of the effects of changing them, but the equations are strictly valid only for small temperature changes. The equations are based on well-established principles (e.g. Edsall and Wyman, 1958Go; Reeves, 1972Go; Burton, 1973Go; Heisler and Neumann, 1980Go; Heisler, 1986aGo), but some may be novel. The model is in line with the general conceptual framework of Reeves (1972Go), but lacks a restrictive focus on imidazole buffering (item 2 of the hypothesis).

This model of what could be cells or extracellular fluid has a constant volume and contains a single, notional, non-bicarbonate buffer having a buffer value ß* and a temperature-dependent pK, denoted pK* as above. Also present is HCO3-, for which the buffer value in open system, ßbic, is 2.303[HCO3-]. For simplicity, the small amounts of CO32- included by Reeves (1972Go) are ignored here. With the principle of electroneutrality in mind, we may define a quantity N, reflecting the concentrations and net charge of all other (non-buffer) ions present, such that, in terms of equivalents:

(3)
where [H+]bound corresponds to the protonated form of the buffer. The term {Delta}N, when positive, can, for example, match movements of HCO3- out of the compartment or the production of dissociating lactic acid within it (see below). Effectively, HCO3- movements in one direction are equivalent to movements of CO32- or OH- in the same direction or of H+ in the opposite direction. All are conveniently described as movements of `base equivalents', so that an increase in N can equate to a loss of base equivalents. Such changes may be part of the acid—base adjustment or unrelated responses to temperature change. In accordance with the definition of buffer value, buffer ionization depends on the difference between pK* and pH (Heisler, 1986aGo), such that:

(4)
From equations 3 and 4,

(5)
The total CO2 content, CCO2, equals ([HCO3-]+[CO2]), where:

(6)
where S is the solubility coefficient of CO2. Therefore,

(7)
Given Equation 6, one form of the Henderson—Hasselbalch equation is:

(8)
From this, since {Delta}log[HCO3-] is {Delta}[HCO3-]/(2.303[HCO3-]) and 2.303[HCO3-] is ßbic,

(9)
It is often useful to treat {Delta}pK1'/{Delta}T-{Delta}logS/{Delta}T as a single term. It happens to be nearly constant at 0.0053°C-1 between 0 and 40°C (as calculated from data on S and pK1' given by Reeves, 1976bGo) and this value is assumed in all calculations below. Two equations follow from Equations 5 and 9:

(10)
and

(11)
Another form of the Henderson—Hasselbalch equation is:

(12)
At typical extracellular and intracellular pH values, the ratio of CCO2 to [HCO3-] is little greater than 1, so that the two may be regarded as interchangable for modelling purposes. Moreover, since PCO2 generally rises with increasing temperature (see below) and S falls, their product, [CO2], tends to be less temperature-sensitive than either, so that {Delta}[CO2]/{Delta}T is often a minor term in Equation 7.

The equations may be adapted to a model of several compartments, each of fractional water volume (or mass) V, such that {Sigma}V=1. Then both {Delta}N/{Delta}T and {Delta}[HCO3-]/{Delta}T may be averaged for the whole, as {Sigma}{V{Delta}N/{Delta}T} and {Sigma}{V{Delta}[HCO3-]/{Delta}T}, respectively. While the terms ß*, ßbic, {Delta}pK*/{Delta}T and {Delta}pH/{Delta}T may differ from one compartment to another, the terms {Delta}pK1'/{Delta}T-{Delta}logS/{Delta}T and {Delta}logPCO2/{Delta}T, in the absence of information to the contrary, may be taken as identical in each. Equations 5 and 9 become:

(13)
and

(14)
Equations 10 and 11 may be adapted likewise.

Initial exploration of the model

Given the set of model equations, what general conclusions may be drawn, using a minimum of specific experimental data and applying other postulates of the alphastat hypotheses? Let us start by modelling either particular cells, or all cells collectively, as a single compartment. The alphastat hypotheses (items 7, 9 and 10) suggest that we consider the special case in which values of both {Delta}N/{Delta}T and {Delta}CCO2/{Delta}T are zero. With the value of {Delta}[HCO3-]/{Delta}T taken as similar to that of {Delta}CCO2/{Delta}T (i.e. zero) in accordance with Equation 12, Equation 11 then reduces to the approximation:

(15)
Then, from Equation 10, {Delta}pH/{Delta}T equals {Delta}pK*/{Delta}T. This accords with the imidazole hypotheses provided that {Delta}pK*/{Delta}T reflects only imidazole buffering. From Equation 15, with {Delta}pK*/{Delta}T lying, say, between -0.011 and -0.018°C-1 (see Introduction), {Delta}logPCO2/{Delta}T is 0.016-0.023°C-1. Inasmuch as losses and gains of CO2 or base equivalents are best minimized, the model thus suggests an appropriate relationship between PCO2 and temperature. These values of {Delta}logPCO2/{Delta}T can also be expressed as Q10 values of 1.45-1.70. These lie within the much wider range for arterial blood tabulated by Heisler (1986cGo), namely 1.17-1.63 in reptiles, 1.08-1.89 in amphibians and 1.04-1.95 in fish. (This overall range corresponds to values of {Delta}logPCO2/{Delta}T of 0.0017-0.029°C-1.) Matches between model and reality cannot establish the correctness of the assumptions for any species, but they do suggest that these could sometimes be about right.

The possibility that water breathers can increase PCO2 as temperature rises has sometimes been dismissed, and muscle PCO2 in Z. viviparus may actually fall with rising temperature (van Dijk et al., 1997Go). In any case, the characteristic fall in plasma pH is partly due, in most fish species studied, to a fall in plasma [HCO3-]. This is discussed below, as also the temperature dependence of PCO2 in reptiles and amphibians.

Item (6) of the hypothesis, reworded in model terms, states that regulation of PCO2 to maintain constant buffer ionization in one compartment maintains buffer ionization in all others. In the imidazole alphastat scheme all values of {Delta}N{Delta}T are zero, while the value of {Delta}pH/{Delta}T in each compartment equals that of {Delta}pK*/{Delta}T. Equation 10 then reduces to Equation 15 (not now an approximation). If values of {Delta}logPCO2/{Delta}T are to be the same in each compartment, so must those of {Delta}pK*/{Delta}T; they are so in the idealized alphastat model since all buffering is due there to similar imidazole compounds. It is noteworthy that buffer values, i.e. {Delta}bic and ß*, do not appear in Equation 15, and are thus irrelevant. In reality, the stated conditions do not generally apply. Nevertheless, regulation of blood PCO2 may well help to regulate pH appropriately in different cell types. Indeed, each type is likely to be adapted — whether in evolution or from day to day — to the characteristic temperature dependence of PCO2 to which it is exposed.

The model suggests an advantage, unrelated to pH optima, for the general fall in intracellular pH with rising temperature: were it not to happen, shifts of base equivalents amongst compartments ({Delta}N) could often be excessive. Modelling cells of a hypothetical fish as a single compartment, suppose that {Delta}pH/{Delta}T were actually zero. Suppose also that {Delta}logPCO2/{Delta}T were 0-0.03°C-1 (encompassing the range of actual averages given above), that intracellular ß* and ßbic were, say, 50 and 4 mequiv kg-1 water pH unit-1, respectively, and that the value of {Delta}pK*/{Delta}T were, say, -0.015°C-1. Then, from Equation 10, {Delta}N/{Delta}T would lie between -0.73 and -0.85 mequiv kg-1 water °C-1, implying a movement of base equivalents into the cells. (The most influential term is ß*{Delta}pK*/{Delta}T.) Suppose now that this shift applied to all cells in the body, that the ratio of extracellular to intracellular water were 0.4 as in the channel catfish Ictalurus punctatus (Cameron, 1980Go) and that the extracellular [HCO3-] were initially, say, 7 mmol kg-1 water. A temperature rise of 4°C would then suffice to deplete the extracellular fluid of all its HCO3-. This is so unrealistic that the postulated constancy of cell pH has to be wrong and the mean intracellular value of {Delta}pH/{Delta}T must be negative. As discussed below, the lesser reduction in plasma [HCO3-] that occurs with increasing temperature in many fish is a separate issue. The argument applies less forcefully to air-breathing vertebrates, in which extracellular [HCO3-] is higher, but, as already noted, values of {Delta}pH/{Delta}T within cells are typically negative even in these animals. As to individual tissues, these may gain or lose base equivalents on warming (e.g. Reeves, 1972Go; Heisler and Neumann, 1980Go; Heisler, 1986cGo; Stinner and Hartzler, 2000Go).

Equation 10 may be rearranged as:

(16)
This includes terms for the active regulation of {Delta}pH/{Delta}T through adjustments in PCO2 and in N, but it is the passive mechanism of non-bicarbonate buffering that acts more promptly. This is represented by the term (ß*bic)({Delta}pK*/{Delta}T-{Delta}pH/{Delta}T), which may be positive or negative. In the special case that {Delta}pH/{Delta}T equals {Delta}pK*/{Delta}T in the steady state, non-bicarbonate buffering makes no ultimate contribution to homeostasis, regardless of how much buffer is present. This contrasts with the determination of {Delta}pH/{Delta}T by {Delta}pK*/{Delta}T in CO2-free buffer solutions.

In extracellular fluids, including cerebrospinal fluid, the buffer value, ß*, is often much less than ßbic (see Introduction). Then, regardless of the value of {Delta}pK*/{Delta}T-{Delta}pH/{Delta}T, the temperature dependence of pH is governed mainly by the values of {Delta}logPCO2/{Delta}T and ({Delta}N/{Delta}T)/ßbic. Further modelling of the extracellular compartment in isolation is unrewarding, especially when compared with previous treatments of true plasma (e.g. Reeves, 1972Go; Rodeau and Malan, 1979Go).

Applying the one-compartment model to fish cells
In its ideal form, the alphastat scheme has proved particularly inappropriate to the water-breathing fish that have been studied, for temperature changes in these produce substantial shifts of base equivalents into and out of both cells and body, as well as changes in plasma [HCO3-]. Here we start by modelling all the cells of the body as a single compartment. As in the source papers, it is assumed that {Delta}N/{Delta}T relates only to shifts of base equivalents across cell membranes, and not to metabolic changes within cells. The first two examples illustrate how old data may be usefully approached in new ways.

The adult dogfish, S. stellaris, is notable for its high value of {Delta}logPCO2/{Delta}T in arterial blood, i.e. 0.029°C-1 on average as compared with 0.0017°C-1 in juveniles (Heisler et al., 1980Go). Warming leads to net loss of base equivalents from the cells collectively, and also from the whole body. Unusually, extracellular [HCO3-] rises (Heisler et al., 1980Go; Heisler, 1984Go). In modelling the cells, we may take the following representative values: {Delta}logPCO2/{Delta}T, 0.029°C-1; ß*, 45 mequiv kg-1 cell water pH unit-1; {Delta}N/{Delta}T, 0.105 mequiv kg-1 cell water °C-1 (Heisler and Neumann, 1980Go; Heisler, 1984Go). Based on an intracellular [HCO3-] of "about 1 mmol l-1" (Heisler and Neumann, 1980Go), ßbic is 2.3 mequiv kg-1 cell water pH unit-1. The value of {Delta}pK*/{Delta}T-{Delta}pH/{Delta}T may now be estimated from Equation 10. If the value of {Delta}pK*/{Delta}T is taken to be between -0.011 and -0.018°C-1 (see Introduction), that of {Delta}pK*/{Delta}T-{Delta}pH/{Delta}T must be 0.0025-0.0028°C-1. For white muscle, which makes up most of the fish, the corresponding difference is also positive, approximately 0.001°C-1, since {Delta}pK*/{Delta}T is approx. -0.017°C-1 (see Introduction) and {Delta}pH/{Delta}T is approx. -0.018°C-1 (Heisler et al., 1980Go). The estimated difference of 0.0025-0.0028°C-1 for all cells collectively, multiplied by ß*, implies that non-bicarbonate buffering generates base equivalents at 0.11-0.13 mequiv kg-1 cell water °C-1. The value chosen for whole-body ß* is based somewhat arbitrarily on measurements on white, red and cardiac muscle, all possibly raised artefactually by inorganic phosphate (see Introduction). If a lower value is used, say 30 mequiv kg-1 cell water pH unit-1, the estimate of {Delta}pK*/{Delta}T-{Delta}pH/{Delta}T becomes 0.0037-0.0042°C-1 and the quantities of base equivalents generated on warming are 2.4% lower.

The mean value of {Delta}pK*/{Delta}T-{Delta}pH/{Delta}T may also be estimated for the intracellular compartment of I. punctatus. Data of Cameron and Kormanik (1982Go) suggest the following representative values. For the whole body: {Delta}logPCO2/{Delta}T, 0.0164°C-1; {Sigma}{V{Delta}N/{Delta}T}, 0.056 mequiv kg-1 cell water °C-1. For all cells: V, 0.726; {Delta}pH/{Delta}T, -0.0148°C-1; {Delta}[HCO3-]/{Delta}T, -0.028 mmol kg-1 water °C-1; ß*, 35 mequiv kg-1 water pH unit-1. For the extracellular fluid: V, 0.274; {Delta}pH/{Delta}T, -0.0141°C-1; {Delta}[HCO3-]/{Delta}T, -0.097 mmol kg-1 water °C-1; ß*, say 4 mequiv kg-1 water pH unit-1 (chosen as slightly below the value of approximately 5.8 for blood with zero haematocrit). The extracellular value of {Delta}pK*/{Delta}T is unknown, but not very critical. If it is taken as, say, -0.013 to -0.019°C-1 (see Introduction), then, from Equation 13, the intracellular value of {Delta}pK*/{Delta}T-{Delta}pH/{Delta}T is 0.0003-0.0006°C-1, again positive. This seems small enough to suggest item (4) of the alphastat scheme and implies very little generation of HCO3- by buffering. From Equation 5, the value of {Delta}N/{Delta}T for the cells is 0.039-0.048 mequiv kg-1 water °C-1.

With increasing temperature, there is a net loss of base equivalents from the cells of both these species, and a rise in PCO2. These effects can be seen as alternative ways of lowering cell pH. According to Equation 10, for constant intracellular values of ßbic, ß*, {Delta}pK*/{Delta}T and {Delta}pH/{Delta}T, reduction in the value of {Delta}N/{Delta}T from x to zero in a model fish would require that the value of {Delta}logPCO2/{Delta}T be raised by xbic. In S. stellaris the value of {Delta}logPCO2/{Delta}T would thus need to be approximately 0.075°C-1. Such a high, perhaps unattainable, value does not explain why the shifts in base equivalents occur in the real fish, since, with only minor changes in {Delta}pK*/{Delta}T or {Delta}pH/{Delta}T, {Delta}N/{Delta}T could be zero even at constant PCO2.

The air-breathing swamp eel Synbranchus marmoratus contrasts with these two species in that, collectively, the cells take up base equivalents on warming, i.e. about 0.25 mequiv kg-1 cell water °C-1 (Heisler, 1984Go). Here, therefore, the value of {Delta}pK*/{Delta}T must be more negative than that of {Delta}pH/{Delta}T. Indeed, values of {Delta}pH/{Delta}T are only -0.009 and -0.003°C-1, respectively, in white skeletal muscle and heart.

White skeletal muscle of Z. viviparus is of interest for its high value of {Delta}N/{Delta}T (van Dijk et al., 1997Go). From the mean value of {Delta}[HCO3-]/{Delta}T, i.e. -0.27 mmol kg-1 cell water °C-1, and from estimates, already noted, of {Delta}pH/{Delta}T in vivo and of {Delta}pK*/{Delta}T and ß*, the value of {Delta}N/{Delta}T is calculated from Equation 5 as 0.58 mequiv kg-1 cell water °C-1. If all the cells were like this, large temperature changes would have major implications for extracellular homeostasis.

Applying the model to tetrapods: the protein titration hypothesis of Stinner et al. (1998Go)
In detailed studies of C. constrictor, Stinner and Wardle (1988Go) and Stinner et al. (1998Go) found an increase in whole-body CO2 stores with cooling, and with it increases in both CCO2 and pH in arterial plasma and skeletal muscle. Little evidence was found for changes in either lactate or the balance of inorganic anions and cations that would suggest shifts of base equivalents. It was concluded that changes in whole-body CO2 stores result from changes in protein ionization coupled with ventilatory regulation of PCO2, such that the overall value of {Delta}pKprot/{Delta}T is more negative than that of {Delta}pH/{Delta}T. Thus there is titration of proteins by carbonic acid (along with other non-bicarbonate buffers), rather than a maintenance of their overall ionization state as in item (4) of the hypothesis.

Stinner et al. (1998Go) extended this idea to other reptiles and amphibians. Whole-body CO2 stores increased with cooling in all 13 species studied (Stinner and Wardle, 1988Go; Stinner et al., 1994Go, 1998Go). The changes took many hours. Mean values of {Delta}CCO2/{Delta}T ranged from -0.02 mmol kg-1 body mass °C-1 in R. catesbeiana, to -0.21 mmol kg-1 body mass °C-1 in the tortoise, Testudo graeca. Presumably the range would be even greater if expressed in terms of body water. Only in the bullfrogs do the results seem close to the alphastat prediction of constant tissue CO2 content.

The further analysis by Stinner et al. (1998Go) may be described in terms of the one-compartment model, in which Equation 7 shows the determinants of {Delta}CCO2/{Delta}T. The term {Delta}N/{Delta}T is regarded as negligible on the basis of the findings for C. constrictor. As already noted, the term {Delta}[CO2]/{Delta}T is also trivial here. Thus Equation 7 reduces to:

(17)
Stinner et al. (1998Go) took the whole-body value of {Delta}pH/{Delta}T as approximating that for arterial plasma and found a linear relationship between that and whole-body {Delta}CCO2/{Delta}T (10 species; r=-0.93). The values of {Delta}pH/{Delta}T are mostly taken from other studies over similar ranges of temperature (Howell et al., 1970Go; Jackson et al., 1974Go; Malan et al., 1976Go; Bickler, 1981Go; Wood et al., 1981Go; Nicol et al., 1983Go). The equation of the regression line is:

(18)
In accordance with Equation 17, this suggested for the whole body a mean non-bicarbonate buffer value, ß*, of 8.24 mequiv kg-1 body mass pH unit-1 and a mean value for {Delta}pK*/{Delta}T of -0.022°C-1.

As Stinner et al. (1998Go) pointed out, this value of {Delta}pK*/{Delta}T is reasonable for some small imidazole compounds. However, the real value is probably no more negative than -0.018°C-1 (see Introduction). As for that whole-body value of ß*, it may be re-expressed in terms of body water using a representative body water content of, say, 76% (Deyrup, 1964Go; Bentley, 1976Go); it then becomes 10.8 mequiv kg-1 water pH unit-1. This is little above the 8.1 mequiv l-1 pH unit-1 calculated for plasma of C. constrictor, despite the greater contribution of the cells, where ß* is presumably much higher (see Introduction). It therefore seems improbably low. Next, the assumption that the whole-body value of {Delta}pH/{Delta}T approximates that for arterial plasma may be inappropriate, since values of {Delta}pH/{Delta}T in C. constrictor averaged -0.009°C-1 in muscle and -0.0028°C-1 in arterial plasma. (Modelling of the sort to be described next, but starting with Equations 13 and 18, also shows the assumption to be implausible.)

The data may be better modelled by treating the body water as two compartments, intracellular and extracellular, and taking account of data on PCO2. Values of {Delta}logPCO2/{Delta}T are assumed to be the same in both compartments, both for simplicity and because the average differences in PCO2 between blood and cells in these air breathers are likely to be small (Burton, 2001Go). Again {Sigma}{V{Delta}[HCO3-]/{Delta}T} is taken as approximating {Delta}CCO2/{Delta}T. Equation 18 is assumed to apply exactly. Plausible values, representing all species collectively, are allotted to other parameters. The water content of the body is again taken as 76%.

Fig. 1 shows the correlations between {Delta}pH/{Delta}T and {Delta}logPCO2/{Delta}T in arterial plasma and between {Delta}CCO2/{Delta}T and {Delta}logPCO2/{Delta}T for the whole body. The nine species fall into two groups and the mean values for each group of {Delta}pH/{Delta}T and {Delta}logPCO2/{Delta}T are shown by crosses marked A and B. For group A they are, respectively, -0.015°C-1 and 0.018°C-1. For group B they are, respectively, -0.004°C-1 and 0.008°C-1. Values of {Delta}CCO2/{Delta}T corresponding to groups A and B, calculated from Equation 18, are -0.058 and -0.148 mmol kg-1 body mass °C-1, or -0.076 and -0.195 mmol kg-1 water °C-1, respectively.



View larger version (12K):
[in this window]
[in a new window]
 
Fig. 1. The temperature dependence of arterial pH and PCO2 and of whole-body CO2 stores (CCO2) in tetrapods. Mean values of {Delta}pH/{Delta}T and of {Delta}CCO2/{Delta}T (mmol kg-1 body mass °C-1) are as tabulated by Stinner et al. (1998Go). Mean values of {Delta}logPCO2/{Delta}T, in °C-1, are from the following sources (where there is more than one per species, the means and ranges of individual means are shown): filled circles, Coluber constrictor (Stinner et al., 1998Go), Varanus exanthematicus (Wood et al., 1981Go); open circles, Trachemys scripta (Jackson et al., 1974Go; Jackson and Kagan, 1976; Hitzig, 1982Go), Alligator mississippiensis (Davies, 1978Go), Chrysemys picta (Nicol et al., 1983Go; Glass et al., 1985Go), Chelydra serpentina (Howell et al., 1970Go), Dipsosaurus dorsalis (Bickler, 1981Go); filled triangles, Bufo marinus (Howell et al., 1970Go), Rana catesbeiana (Howell et al., 1970Go; Reeves, 1972Go; Mackenzie and Jackson, 1976). As described in the text, the crosses A and B show representative values for the two distinct groups of species (i.e. filled circles, and open circles together with filled triangles).

 

For the extracellular and intracellular fluids, respectively, the values of V are taken as 0.4 and 0.6 and the bicarbonate buffer values, ßbic, are taken as 60 and 24 mequiv kg-1 water pH unit-1. From these parameters and the data of the previous paragraph, the value of {Delta}pH/{Delta}T for the cells is calculated using Equation 14. For both sets of data it is -0.0141°C-1. The important point here is not its exact value, which depends on the chosen parameters, but the fact that the mean values for the two groups of species are plausibly modelled as similar. This seems a reasonable postulate (despite differences within groups) if optimum cell function depends on the relationship between intracellular pH and temperature.

The parameters ß* and {Delta}pK*/{Delta}T for extracellular fluid are now allotted plausible values, say 5 mequiv kg-1 water pH unit-1 and -0.018°C-1, respectively (see Introduction). Then the extracellular values of {Delta}N/{Delta}T, calculated from Equation 10, are 0.123 mequiv kg-1 water °C-1 for group A and 0.008 mequiv kg-1 water °C-1 for group B. Warming therefore leads to a loss of base equivalents from the extracellular fluid. These calculations may be repeated for the whole body using Equation 13, with values for ß* and {Delta}pK*/{Delta}T in the intracellular fluid taken, say, as 25 mequiv kg-1 water pH unit-1 and -0.0130°C-1, respectively. Then the whole-body value of {Delta}N/{Delta}T, i.e. {Sigma}{V{Delta}N/{Delta}T}, is 0.087 mequiv kg-1 water °C-1 for group A and 0.183 mequiv kg-1 water °C-1 for group B. These are almost equal, but opposite in sign, to the respective values of {Delta}CCO2/{Delta}T given above. (That this is about equally true of the two group means was arbitrarily achieved by adjusting the value of {Delta}pK*/{Delta}T.) For groups A and B, the sums {Sigma}{V{Delta}N/{Delta}T}+{Delta}CCO2/{Delta}T are, respectively, +0.011 and -0.012 mequiv kg-1 water °C-1. These small differences correspond to the titration of proteins and other buffers (Equation 13).

The model compares and integrates data from two groups of species, but the diagnoal lines in Fig. 1 can also represent the changes in a single hypothetical individual as the value of {Sigma}{VN} alters after temperature changes. After warming, there is a loss of gaseous CO2 from its body and, as modelled, this loss is nearly matched by a loss of base equivalents. These come partly from the cells, with the CO2 generated from HCO3- and H+ ions (and almost entirely so in the case of group B). The reduction in [HCO3-] in the cells is matched by a fall in PCO2 that keeps the value of {Delta}pH/{Delta}T constant. (For the PCO2 to fall even as CO2 is generated from HCO3-, ventilatory adjustments to PCO2 would have to be rapid; the many hours needed to achieve a steady state would thus reflect slow changes in N rather than slow gas exchange.) If the value of {Sigma}{V{Delta}N/{Delta}T} were zero, {Delta}CCO2/{Delta}T would be positive instead of negative, i.e. 0.036 mmol kg-1 water °C-1 (calculated from Equations 13 and 18).

No one set of parameters can be right for all species, and each species or individual should ideally be modelled with its own set. Moreover, data for real cells are generally for particular muscle tissue rather than for the whole intracellular compartment. The chosen parameters are broadly in line with data given in the Introduction, but the constant intracellular value of {Delta}pH/{Delta}T in the model (-0.0141°C-1) is more negative than the values of -0.009, -0.012 and -0.007°C-1 measured in skeletal muscle of C. constrictor, R. catesbeiana and B. marinus, respectively (Stinner et al., 1998Go; Stinner and Hartzler, 2000Go). It is closer to the mean whole-body intracellular value (-0.0151°C-1) obtained by Bickler (1982Go) in the lizard Dipsosaurus dorsalis, itself more negative than his values for skeletal and cardiac muscle (-0.0098 and -0.0104°C-1, respectively).

According to the model, {Delta}[HCO3-]/{Delta}T for the extracellular fluid is negative, having values of -0.14 and -0.08 mmol l-1 °C-1, respectively, for groups A and B (calculated from Equations 5 or 9). Some values determined for real arterial plasma in these tetrapods are similar in sign and magnitude (Wood et al., 1981Go; Stinner and Wardle, 1988Go; Stinner et al., 1998Go), but others do not differ significantly from zero (Jackson et al., 1974Go; Bickler, 1981Go; Nicol et al., 1983Go; Stinner et al., 1994Go). Shifts of base equivalents ({Delta}N/{Delta}T) between compartments are hard to quantify experimentally, because accurate analyses are needed for all ions present. The shifts seem insignificant in D. dorsalis (Bickler, 1984Go) and, although they do occur in B. marinus and R. catesbeiana, a consistent trend is not evident (Stinner and Hartzler, 2000Go). Neither these discrepancies and uncertainties, nor the arbitrariness of some model parameters, invalidate the semi-quantitative conclusions summarized next.

Two important conclusions have emerged. Firstly, it is shown that mean intracellular values of {Delta}pH/{Delta}T could be similar in the two groups of species. Secondly, the net loss of gaseous CO2 following a rise in temperature could be due largely to titration of HCO3- as base equivalents are lost from cells and body (or proton equivalents gained). The latter idea is absent from the model of Stinner et al. (1998Go), but was originally suggested by Stinner (1982Go) for the snake Pituophis melanoleucus. Bickler (1984Go) did not find evidence for a major role of excretion in the acid—base responses of D. dorsalis to temperature. Moreover, in none of the species can the loss of base equivalents be due mainly to excretion of HCO3- since the accompanying reduction in whole-body CO2 stores is measured as gaseous CO2. It is therefore more likely that the whole-body gains and losses of base equivalents involve metabolic adjustments to intracellular concentrations of organic ions (see below). Because {Delta}pH/{Delta}T does not alter, these would not be homeostatic for pH. A major temperature-dependent process modifying N, and best developed in species where CCO2 changes most, should now be sought. Although the model is made consistent with the relationship of Equation 18, that remains unexplained.

Discussion

The equations provide a convenient approach to the effects of temperature on pH, PCO2, [HCO3-] and buffer ionization. They can be used both to explore the effects of varying buffer properties etc. in hypothetical animals and to complement previous analyses of experimental data. The tetrapod model unites various facts, ideas and uncertainties in what may be less a true description than a step towards better understanding. Indeed there are yet more uncertainties involved and some of these are discussed under the next heading. Also discussed below are other parts of the imidazole alphastat scheme (items 5, 11 and 12) that do not relate directly to the model, and three phenomena that are excluded from the strict alphastat scheme. These are metabolic adjustments to non-buffer ions within cells, and movements of base equivalents both amongst cells of different types and between body and environment.

Some limitations of the model in interpreting measurements on real animals
In the model `total CO2' consists just of the dissolved gas and HCO3- ions, and their reactions are treated in terms of a straightforward apparent equilibrium constant, pK1'. However, this has been found in many studies to decrease with increasing pH (e.g. Dill et al., 1937Go; Boutilier et al., 1985Go; Heisler, 1986aGo). This effect can markedly influence estimates of intracellular [HCO3-] when this is calculated from PCO2 and pH (Reeves, 1976aGo). Unfortunately, the influence of pH on pK1' varies from study to study and a relationship quantified for arterial plasma may be wrong for cells, especially when some of the cell water is `bound' (Garner and Burg, 1994Go). Uncertainties regarding {Delta}pK1'/{Delta}T and {Delta}[HCO3-]/{Delta}T are much less. This effect of pH is not fully understood. `CO2' as measured gasometrically exists not only as free HCO3- and dissolved gas (plus minute amounts of carbonic acid), but as CO32- (generally in small amounts), as carbamate (barely studied outside of erythrocytes), possibly as the compound H2CO3.HCO3- (Covington et al., 1981Go) and as ion pairs of HCO3- with cations such as Na+, Mg2+ and Ca2+. Boutilier et al. (1985Go) and Burton (1987Go) discuss these and other uncertainties in calculating `[HCO3-]' from pH and PCO2. According to formulae given by Heisler (1986aGo), values of {Delta}pK1'/{Delta}T-{Delta}logS/{Delta}T for solutions resembling protein-free plasma are about 0.0053°C-1 for 0-25°C, as above, and about 0.0069°C-1 for 25-35°C. As discussed in the next paragraph, the intracellular PCO2 cannot be assumed to be exactly that of accessible extracellular fluids.

The single-compartment model is homogeneous, unlike both real extracellular space and real cells. Regarding PCO2, this is generally higher in venous than arterial plasma and higher still in interstitial fluid (Pörtner and Sartoris, 1999Go), and cells vary in their relationships to blood vessels. As modelled by Burton (2001Go), the discrepancy between arterial and mean whole-body interstitial or cellular PCO2 varies inversely with arterial PCO2 and is therefore greatest in water-breathing fish. How far the discrepancy varies with temperature is unclear, because it depends also on respiratory quotient, oxygen tensions, the relative solubilities of the two gases, and the possible disequilibrium of CO2 in blood. As modelled for the whole body, average interstitial and cellular PCO2 in some fish can be more than twice the arterial PCO2. For real cells CCO2 has sometimes been calculated from cell pH and arterial PCO2 by the method of Cameron (1980Go); for fish especially, the results could be much too low. Equation 14 is based on the assumption that {Delta}logPCO2/{Delta}T, but not necessarily PCO2, is the same in all compartments.

Cytoplasm is heterogeneous too. Much of a cell may be taken up with acidic organelles or the very alkaline mitochondrial matrix. In addition, local variations in net fixed charge density on proteins and membranes must cause inhomogeneities of pH and [HCO3-]. Estimates of cell pH made using DMO (5,5-dimethyloxazolidine-2,4-dione) yield values that approximate to averages for the whole cell contents, but, more exactly, what is averaged is 10-pH (Waddell and Bates, 1969Go). There is little quantitative information on mitochondrial pH in vivo and on its temperature sensitivity in ectotherms. However, Moyes et al. (1988Go) have studied mitochondria isolated from red muscle of the carp, Cyprinus carpio: provided extramitochondrial pH varied as in vivo, the transmembrane pH gradient remained constant. If this gradient is generally insensitive to temperature in ectothermic vertebrates, then values of {Delta}pH/{Delta}T in cells obtained with DMO should reflect cytosolic values. Cell [HCO3-] may be calculated from pH and PCO2. With pH values obtained by the DMO method, the resulting [HCO3-] averaged over all subcompartments, each with its fractional volume V and concentration [HCO3-], equals 1/{{Sigma}(V/[HCO3-])}, where {Sigma}V=1. Given, for example, two subcompartments of equal volume differing in pH by 0.3, the true mean [HCO3-] is 12% higher than that calculated from the pH measured by DMO. Pörtner and Sartoris (1999Go) give a detailed analysis of the effects of cytoplasmic heterogeneity on pH measurements and calculations of [HCO3-].

The model compartments are of constant vol