Influence of initial molecular substance on the diffusion flux across cell membranes

The influence of initial placement of molecular or ion substance is investigated on the diffusion fluxes across the cell membrane. The diffusion fluxes and recovery curves are obtained by considering both the singlespot and double-spot concentrations inside the cell membrane. The results show that the additional concentration inside the membrane reduces the net fluxes at the cell interior as well as the exterior. In addition, it is found that the change in diffusion flux at the two outer walls of the membrane by the two-spot concentrations in the cell membrane is weaker than that of the single-spot concentration at the center. The variation of the influence of initial locations of the molecular concentrations inside the cell membrane on the diffusion fluxes is also discussed. This result can be applied to the diffusion process in avascular collagenous tissues.


INTRODUCTION
A diffusion process has been one of the most interesting processes in physical sciences since the diffusion provides physical insights into the flux variation in various scientific areas, such as biology, chemistry, mathematics, and physics.In biological fields, the diffusion process is used in animal dispersal, bacteria motion, cell movement, chemical diffusion such as drug release in membrane and tissues, and plant electrophysiology [1][2][3][4][5][6][7][8].Especially, the theory of transport into or out of a cell by the intra-cellular diffusion has received considerable attention since this process has wide applications in biological sciences.
It has been shown that the diffusion in the human body is noteworthy over small distance scales on the order of 1 100 m µ − and the diffusion plays an important role in the metabolic activity [9].It has been also shown that the diffusion phenomena in cells would be strongly biased by the influence of environments.Hence, it is expected that the diffusion in a cell membrane involving double lipid layers would be significantly influenced by the location of the initial location of a molecular solute.Thus, in this paper we investigate the influence of initial placement of molecular substance on the diffusion fluxes at outer walls of the cell membrane by using the theoretical analysis with the boundary conditions since the diffusion processes are significant for the characterization of cellular and physiological processes and for the operation of biomedical machines [10], such as the heart-lung bypass device, kidney dialysis device, and membrane oxygenator.The diffusion fluxes and recovery curves in the cell membrane are obtained by considering the single-spot and double-spot concentrations, respectively, inside the cell membrane.The influence of initial location of the concentration inside the cell membrane on the net diffusion fluxes at the outer walls of the cell membrane is also discussed.
This paper is composed as follows.In Section 2, the diffusion equation in the cell membrane with the detailed initial and boundary conditions and the Danckwert's method for the diffusion-reaction equation are discussed.In Section 3, we obtain the solutions of the diffusion equation for the single-spot and double-spot concentrations inside the cell membrane and the recovery curves for both single-spot and double-spot concentrations.In Section 4, we obtain the net flues at the wall boundary for the single-spot and double-spot concentrations inside the cell membrane.In Section 5, the influence of initial placement of molecular substance on the diffusion fluxes OPEN ACCESS at outer walls of the cell membrane is discussed.In addition, we discuss the variations of the concentrations and recovery curves inside the cell membrane.Finally, the conclusions are given in Section 6.

DIFFUSION IN THE CELL MEMBRANE
From Fick's first and second equations for the flux ( ) , J x t , the diffusion equation [5] for the concentration ( ) where x is the position in one-dimension, t is the time, and D is the diffusion coefficient.It is known that a typical value of the diffusion coefficient for ions such as K + , Cl − , and Na + is about 6 2 2.5 10 cm s − × [11].For the sake of the simplicity, we consider the one-dimensional diffusion problem for the investigation of the diffusion fluxes across cell membranes.In Equation ( 1), we retain only the diffusive flux term, but neglect the electrical drift effect since the influence of initial placement of molecular or ion substance on the diffusion fluxes at outer walls of the cell membrane is the main purpose of this work.Hence, we assume that the cell membrane has been localized in the domain of 0 x L < < and the molecular concentrations at the domain boundaries as the boundary conditions are ( ) Two cases for the initial intracellular con- centration as the initial conditions shall be considered: for the single-spot concentration case and ( ) ( ) ( ) ( ) ( ) for the double-spot concentration case, where ( ) x δ is the Dirac delta function and a and b are initial locations of the molecular concentrations.If the reaction phenomenon is involved in the diffusion process, the righthand side of Equation (1) would have an extra term proportional to the concentration: , where k is a constant.In this case, the diffusion-reaction equation in the cell membrane is given by where ( ) , kC x t represents the rate of removal of dif- fusing substance.It has been also shown that the diffusion-reaction equation is useful for the mathematical modeling of glioma growth and also for the investigation on the spreading of brain tumors [12] and ( ) , kC x t term corresponds to the death of cells in the tumor growth model.It has been well known that the general solution of Equation (2) can also be obtained the following transformation based on Danckwert's method [13,14]: where ( ) 1 , C x t is the solution of Equation (1) with the same boundaries.As a result, the standard diffusion equation, Equation (1), is used to examine the change in diffusion fluxes at outer walls of the cell membrane depended the initial molecular substance.In the following Section III, the solutions, ( ) , C x t for the single-spot and for the double-spot concentrations inside the cell membrane are discussed.

SINGLE-AND DOUBLE-SPOT CONCENTRATIONS
Using the separation of variables method, the solution of the diffusion equation for the domain 0 x L < < would be represented by where ( ) is the separation constant determined by the boundary condition at x L = , n A are coefficients to be determined by the initial condition at 0 t = , and the first two terms ( )( ) represents the solution when the separation constant is zero.Since we are interested in the effect of the initial molecular concentration ( ) ,0 C x inside the cell membrane, the coefficients n A are determined by ( ) ( ) where ( ) for the single-spot concentration and for the double-spot concentrations with 0 a L < < and 0 b L < < .The consideration of the initial concentration inside the cell membrane would be also quite useful to understand the behavior of the tumor growth since it has been known that the spatial arrangement effect is crucial for the spatial spread of the cancerous cells and also for the appearance of the tumor in the human body [12].Hence, the general solutions for the diffusion equations with the single-and double-spot concentrations are, respectively, found to be the following forms: Once the concentration function ( ) is obtained, the measurement of the concentration fluctuations would be determined by the signal correlation function where ( ) ( ) ( ) , A is a factor related to the experimental detection, and ( ) x is the intensity of the laser beam [15].A detailed investigation on the correlation spectroscopy for the influence of initial placement of molecular or ion substance on the concentration fluctuations will be treated elsewhere.The recovery curves for the single-spot

NET DIFFUSION FLUX IN THE CELL MEMBRANE
The diffusion fluxes at the wall boundaries 0 x = , i.e., 0 J , and x L = , i.e., L J , for the case of the single-spot concentration are, respectively, given by ( ) If we retain on the first term in the series, the net flux ( ) for the single-spot concentration is then found to be ( ) ( ) As it is seen in Equation ( 12), the net flux in the cell membrane is completely determined by the location ( ) sin a L π and amount a C of the initial molecular concentration inside the cell membrane.Likewise, the diffusion fluxes at the wall boundaries 0 x = and x L = for the case of the double-spot concentrations are, respectively, obtained by ( ) If we also retain on the first term as we did in Equation (12), the net flux ( ) for the twospot concentrations is then obtained as ( ) As shown in Equation ( 15), the net flux in the cell membrane is also completely controlled by the locations ( ) Since the ion diffusion in biological membranes is one of the main transport processes [16], the results of S J ∆ and D J ∆ , Equations ( 12) and ( 15), in this model analy- sis would provide useful information on the electric current and the membrane potential.Recently, the diffusional anisotropy has been investigated in avascular collagenous tissues since the diffusion in articular cartilage is the main molecular transport process [17].Hence, the effects of the initial placement of molecular substance on the two-dimensional anisotropic diffusion process will be treated elsewhere.

INFLUENCE OF INITIAL ION PLACEMENT
In this section, the graph plotting for the 2D and 3D plots are obtained by the technical computing software Mathematica. .As it is seen, it is found that two non-flux positions for S J ∆ can be existed when the initial molecular concentration inside the cell membrane is equal to the larger one between two molecular concentrations at the domain boundaries.We have also found that that two non-flux positions are placed near the wall boundaries.In addition, it is found that the net flux decreases with an increase of the diffusion time t.  are disappeared inside the cell membrane when the initial molecular concentration inside the cell membrane is smaller than two molecular concentrations at the wall boundaries.In addition, the net flux S J ∆ is found to be quite small at the center of the cell membrane.
Figure 3 represents the scaled recovery curve ( ) for the initial single-spot concentration  From this figure, it is found that the recovery curve S N would have the maximum value when the initial single-spot concentration is placed near the right wall boundary and decreases with an increase of the diffusion time t .However, we have found that the recovery curve S N would have the minimum value when the initial single-spot concentration is placed near the left wall boundary and increases with an increase of the diffusion time t .It is also found that that the recovery curve S N decreases monotonically with an increase of the diffusion time when the initial single-spot concentration is placed at the center of the cell membrane.
Figure 4 shows the scaled recovery curve S N for the initial single-spot concentration case as a function of the scaled diffusion time t for various values of the initial location a when As shown in this figure, it is found that the recovery curve S N would have the minimum value when the initial single-spot concentration is placed near the left wall boundary at the center of the membrane and increases with an increase of the diffusion time t .However, we have found that the recovery curve S N has the maximum value when the initial single-spot concentration is placed near the right wall boundary and decreases with an increase of the diffusion time t .for the initial double-spot concentration case is found to be also quite small at the center of the cell membrane.  .From this figure, we have found that the recovery curve D N has no maximum and increases and finally saturates with increasing diffusion time t when the two locations a and b are away from the center of the membrane.However, it is found that the recovery curve D N has the maximum value at the small diffusion time and saturates with an increase of the diffusion time t .

CONCLUSIONS
In this work, we inspected the influence of initial placement of molecular or ion substance on the diffusion fluxes at outer walls of the cell membrane.The recovery curves in the cell membrane were obtained for the single-  spot and double-spot concentrations, respectively.The net diffusion fluxes were obtained by considering both single-spot and double-spot concentrations across the cell membrane.For the initial single-spot concentration case, it was found that two non-flux positions would exist near the wall boundaries if the initial molecular concentration inside the cell membrane were equal to the larger one between the two molecular concentrations at the domain boundaries.It is also found that the two non-flux positions were placed near the center of the cell membrane when the initial molecular concentration inside the cell membrane was equal to the smaller one between the two molecular concentrations at the wall boundaries.In addition, we have found that the non-flux positions disappeared inside the cell membrane when the initial molecular concentration was smaller than two molecular concentrations at the wall boundaries.The molecular concentration S C could be increased more than 0 C when the initial single-spot concentration was placed near the left wall boundary and . However, we have found that the molecular concentration monotonically decreases with an increase of the position when the initial single-spot concentration is placed at the center or near the right wall boundary.Moreover, we have found that the molecular concentration can be smaller than L C when the initial single-spot concentration was placed near the left wall boundary.Thus, the recovery curve would have the maximum vale when .It is also found that the net flux had the maximum value at the center of the cell membrane.We found two non-flux positions could also be existed near the center of the cell membrane when .However, we have found that the molecular concentration decreased with an increase of the distance x regardless of the locations of the ini- tial double-spot concentrations when . It is found that the recovery curves had maximum values at very short diffusion times and the recovery curve increased as the two positions were getting close to the center of the cell membrane when . It is found that the recovery curve at first had no maximum, but increased and finally saturated with increasing diffusion time when the two initial locations were away from the center of the membrane.It had the maximum value at the small diffusion time and saturated with an increase of the diffusion time when .From this work, we have found that the influence of initial locations and amounts of the molecular substance plays a significant role on the diffusion fluxes across the cell membrane.Since the results are obtained in the analytic expressions with the appropriate physical conditions, the equations in this work would be quite reliable for understanding the diffusion process in biological systems.Hence, the results of this work would be useful for understanding the characteristics and properties of the intracellular diffusion process as well as the spread of the cancerous cell and appearance of the tumor.The results in this work can also be applied to the enhancement of the diffusion process in articular cartilage.Moreover, these results would be applied to the biomedical machines related to the biological diffusion processes including the initial ion concentrations inside the wall boundaries.
are then, respectively, obtained as follows:( a C and b C of the initial molecular concentrations inside the cell membrane.

Figure 1
represents the net flux S J ∆ for the initial single-spot concentration case as a function of the scaled position ( )

Figure 1 .
Figure 1.The net flux Δ S J for the initial single-spot concentration case as a function of the scaled position ( ) a a L ≡

Figure 2
Figure 2 shows the net flux

Figure 2 .
Figure 2. The net flux Δ S J for the initial single-spot concentration case as a function of the scaled position a when 0 = 2 C , = 1 L C , and = 0.8 a C , i.e.,

Figure 3 .
Figure 3.The scaled recovery curve

Figure 5
represents the surface plot of the net flux D J ∆ for the initial double-spot concentration case as a function of the scaled positions() can see from this figure, it is found that two non-flux posi-

Figure 4 .
Figure 4.The scaled recovery curve S N for the initial singlespot concentration case as a function of the scaled diffusion time t when 0 = 2 C , = 1 L C , and = 1 a C , i.e.,

Figure 6 Figure 5 .
Figure 5.The surface plot of the net flux Δ D J for the initial double-spot concentration case as a function of the scaled positions ( ) a a L ≡ and ( ) b b L ≡ when = 0.1 t , 0 = 2 C ,

Figure 6 .
Figure 6.The surface plot of the net flux Δ D J for the initial double-spot concentration case as a function of the scaled positions a and b when = 0.1 t , 0 = 2 C , = 1 L C , = 0.4 a C ,

Figure 7
shows the scaled recovery curve D N for the initial double-spot concentration case as a function of the scaled diffusion time t for various values of the scaled positions a and b when is seen in this figure, the recovery curves D N have maximum values at very short diffusion times t .It is also found that the maximum time of the recovery curve D N has been increased as the two positions are getting close to the center of the cell membrane.

Figure 8
represents the scaled recovery curve D N for the initial double-spot concentration case as a function of the scaled diffusion time t for various values of the scaled positions a and b when

Figure 7 .
Figure 7.The scaled recovery curve D N for the initial double-spot concentration case as a function of the scaled diffusion time t when = 0.1 t , 0 = 2 C , = 1 L C , = 1 a C , and

Figure 8 .
Figure 8.The scaled recovery curve D N for the initial double-spot concentration case as a function of the scaled diffusion time t when = 0.1 t , 0 = 2 C , = 1 L C , = 0.5 a C , and the initial single-spot concentration was placed near the right wall boundary and the recovery curve decreases with an increase of the diffusion time.The minimum value would be when single-spot concentration was placed near the left wall boundary at the center of the membrane and increases with an increase of the diffusion time.For the initial double-spot concentration case, it was found that two non-flux positions were located near the wall boundaries when flux was still quite small at the center of the cell membrane.However, we have found that the non-flux positions were disappeared inside the cell membrane when the net flux for the initial double-spot concentration case was found to be also quite small at the center of the cell membrane.We have found that the molecular concentration D C had the maximum value when two positions were near the center of the cell membrane when