Research Article  Open Access
S. Băcăiţă, C. Urîtu, M. Popa, A. Uliniuc, C. Peptu, M. Agop, "Drug Release Kinetics from Polymer Matrix through Fractal Approximation of Motion", Smart Materials Research, vol. 2012, Article ID 264609, 8 pages, 2012. https://doi.org/10.1155/2012/264609
Drug Release Kinetics from Polymer Matrix through Fractal Approximation of Motion
Abstract
The present paper analyzes the process of drug release from polymer matrix. This process has been considered as fractal polymer process. Since complexity of physical processes is replaced by fractality, the paper studies the process through fractal approach. In drug dynamics, fractal “diffusion” equation can be obtained through fractal approximation of motion. All experimental release curves have been best demonstrated by Weibull relation (which was, in its turn, also demonstrated). Weibull parameters are related to the fractal dimension of drug release kinetics from a polymer matrix. Such a dimension can characterize and measure the complexity of the system. In the abovementioned context, some experimental results of our researchers are presented and analyzed by comparing them with Peppas relation, a basic law in the description of drug release kinetics. Consequently, experimental data for Weibull relation are better correlated with certain resulting factors. At the same time, some conclusions regarding the phenomena involved in the process are considered as being based on the approach.
1. Introduction
The adequate description of real, natural, and artificial objects is restricted by the mere use of Euclidean geometry, that is, the description of integerdimensioned objects. In such perspective, this happens because numerous objects with noninteger dimensions, such as plants, galaxies, population patterns, and crystal growth, are left beyond analytical purpose. The characteristics of such objects can be described by means of fractal geometry [1, 2]. Natural and synthetic polymers considered as fractal objects are also included in the above category. Their main structural unit, the macromolecular coil, is known to be a fractal with typical fractal behavior [3–8].
The main idea of the paper focuses on this type of processes in which polymer fractality is responsible for drug release from various polymer matrixes. In what follows, we shall produce arguments to sustain the above statement.
The experimental drug release kinetics indicate that drug loaded polymer matrix structures are thermodynamically unstable and evolve towards equilibrium. Specific parameters for each structure (drug type, incorporated drug dose(s), types and amounts of excipients, preparation technique, environmental conditions during drug release, and geometry and dimensions of drug delivery system) will lead to different evolutive trajectories, considered as consequences of internal collective processes. Thermodynamically nonbalanced processes generate the formation of fractal structures [5, 9]. That is why fractal structures appear in drug release environment, too.
Drug release phenomena (water penetration into the device, drug and excipient dissolution, phase transitions, drug and/or polymer degradation, polymer swelling, physical drugexcipient interactions, and chemical reactions of drug, excipients, and/or water) are complex processes. Nevertheless, mathematical expressions pharmaceutics uses for the description of drug release kinetics are rather simple, namely Higuchi, RitgerPeppas, PeppasSahlin, and Alfrey power laws [10].The structure whose dynamics is ruled by power laws tends to critically organize itself [11], as fractals do.
Thus, the question of whether or not fractal analysis should be used in the description of the structure and evolution of such a structure, that is, drugloaded polymer matrix, is not a matter of scientific choice, but a consequence of a correct approach.
This is the reason why we will analyze this process considering that fractality can replace the complexity of physical processes. Moreover, it will be no longer necessary to use the entire classical “arsenal” of quantities of standard physics (differentiable physics). Physical systems will behave as special “fluids” that lack interaction. Thereby, we introduced fractal approximation of motion in the study of the complex dynamics of physical systems, taking into account that drug trajectories are continuous, but nondifferentiable curves, called fractal curves. Interactions are dealt with in the abovementioned manner by a physical model, which is the scale relativity theory (SRT) [12, 13]. Accordingly, we will obtain a fractal “diffusion” equation and deduce Weibull law that, in comparison with power laws, better describes all experimental release curves [14], thus eliminating all possible critical opinions regarding the absence of kinetics in the use of Weibull’s law or the nonphysical nature of its parameters [15].
This paper is structured as follows: Section 2 presents general notions on the fractal approximation of motion through scale relativity theory. Fractal “diffusion” equations, both for Fick type “diffusion” and for anomalous type “diffusion,” as well as Weibull relation are exhibited in Section 3. Experimental results of our theoretical model are presented in Section 4. Section 5 comprises the conclusions.
2. Fractal Approximation of Motion Applied in Drug Release Process
Drugloaded polymer matrix can be included in the category of complex dynamic systems which display nonlinear behavior, selfsimilarity, and selfstructure. At the same time, these systems intensely fluctuate at all possible scales [9–21]. Since fractality appears to be the universal feature of such systems, it is essential to construct fractal physics [13, 18] when describing such evolution.
2.1. Nondifferentiability Consequences in Drug Release Process
Drug release process is a very complex phenomenon, with many significant variables in the evolution of the system. That is why we assume that this type of process takes place on continuous, but nondifferentiable curves (fractal curves). Nondifferentiability implies the following [13, 18].(i)A continuous and nondifferentiable curve (or almost nowhere differentiable) is explicitly scale dependent, and its length tends to infinity, when the scale interval tends to zero. In other words, a continuous and nondifferentiable space is a fractal one, in that general meaning Mandelbrot gave to the concept [11].(ii)There is an infinity of fractal curves (geodesics) that connects its points (or starting from any point) and this is applicable for all scales.(iii)Local differential time reflection invariance breaking: the time derivative of relative variation for the released drug can be written twice:
Both definitions are equivalent for the differentiable case. As regards nondifferentiable situations, these definitions are no longer valid, since limits cannot anymore be defined. Within fractal theory, physics is related to the behavior of the function during the “zoom” operation on the time resolution δt, here identified with the differential element dt (substitution principle), which is considered as an independent variable. The standard concentration field (t) is therefore replaced by fractal concentration field (t,dt), explicitly dependent on the time resolution interval. The characteristic derivative of this interval is undefined only at that unnoticed limit which is described by a fractal function. Consequently, this will lead to the definition of the two derivatives of the fractal concentration field as explicit functions of and dt variables: “+” corresponds to the forward process and “−” to the backward one.(iv)The differential of the coordinates, , can be decomposed as follows: where is the “classical part,” and is the “fractal part,” depending on resolution. The differential of the “fractal part” components , , satisfies the relation (the fractal equation) [22–24] where are constant coefficients, and is the fractal dimension. For fractal dimension, we can use any definition (Kolmogorov, Hausdorff [1, 13, 25–27], etc.), but once a certain definition is admitted, it should be kept up to the end of analysis.(v)The local differential time reflection invariance is recovered by combining the two derivatives, and , in the complex operator When this operator is applied to the “position vector,” a complex speed yields with
The real part of the complex speed represents the standard classical speed, which is independent of resolution, while the imaginary part is a new quantity arising from fractality, which is resolution dependent.(vi)The average values of the quantities must be considered in the sense of a generalized statistical fluid. Particularly, the average of is with
2.2. Covariant Total Derivative in Drug Release Process
Let us now admit that the curve (continuous, but nondifferentiable), describing the drug release, is immersed in a threedimensional space, and that of components is the position vector of a point on the curve. Let us also consider the fractal concentration field and expand its total differential up to second order
The relations (10) are manifold valid in any point of the space and, implicitly, for the points on the fractal curve which we have selected in relations (10). Hence, the forward and backward average values of this relation take the form
We formulate the following stipulation [22–24]: the mean value of the fractal concentration field and its derivatives coincide with themselves and the differentials and are independent. Therefore, the average of their products coincides with the product of the averages. Thus, using (3) with the property (9), (11) becomes
Even the average value of the fractal coordinate is null (see (9)), whilst for the higher order of the fractal coordinate average, the situation can be different. Firstly, let us focus on the mean . If , this average is zero due the independence of and . So, using (4), we can write
Then, if we divide by and neglect the terms which contain differential factors (for details on the method, see [16, 17]), (12) may be written under the form
These relations also allow us to define the operator
Under these circumstances, let us calculate . Taking into account (15), (5), and (6), we obtain
This relation also allows the definition of the fractal operator
Particularly, by choosing we obtain the following results:(i)the relation (13) becomes (ii)the fractal operator (17) takes the form [22–24](iii)Nottale’s model is formulated in the fractal dimension , that is, for movements on Peano curves and for Wiener’s stochastic processes [1, 13, 25–27].
3. Generalized Principle of Scale Covariance: Fractal “Diffusion” Type Equation
We now apply the principle of scale covariance and postulate that the transition from classical (differentiable) to “fractal” physics can be implemented by means of replacing the standard time derivative operator with the complex one (this result is a generalization of Notalle’s scale covariance principle as in [13]). As a consequence, we can now write the fractal “diffusion” type equation in its covariant form
This means that at any point of a fractal curve, the local temporal term , the nonlinearly “convective” term , and the dissipative one equilibrate.
3.1. Standard “Diffusion” Type Equation
Separating the real part from the imaginary one in (21), that is, we can add these two equations and obtain a generalized “diffusion” type law in the form
The standard “diffusion” law, that is, results from (23) the following assertions:(i)diffusion paths are Peano type fractal curves. In other words, the fractal dimension of the fractal curves is . Moreover, the average values (19) are defined through Wiener’s stochastic processes [1, 13, 25–27], that is,(ii)the movements at differentiable and nondifferentiable scales are synchronous, that is, ,(iii)the structure coefficient , that characterizes fractalnonfractal transition, identifies with diffusion coefficient, that is,
This kind of “diffusion” is analyzed in [28, 29], a situation for which Papadopoulou et al. [30] calculated the value of the constant and obtained , equivalent with .
3.2. Anomalous “Diffusion” Type Equation: Weibull Relation
The anomalous diffusion law results from (23) on the following assumptions:(i)diffusion paths are fractal curves with fractal dimension ,(ii)time resolution δt is identified with the differential element dt, that is, in this case, the substitution principle can be also applied,(iii)movements at differentiable and nondifferentiable scales are synchronous, that is, .
Then, (23) can be written as follows:
For the onedimensional case, applying the variable separation method [31] with standard initial and boundary conditions implies where is a system characteristic length, and is a separation constant, dependent on diffusion order .
Accepting the viability of the substitution principle (see Section 2.1), from (30), through integration, it follows that
According to certain results of the fractional integraldifferential calculations [32, 33], (31) becomes
Moreover, (32) can be written under the form
The timedependent relative variation of concentrations is defined as where and are cumulative amounts of drug released at time and infinite time.
From (33) and (34), it follows that equation similar to Weibull relation , where and represent characteristic constants for each system that can be defined by
We can notice that both and constants are functions of fractal dimension of the curves on which drug release mechanism takes place. This represents the complexity and nonlinearity dynamics of the system. Moreover, constant is also dependent on the “diffusion” order n.
Let us note that the fractal processes [1, 13, 25–27] given by (21) with are considered as being “anomalous diffusion” (subdiffusion for and superdiffusion for ). The “FokkerPlanck equations” for anomalous diffusion do not usually have the form of an ordinary diffusion equation. Indeed, it is well known that the “FokkerPlanck equations” for anomalous diffusion take the form of fractional derivative equations, called fractional FokkerPlanck equations [25–27].
This new approach will be presently applied in the analysis of one polymerdrug type system that has been obtained by our researchers.
4. Experimental Results and Discussion
This paragraph presents some experimental results obtained by our researchers, for different kind of polymer matrices.
4.1. ChloramphenicolLoaded GELPVA Microparticles
Gelatin and poly(vinyl alcohol) (GELPVA) microparticles crosslinked with glutaraldehyde (GA), for samples with different amount of crosslinking agent (2%, 6%, 8%, 10%, …, the sample code indicates the crosslinking amount: for example, GA2 represents a sample with 2% cross linking amount) loaded with chloramphenicol were tested (details regarding materials and experimental protocol can be found in [34]).
The experimental points of release kinetics are shown in Figure 1.
4.2. LevofloxacinLoaded StarchBased Hydrogels
Hydrogels with controlled hydrophilicity were obtained from new amphiphilic starchpolyester graft copolymer and unmodified starch with fatty acid chains, crosslinked with citric acid, at different crosslinking time (4 h, 5 h, 6 h, … the sample code indicates the crosslinking time: for example, HS4 represents a sample with 4 h crosslinking time). The samples were loaded with levofloxacin, and their release kinetics was studied (details regarding materials and experimental protocol can be found in [35]).
The experimental points of release kinetics are shown in Figure 2.
4.3. AdrenalineLoaded GELCHI Microparticles
Gelatin and chitosan (GELCHI) microparticles were prepared within twostep crosslinking process in an emulsionphase separation system crosslinked with glutaraldehyde (GA). Samples were prepared with different gelatin: chitosangelatin amount ratios of (10%–90%, 30%–70%, 50%50%, …, the sample code indicates the gelatin/chitosan ratio: for example, GC 10–90 represents a 10% gelatin and 90% chitosan sample).
The samples were loaded with adrenaline, and their release kinetics was studied (details regarding materials and experimental protocol can be found in [36]).
The experimental points of release kinetics are shown in Figure 3.
4.4. Experimental Results Analysis and Discussions
As a result of different experimental protocols, all the above drug carriers represent differentiated polymer matrices. Despite this, they have qualitatively similar time behavior. Therefore, qualitative analysis should be performed according to the loaded drug amount that is directly proportional to the released drug amount.
These results have been analyzed matching experimental data with Peppas and Weibull laws, as it has been demonstrated in the previous paragraph. As a result, for each of the samples (see Table 1), we obtained Peppas and Weibull parameters, the correlation factors, and release kinetics fractal dimension.

The first observation refers to the correlation coefficient between the experimental curves and Weibull fitted curve, which proved to be a very good one, even better than in the case of Peppas curve. This means that the entire release process can be better described by Weibull law as compared to Peppas law, thus showing the wide applicability of the former (the best values for correlation factors are italicized, in Table 1). Experimental Weibull curves for HS and GA samples are graphically illustrated in Figure 4.
(a)
(b)
At the same time, the values for n, namely no. 0,5 value, show that in all these cases the diffusion mechanism is a nonFickian one, while diffusion does not prevail. There are also some other phenomena with relevant contribution: the physical interaction of drugpolymer matrix with release environment, chemical reactions, and drug/polymer degradation.
Consequently, the complexity of the phenomena manifestly generates a complex trajectory for drug particles. It is known that the complex trajectory can be measured through fractal dimension, namely, in this particular case, the fractal dimension of the release curve. This is confirmed by fractal dimension values, according to (36b). The values between 1 and 3 generally correspond to the usually accepted values in the case of fractal processes [37]; higher values denote either fractal dimension should be redefined as a function of structure “classes,” or the drug release process is a complex one, with many freedom degrees in the phase space [38]. Another observation refers to samples with (“subdiffusion”) or with , (superdiffusion) (see paragraph 3.). This classification is in compliance to the experimental observation according to which these samples, unlike others, exhibit “faster” diffusion, with higher diffusion rate.
At the same time, the value for fractal dimension of the adrenaline release curve for GC 30–70 sample is the highest, reflecting the high degree of complexity. Consequently, the behavior of such systems should be differently described. Such an approach will be subject of future research. GC samples exhibit strong variations on the first part of the release curve. For two of them, the attempt to match with Peppas and Weibull equations failed, the parameters having negative values that led to negative fractal dimensions (written in bold, in Table 1). These are values that can point to some phase transitions in multiphase complex space [39].
5. Conclusions
In this paper, we have replaced the complexity of the physical processes that generates drug release from polymer matrix with fractality. That is why the whole classical “arsenal” of standard physics quantities (differentiable physics) was no longer used. Consequently, fractal approximation of motion was introduced in the study of the complex physicchemical dynamics of the system. By means of fractional calculus, the fractal “diffusion” equation gives rise to Weibull relation, a statistical distribution function of wide applicability, including drug release studies. Our approach has considered all simultaneously involved phenomena that are equivalent to complexity and fractality. Consequently, a physical basis to this equation and for its parameters has been created. The parameters are functions of fractal dimension of the curves on which drug release mechanism takes place. This dimension measures the complex and nonlinear dynamics of the system that depends on diffusion order.
This type of approach is a viable one, since the analysis of our experimental results shows that experimental curves can match with some very good correlation factors, through Weibull laws, better than in the case of power type laws. At the same time, the fractal dimension of the drug release curve depends on the complexity of the system and can be considered as control parameter for the nonlinear behavior of it.
Acknowledgments
This paper was supported by PERFORMERA (“Postdoctoral Performance for Integration in the European Research Area” (ID57649)) project, financed by the European Social Fund and the Romanian Government. Thanks to the referent for his pertinent observations.
References
 B. B. Mandelbrot, The Fractal Geometry of Nature, Freeman, San Francisco, Calif, USA, 1983.
 D. Stauffer and H. E. Stanle, From Newton to Mandelbrot, Academic Press, New York, NY, USA, 1996.
 V. U. Novikov and G. V. Kozlov, “Structure and properties of polymers in terms of the fractal approach,” Russian Chemical Reviews, vol. 69, no. 6, pp. 523–549, 2000. View at: Google Scholar
 G. V. Kozlov and G. E. Zaikov, Fractals and Local Order in Polymeric Materials, Nova Science Publishers, New York, NY, USA, 2001.
 E. Hornbogen, “Fractals in microstructure of metals,” International Materials Reviews, vol. 34, no. 6, pp. 277–296, 1989. View at: Google Scholar
 M. E. Cates, “Statics and dynamics of polymeric fractals,” Physical Review Letters, vol. 53, no. 9, pp. 926–929, 1984. View at: Publisher Site  Google Scholar
 J. Havlin and D. BenAvraham, “Fractal dimensionality of polymer chains,” Journal of Physics A, vol. 15, pp. L311–L316, 1982. View at: Google Scholar
 M. Muthukumar, “Dynamics of polymeric fractals,” The Journal of Chemical Physics, vol. 83, no. 6, pp. 3161–3168, 1985. View at: Google Scholar
 M. H. Bessendorf, “Stochastic and fractal analysis of fracture trajectories,” International Journal of Engineering Science, vol. 25, no. 6, pp. 667–672, 1987. View at: Google Scholar
 Y. Fu and W. J. Kao, “Drug release kinetics and transport mechanisms of nondegradable and degradable polymeric delivery systems,” Expert Opinion on Drug Delivery, vol. 7, no. 4, pp. 429–444, 2010. View at: Publisher Site  Google Scholar
 S. Popescu, Actual Issues in the Physics of SelfStructured Systems, Tehnopress, Iasi, Romania, 2003.
 L. Nottale, “Fractals and the quantum theory of spacetime,” International Journal of Modern Physics A, vol. 4, pp. 5047–5117, 1989. View at: Google Scholar
 L. Nottale, Fractal SpaceTime and Microphysics: Towards a Theory of Scale Relativity, World Scientific, Singapore, 1993.
 K. Kosmidis, P. Argyrakis, and P. Macheras, “Fractal kinetics in drug release from finite fractal matrices,” Journal of Chemical Physics, vol. 119, no. 12, pp. 6373–6377, 2003. View at: Publisher Site  Google Scholar
 P. Costa and J. M. Sousa Lobo, “Modeling and comparison of dissolution profiles,” European Journal of Pharmaceutical Sciences, vol. 13, no. 2, pp. 123–133, 2001. View at: Publisher Site  Google Scholar
 S. Popescu, “Turing structures in dc gas discharges,” Europhysics Letters, vol. 73, no. 2, pp. 190–196, 2006. View at: Publisher Site  Google Scholar
 A. J. Lichtenberg and M. A. Lieberman, Regular and Stochastic Motion, Springer, New York, NY, USA, 1983.
 L. Nottale, Scale Relativity and Fractal SpaceTime—A New Approach to Unifying Relativity and Quantum Mechanics, Imperial College Press, London, UK, 2011.
 R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, MacGrawHill, New York, NY, USA, 1965.
 D. G. Dimitriu, “Physical processes related to the onset of lowfrequency instabilities in magnetized plasma,” Czechoslovak Journal of Physics, vol. 54, supplement, no. C468, p. C474, 2004. View at: Google Scholar
 D. G. Dimitriu, V. Ignatescu, C. Ioniţǎ, E. Lozneanu, M. Sanduloviciu, and R. W. Schrittwieser, “The influence of electron impact ionisations on low frequency instabilities in a magnetised plasma,” International Journal of Mass Spectrometry, vol. 223224, pp. 141–158, 2003. View at: Publisher Site  Google Scholar
 M. Agop, N. Forna, I. CasianBotez, and C. Bejenariu, “New theoretical approach of the physical processes in nanostructures,” Journal of Computational and Theoretical Nanoscience, vol. 5, no. 4, pp. 483–489, 2008. View at: Google Scholar
 I. CasianBotez, M. Agop, P. Nica, V. Paun, and G. V. Munceleanu, “Conductive and convective types behaviors at nanotime scales,” Journal of Computational and Theoretical Nanoscience, vol. 7, no. 11, pp. 2271–2280, 2010. View at: Publisher Site  Google Scholar
 G. V. Munceleanu, V. P. Paun, I. CasianBotez, and M. Agop, “The microscopicmacroscopic scale transformation through a chaos scenario in the fractal spacetime theory,” International Journal of Bifurcation and Chaos, vol. 21, no. 2, pp. 603–618, 2011. View at: Publisher Site  Google Scholar
 J. F. Gouyet, Physique et Structures Fractals, Masson, Paris, France, 1992.
 M. S. El Naschie, O. E. Rössler, and I. Prigogine, Quantum Mechanics, Diffusion and Chaotic Fractals, Elsevier, Oxford, UK, 1995.
 P. Weibel, G. Ord, and O. E. Rösler, Space Time Physics and Fractality, Springer, New York, NY, USA, 2005.
 K. C. Sung, P. R. Nixon, J. W. Skoug et al., “Effect of formulation variables on drug and polymer release from HPMCbased matrix tablets,” International Journal of Pharmaceutics, vol. 142, no. 1, pp. 53–60, 1996. View at: Publisher Site  Google Scholar
 J. Siepmann, A. Streubel, and N. A. Peppas, “Understanding and predicting drug delivery from hydrophilic matrix tablets using the “sequential layer” model,” Pharmaceutical Research, vol. 19, no. 3, pp. 306–314, 2002. View at: Publisher Site  Google Scholar
 V. Papadopoulou, K. Kosmidis, M. Vlachou, and P. Macheras, “On the use of the Weibull function for the discernment of drug release mechanisms,” International Journal of Pharmaceutics, vol. 309, no. 12, pp. 44–50, 2006. View at: Publisher Site  Google Scholar
 L. Jude, Mathematics Physics Equations. Theory and Applications, Matrix Rom Publishing, Romania, 2010.
 K. B. Oldham and J. Spanier, The Fractional Calculus: Theory and Applications of Differential and Integration to Arbitrary Order, Dover Publications, New York, NY, USA, 2006.
 A. A. Kilbas, H. M. Srivastava, and J. J. Trujilto, Theory and Applications of Fractional Differential Equations, Elsevier, Armsterdam, The Netherlands, 2006.
 C. A. Peptu, A. Perichaud, and M. Popa, “Hydrogel microspheres based on environmentally friendly polymers with potential biomedical applications,” Environmental Engineering and Management Journal, vol. 10, no. 5, pp. 717–727, 2011. View at: Google Scholar
 A. Uliniuc, Chemical modifications of polysaccharides and their hydrogels through “clickchemistry” procedure, Ph.D. thesis, “Ghe. Asachi” Technical University of Iasi, Iasi, Romania, 2011.
 C. A. Peptu, G. Buhus, M. Popa, A. Perichaud, and D. Costin, “Double crosslinked chitosan′"Gelatin particulate systems for ophthalmic applications,” Journal of Bioactive and Compatible Polymers, vol. 25, no. 1, pp. 98–116, 2010. View at: Publisher Site  Google Scholar
 P. S. Addison, Fractals and Chaos—An Illustrated Course, Institute of Physics Publishing, Londra, Regno Unito, 1997.
 A. J. Lichtenberg, PhaseSpace Dynamics of Particle, John Wiley & Sons, New York, NY, USA, 1969.
 M. H. Jensen, G. Paladin, and A. Vulpiani, “Random fractals, phase transitions, and negative dimension spectra,” Physical Review E, vol. 50, no. 6, pp. 4352–4356, 1994. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 S. Băcăiţă et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.