Solution of non-linear Poisson-Boltzmann equation for erythrocyte membrane

Solution of non-linear Poisson-Boltzmann equation for erythrocyte membrane

Frederico A. O. Cruz, Fernando S. D. S. Vilhena, and Célia M. Cortez

Universidade do Estado do Rio de Janeiro

Departamento de Ciências Fisiológicas

Av. Prof Manuel de Abreu, 48

20550-170, Rio de Janeiro, RJ - Brazil

Received 25 October, 1999

The space-dependent electric potential equations were determined by solution of non-linear Poisson-Boltzmann equation, which decribes an erythrocyte membrane model which takes into account the charges of molecules associated with lipidic bilayer (glycocalyx molecules, outer side, and spectrin, inner side). It was assumed that on both sides of the membrane the charges are homogeneously distributed. Therefore, the potential was dependent on the z-coordinate only. Thus, we obtained the potential profile curve. The results are compared to previous ones obtained from more simple membrane model (Bioelectroch. Bioenerg. 1993,

32

:305-315).

I Introduction

It is well known that the passive transport of charged particles across biologic membranes is strongly affected by electrostatic forces and that the flux rates are dependent on intramembrane field forces [1-5]. But realistic values of membrane potential from experimental searches are obtained only for transmembrane potential (i.e., difference of potential between two points distant from outer and inner surfaces) [6]. The determination of potential profile, which describes the potential behaviour on all adjacent points along perpendicular axis to membrane, is performed based on theoretical models [6, 7, 8]. This is a way to estimate the potential in the vicinity of the membrane. However, informations from experimental searches are essential to formulate such mathematical models which allows the calculation of the membrane potential in function of space. Thus, several models have been suggested for biologic membrane and several electrostatic theories have been developped [1, 2, 6, 9]. Some modern concepts, such as interionic correlations and surface conductivity, were already introduced within of the Gouy-Chapman theory, being its range of validity still discussed [10]. According to Cevc [6], the effects due to surface irregularities and to interfacial hydration decrease the significance of interionic correlation effects in several membranes.

In previous work [1], we associated some recent knowledge on membrane electrostatics with classical Gouy-Chapmann model to describe erythrocyte membrane. For this model we performed a systematic study on the effects of ionic strength change on the electric potential profile. In the present work, we determine the potential profile for a complex membrane model, which takes into account electric charges of glycocalyx and spectrin layers. Donath and Patushenko [11] already have taken into account space charges of the glycocalyx for calculations of the outer surface potential. Heinrich et al [2] used space distributions of charges of glycocalyx, spectrin layer and hemoglobin for numerical calculation of potential profiles. Here, we determined the analytical solution for non-linear Poisson-Boltzmann equation, which includes parameters due to fixed charges on glycocalyx and spectrin. The results are compared with previous results obtained for a simple model [1]. Analytical solution of the electric potential in neighboring of the membrane is of special importance for studies involving more complex mathematics, such as the study of mechanical stability of membrane by hydrodynamic analysis [12-16].

II Adopted model

As already known, the erythrocyte membrane is constituted by a double bilayer of phospholipids, cholesterol and proteins [17], being asymmetric by distribution of phospholipids between the two lipid monolayers [18]. Molecules of phosphatidilserine and phosphatidylethanolalanine are preferentially located in the inner monolayer and phosphatidylcholine and sphingomyelin are founded in the outer monolayer [19]. In the outer surface there is an additional layer rich in polysacharides, which are associated with the bilayer molecules, this is called the membrane glycocalyx layer. On the inner surface there is a protein trap that supports the membrane, the spectrin [20]. In a solution of NaCl at pH 7.0, erythrocytes exhibit a negative electrophoretic mobility. This is especially due to charged molecules of the glycocalyx layer [6]. In the same form, spectrin has a net negative charge at physiological pH [21]. In addition, all membranes are embedded in diffuse inorganic salt double layers [6].

Based on these knowledge, we adopted for our study the model shown in Fig. 1(a) and (b). The system is formed by five regions: extracellular bulk region (phase 1), glycocalyx layer (phase g), lipidic bilayer (phase f), spectrin layer (phase p) and cytoplasmatic region (phase 2). Within glycocalyx and spectrin layer there are space charge distributions denoted respectively by r_g and r_p. The charge densities of these two later layers are considered almost constant in the z direction. The lipidic bilayer has dielectric constant e_f and thickness h (extends from z=-h/2 to z=h/2), and it is limited by surfaces S₁ and S₂, as can be seen in Fig. 1(b). At these points the glycocalyx/bilayer and spectrin/bilayer interfaces are located, respectively. Similar to the Heinrich model [1], the inner surface charges are due only to the head group of phospholipids. So, there are electric charge distributions on S1 and S2 and their densities are denoted by Q_S1 and Q_S2 (Fig. 1(b)). The origin of the Cartesian coordinate axis is located at the bilayer centre. The electric potential depends only of the z coordinate, because charge distribution is considered constant along x and y axis. In aqueous phases, the electrochemical potentials, m_i, of each solute ion g are considered constant. The charge density r in the bilayer film is taken as zero (r_i(+) = r _i(-) = 0). S_g is the phase 1/glycocalyx interface and S_p is spectrin/phase 2 interface.

Figure 1.
(a) Adopetd membrane model. The erythrocyte membrane is constituted by a lipidprotein bilayer, a outer layer rich in polysacharides, the glycocalyx, and a inner layer of spectrin. The bilayer surfaces located at positions z = ±h/2. The glicocalyx extends from -h/2 to -h_g and the spectrin layer extends from h/2 to h_p. (b) Phases 1 and 2 correspond to extracellular medium and cytoplasm medium, respectively. f₁(z) represents electric potential along z axis into phase 1, being f₀₁ the potential for z ® -¥ and f_Sg the potential on z = -h_g (on S_g interface). In phase 2, we have f₂(z) as the electric potential along z, f₀₂ is the potential for z ® ¥ and f_Sp is the potential on z = -h_p (on S_p interface). Phases g and p correspond glycocalyx and spectrin layer, respectively, being f_g(z) and f_p(z) the electric potentials along z axis in these phases. S₁ and S₂ are limit surface of bilayer and f_S1 and f_S2 are potentials on these interfaces. Q_S1 and Q_S2 are outer and inner surface charge density of bilayer.

III Calculation of potential profile

The Poisson equation applied for the five considered phases is

and it solution give us the electric potential f(x, y, z) at any point in each one. The considered boundary are as follows (see Fig. 1 (b)):

• (I) When

  z ® -¥ and

  z ® ¥, (a) the electric potential tends to limiting values f

  ₀₁ and f

  ₀₂, respectively, and (b) the corresponding ionic concentration tends to limiting values h

  _g01 and h

  _g02. In these regions one has the electroneutrality conditions: S

  _g

  Z

  _g

  e h

  _g01 = 0 and S

  _g

  Z

  _g

  e h

  _g02 = 0.

• (II) At

  z = -

  h

  _g,

  z = ±

  h/2 and

  z =

  h

  _p, the electric potential is continuous through the surfaces and takes the values f

  _Sg at the surfaces

  S

  _g and f

  _Sp at surface

  S

  _p, f

  _S1 and f

  _S2 at

  S

  ₂.

• (III) In these surfaces the condition of discontinuity of the electric displacement vector is valid.

Solution within phases 1 and 2. The potential equation for 1 and 2 phases already was calculated in previous paper [1]. In this regions we have only aqueous phase. Using the Boltzmann equilibrium, symmetric electrolytes and electroneutrality condition in both phases we obtained,

k_i is related to the of Debye constant, f_Sj the surface potential on S_j (j = g, p), f_oi is the bulk potential in phase i, b = Ze/KT, K is Boltzmann constant, Ze is molar charge. The reference state is assumed to be stationary.

IV Solution within glycocalyx and spectrin layers

Within the Glycocalyx and spectrin layers the Poisson equation is non-linear due to addition of the term related to fixed charges on these molecular nets [2].

where r_j is the electrolytic charge density and r_fj refer to fixed charges into glycocalyx (r_fg) or spectrin (r_fp).

Using a Boltzmann distribution for the ions in aqueous phases and considering charge homogeneity on directions x and y, we have for each type of ions:

where Z₀e is the molar electric charge h_0j the corresponding molar concentration of ions in solution within the phase j, and e_j dielectric constant. In simplified form we have,

where

The first order solution of equation 6 give us

Applying the boundary conditions II on S_g e S_p, we have

where E_Sj = ¶f_j/¶z, at S_j.

Thus, we can write

being

The solution of equation (13) is not trivial, hence simplifications was required. Using Taylor's expansion, we obtained that

which is a infinite series, where a_j = (l/c_j)^1/2.

Taking into account experimental data from literature for erythrocyte parameter values, and previous results [1], we studied the series terms. We found that, for u_j < 200mV, the terms decrease significantly as the u_j-exponent increases. So the higher order terms do not contribute appreciably and the equation can be reduced to

where b_j = (l/2) g_ja_j³.

Solving the first order differential equation and considering electric potential continuity (condition II) at surfaces S₁ and S₂,

being

So, we can write: for glycocalyx region:

where

for spectrin region:

where

Thus, from equations 18 and 19 we obtain, respectively,

and

V Solution within bilayer

Considering that in lipidic bilayer charge density can be vanished, we have

being the solution of a linear function family.

Taking into account the vector electric field displacement discontinuity (condition III) on the S₁ and S₂, and making some algebraic manipulations, we have:

being

VI Results

Applying the referent values of the erythrocyte membrane (table 1) in results obtained here and using the interpolation method, we could trace the potential profile curve (Fig. 2). Fig. 3 shows the potential profile for the simple model based on the Gouy-Chapman theory of difuse layers, presented in anterior work [1].

Table 1.
Numerical values os the parameters used of the calculations

Figure 2.
Electric Potential Profile of the erythrocyte membrane calculated by use of the potential equations from analytical solution of nonlinear Poisson-Boltzmann equation. Numerical values adopted are shown in table 1. (Ionic strength = 172 mM).

Figure 3.
Electric Potential Profile of the erythrocyte membrane calculated by use of the potential equations from analytical solution of the linear Poisson-Boltzmann equation. Q_S1 = 1.41 × 10^-3 A.s/m², Q_S2 = -9 × 10^-2 A.s/m², h = 10 nm. (Ionic strength = 172 mM).

VII Discussion and conclusions

The inclusion of the volumetric charge distribution into Poisson's equation (equation 1) led us to a Poisson-Boltzmann equation one with non linear characteristics and not trivial analytical solution. Thus, after search into possibilities we chose Taylor's series expansion to simplify its second integration (equation 13). However, we observed that the series (equation 15) was an infinit series, and did not converge for any potential value. Thus, it had required a detailed study to find the intervals where it became convergent, and if it was possible into expected interval of potential values, taking into account numerical values of the involved potential equation parameters. So, we investigated the contribution of each series term and we verified convergence for numerical values of erythrocyte parameters. Indeed, the values of the terms decrease significantly with increasing of u_j-exponent, considering a potential interval between 0-200 mV, or 0 £ (u_i-u_j) £ 200 mV. Thus, it became possible to obtain a analytical result for potential, which can be easely used for posterior studies involving more complex calcules. The potential profile curve shown in the Fig. 2 was traced based on these obtained equations.

Using numerical integration, Heinrich et al [2] obtained the potential profile for their membrane model, which considered glycocalyx and spectrin charges hemoglobin charges. They studied the influence of ionic strength on this profile. Their results showed that, for all osmotic states studied, a negative inner surface potential (S₂) of more than 60 mV was calculated. For a ionic strength of 145 mM the potential on the S₂ was about -75 mV, and on S₁ (outer bilayer surface) was about -10 mV. They also shown that at a decrease of ionic strength causes a increase of negativity of potential on S₂. Our potential profile (Fig. 2), which was obtained from analitical solution of Poisson-Boltzmann equation, presents the same shape than those shown by Heinrich et al [2]. We considered in our calculation a ionic strength of 172 mM.

In a previous paper [1], we adopted a simple model, which was based on the Gouy-Chapman theory of diffuse layers, in which space distributions of fixed charge of glycocalyx and spectrin layers was neglected. Comparing the curve of Fig. 2 with that of previous results (Fig. 3), we observed that the profile suffers some change when we took into account these additional charges. Within layers where there are fixed charges, the potential varies less pronouncedly than when these space distribution of charges are neglected. In this case (Fig. 3), the observed potential varies exponentially in both two sides, and the linear part of the inner curve presents a pronounced inclination. As can be seen in Fig. 2, the inclination of potential profile in the neighboring of the bilayer is smaller than that shown in Fig. 3 in the same place. In the simple model [1], we adopted a distance of 10 nm between the two parallel planes which limit the membrane region. So, the outer surface was located approximately at half the distance between the bilayer surface and the more extreme point of glycocalyx layer. All charges belonging to this membrane side (bilayer surface and neighbouring fixed charges) were taken into account as surface charge. Thus, in that case, we circumvented space charges by treating them as surface charges becoming trivial the solution of the Poisson-Boltzmann equation by analytical means. By confront of these two figures, it can be observed the importance of glycocalyx charges on potential within the outer membrane surface, as well as the potential within spectrin layer at the inner surface is mainly determined by the charges in this region.

Applying the analytical solution of the electric potential of the membrane in the hydrodynamic analysis, we can study the influence of these charges on the mechanical stability of the membrane. As can be seen in Fig. 2, in spectrin region there is an important electric field, and its effect on transport phenomena through membrane may be imperative, as well as on regulation of other membrane functions.

• [1] C. Cortez Maghelly and P.M. Bisch. Bioelectroch. and Bioenerg 32, 305 (1993).
• [2] R. Heinrich, M. Gaestel and R. Glaser. J. Theor. Biol. 96, 211 (1982).
• [3] S.J. Kennedy. J. Membrane Biol. 42, 265 (1978).
• [4] R. Glaser, I. Bernhardt and E. Donath. J. Electroanl. Chem. 116, 281 (1980).
• [5] G. Schwarz. J. membrane Biol. 43, 127 (1978).
• [6] G. Cevc. Biochim. Biophys. Acta 1031, 311 (1990).
• [7] J. Baker, W.T. Coakley and D. Gallez. Eur. Biophys. 22, 53 (1993).
• [8] H.P. Schnebll, C. Roeder and L. Tarcsay. Experimental Cell Research 98, 273 (1976).
• [9] C.C. Wang, L.J. Bruner. Biophys. J. 24, 749 (1978).
• [10] E. Donath and A. Voigt. Biophys. J. 49, 493 (1989).
• [11] E. Donath and V. Pastushenko. J. Electroanal. Chem. 104, 543 (1979).
• [12] M. Prevost, P.M. Bisch and A. Sanfeld. J. Colloid Interface Sci. 88, 353 (1982).
• [13] P.M. Bisch, H. Wendel and D. Gallez. J. Colloid Interface Sci. 92, 105 (1983).
• [14] W.T. Coakley and D. Gallez. Interfacial Stability. In: Biophysics of the Cell Surface (Glaser, R. & Gingell, D., eds), p.287. Berlin (1990) Springer Verlag.
• [15] C. Cortez Maghelly and P.M. Bisch. J. Theor. Biology 176, 325 (1995).
• [16] C. Cortez-Maghelly and P.M. Bisch. Braz. J. Phys. 26, 604 (1996).
• [17] E.M. Manevich, K.M. Labin and A.I. Archakov. Biochim. Biophys. Acta, 815, 455 (1985).
• [18] P. Gascard, D. Tran, M. Sauvage, J.C. Sulpice, K. Fukami, T. Takenawa, M. Claret and F. Girand. Biophys. Chem. 1069, 27 (1991).
• [19] P.J. Raval and D. Allan. Biochim. Biophys. Acta. 722, 192 (1982).
• [20] G.L Nicolson, V.T. Marchesi and S.J. Singer J. Cell. Biol. 51, 265 (1971).
• [21] A. Elgsaeter, D.M. Shotton and D. Branton. Nature 307, 378 (1984).
• [22] E. Donath and D. Lerche. J. electroanal. Chem. 116, 41 (1980).
• [23] M. Blank, R.G. King, L. Soo, Re. Abbott and S. Chim. J. Coll. Interface Sci. 69, 67 (1979).
• [24] D.S. Parson. Biological Membranes Oxford (1975) clarendon Press.
• [25] G.M.W. Cook, D.H. Heard and G.V.F. Seaman. Nature. 191, 4449 (1961).
• [26] R. Schauer, A.P. Cornfield, M. Wember and D. Danon, Physiol. Chem. 356, 1727 (1975).

Publication Dates

• Publication in this collection
  07 Jan 2002
• Date of issue
  June 2000

History

• Received
  25 Oct 1999