Cellular information dynamics through transmembrane flow of ions.

We propose cells generate large transmembrane ion gradients to form information circuits that detect, process, and respond to environmental perturbations or signals. In this model, the specialized gates of transmembrane ion channels function as information detectors that communicate to the cell through rapid and (usually) local pulses of ions. Information in the ion "puffs" is received and processed by the cell through resulting changes in charge density and/or mobile cation (and/or anion) concentrations alter the localization and function of peripheral membrane proteins. The subsequent changes in protein binding to the membrane or activation of K+, Ca2+ or Mg2+-dependent enzymes then constitute a cellular response to the perturbation. To test this hypothesis we analyzed ion-based signal transmission as a communication channel operating with coded inputs and decoded outputs. By minimizing the Kullback-Leibler cross entropy [Formula: see text] between concentrations of the ion species inside [Formula: see text] and outside [Formula: see text] the cell membrane, we find signal transmission through transmembrane ion flow forms an optimal Shannon information channel that minimizes information loss and maximizes transmission speed. We demonstrate the ion dynamics in neuronal action potentials described by Hodgkin and Huxley (including the equations themselves) represent a special case of these general information principles.

Introduction

Investigation of information dynamics in living systems typically focus on the role of the genome in transgenerational flow of information and, via epigenetics, variations in gene transcription in response to internal or external signals. Here we propose an additional mechanism of information dynamics that is based on ion fluxes at the cell membrane. As outlined below, we posit these dynamics are localized around the cell membrane and are independent of additional input from the genome, thus allowing both rapid and local responses to environmental perturbations or signals.

In brief, our hypothesis begins with the observation that virtually all cells generate transmembrane ion gradients produced by ATP-dependent membrane pumps. These gradients form high information[1,2] states (e.g., large Shannon entropy) and possess high potential energy (similar to a coiled spring) initial conditions (Fig. 1A). External perturbations are detected by specific gates in the transmembrane ion channels. Information triggered by the external cue is then communicated to the cytoplasm through a local flow of ions that is controlled by opening the gates in these channels (Fig. 1B). The information from the “pulse” of ions (similar to Ca2+ puffs[3]) around one or more membrane channels is then translated into signaling circuits by changing local charge density and/or ion concentration. We note that the ion could flow into the cells creating the equivalent of the Ca2+ puffs or out of the cell generating a “negative puff” caused by absent ions (Fig. 1B).

Figure 1

Summary of ion-based information dynamics around the cell membrane (A). Baseline state. The initial conditions of the cell membrane produced through energy-dependent ion pumps. Concentrations of Na+ and Cl− are much higher in the extracellular space (left) than in the intracellular space (right) and vice versa for K+. Eukaryotic cells can use up to 40% of their energy budget to produce ion gradients. Specific transmembrane ion channels are present with their gates closed, thus allowing limited ion flow along concentration gradients. Peripheral membrane proteins included structural elements that support the membrane (proteins A) and enzymes (proteins B). Negative charges on enzymes (proteins B) are shielded by K+ ions causing the enzyme to be inactive. At the inner leaflet of the cell membrane primarily K+ ions shield proteins from the negative chargers of the membrane. As a result the positively charged protein A cannot bind. Note that K+ ions are about twice the  volume of Na+ ions. In our model, this subtle difference is posited to be sufficient to alter the shape of the active site and, therefore the activity of some enzymes (as shown in enzyme B). (B) Ion flux due to perturbation. An external cue causes the gate in the K+ channel to open such that K+ ions flow along concentration gradients out of the cell. This results in decreased shielding of the negative charges in the inner membrane leaflet so that positive charges on protein A can form electrostatic bonds. The binding of the structural protein to the lipid membrane stabilizes it as a response to the perturbation. Further, loss of K+ ions causes the enzyme (protein B) to change in both position and shape so that it binds to its substrate to catalyze the reaction S → P1 + P2. Note that biological enzymes can increase the reaction rate by as much as 15 orders of magnitude allowing the response to be both local and rapid. (C) Ion-Flux Induces response to perturbation. Activation of the enzyme B results in rapid increase in the local concentrations of the reaction products P1 and P2. In this model, P2 blocks the actions of the perturbation closing the gate on the K channel. The other product, P1, opens the Na channel allowing rapid flow of cations, which restore the shielding on the inner leaf of the cell membrane, displacing the A proteins from the inner leaflet of the cell membrane back into the cytoplasm.

The change in ions shielding negative charges on the inner leaf of the cell membrane can alter the binding of local peripheral membrane proteins[4] thus altering the membrane structure and function (Fig. 1C). Similarly, opening of the potassium channel receptor by ligand binding to the Epidermal Growth Factor Receptor[5,6] results in loss of K+ shielding of negative charges on the inner leaf of the cell membrane allowing “recruitment”[7,8] of positively charged RAF proteins to the GTP-bound RAS proteins (which are also negatively charged).

An additional cellular response is elicited because the functions of many enzymes[9] is dependent on local cation concentrations (e.g. Na+, Mg2+, Ca2+,or K+). The subsequent change in localization, activation or repression of local enzymes can produce an additional response to the local perturbation (Fig. 1C). When a successful response to the perturbation is complete, the ion gates close and the initial conditions (i.e. large transmembrane ion gradients) are restored to respond to any subsequent environmental changes.

These ion dynamics are similar to those observed in propagating action potentials of neurons[10]. However, we note that the steep transmembrane ion gradients that permit rapid flow of ions through the transmembrane channels of neurons are widely present in eukaryotes. Hence, we hypothesize that the transmembrane ion flows observed and modeled by Hodgkin and Huxley represent a special case of a broader phenomenon in which cells use transmembrane ion gradients to access, process and respond to information cues in the environment.

Here we assessed the potential of pulses of ions through cell membrane channels as a mode of information transmission, and define a general mathematical expression of these dynamics. We demonstrate that communication through ion fluxes represents an optimal mechanism for information transmission that both minimizes information loss and maximizes the speed of information transmission. Furthermore, our analyses support the notion that the ion dynamics observed in the Hodgkin-Huxley experiments and, indeed, the Hodgkin-Huxley equations themselves, represent a special case of these broader principles that are valid for a wide range of eukaryotic cells.

Methods

Changing ion concentrations around the cell membrane

Consider a problem of determining, over a time interval , the generally varying change in the intracellular ion concentrations of N types of ions that result from flow of ions through channels in the cell membrane. Examples might be the ions Na+, K+, Cl−, Ca+2, etc. Thus, are the probability density functions over the time interval and so can also represent a rate or flow. We also use this terminology. Thus each type i of ion obeys normalizationAlso, definite boundary values are taken to hold. What are the unknown fluctuations of ion probability density functions ?

Transmembrane ion flows as information

Constant extracellular charge density

Here we assume that, over the time domains in which a model functions, the probability densities of the ions outside of the membrane are fixed. That is, ATP-dependent membrane ions pumps will use energy to move ions into or out of the cell. This will have the effect of changing the intracellular ion concentrations but not those of extracellular ions , which are approximately constantin comparison to the generally varying probability densities of ions within the cell. Note that this constancy, or uniformity, is mathematically equivalent to the assumption of an ideal enzyme of uniform density, as is often used to characterize the Arrhenius law[11].

Constancy (2) arises as follows:

In the case of a single-cell organism the extracellular space is effectively the ocean, wherein ion concentrations do not change in biologically relevant time. Or

For a multicellular organism, the extracellular space is in equilibrium with the interstitial fluid of the entire organisms (including, for example, the blood). Therefore the ion concentrations in the stable physiological state of a healthy individual should likewise tend to be constant. While there may be a small local changes in concentration exterior to the cell due to ion flow across the membrane, this will be quickly dissipated such that the extracellular ion space remains constant, serving as a source or sink for intracellular ions.

With the assumption that we can therefore assume that changes in intracellular ion concentrations are solely dependent on functions within the cell membrane. That is, membrane ion pumps generate transmembrane gradients by varying the intracellular ion concentrations relative to fixed extracellular concentrations. Further, when the channel gates open, the rapid flow of ions (measured to be about 105 ions/second) of ions allows . We propose this ion-based communication channel represents a method for information transfer that both minimizes loss and maximizes speed..

Criterion of minimum information loss

We first seek a general solution for intracellular ion concentration changes that is consistent with the constant concentrations outside the cell. We first useful to consider how informative flow rates occur.

We start with the assumption that ATP-dependent membrane pumps transfer ions across the cell membrane. Because of these ion-specific pumps, the intracellular concentrations of Na+ and Cl− are much lower than those in the extracellular space, and vice versa for the concentrations of K+ (Fig. 1). The non-equilibrium condition is manifested by an often large transmembrane potential V. Notably, the asymmetric distribution of ions produced by these pumps (compared to random distribution in the absence of those pumps) encodes on the order of 1014 bits of Shannon[12] information[1,2,13,14].

Suppose the flow of some number of ions from outside to inside the cell membrane can signal the occurrence of an environmental event. For this to occur, the event must trigger the gates in one or more of the ion-specific membrane channels. These gates permit ion flow along the pre-existing concentration gradients between the inside and outside of the cell. Several hundred types of channel gates are encoded in the human genome, permitting responses to a potentially large range of external perturbations. The role of ion flow in cell functions is most well-studied in neurons. To propagate nerve impulses, the channel gates along the axis of the neuron respond to voltage fluctuations cause by the local arrival of the action potential. When the gates open, the transmembrane gradient permits rapid ion flow along steep concentration gradients producing the action potential described by the Hodgkin Huxley equations.

We assume that evolution selects for optimal means for addressing risk and capitalizing on opportunities in the environment. We propose that this optimization principle selects for information dynamics at the membrane, which conveys maximum Shannon information regarding the occurrence times of these events. Of course the system is passive so that, in general, there can only be a loss of such information relative to the constant ion density probability functions outside. Accordingly, the maximal Shannon information is conveyed when the resulting intracellular ion probability density is minimally different the of Eq. (2). As noted below, this also defines how fast the information from an environmental source can be decoded by the cell surface.

Acceleration of transmembrane ion flows through gated channels

A reasonable assumption is that an environmental change is sensed at the cell membrane. This could represent a change in a physical interaction with an adjacent cell or other surface or a change in the local concentration of some substance (e.g., substrate molecules in the case of a free swimming eukaryote or, in the case of neurons, a neurotransmitter). How is this information communicated into the cell? We propose that eukaryotic cells have efficient mechanisms for sensing and reacting to environmental change through gated ion channels on the cell membrane. ATP-dependent membrane ion pumps typically generate large transmembrane gradients but the flow of ions along such gradients is ordinarily very slow (“leakage current”) due to high resistance in nonpolar regions of the cell membrane. Opening of channel gates by an external perturbation has the effect of greatly expanding the transmembrane flow of the specific ion that passes through that channel. Indeed, it has been noted that the actions of the channel in the membrane are analogous to a catalyst[15] for diffusion, in the sense that the acceleration of the ion flow are qualitatively and quantitatively similar to enzyme-induced acceleration of a chemical reaction rate.

What can we learn by considering such an analogy? As will be shown, it leads to derivation of the Hodgkin-Huxley equations that govern the propagation of ions across the neural cell membrane. Moreover, the derivation is sufficiently general to apply to ion flow through gates in a wide class of cell membranes.

In general the passage of a given type of ion, say K+, is restricted to a specific channel just as each of the key ions (e.g., Na+, K+, Cl−, Ca2+, Mg2+) can pass through one and only one type of channel[16]. Each such ion then represents a distinct component of the transmembrane information flow. The channels contain gates that are ordinarily closed, thus preventing ion flow. However, this still allows a “leakage” current, which is the usually small flow of ions through the membrane. Once a transmembrane ion gradient is established through ATP-dependent membrane ion pumps, maintenance of that gradient requires only that the pumps replace the small loss of ions due to leakage current. A detailed description of the gate channels can be found in[16].

Finally, the passage of a given type of ion, say K+, requires that more than one gate in the channel opens. For example, the K+ channel allows ion flow only when all 3 channel gates are open. These appear to form a ring shaped gap through which the required K+ ion passes[17,18].

Principle of minimum Kullback-Leibler divergence

The cell membrane provides a ‘window’ to the world outside for each cell. Some signals that arrive at the membrane, such as growth factors, require transmission through well recognized cytoplasmic signaling pathways to the nucleus for further processing and action. However, others, particularly life-threatening perturbations, may require an immediate response that is mediated entirely by the molecular machinery in and around the cell membrane.

The rapid transmembrane flow of ions within one or more channels allows cells to receive, locally process and respond to information in the environment. Thus, ions of type i = 1,.., N may pass through the cell membrane, each potentially signaling vital information to the cell. Hence the changing concentration, and correct identification, of the various ion types is vital to cell survival and proliferation. We shall assume that natural selection, acting through geological time, has optimized the speed and accuracy of identifying i on average. This implies that the intracellular probability functions are formed, coded, and decoded at the membrane according to a well-defined optimization principle. We use minimum average code length as the criterion[19], as follows (in brief, see Eqs (3), (4)).

Information flow optimization by maximally efficient code

Claude Shannon[12,20] posed the following question of optimization: If we are given a probability on event times t for events i, what code lengths , i = 1, .., N achieve their minimum expected values over all times t? This defines a maximally efficient code. Shannon and Fano showed that if symbols (say letters of an alphabet) , i = 1, .., N occur with respective frequencies these are coded with maximum efficiency (minimum average code length) if each symbol i is coded with a length . Thus, the most frequently occurring symbols would be coded with minimum lengths l. Or, minimum code lengths are required, and this is independent of time. This is intuitive for maximally fast response. Thus, the information dynamics at the cell membrane in which the external ion concentrations (expressed as probabilities ) are actually constants , manifest the conditions for optimum coding.

However, the same constancy does not hold in practice for data. Their occurrence rates randomly differ rates from the ideal rates . (Think of the letter rates on a single page of print). Hence, if these are coded according to Shannon and Fano’s scheme, the result is a wasteful excess of code lengths by amountat each time. The expected value of the excess code length over all time for each ion type i is then . This average excess code length translates into excessive total average communication time required for responding to the environmental event being reported on by the ion. In turn, this implies excess reaction time for coping with that event. Excess reaction time to an external threat, for example, can be deadly for the cell. Hence we assume that natural selection minimizes reaction time.

The role of multiple channels

The above models focus on a single specific channel and channel state. This is clearly an oversimplification as it is likely that ion flows into mammalian cells are the result of an ensemble of single channels that superimpose at various times with various weights to form the net channel state[21]. This may vary with time both stochastically and with subtle changes in the environmental condition. These ensemble dynamics leading to a net channel can be accommodated using KL principle because it is a linear differential equation. Thus, the net channel state P (t) is, by hypothesis, the statistically average state. Therefore any linearly weighted superposition of single channel outputs will also satisfy our KL principle so that the net state P (t) obeys KL principle. However, the detailed shape of the curve P (t) will depend upon what the times , n = 1, …, 4 are in the various channel states. Nevertheless, even allowing for a quite complex ensemble of independent channels, the total ion flux gives an output activation wave obeying the KL principle.

Role of natural selection

The average such excess is assumed to be minimized, as an expression of natural selection (see above). The average is evaluated by simply weighting the r.h.s. of (3) with and integrating over all t, as

In words, Eq. (4) states the average amount of wasted code length over all ion occurrences i equals the relative entropy or the ‘Kullback–Leibler divergence[22]. This is effectively a distance between probabilities and , . For example, at any i and t for which probability , there is contribution to in Eq. (4). There is thus zero distance between such pairs. This implies that if probabilities are formed, coded, and decoded at the membrane according to principle (4) the external environmental events will be identified by the cell membrane with minimum error and maximum speed. We propose that these optimization principles, because the allow maximally fast and accurate responses to external perturbations,  are fundamental Darwinian selection forces in the evolution of all cells.

Quite aside from the above requirement by natural evolution, the minimization requirement in Eq. (4) is also suggested by the dual facts that (i) the Kullback–Leibler divergence Eq. (4) has the form of the Shannon information [12], and (ii) the transmission of the ion through the membrane is an irreversible thermodynamic process which, therefore, obeys coarse graining[23-26]. Such processes can only lose information, and so we again arrive at principle (4) which states that optimal function requires this loss of information to be a minimum[27].

We note that the Kullback–Leibler principle Eq. (4) is being used as a principle of thermodynamics. Such use allows microscale effects, for example occurrence rates , to be found without detailed knowledge of that microscale, such as the mathematical dependences of the fields of potential (e.g., inverse square) that are provided by ions. Thermodynamics is thereby seen to provide a basis for the phenomenon of natural selection.

Convex nature of

The minimum in principle (4) is always an extreme value since is a convex functional. This is fortunate, as finding its extreme value is then elementary (see examples below).

In summary to this point, the solution to principle (4) is found under the special condition Eq. (2) of a constant . Also the solution obeys:

(i) minimum departure from that of the constant ion rate (of Eq. (2)) exterior to the cell, in the sense of carrying a level of Shannon information that departs minimally from that of the constant , and

(ii) optimally efficient information transfer of environmental messages (such that they are decoded with maximum speed within the membrane).

Tendency for solution to obey

Note that in Eq. (4) that at all t for which since there is a very small contribution to . Conversely, for all t for which there is large contribution to . Thus, mathematically, is a type of logarithmic measure of the ‘distance’ between corresponding and over all . (Note: isn’t a formal ‘distance’ measure as it isn’t symmetric in exchange of corresponding and , but this is not important for our purposes, as will become apparent).

Use of prior knowledge

The minimization in principle (4) is always found for a given state of prior knowledge, usually in the form of equality constraints on the . These are implemented by the simple weighted addition of the constraints to principle (4), where each is multiplied by a corresponding constant Lagrange multiplier constant. These constants are found by requiring them to obey the known constraints (see below).

Data Availability

No datasets were generated or analyzed during the current study.

Results

Application of the principle to membrane ion channels

The membrane is regarded to generally have some level of transmembrane potential V that can change with time t. Hence principle (4) is applied for some arbitrary values of V . In addition, the minimum in Eq. (4) is attained in the presence of known constraints obeyed by the unknown laws . In general, these obey two, trivial normalization constraints, which are obeyed in all problems plus other, more substantive constraints such as mean passage time through the membrane. These effectively define the particular problem at hand.

We next compared results with known experimental observations, specifically of those of ions that pass through membrane ion channels in the membrane of a neuron – as described by the classical Hodgkin-Huxley equation. Hence, let the law obey substantive constraints as follows. As with Hodgkin and Huxley[10] let the ions pass through the cell membrane over fixed average times ,. These , were empirically observed by Hodgkin and Huxley. It is then seen that the particular values of are substantive in strongly influencing the output rates obeying principle (4).

Hodgkin and Huxley investigated, in particular, the case of N = 3 types of ions, specifically Na+, K+ and Cl−. They observed average path times for these types of ions. Hence, these are used as constraints on the time behavior determined by principle (4). As is usual, these constraints are entered via Lagrange constraint multipliers . Constancy effect Eq. (1) and normalization Eq. (2) are also so used. The resulting constrained principle (4) isWe will see that this principle gives rise to the known dynamical equations of Hodgkin and Huxley for occurrence rates at the assumed potential value V.

Derivation

We note the absence of any derivative terms in the integrands of Eq. (5). Then its Euler-Lagrange solution simplifies to that obtained by simply setting partial derivative of its Lagrangian (total integrand) to zero. This gives a requirementat the arbitrary time t. Solving for ,Boundary values were likewise observed by H and H. These should be obeyed by solution (7). Enforcing these in (7) gives the following solutions.Evaluating this at t = 0 givesSolving,Putting this in (8) then gives

Evaluating Lagrange constants

In principle (5), the mean time constraint and normalization constraintgiving the resultThen result (9a) becomesThese are equivalent to the Hodgkin and Huxley (H and H) equations numbered (8), (17) and (18) in[10]. Thus, the dynamics described by the “Hodgkin and Huxley” equations apply to a wider scope of cell membranes than merely neurons. As a check, we compare Eq. (10) with Eq. (8) of H and H. In their notation this isThus our notation for the ith rate in Eq. (10) becomes simply n for H and H, with their subscript denoting the time. Hodgkin and Huxley take their end time T to be very large, effectively . Thus our notation translates into in their Eq. (8) [Eq. (11) preceding]. Also as in Hodgkin and Huxley[14], a ‘resting state’ is defined by V = 0. Let this occur at time t = 0. Then in Eq. (10) becomes in their Eq. (8) [same as preceding Eq. (11)], although the argument 0 in explicitly represents their case of zero potential.

Some consequences of flow solution

Note that if all time constants in Eq. (10) are very large this gives no environmental information significantly enters the cell. However, such large have their use. As noted above, passage of any ion across the membrane requires that a specific number of gates open through which the ion then passes. For example, passage of a K+ ion through a membrane channel will not occur until all 3 gates open. The latter probability is of approximate size . Alternatively, the passage of an Na+ ion requires that 2 gates in the ion channel open, which then has probability. As the square of a probability tends to be greater than the cube of a comparable one, the Na+ ions are more probable to pass through the channel than are 3 K+ ions. Consequently, a cell tends to efflux Na+ before influxing K+. Further, there is a need for this time gap between the two since, otherwise, the passages of the two ions might so overlap in time such that they collide with one another. Thus, it is important for K+ ions to be delayed in their passage through the membrane. Collectively, the flow solutions both accelerate and slow selected ions to optimize decoding times and avoid collisions.

Ion transition rate constants α, β

Probabilities must be further defined. This involves the vital role played by the membrane potential V.

Ion transitions results from opening of channel gates that are governed by a perturbation in the transmembrane voltage V. However, as in the preceding, at one assumes no perturbation is applied to the membrane so voltage Also, let represent, respectively, the rate of ion transfer at this voltage from outside the cell membrane to inside and from inside to outside. This represents slow “leakage”[10] current with a probabilityNext, we have assumed that a fixed, finite voltage V is applied to the membrane due to arrival of a propagating wave down the axon of the nerve. Then rates instantly take on fixed values (1), appropriate to the value of V over time t that represents the course of the action potential. This allows us to suppress all V notation (as in ref.[10]). Thus in particularEndpoint values (12) and (13) supplement the interior H-H solution Eq. (10). Finally, by the above definition of with V fixed and finite, it must be that the out-to-in ion rate obeys and the in-to-out rate obeys Adding these givesThis shows how the rates determine the rapidity of exponential change of ion rate in Eqs (10).

These rates obviously play a similar role in our solution Eq. (10), supplemented by Eq. (14), to that of the activation energy in the Arrhenius law[13]. They both govern the rapidity of an exponential falloff.

Discussion

In summary, we propose the transmembrane ion gradients that permit neurons to form and propagate an action potential also play a critical role in the function and survival of non-neuronal cells. We propose this difference between extra- and intra-cellular ion concentrations produces initial conditions that allow reception, communication, processing and response to environmental information (Fig. 1). The information is detected by specialized gates within transmembrane ion channels that, by opening, allow transmission of the information into the cell in the form of an ion puff. By changing the structure of the proteins that form the gates, they can be made to be sensitive to a wide range of environmental factors. Furthermore, the information encoded in the flow of ions is received and processed in the cell through a local change of ion density and cation/anion concentrations around the channel. This may reduce (or increase) the mobile ion screening of negative charges in the inner leaf of the cell membrane and by increasing (or decreasing) the Debye length[28,29], alter the attachment of peripheral membrane proteins (Fig. 1B). This can alter the structure of the membrane to respond to the perturbation. Furthermore, the function of a large number of enzymes is dependent on local is dependent of the local ion concentrations. For example, hexokinase[30,31] is dependent on local Mg2+ concentrations and pyruvate kinase[9] is dependent on K+ concentrations. In this way, the ion flux through K+ or Mg2+  [32] channels can promote or suppress reactions to rapidly respond to a perturbation.

Since the discovery of the structure and function of DNA, investigations of biological information dynamics have focused almost exclusively on gene replication and expression, and on translation. These include transgenerational transmission of information through chromosomes, epigenetic control and time dependent changes in transcription of genetic information in response to local environmental changes. However, in prior theoretical studies we have noted that the transmembrane ion gradients potentially may encode in excess of 1014 bits of Shannon information[1,13,14]. Further, the membrane ion protein pumps that generate these gradients can consume up to 40% of a cell’s energy budget[33]. Why does a cell invest such a large fraction of its resources to the formation of a transmembrane ion gradient? Indeed, at first glance, the large expenditure of energy to achieve the seemingly trivial task of replacing intracellular Na+ by K+ (or replacing Ca2+ with Mg2+), seems at best puzzling.

Here we propose a specific mechanism by which the transmembrane ion gradients, and their associated information dynamics, can detect, process and respond to environmental perturbations on a very rapid time scale without any direct input from the genes. That is, while the genome clearly contributes to these dynamics by encoding and translating the information in the associated proteins, the ion-based dynamics at the membrane can proceed locally and rapidly without additional input from the nucleus. As noted above, this model of information dynamics optimizes both speed and accuracy of information transmission. Thus, we propose that the large resource expenditure necessary for generation of the ion gradient is an indication of the critical role of these information dynamics in cell survival and function.

In general, our results support a more distributive model of biological information. Thus, information dynamics at the cell membrane are deeply linked with heritable information which contributes the molecular machinery that generates the gradients and permits detection, processing and response to environmental signals and perturbations. Similarly while we view the membrane information dynamics to occur independently on fast and local scales, we note that ion fluxes can generate molecular signals including lipid second messengers[34] that communicate with the nucleus and other cellular organelles. Furthermore, ion flows may facilitate direct communication of cellular information to the nucleus. For example, ligand binding to EGFR is associated with gate opening in K+ channels[5,35]. This loss of K+ screening for the negative charges in the inner leaf of the membrane (Fig. 1A,B) may facilitate transmission of signal along the MAPK pathway by allowing coulomb interactions between positively charged Raf protein (pK of 9.3) or its associated scaffolding proteins to the GTP-activated Ras protein. Indeed, it is often stated that during MAPK signaling that Ras “recruits”[7] Raf to the membrane. We note that the term “recruits” has no physical meaning but rather represent empirical observations that Raf moves very rapidly to the location to the membrane site containing an activated Ras. We propose the ion dynamics summarized in Fig. 1 could account for this rapid and specific movement.