The surface charge of a cell lipid membrane - arXiv

The surface charge of a cell lipid membrane

M. Pekker 1 and M. N. Shneider 2

1 MMSolution, 6808 Walker Str., Philadelphia, PA 19135, USA 2Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, NJ 08544, USA

E-mail: [email protected] and [email protected]

Abstract In this paper the problem of surface charge of the lipid membrane immersed in the physiological solution is considered. It is shown that both side of the bilayer phospholipid membrane surface are negatively charged. A self-consistent model of the potential in solution is developed, and a stationary charge density on the membrane surface is found. It is shown that the ions of the surface charge are in a relatively deep (as compared to kBT) potential wells, which are localized near the dipole heads of phospholipid membrane. It makes impossible for ions to slip along the membrane surface. Simple experiments for verifying the correctness of the considered model are proposed. A developed approach can be used for estimations of the surface charges on the outer and inner membrane of the cell.

Introduction

The question about the distribution of electric charge near the cell membranes is the key for many problems associated with the interaction of cells with external electromagnetic fields [1-6]. For example, if ions are not electrically neutral, the areas are not strongly bound with the membrane and can move freely along it, the electromagnetic field has a weak influence on the membrane. If the ions are firmly bound to the membrane, the electrical component of electromagnetic fields could lead to the membrane deformation [7, 8].

The surface charge of the cell membrane has been studied in many experimental works. We will not consider all works and may refer only to [9,10], where the experimental data of the surface charge of cell membranes obtained by different methods are given. A range of results is broad enough: 2/002.03.0 mCm −=σ . The problem of spatial distribution of charge in the vicinity of the biological membranes surface has been considered in many теоретических papers, see eg [1-5]. In all these works, the near-surface potential of the membrane was considered under the Gouy-Chapman theory [11, 12] or its later modification by Stern [13], in which the charge on the membrane surfaces is considered to be given. In these theories, the membrane was considered as a continuous dielectric, without taking into account its fine structure, and a surface charge was determined on the basis of the electrochemical properties of the dielectric surface (see, eg [14, 15]).

In the classical Gouy-Chapman-Stern theory the surface charge ions can freely slide along the surface. Because, it is assumed in this theory, that the interaction of ions with the surface occurs due to electrostatic

forces witsurface. Inform a mo(Ris.1b), iinteractionthat the lovibrations

Fig. 1a. Athe inside

Fig. 1b. distance e

In our papthat the suthe dipoleion bindinallowed tdetermine

The memimpact ofchannels.

th the inducen the case of aosaic (lattice)i.e. sliding thens of the weaongitudinal es of the memb

A simplified and the outsi

A fragment equal to 0.15 o

per, we take iurface of the ce heads of theng energy U the constructie the average

mbrane surfacf weak microw

ed charges ana phospholipi) structure (Fem along the ak microwavlectric field o

brane, becaus

mosaic modeide of the mem

of a two-dimof the lattice p

into account cell membran

e membrane pwith the me

ion of a selfcharge densit

ce charge thewave electrom

d absorption id membrane,

Fig.1a) [16]. Asurface of the field with tof the microwe the surface

el of a phosphmbrane.

mensional disperiod (a desc

the fine strucne is charged phospholipidsembrane surff-consistent thty of sitting o

eory presentedmagnetic field

forces, which, the electrolyAs a result, the membrane the cell memwave acting ocharge ions a

holipid memb

stribution patcription of the

cture of the binegatively w

s. These ions face is much heory of Goun the membra

d in this papds on the surf

h depends onyte ions interahe surface chis suppressed

mbranes, propoon the surfacare rigidly bon

brane [16]. Sp

ttern of potene computation

ilayer phosphwith ions, trapp

are strongly bgreater than

uy-Chapman ane.

per supports face density r

nly on the disact with phosp

harge ions ared. This fact is osed in [17], e charge caunded to the su

pheres represe

ntial along mnal model is g

holipid membped in the pobound to the

n the thermal and Stern th

the idea [17redistribution

stance of the pholipid dipo

e in a periodiccrucial in thewhere it wa

uses forced lourface and can

ent the dipole

membrane surgiven in Parts

branes [16]. Iotential wells

membrane suenergy kBT.

heories, and

] of electromn of transmem

ion to the oles, which c potential e theory of s assumed

ongitudinal nnot slip.

e heads on

rface, at a s 1 and 2).

t is shown formed by urface: the

That fact helped to

mechanical mbrane ion

We emphasize that this article does not discuss the passage of ions through the transmembrane ion channels or pores. We restrict ourselves to considering only the surface processes in the lipid membrane placed in the salt-water solution. The main ions in the physiological solution аrе: Na+, K+, Са2+, Mg2+ and Сl-. The total concentration of these:

26108.1 ⋅≈in m-3 (300mmL/Mole) [21], но для простоты мы рассмотрим только случай, когда the bilayer phospholipid membranes immersed in water solution of NaCl.

1. A model of the phospholipid membranes.

The dielectric permittivities of the membrane and water are correspondingly: εm ≈2 and 80≈wε [18, 19]. It

would seem that in order to evaluate the force acting on the ion of charge q near the membrane surface the formula [20] can be applied:

( )mw

mwq h

qFεεεεε+−

=0

2

2

4 (1)

Here ε0 is the permittivity of free space, and h is the distance to the interface between dielectrics. However, the paradoxical result follows from the formula (1) that the ion cannot get from the dielectric with the higher dielectric permittivity into the dielectric with the lower dielectric permittivity, because the force acting on it increases inversely proportional to the square of the distance to the boundary between dielectrics. This statement contradicts the experimental facts and theoretical models of Hodgkin and Huxley’s kind, based on the phospholipid membrane permeability to ions. Usually, in theories of the interaction of ions with the surface of dielectrics, the minimum value of h is selected as the size of the ion, or the Debye radius, if the insulator’s surface is saturated with ions [15]. However, this does not resolve the problem of phospholipid membranes "impenetrability". In fact, the problem lies in that the formula (1) is derived in the approximation of continuous medium with the dipoles of the infinitely small size, without taking into account the real structure of the surface layer of the membrane.

It is known that the phospholipid molecules of cell membrane are forming a mosaic (matrix) structure in which dipole heads are directed towards the liquid (positively charged head faces outward membrane) [16]. The average surface area per molecule of the lipid is 5.0≈ nm2, the length of the polar head is 15.0~ − nm, the radius of the head is 3.02.0~ − nm, and the distance between the hydrophilic heads of the membrane is in the range of 5-7 nm [21, 22]. The dipole moment of the phospholipid head is 18.5-25 D [23] (1 D =3.34·10−30 C.m), i.e. more than 10 times greater than the dipole moment of water molecules. On the basis of geometrical dimensions of the cell membrane and the size of water molecules, it can be concluded that the free space between the head does not exceed the size of a water molecule (~ 0.2nm). That is, the membrane, interacting with the ions of surrounding liquid, cannot be considered as a dielectric medium with an infinitesimal dipoles size.

These facts allow the consideration of the following simplified model of the ion interaction with the membrane: 1. The membrane represents a matrix (Fig. 2) with a mesh size aa × . In the nodes of cells the dipoles are

located; the dipole charge is q; the distance between the charges (the dipole length), d; the distance between the dipoles along the axis z, l.

2. An ion is a classical particle and cannot approach a dipole at a distance less than the sum of the radii of the head and the size of the ion. It is important that the ion can approach the membrane dipole heads close

enouginterathan ttwo wwatermemb

3. We wconsid

Fig. 2. Athe mesh the dipole

The main 8.1 ⋅≈in

phospholi

2. The po

The expresurface of

a

gh that wouldaction of ions the dipole mowater moleculr molecules cbrane. want to emphder in this wo

A simplified msize of the m

es along the z

ions in the p2610 m-3 (30

ipid membran

otential near

ession for thef the membran

a

d be "capturedwith the surfa

oment of the les, the intera

can be neglec

hasize one moork, for us it is

mosaic modelmatrix, d is the

axis.

physiological 00mmL/Molenes are immer

the surface o

e membrane ne (Fig. 3):

y

d" by the poteface of the diephospholipid

action betweected, as comp

ore time that s important to

l of a membrae distance bet

solution аrе: ) [21], but frsed into the w

of the membr

potential in a

ential well. Telectric [15]. Sd head molecuen the ions lopared with th

ionic permeao know how s

ane: the memween the cha

Na+, K+, Са2

for simplicitywater solution

rane

a cylindrical

x

This is a standSince the dipoule and near ocated near the interaction

ability membstrong ions at

mbrane is a marges of the di

+, Mg2+ and Сy we considen of NaCl.

coordinate sy

dard assumptiole moment othe head can

he surface of of the ions w

brane is a septtached to surf

matrix with dipipole head, l

Сl-. The total er only the

ystem ( zr ,,θ

ion in the theof water is 10

not be moref the membranwith the dipo

parate problemface of memb

poles in the nis the distanc

concentrationcase when th

)z at a point

eory of the times less than one-ne and the oles of the

m and not brane.

nodes. a is ce between

n of these: he bilayer

above the

( ) ( )

( )

( ) ( )

( ) ( )

( ) ( ) ⎟⎟⎟

⎠

⎞

+++−⋅−++

++−⋅−+−

++−⋅−+−

⎜⎜⎜

⎝

⎛

+−⋅−+

==

∑ ∑

∑ ∑

−= −=

−= −=

2,,

2,

2

2,,

2,

2

2,,

2,

2

2,,

2,

20

,

cos2

1

cos2

1

cos2

1

cos21

4

,,,,

dlzrrrr

lzrrrr

dzrrrr

zrrrrq

zrUzrU

jijiji

jijiji

jijiji

n

ni

n

nj jijiji

n

ni

n

njji

θθ

θθ

θθ

θθπε

θθ

(2)

Here, the radial and angular coordinates jir , , ji ,θ correspond to the position of the dipole in the node with the

numbers (i, j), q is the dipole charge, and 0ε is permittivity of free space. The value of n in (2) is chosen such

that the potential near the membrane does not depend on the size of the matrix n. The count goes from the

node (dipole head) (Fig. 3): ( ) 2/122 yxr += , ( )xy /arctan=θ .

We expand the potential U in a Fourier series byθ . Taking into account the rotational symmetry on 2/π rotations for the accepted square matrix of the dipole heads, we obtain:

( ) ( ) ( ) ( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( ) ( ) ( ) ( )∫∫

∫

∑∑

==

==

++=

≠≠

==

2/

00

2/

00

2/

000

110

4sin,,4, 4cos,,4,

0, ,,2,

4sin,4cos,,,,

ππ

π

θθθπ

θθθπ

θθπ

θθθ

djzrUzrVdjzrUzrU

zrVdzrUzrU

jzrVjzrUzrUzrU

jj

J

jp

J

jj

(3)

Fig.3. Spalattice site

In the framlower sidethe effectdependennot exceed For exampthe follow Fig. 4. sh

TkU B/ o

atial lattice oes correspond

mework of the of membrant of dipoles ce of U on θd 7%, so we w

ple, let considwing set of par

hows contour

on the axis z.

f dipoles in bd to the dipole

he model presene (Fig. 2) on in the upperθ , in the regiwill hencefort

der the case, rameters: =a

plot of reduc

bilayer lipid me heads of pho

ented above, the potential

r side of theon 2/ =≤ arth neglect the

when bilayer7.0= nm, =d

ced TkU B/

membrane. Poospholipid mo

estimates sho distributionU

e membrane 35.0= nm th

em.

r phospholipid5.0= nm, =l

for a singly c

otential is calolecules.

ow that the inU at 0>z con the poten

he contributio

d membrane i8= nm, and q

charged ion a

lculated at th

nfluence of thecan be neglectntial at z −<

on of above th

is immersed i19106.1 −⋅=q

at T = 300K,

he point P. El

e dipoles located. The same( )dl +− . As

he zero harmo

in the NaCl sC [16, 21].

Fig. 5 the de

lipsoids at

ated on the e is true to far as the onics does

solution, at

ependence

Fig. 4. Contour plot TkU B/ , corresponding to a half of the cell shown in Fig. 1b.

0.25 0.50 0.75 1.000.1

1

10

100

4

32

U/k

BT

z/a (a=0.7nm)

r/a=0 r/a=0.14 r/a=0.21 r/a=0.2975

1

Fig. 5. The parameter TkU B/ dependence on z. Curve 1 corresponds to ar / = 0, 2 – ar / = 0.14, 3 - ar / =

0.21, and 4 – ar / =0.30.

It is seen from Fig, 5 that the value TkU B/ decreases exponentially with the distance from the membrane,

wherein a distance 5.0/ =az , the ratio TkU B/ is about 10. That is still sufficiently deep potential well, but

already at 1/ =az , the influence of the membrane on the ions can be neglected.

11

1

55

5

5

10

10

10

20

20

20

3030

30

40

4060

6080100

r/a

z/a

0 0.05 0.1 0.15 0.2 0.25 0.3

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

1.1

0

10

20

30

40

50

60

70

80

90

100

An effective potential of the ion across the membrane can be determined. Let the ion be a sphere of radius ir ,

and hr is the distance from the positive charge of the dipole to the head of the ion (Fig. 6).

Fig. 6 The negative ion located near the head of the phospholipid molecule dipole. Ion is a charged sphere of radius ir , and hr is the distance from the positive charge of the dipole to the head of the ion.

In this case, the effective reduced ion potential,

( )( )

∫ ∫+ −+−

=ih

h

ihirr

r

zrrr

iBB

ef rdrdzzrUdzrTkTk

U 2

03

22

,4

32π

π (4)

The dependence of the effective potential of the ion −Cl in the potential well on the size of the head hr is

shown in Figure 6. The chlorine ion radius ir = 0.18 nm, and the atomic weight is 4.35=iM .

It is seen that for the typical size of the dipole head 2.0=hr nm ( 29.0/ =arh ) and the ion radius 18.0=irnm, the reduced potential depth is 10/ ≈TkU Bef . Thus, the trapped ion is quite firmly bound to the

membrane. It should be noted that for electrolytes of a different composition of negative ions, for example in the axon, where the major negative ions are anion groups of macromolecules and phosphates, estimations of the potential will differ from that shown in Fig.7.

ir

hr

z

Ion

Dipole head

Fig. 7. Thion ir = 0.

3. Evalua

Foras the miapproach

We Since the of mesh ssize, we cmembrane

Fig. 8. A the ions fr

Let us estthe memb

he dependence.18 nm (Cl-).

ation of the m

simplicity, wnimal distancin models of assume that field of the m

size is negligican relate the e/liquid boun

scheme of throm the poten

timate the fluxbrane:

0.150

5

10

15

20

25

30

35

U/k

T B

e of the effec

membrane su

we assume thece that the iothe ions interthe space ch

membrane (wibly small, ancharge on th

ndary (Fig. 8).

he ions flows ntial wells of t

x of ions falli

0.20 0.25

tive potential

urface charge

e ions as a poon can approaraction with thharge on the m

without takingnd the distancee membrane .

near the memthe depth U.

ing on the me

0.30 0.35

rh/a (a=0.5n

l of the ion on

e.

int charges, bach to the mhe dielectric smembrane do

g into accounte between thesurface to the

mbrane. The

embrane. Assu

5 0.40 0.4

nm)

n the size of th

but will take imembrane. As

surfaces. oes not affectt the ionic chae free ions ine charges in s

ions flux to t

uming the Bo

45 0.50

the head at the

into account t indicated ab

t its structurearge on it) at

n the solution solution by eq

the membran

oltzmann dist

e assumed rad

the finite sizebove, this is a

e, i.e. the dipt a distance ofis larger than

quating the fl

ne and “evapo

tribution of io

dius of the

e of the ion a standard

pole heads. f the order n the mesh lows at the

oration” of

ons outside

⎟⎟⎠

⎞⎜⎜⎝

⎛−= ∞ Tk

qnnB

ϕexp , (5)

where ∞n - is ion density at infinity. Assuming that the ions are in thermal equilibrium with the water molecules, we obtain the following estimate for the ion flux on the surface of the membrane:

⎟⎟⎠

⎞⎜⎜⎝

⎛−⎟⎟

⎠

⎞⎜⎜⎝

⎛= ∞ Tk

qM

MvnP

Bi

OHOHin

ϕπ

exp2

12/1

2

2. (6)

Here 2/1

2

2 ⎟⎟⎠

⎞⎜⎜⎝

⎛

i

OHOH M

Mv is the averaged ion velocities in solution, and iOH MM ,

2 are the masses of the water

molecule and the ion, correspondingly.

On the other hand, an estimate of the ion flux “evaporating” from the surface of the membrane is as follows. The flow of ions from the surface of the membrane can be considered as a process of evaporation of the liquid with a work function equal to the potential well U. The ion in the well is exposed to random impacts by water molecules, while the ion’s presence time in the potential well can be estimated as:

( )( )1/exp2

2 0, −⎟⎟⎠

⎞⎜⎜⎝

⎛≈ kTU

MM

OH

iHii ττ . (7)

Here ( )OHiHi MM22

/0,τ is the characteristic time of energy exchange between the ion and the water

molecules. In (7) we have taken into account that at 1/ <<kTU the ions are reflected from the membrane, i.e., the time delay is zero. For estimates 0, 2Hiτ can be assumed equal to the collision time of water molecules

in a solution TOHHi v/22 0, λτ = , 9102.0

2

−⋅=OHλ nm is the mean free path of a molecule of water, and Tv is

the thermal velocity of water molecules. Accordingly, the number of the ions leaving the unit surface of the membrane per unit time is equal to:

( ) 1/exp11 2

2

2

20

2 −⎟⎟⎠

⎞⎜⎜⎝

⎛⋅

==TkUM

Ma

vNNa

PBi

OH

OH

Hi

i

iout λτ

. (8)

Here 2/1 a is the number of potential wells per unit membrane surface; iN is the relative population of the

potential wells with ions. Equating the incident flux of the ions on the membrane to the "evaporating" flow from it, we obtain the relative population of the potential wells:

⎟⎟⎠

⎞⎜⎜⎝

⎛−⎟⎟

⎠

⎞⎜⎜⎝

⎛⎟⎟⎠

⎞⎜⎜⎝

⎛−⎟

⎟⎠

⎞⎜⎜⎝

⎛⋅=

− ∞ 1expexp2

11

0

2/1

2

2

2 TkU

Tke

MMan

NN

BBOH

iOH

i

i ϕλπ

. (9)

In (8) we took into account that the probability of the incident ion being captured for in a potential well is proportional to iN−1 .

iτ

Now we will find the distribution of electric potential in the liquid. Following the Debye approximation for the liquid electrolyte [11]:

( )Dx λϕϕ /exp0 −= (10)

Where 20

2 qnTkB

D∞

=εε

λ is the Debye length.

The field on the membrane surface:

00 2εεσϕ

=∂∂

=xx (11)

Where σ is the surface charge density. For the case when the upper and lower surfaces of the membrane are equally charged 2/2 aqNi=σ , we obtain the potential at the membrane surface:

20

20

0 2 aN

nTk

aNq iBiD ⋅=⋅

=∞εεεε

λϕ (12)

Substituting (12) into (9), we obtain the equation for the value of iN :

( )

Tknaq

TkUN

MMan

NN

B

Bi

OH

iOH

i

i

∞

∞

=

⎟⎟⎠

⎞⎜⎜⎝

⎛−⎟⎟

⎠

⎞⎜⎜⎝

⎛⋅−⎟

⎟⎠

⎞⎜⎜⎝

⎛⋅=

−

02

2/1

2

21

1expexp2

11

2

2

εεη

ηλπ

(13)

For parameter values that are typical for a physiological solution: 26109.0 ⋅=∞n m-3, 182 105.0 −⋅=a m2,

KT 300= , 2.02=OHλ nm, for the chloride ion ( ) 41.1/ 2/1

2≈OHi MM we obtain from (13) 9.13≈η , and,

therefore,

( )( ) ( )iBi

i NTkUN

N⋅−⋅−⋅=

−9.13exp1/exp025.0

1 (14)

The calculated dependence of iN on the reduced potential Tk

U

B

is shown in Figure 9.

0 2 4 6 8 10 12 14 16 18 20 22 240.00

0.04

0.08

0.12

0.16

0.20

0.24

0.28

0.32

σ=[C

/m2 ] ~

qN

i

U/kBT

1 23

4

5

6

Fig. 9. Dependencies of the relative occupancy of the membrane iN and the surface charge density σ on the

reduced potential TkU B/ . Curve 1 corresponds to 26105⋅=∞n , 2 – 2610=∞n , 3 – 25105⋅=∞n , 4 – 2510=∞n , 5 – 24105⋅=∞n , 6 – 2410=∞n .

The maximum surface charge density in Fig. 9 corresponds to the case where all of the cells (Fig. 1b) are full,

i.e. one negative ion is in every cell. It is seen that, even for small values TkU B/ , the relative occupancy of the potential wells with charges is large enough and the charge density of the charges, which are tightly bounded to the membrane, can reach hundredths of coulomb per square meter. It should be noted that in this estimation we did not consider the interaction between the ions bonded to the membrane, which can certainly play an important role for a more accurate calculation of the occupancy of the potential wells.

Let us estimate the charge of the free ions of one sign in the non-quasineutral region in the vicinity of the membrane surface.

2/013.0 mCrnqq Dfree ≈⋅⋅≈ ∞ (15)

For typical values of the potential 155/ −≈TkU B , which correspond to the size of the head 3.02.0~ −hr

nm (Fig.9), a bounded charge sitting on the membrane surface is much larger than the charge associated with unbounded ions.

4. Discussion

In this article, we examined the behavior of the membrane surface with the dipole heads placed in an electrolyte solution. However, despite of this very simple approach, the results are directly related to the processes in biological lipid membranes, such as the membranes of the axon. Because the negative ions inside of the axon differ from those located on the outside (inside the main negative ions are anion groups of macromolecules and phosphates, but chloride ions outside), it is expected that the surface charge density of ions sitting on the membrane is different. Figure 10 shows the potential distribution inside and outside of the

membranecharges.

Fig. 10. membrane

A few wo

parameter

and kU /of the melongitudinthe surfacdetermine

In [24,25molecularnot illustrcalculatio

Conclusio

It is show

1. Tfi

2. Tgr

3. Tph

4. Ash

5. Tbofi

e where the

Potential dise have the sam

ords about th

rs that determ

TkB . Accordinembrane, the cnal displacemce charge dene the surface c

5], the distribr dynamics, hrated. Therefns performed

ons

wn in this pape

The effective birmly bounded

The value of threater than 0.

The surface chhospholipid m

A qualitative shown that the

This work supound with it,ield, interacti

inner and ou

stribution insme (A) and di

he experimen

mine the relat

ngly, by chancurves, simila

ments of the mnsity (Fig. 1charge of the

bution of ionhowever, in thfore, it wouldd by using the

er that:

binding energd to the surfache potential d7 nm the influharge ions aremembrane, whself-consistene density of thpports the ass, which undeing with the

uter sides of

ide and outsifferent (B) su

ntal verificati

tive occupanc

nging the denar to those sh

membrane in 0). Accordinmembrane an

ns near the shose articlesd be interestie molecular dy

gy of the ion wce of the mem

decays exponeuence on the e in relativelyhich prevents

nt theory of thhe bound charumption that rlies the worions, transfe

the membran

side the memurface charge

on of the pro

cy of the mem

nsity of NaCl hown in Fig.8the microwavgly, by channd to compare

surface of a the binding

ing to compaynamics.

with the memmbrane. entially with dions of the ely deep potents them from slhe potential dirges on the me

the ions locark [17]. Therers energy an

ne have the s

mbrane whenes.

oposed mode

mbrane with

in the solutio8, can be reprve field is re

nging the salie with estima

phospholipienergy of theare the corre

mbrane is of o

distance fromlectrolyte by ttial wells, locliding along tistribution neembrane can ated on the erefore, the elend momentu

same (A) and

n the inner a

el. As follow

negative ion

on and measuroduced. It issonant and isinity of the w

ates obtained a

d membranee ions with psponding res

order of sever

m the membrathe membranecalized near tthe surface ofear the membrreach hundre

excitable axonectric compon

um directly to

d different (B

and outer sid

ws from (13),

s iN , are ∝η

uring the surfa shown in [1s linearly depwater, it is pabove.

e was calculaphospholipid hsults of our m

ral kBT. i.e., th

ane, so that at e can be neglethe dipole hef the membranrane is consid

edths of C/m2

n membrane nent of the mo the membr

B) surface

des of the

, the main

∞∝ N/1

ace charge 7] that the

pendent on possible to

ated using heads was model and

he ions are

a distance ected.

eads of the ne. dered. It is . are tightly

microwave rane. This

intr

It should bfor differe

Fig. 11. membrane

Referenc

[1] J. Craielectrostatvesicles. B [2] Mohalipid mem [3] Grego [4] Stuart Chern, 18 [5] Ross FBiophysic [6] V. Pikactivity anNeural En

[7] P. M. Electroma

[8] P. M. Biophysic

[9] N. Lakmembrane

nteraction leadransmembrane

be noted thatent types of b

Membrane de; 1σ , 2σ are

es

ig Franklin, Dtics: synthesiBiophysical J

mmed Awadmbrane. INTE

or Cevc (1990

McLaughlin 8:1I3-36.

F. Flewelling cal Journal 49

kov, X. Arakand plasma meng. 7, 045003

Krasil’nikov.agnetic Field.

Krasil’nikov.cs (Russia), 4

kshminarayane. J. Membran

ds to forced me protein ioni

the proposediological mem

deformations e the inner an

David S. Cafisis and voltagJ. 64, 642-653

d Ali Khalid. CH, Open Sc

0). Membrane

(1989). The

AND Wayne9, February, 5

aki, M. Harrinembrane prop. (2010)

. (1999). Reso Biophysics

. (2001). Effe6 (3), 437.

naiah and K. Mne Biol. 23, 2

mechanical viic channels.

d model of thmbranes, takin

caused by tnd outer surfac

so, Ross F. Fle dependent

3.

(2013). Memcience, Open m

e electrostatics

electrostatic p

e L. Hubbell (531-540.

ngton, S. E. Frperties with lo

onance Intera(Russia), 44 (

ect of the Surf

Murayama. (1279-292.

ibration of the

e potential ding into accou

the microwavce charge den

lewelling, andpartitioning o

mbrane electrominds, http://

s. Biochimica

properties of m

(1986). Hydro

raser, and P. How-power mil

action of Surfa(6), 1041.

face Charge o

1975). Estima

e membrane a

istribution neaunt their chara

ve electric fnsities.

d Wayne L. Hof negative h

ochemistry: e/dx.doi.org/10

a et Biophysic

membrane. A

ophobic ion in

H. Siegel. (20llimeter wave

face-Charged

on the Strain D

ation of surfa

and, as a resu

ar the membracteristics.

field parallel

Hubbell (1993hydrophobic

electrochemic0.5772/55507

ca Acta, 1031

Annu. Rev. Bi

nteraction wi

010). Modulaes in organoty

Lipid Vesicle

Dynamics of

ce charges in

ult, to a redistr

rane can be ex

to the surfa

3). Probes of mion spin labe

cal processes 7

1-3, 311-382.

iophys. Bioph

th membrane

ation of neuroypic cortical s

es with the M

Lipid Membr

n some biolog

ribution of

xtended to

ace of the

membrane els in lipid

in bilayer

hys.

e.

onal lices. J.

Microwave

ranes,

ical

[10] Thomas Heimburg. (2009). Physical Properties of Biological Membranes. arXiv:0902.2454v2 [physics.bio-ph] 16 Feb 2009.

[11] Gouy, G. (1909). (1910). Comt.Rend. 149, 654 (1909), (1910). J.Phys. 4, 9, 457. [12] Chapman, D.L. (1913). Phil.Mag., 6, 25, 475. [13] Stern, O. Z. (1924). Electrochem, 30, 508. [14] A. V. Delgado, F. Gonzalez-Caballero and et al. (2007) Measurement and interpretation of electrokinetic phenomena. IUPACK Pure Appl. Chem., Vol. 77, No. 10, pp. 1753–1805, 2005, DOI: 10.1351/pac200577101753 (J Colloid Interface Sci. 2007 309, 194-224). [15] S M Korobeynikov, A V Melekhov, G G Furin,V P Charalambakos and D P Agoris, (2002), Mechanism of surface charge creation due to image forces. J. Phys. D: Appl. Phys. 35, 1193–1196 [16] Bruce Alberts, Alexander Johnson, Julian Lewis, Marin Raff, Keith Roberts and Peter Walter. (2002). Molecular Biology of the Cell, 4-th edition. New York, Garland Science.

[17] M.N. Shneider, M. Pekker, (2013) “Non-Thermal Mechanism of Weak Microwave Fields Influenc on Nerve Fiber”, Journal of applied physics, 114, 104701 [18] G. Gramse, A. Dols-Perez, M.A. Edwards, L. Fumagalli and G. Gomila. (2013). Nanoscale Measurement of the Dielectric Constant of Supported Lipid Bilayers in Aqueous Solutions with Electrostatic Force Microsco Biophysical Journal, Volume 104, Issue 6, 1257-1262, doi:10.1016/j.bpj.2013.02.011 [19] M. Uematsu and E.U. Franck (1980). Static dielectric constant of water and steam. J. Phys. Chem. Data, Vol. 9, No. 4, 1291. [20] L.D. Landau & E.M. Lifshitz. (1960). Electrodynamics of Continuous Media (Volume 8 of A Course of Theoretical Physics ) Pergamon Press. [21] John F. Nagle and Stephanie Tristram-Nagle. (2000). Structure of lipid bilayers. Biochimica et Biophysica Acta 1469, 159-195. [22] Leonor Saiz and Michael L. Klein (2001). Electrostatic interactions in a neutral model phospholipid bilayer by molecular dynamics simulations. J. Chem. Phys. Vol. 116, N 7, 3052. [23] Alireza Mashaghi, P. Partovi-Azar and et al., (2012). Hydration strongly affects the molecular and electronic structure of membrane phospholipids. The Journal of Chemical Physics 136, 114709. [24] D. Walther, P. Kuzmin, E. Donath. (1996) Brownian dynamics simulation of the lateral distribution of charged membrane components. EUROPEAN BIOPHYSICS JOURNAL WITH BIOPHYSICS LETTERS Volume: 24 Issue: 3 Pages: 125-135. [25] Robert Vácha , Max L. Berkowitz, Pavel Jungwirth. Molecular Model of a Cell Plasma Membrane With an Asymmetric Multicomponent Composition: Water Permeation and Ion Effects, http://dx.doi.org/10.1016/j.bpj.2009.03.010