Fractional derivatives in the transport of drugs across biological materials and human skin
Introduction
The penetration of drugs through intact skin can be considered as a process of molecular diffusion across a composite multilayer membrane, whose principal barrier to diffusion is localized within the stratum corneum (SC). As a matter of fact, skin is a multilayered membrane composed by the stratum corneum that represents the outermost layer and by the epidermis and the dermis, which of them characterized by different structural and physico-chemical properties. While epidermis and dermis exert a barrier towards lipophilic compounds, stratum corneum, consisting of several layers of hydrated and keratinized cells, suspended in an extracellular lipid matrix, acts as a rather impermeable layer to most of drugs and other hydrophilic substances.
This biological system, because of its structural organization, is particularly suitable to be described, as far as the drug permeation is concerned, by a diffusion coefficient able to take into account the intrinsic heterogeneity of the system. This can be conveniently taken into account by a diffusion equation with fractional derivative that, because of its intrinsic nature, is constructed by adding to the initial value the successive weighted increments over time.
We have recently proposed that in highly heterogeneous systems, and particularly in the case of the skin [1], the overall transport process could be described by the second Fick law modified by introducing a memory formalism (diffusion with memory) [2], [3], [4], [5]. This approach has been employed to describe the concentration profile inside a biological membrane when a sudden change of the concentration in the medium bathing one of its face is applied.
The diffusion with memory based on a fractional derivative approach has been already successfully used in different fields, ranging from electromagnetism [6] to fractal media [7]. In particular, fractional derivatives have been employed to model the rheological properties of solids [8], in the constitutive equations of polarizable media in the time domain [6] and in chaotic systems [9] and, more recently, in the dynamics of glass-forming materials [10] and in demography and economy [11], [12].
In this note, we apply the diffusion with memory formalism to the penetration of low molecular weight solutes through a composite membrane modeled as a bilayer structure consisting of a layer of thickness adjoining a semi-infinite medium, each of them characterized by different diffusion parameters.
The introduction of memory into the classical diffusion equations (Fick’s equations) is carried out by introducing, besides the classical Caputo definition of the fractional derivative, a new definition of the fractional derivative recently proposed by Caputo and Fabrizio [13] with a different kernel, removing the singularity present in the classical Caputo fractional derivative.
We compare the model predictions based on these two different definitions with some experimental results of drug penetration into skin in vitro we have previously analyzed in the framework of the classical Caputo definition [6]. Moreover, this analysis was extended to other diffusion processes, such as permeation of water across the human skin and the diffusion of an antiviral agent through the skin of male rats. While both the diffusion models are able to give a satisfactory agreement with the experimental data, the model based on the new definition of fractional derivative proposed by Caputo and Fabrizio [13], on the basis of the simple residual analysis furnishes a little bit better agreement.
Section snippets
Definition of fractional derivatives
The usual Caputo fractional derivative of order is defined as [14], [15] with and the Euler Gamma function. It is worth noting how the memory function captures the past. What the fractional derivative memory functions is remembering is the past values of the function, which implies that the function is constructed by adding to the initial value the successive weighted increments over time. The increments per unit time are represented by
The diffusion model: extension of the classical Fickian model
The classical constitutive equations involved in the study of diffusion, based on the Fick’s laws read where is the flux and the solute concentration at position and time . is the usual diffusion coefficient.
By introducing the memory formalism through the fractional derivatives, the above stated equations become where now is the diffusion coefficient (with memory). Eqs. (7), (8)
Comparison with experiments
The first example we are dealing with is shown in Fig. 6, where the concentration of a model compound, 4-cyanophenol, dissolved in a saturated aqueous solution is plotted as a function of the position within the stratum corneum, after 15 min exposure (experimental data from Ref. [16]). The concentration of the chemical was obtained using quantitative attenuated total-reflectance Fourier transform infrared [AFR-FIR] spectroscopy [17] in sequentially tape-stripped layers of the stratum corneum
Conclusions
The comparison of the results previously obtained with the usual Caputo fractional derivative with those obtained with the new definition (Caputo and Fabrizio definition) clearly shows that, at least in some cases, the new definition gives better results. However, it may be possible the different phenomena have a memory represented by different mathematical memory formalisms. The two different fractional derivatives now available, together with more experimental data analysis, may help in
References (20)
- et al.
Memory formalism in the passive diffusion across highly heterogeneous systems
J. Membr. Sci.
(2005) - et al.
Time and spatial concentration profile inside a membrane by means of a memory formalism
Physica A
(2008) - et al.
Diffusion with memory in two cases of biological interestm
J. Theoret. Biol.
(2008) - et al.
Introduction to transfer motion in fractal media: the geometry of kinetics
Solid State Ion.
(1983) - et al.
Fractional calculus and continuous-time finance: Ii. The waiting time distribution
Physica A
(2000) - et al.
Piroxicam delivery into human stratum corneum in vivo: iontophoresis versus passive diffusion
J. Control. Release
(2001) - et al.
Acyclovir permeation through rat skin: mathematical modelling and in vitro experiments
Int. J. Pharm.
(2003) - et al.
Diffusion modelling of percutaneous adsorption kinetics: 4. Effect of slow equilibrium in the stratum corneum on adsorption and desorption kinetics
J. Pharm. Sci.
(2009) - M. Caputo, C. Cametti, Permeation of water through the stratum corneum. a mathematical model based on diffusion with...
- et al.
The memory formalism in the diffusion of drugs through skin membrane
J. Phys. D: Appl. Phys.
(2009)
Cited by (12)
Unsteady rotating electroosmotic flow with time-fractional Caputo-Fabrizio derivative
2020, Applied Mathematics LettersCitation Excerpt :Numerical analysis for EO flow of fractional Maxwell fluids was conducted by Yang et al. [24]. In addition, for the transport problem of drugs through human skin, Caputo and Cametti [25] compared theoretical models of classical Fick’s law, usual Caputo fractional model and new definitions of Caputo–Fabrizio derivative [26] with experimental data. They found a closer agreement is obtained in the case of model based upon the new definition of fractional derivative.
Modeling non-Darcian flow and solute transport in porous media with the Caputo–Fabrizio derivative
2019, Applied Mathematical ModellingConformable derivative approach to anomalous diffusion
2018, Physica A: Statistical Mechanics and its ApplicationsCitation Excerpt :The fractional derivative generalized the classical derivative of integer order to a differential operator of arbitrary order. Due to its non-local behavior which yields the memory effects and history dependence, in the past decades, the fractional derivative has gained remarkable applications in various areas of science and engineering [11,12] such as the time-dependent characteristics of rocks [13,14] and composite materials [15], fluid mechanics [16], nanotechnology [17], control system [18], electrical circuits [19], thermo-elasticity [20], nuclear reactor dynamics [21] and biology [22]. Moreover, the fractional derivative approach to anomalous diffusion can be found in many literatures. [23]
In vitro and in vivo characterization of pharmaceutical topical nanocarriers containing anticancer drugs for skin cancer treatment
2018, Lipid Nanocarriers for Drug TargetingA NEW APPROACH TO FRACTIONAL DIFFERENTIAL EQUATIONS
2023, Thermal Science