Medical, Pharma, Engineering, Science, Technology and Business

**Abd-el-Malek MB ^{1,2*} and Amin AM^{1}**

^{1}Department of Engineering Mathematics and Physics, Alexandria University, Alexandria 21544, Egypt

^{2}Department of Mathematics and Actuarial Science, The American University in Cairo, New Cairo 11835, Egypt

- Corresponding Author:
- Abd-el-Malek MB

Department of Engineering Mathematics and Physics

Faculty of Engineering, Alexandria University

Alexandria 21544, Egypt,

**Tel:**01000066496

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

**Received date:** August 10, 2015; **Accepted date:** November 20, 2015; **Published date:** November 25, 2015

**Citation:** Abd-el-Malek MB, Amin AM (2015) Lie Group Analysis for Solving the Problem of Diffusion of Drugs across a Biological Membrane. J Generalized Lie Theory Appl 9:236. doi:10.4172/1736-4337.1000236

**Copyright:** © 2015 Abd-el-Malek MB, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Generalized Lie Theory and Applications

The Lie group method is applied to study the diffusion process of drugs across a biological membrane which tends to partially absorb the drug. For the diffusion coefficient, we considered two cases. The Lie group analysis is based on reducing the number of independent variables by one, and consequently the mathematical model described by nonlinear partial differential equation to, covers the diffusion process with the boundary and initial conditions, and is transformed into an ordinary differential equation with the corresponding conditions. The obtained nonlinear ordinary differential equation is solved numerically using the 4th and 5th Runge Kutta method, and the results are illustrated graphically and in tables too.

Biological membrane; Diffusion of drugs; Concentration of drug; Absorption

The diffusion of drugs across a composite structure such as a biological membrane is one of the interesting types of diffusion studies and is of great importance in the pharmaceutical sciences [1-4]. In 2002, Abd-el-Malek et al. [5] considered the problem of diffusion of a drug through a thin membrane to define the concentration of drug in the membrane in addition to that in the donor and the receiver cells, by applying the group theoretic approach which is a special case of the Lie group analysis and leads to limited results.

Several experiments were carried out to study such diffusion process, one of such experiments was carried out by Hoogervorst et al. [6] and Spacek and Kubin [7].

Spoelstra and Van Wyk [8] described the process as it is modeled by two cells of the same volume, called donor and the receiver cells, which are separated by a thin membrane. In the donor cell a high concentration of a drug in a saline solution is placed and the receiver cell contains only a saline solution. Initially, the membrane is free from the drug, the two cells are stirred continually, and the drug starts to diffuse through the membrane. The concentration of drug in the receiver cell and also in the membrane begins to increase. They solved the problem numerically using the finite-difference method only to evaluate the values of the parameters of the model and to define the concentration in the donor and the receiver cells.

Here we apply the Lie group analysis, which is a powerful and fundamental tool that provides invariant solutions to boundary value problems. It is applicable to both linear and nonlinear differential equations. The mathematical technique we use is one-parameter group transformation. The Lie group analysis expresses the infinitesimals of a group in terms of one or more functions, called infinitesimal functions, each of which depends on independent and dependent variables. The procedure for finding the infinitesimals is then reduced to finding the auxiliary equation. This can be obtained by solving set of partial differential equations that arise as a result of invoking invariance of the partial differential equation and its auxiliary conditions. Hence, the Lie group analysis can be used to solve a wider variety of nonlinear problems.

**Mathematical formulation**

The following assumptions are considered to construct our model:

(1) The rate of absorption at any point in the membrane is proportional to the concentration of the drug at the point,

(2) The x-axis of the coordinate system is taken in the direction of the normal on the membrane,

(3) No diffusion takes place through the edges of the membrane,

(4) The concentration in the membrane is considered as a function of time and of the x-coordinate, and

(5) The rate of transfer of diffusion in the membrane, per unit area, is proportional to the gradient of the concentration across the membrane.

Let *C(x, t)* be the concentration of the drug at a distance *"x"* from one face of the membrane and at an instant of time *"t"* and *P(x, t)* be the function which determines the coefficient of diffusion. The concentrations of drug in the donor and the receiver cells are *D(t)* and *R(t)*, respectively.

Consider a membrane of unit thickness and if *"Q"* is the coefficient of partition, *"q"* is the coefficient which determines the rate of absorption of the drug in the membrane.

The equation which governs the diffusion process may be written as

, (2.1)

with the following conditions.

(i) Initial condition:

(2.2)

(ii) Boundary conditions:

(2.3)

(2.4)

Where

(2.5)

(2.6)

**Lie symmetry group method**

According to Lie’s method [9-11], the infinitesimal generator of the symmetry group admitted by (2.1) is given by

(2.7)

where

*i* = 1, 2, *j* = 1, 2, and *x ^{i}* is independent variable (

(2.8)

in which are infinitesimal functions of the group variables. Then, the corresponding one-parameter Lie group of transformations is given by

If (2.9) is left invariant by the transformation then Equivalently, we can obtain by solving

subject to initial condition

Since the equation which governs the diffusion has at most second – order derivatives, the second prolongation of the generator should be considered in the form:

(2.11)

where

(2.12)

A vector X given by (2.11) is said to be a Lie point symmetry vector field for equations (2.1) with Invariant condition X^{(2)}{2.1} = 0, and solving this system in view of the invariance of the initial and boundary conditions we get the determining equation:

(2.13)

Solving we get

(2.14)

The auxiliary equation will be

(2.15)

According to our initial and boundary condition we found the solution when c_{2} ≠ 0 and all others will be zeros and substitute in (2.15) we get

(2.16)

Solving we get

where *γ* is an arbitrary constant and *F _{1}(x)* and

By substituting (2.17) in (2.1), we get

(2.18)

and if we consider:

*F _{1}(0) = α Q*, and substitute in (2.3) and (2.5) we get

(2.19)

*F _{1}(1) = Q*, and substitute in (2.4) and (2.6) we get

(2.20)

Now consider the following two cases for *F _{2}(x):*

**Case (1): ** *F _{2}(x) = e^{-x}* and substitute in (2.18) we get

(2.21)

With boundary condition

*F _{1}(0) = α Q * (2.22)

*F _{1}(1) = Q * (2.23)

From the results of **Figure 1**, it is clear that the concentration *C(x, t)* is slowly decreasing from the initial border of the membrane to the other.

From **Figure 2**, it is clear that the concentration in the donor cell *D(t)* decreases with time and the concentration in the receiving cell *R(t)* increases time which is compatible with our physical model.

From **Figure 3**, it is clear that the concentration *C(x, t)* in the membrane is increasing with time since our membrane does not contain any drug but at the same time the points near to initial border has higher concentration.

From **Table 1**, it is clear that the concentration *C(x, t)* is slowly decreasing with distance and is increasing with time.

x |
F1(x) |
C (x, 50) |
C (x, 150) |
---|---|---|---|

0.0 | 163.9000 | 0.00204875 | 0.00614625 |

0.1 | 152.0091 | 0.00190011 | 0.00570034 |

0.2 | 140.6348 | 0.00175794 | 0.00527380 |

0.3 | 129.8698 | 0.00162337 | 0.00487012 |

0.4 | 119.8154 | 0.00149769 | 0.00449308 |

0.5 | 110.5853 | 0.00138232 | 0.00414695 |

0.6 | 102.2979 | 0.00127872 | 0.00383617 |

0.7 | 095.1040 | 0.00118880 | 0.00356640 |

0.8 | 089.1688 | 0.00111461 | 0.00334383 |

0.9 | 084.6975 | 0.00105872 | 0.00317616 |

1.0 | 081.9500 | 0.00102438 | 0.00307312 |

**Table 1:*** F _{1}(x), C (x, 50) and C (x, 150)* for

**Case (2):** , and substitute in (2.18) we get

(2.24)

with boundary condition

*F _{1}(0) = α Q *(2.25)

*F _{1}(1) = Q* (2.26)

From **Figure 4**, we notice that the concentration *C(x, t)* is decreasing from the initial border of the membrane to the other border.

In **Figure 5**, we studied effect of time on the concentration; we found that firstly the concentration *C(x, t)* in the membrane increased with time since our membrane does not contain any drug, but at the same time the points near to the initial border has higher concentration.

**Table 2**, illustrates the concentration *C(x, t)* is slowly decreasing with distance and is increasing with time.

x |
F1(x) |
C (x, 50) |
C (x, 150) |
---|---|---|---|

0.0 | 163.9000 | 0.00204875 | 0.00614625 |

0.1 | 151.5883 | 0.00189485 | 0.00568456 |

0.2 | 140.5958 | 0.00175745 | 0.00527234 |

0.3 | 130.6503 | 0.00163313 | 0.00489939 |

0.4 | 121.5661 | 0.00151958 | 0.00455873 |

0.5 | 113.2246 | 0.00141531 | 0.00424592 |

0.6 | 105.5353 | 0.00131919 | 0.00395757 |

0.7 | 098.4788 | 0.00123098 | 0.00369296 |

0.8 | 092.0940 | 0.00115118 | 0.00345352 |

0.9 | 086.5033 | 0.00108129 | 0.00324387 |

1.0 | 081.9500 | 0.00102438 | 0.00307312 |

**Table 2 :*** F _{1}(x), C (x, 50) and C (x, 150)* for

The relative difference between two cases:

We find *R.D*% in three cases typical as below

**Table 3** illustrates that, the small relative variation in concentration for the two cases is explained in **Tables 1** and **2**, therefore the effect of the diffusion coefficient in Eq. (2.1) very small (**Figure 6**).

x |
0 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 | 0.6 | 0.7 | 0.8 | 0.9 | 1 |

R.D.% |
0 | - 0.28 | - 0.03 | 0.6 | 1.44 | 2.33 | 3.07 | 3.43 | 3.18 | 2.09 | 0 |

**Table 3:** Illustrates *R.D%.*

The value of R.D.% increases gradually against *x* and reached to its maximum value at *x* = 0.7, then starts to decrease again and reaches zero at *x* = 1.

Non-linear partial differential equation, has been solved using Lie group analysis, and *C(x, t) *and *P(x, t)* have been evaluated analytically using the parameters of an experimental Study [8]. The results obtained map the state of the membrane in a way similar to the one obtained by Spoelstra and van Wyk [8], using a finite-difference method.

The other transformations, like spiral transform and linear transform, do not work for leaving the differential equation as well as the initial and boundary conditions invariant.

The concentration proportional linear with time and *P(x, t)* which describes the coefficient of diffusion is inversely proportional with time. Value of *R.D.*% depends on the method of numerical solution to equation (2.18) and the value of *F _{1}(x).*

- Steffi H, Claus-Michael L, Ulrich FS (2013) Improved input parameters for diffusion models of skin absorption. Advanced Drug Delivery Reviews 65: 251-264.
- Mehdi A, Maryam K, Monireh A (2006) The study of drug permeation through natural membranes. International Journal of Pharmaceutics 327: 6-11.
- Francesca S, Francesco C, Annalisa A, Silvia F, Paola M (2012) Regenerated keratin membrane to match the
*in vitro*drug diffusion through human epidermis. Results in Pharma Sciences 2:72-78. - Michele C,Cesare C (2009) The memory formalism in the diffusion of drugs through skin membrane. Journal of Physics D: Applied Physics: 42.
- Abd-el-Malek MB, Kassem MM,Meky MLM (2002) Group theoretic approach for solving the problem of diffusion of a drug through a thin membrane. J. Comput. Appl. Math. 140:1-11.
- Hoogervorst CJP, Van Dijk JAPP,Smit JAM (1977) Nonstationary diffusion through membranes. 1. Transient diffusion through a membraneseparating two unequal volumes of well stirred solution. J. Phys. Chem. 82: 1311-1318.
- Spacek P,Kubin M (1967) Diffusion in gels. J. Polymer Sci. C 16: 705-714.
- Spoelstra J, Van Wyk DJ (1987) A method of solution for a non-linear diffusionmodel and forcomputing the parameters in the model. J. Comput. Appl.Math. 20: 379-385.
- Bluman GW, Kumei S (1989) Symmetries and Differential Equations.Applied Mathematical Sciences: 154.
- Olver PJ (1986) Applications of Lie Groups to Differential Equations. Springer.
- Ovsiannikov LV (1982) Group Analysis of Differential Equations. American Mathematical Society: 8.

Select your language of interest to view the total content in your interested language

- Adomian Decomposition Method
- Algebra
- Algebraic Geometry
- Analytical Geometry
- Applied Mathematics
- Axioms
- Balance Law
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Biometrics
- Biostatistics methods
- Clinical Trail
- Combinatorics
- Complex Analysis
- Computational Model
- Convection Diffusion Equations
- Cross-Covariance and Cross-Correlation
- Deformations Theory
- Differential Equations
- Differential Transform Method
- Fourier Analysis
- Fuzzy Boundary Value
- Fuzzy Environments
- Fuzzy Quasi-Metric Space
- Genetic Linkage
- Geometry
- Hamilton Mechanics
- Harmonic Analysis
- Homological Algebra
- Homotopical Algebra
- Hypothesis Testing
- Integrated Analysis
- Integration
- Large-scale Survey Data
- Latin Squares
- Lie Algebra
- Lie Superalgebra
- Lie Theory
- Lie Triple Systems
- Loop Algebra
- Matrix
- Microarray Studies
- Mixed Initial-boundary Value
- Molecular Modelling
- Multivariate-Normal Model
- Noether's theorem
- Non rigid Image Registration
- Nonlinear Differential Equations
- Number Theory
- Numerical Solutions
- Operad Theory
- Physical Mathematics
- Quantum Group
- Quantum Mechanics
- Quantum electrodynamics
- Quasi-Group
- Quasilinear Hyperbolic Systems
- Regressions
- Relativity
- Representation theory
- Riemannian Geometry
- Robust Method
- Semi Analytical-Solution
- Sensitivity Analysis
- Smooth Complexities
- Soft biometrics
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Super Algebras
- Symmetric Spaces
- Theoretical Physics
- Theory of Mathematical Modeling
- Three Dimensional Steady State
- Topologies
- Topology
- mirror symmetry
- vector bundle

- 6th International Conference on
**Biostatistics**and**Bioinformatics**

Nov 13-14, 2017, Atlanta, USA

- Total views:
**8273** - [From(publication date):

December-2015 - Aug 22, 2017] - Breakdown by view type
- HTML page views :
**8086** - PDF downloads :
**187**

Peer Reviewed Journals

International Conferences 2017-18