Next Article in Journal
Mechanical Properties of Polyvinyl Alcohol Fiber-Reinforced Cementitious Composites after High-Temperature Exposure
Next Article in Special Issue
Photocatalytic and Antibacterial Properties of Doped TiO2 Nanopowders Synthesized by Sol−Gel Method
Previous Article in Journal
Hydrogels, Oleogels and Bigels as Edible Coatings of Sardine Fillets and Delivery Systems of Rosemary Extract
Previous Article in Special Issue
A Molecular Description of Hydrogel Forming Polymers for Cement-Based Printing Paste Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multifractal Vision of 5-Fluorouracil Release from Chitosan-Based Matrix

by
Maria-Alexandra Paun
1,2,
Vladimir-Alexandru Paun
3 and
Viorel-Puiu Paun
4,5,*
1
School of Engineering, Swiss Federal Institute of Technology (EPFL), 1015 Lausanne, Switzerland
2
Division Radio Monitoring and Equipment, Section Market Access and Conformity, Federal Office of Communications (OFCOM), 2501 Bienne, Switzerland
3
Five Rescue Research Laboratory, 75004 Paris, France
4
Physics Department, Faculty of Applied Sciences, University Politehnica of Bucharest, 060042 Bucharest, Romania
5
Academy of Romanian Scientists, 050085 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Gels 2022, 8(10), 661; https://doi.org/10.3390/gels8100661
Submission received: 25 September 2022 / Revised: 10 October 2022 / Accepted: 14 October 2022 / Published: 16 October 2022

Abstract

:
A suite of four drug deliverance formulations grounded on 5-fluorouracil enclosed in a chitosan-founded intercellular substance was produced by 3,7-dimethyl-2,6-octadienal with in situ hydrogelation. The formulations have been examined from a morphological and structural point of view by Fourier transform infrared (FTIR) spectroscopy and microscopy with polarized light, respectively. The polarized optical microscopy (POM) pictures of the three representative formulations obtained were investigated by fractal analysis. The fractal dimension and lacunarity of each of them were thus calculated. In this paper, a novel theoretical method for mathematically describing medicament deliverance dynamics in the context of the polymeric medicament constitution limit has been advanced. Assuming that the polymeric drug motion unfolds only on the so-called non-differentiable curves (considered mathematically multifractal curves), it looks like in a one-dimensional hydrodynamic movement within a multifractal formalism, the drug-release physics models are provided by isochronous kinetics, but at a scale of resolution necessarily non-differentiable.

1. Introduction

A topical field of pharmaceutical research is the development of new semi-solid pharmaceutical systems as potential vehicles for the topical release of medicinal substances using excipients of natural origin. The benefits presented by these types of formulations are determined both by the diversity of excipients and by their particularly advantageous specific properties, namely, relative abundance, biocompatibility, biodegradability, non-irritability, innocuousness, and low-cost price. Hydrogels represent an important category of ointment bases, and are frequently used in the formulation of dermatological and cosmetic preparations due to the many advantages they offer. This group also includes biocompatible hydrogels based on natural polymers (e.g., chitosan, pectins, gums) which have been intensively studied in the last decade due to the benefits provided by the specific characteristics of the gel-forming agents.
An initial direction of research was aimed at the development and pharmaco-technical characterization of some hydrogel bases containing a natural hydrophilic polymer, chitosan, as a gelling agent, and different co-solvents able to incorporate and dissolve a large number of medicinal substances insoluble in water. At this stage of the study, the influence of the type and concentration of the co-solvent on the physico-chemical and rheological characteristics of the hydrogel bases with chitosan in different concentrations was evaluated. At the second stage of the experimental research, the possibility of formulating and preparing new biocompatible topical hydrogels based on chitosan together with different natural absorption promoters was investigated.
Chitosan-based pharmaceutical compounds are of great interest in the medicament delivery domain, and are suitable due to their inherent properties of biocompatibility, degradation comportment, and biodegradability, which are advised for current practices in vivo [1,2,3]. In order to further ameliorate chitosan’s capacity to anchor great medicaments’ quantity and to liberate them in a supervised fashion, a lot of tests consisting particularly of chitosan composition crosslinking by different agents have been designed. In this line of conduct, an ecological method was detailed by crosslinking chitosan with eco-friendly products of the monoaldehyde type [4]. The mentioned procedure was found to be a remarkable one, delivering hydrogels whose properties of an aldehyde nature were easily commanded. Therefore, by utilizing the native aldehyde 3,7-dimethyl-2,6-octadienal, well known typically as a right commercial denomination by citral, hydrogels with superior biocompatibility, biodegradability, non-toxic, and non-irritating qualities and with their own antifungal activity were achieved for matrix components for medicament release systems.
Principally, the approach that will be made in the present work, otherwise known as the fractal manner of pharmaco-kinetics (PK) attitude, assumes the utilization of fractal analysis and fractional calculus. Additionally, by increasing the accuracy of performance on the fractal dimension and lacunarity notion, maximum confidence will be achieved when interpreting the images obtained through microscopy. Similarly, it is possible in the context of the analysis of several compartments [5] to successfully describe the polymeric networks’ dynamics and diffusion in porous fractal media. It can be said that, recently, the analysis of several compartments by PK has allowed the modeling of processes such as medication dissolution [6], drug absorption [7], as well as real distribution, together with arranging the constituent parts to achieve the whole [8].
Other things worth mentioning in the introduction are studies that refer to the design of biocompatible chitosan with a high photothermal conversion capability [9], ultrathin fibers and their antibacterial food packaging applications [10] and, last but not least, electrospun structural nanohybrids combining three composites for fast helicopter delivery [11].
In the presented context, special attention was paid to the open problem related to the time development of the kinetics dependent on the major bio-molecular reactions [6]. In continuation, we can say that the compartmental ideal analysis, the newest method for reporting drug delivery dynamics in elaborate systems (operating at a level of fractional derivatives or alternative ordinary procedures utilized in PK), can be depicted, including non-differentiable curves (named multifractal curves).
Afterwards, as a work alternative with a single mathematical variable represented by a non-differentiable rigid function, it is feasible that we trust only in the approximations of these mathematically defined functions, procured by averaging functions at various scales of determination–resolution. In the aftermath, all proposed variables to depict drug delivery processes will continue to work as the mathematical functions are limited, thus dictating a life with non-differentiable character in the case of scale-zero/nil value resolutions differentiable in different circumstances. Ultimately, the theoretical multifractal model is consequently certified by the empirical experimental data associated with the 5-fluorouracil deliverance out of the chitosan-founded matrix.

2. Results and Discussion

Theoretical Model

Typical investigations on drug-release kinetics deliver significant information on chemical systems’ function. To better explain the subtle transport mechanism and the structure versus function connection of a chemical system, it is decisive to diminish the difference between the experimental macroscopic data and the transport comportment at the atomic and molecular structural level. Obviously, the transport mechanisms at the microscopic level will prevail; these are mechanisms that we also have in mind for the model proposed in the article.
In these circumstances of pharmaco-kinetic behavioral consideration, a novel theoretical method for mathematically describing medicament deliverance kinetics, in the context of the polymeric medicament constitution limit, is advanced [12,13]. Presuming that the polymeric medicament dynamics unfolds only on the so-called non-differentiable curves (considered mathematically multifractal curves), it looks like in a one-dimensional hydrodynamic conventionalism of multifractal arguments, the drug-release physics models are provided by kinetics of synchronous type, but at a scale of resolution necessarily non-differentiable.
The one-dimensional multifractal hydrodynamic equation [11,12,13,14,15] is
t V D + V D x V D = x [ 2 λ ( d t ) [ 4 f ( α ) ] 2 x ( x ρ ) ρ ]
t ρ + x ( V D ρ ) = 0
In the communicated expressions above, x is the fractal spatial coordinate (one-dimensional), t is the classical (non-fractal) time which plays the role of an affine parameter within the movement curves [14,15]. The speed VD is the differential velocity (scale resolution dt, independent of time), the function f(α) is the singularity spectrum, α is dependent, and √ρ is a function of the state’s amplitude.
The initial conditions (3) and boundary conditions (4) of the above equations are
V D ( x , t = 0 ) = V 0 ,     ρ ( x , t = 0 ) = 1 π α e x p [ ( x α ) 2 ]
V D ( x = V 0 t ) = V 0 ,     ρ ( x = , t ) = ρ ( x = + , t ) = 0
where V 0 is the initial speed, while α is the Gaussian distribution parameter of the position variable x.
According to the mathematical procedures recommended in bibliographical references [16], differential Equations (1) and (2) accept the following solutions
V D ( x , t , σ , α ) = V 0 α 2 + ( σ α ) 2 x t α 2 + ( σ α ) 2 t 2
ρ ( x , t , σ , α ) = π 1 / 2 ( α 2 + ( σ α ) 2 t 2 ) 1 / 2 exp [ ( x V 0 t ) 2 α 2 + ( σ α ) 2 t 2 ]
where
σ = λ ( d t ) [ 2 f ( α ) ] 1
is the multifractal degree.
According to the stated fractal theory, the non-differentiable velocity V F will have the formula
V F ( x , t , σ , α ) = σ ( x V 0 t ) α 2 + ( σ α ) 2 t 2
By inserting the dimensionless variables in the above expressions
ξ = x V 0 τ 0 ,       η = t τ 0
as well as the dimensionless parameters
μ = σ τ 0 α 2 ,       ϕ = α V 0 τ 0
where τ 0 is named specific time. Formulas (5), (6) and (8) are written as
V V D ( ξ , η , μ ) = V D ( x , t , σ , α ) V 0 = 1 + μ 2 ξ η 1 + μ 2 η 2
ρ ( ξ , η , μ , ϕ ) = π 1 2 α ρ ( x , t , σ , α ) = ( 1 + μ 2 η 2 ) 1 2 e x p [ ( ξ η ) 2 ϕ 2 ( 1 + μ 2 η 2 ) ]
U V F ( ξ , η , μ ) = V F ( x , t , σ , α ) V 0 = μ ( ξ η ) 1 + μ 2 η 2
In Figure 1a,b, 3D plot representations of the U ( ξ , η )   multifractal function are presented for two distinct sets of variable values ξ and η, utilized for drug-release analysis.
In Figure 2, 3D plot representations of the V ( ξ , η )   multifractal function are presented for two distinct sets of variable values ξ and η, utilized for drug-release analysis.
If we now make the ratio between U and V, respectively, U/V, we obtain a homographic expression dependent of ξ and η, as presented below:
U V = μ ( ξ η ) 1 + μ 2 ξ η
In Figure 3, 3D plot representations of the U/V multifractal function, with expressions dependent of the ξ and η variables, are presented.
In the case of the condition validity in which the differential of the ratio U/V is zero (so-called dynamic simultaneity)
d ( U V ) = d U d V V 2 = 0 V = c o n s t U
we have between U and V the relationship V = constU.
This coincides with the prolongation of the Newton first principle to any scale resolution, or similarly, with to the “synchronizations” of drug-release dynamics at a non-differentiable scale. Additionally, it can be mentioned that the varied “mechanisms” implicated in drug-release activity may be mimicked by doubling the period, or a manifested quasi-periodicity and the process’ intermittent solicitations [15].
Following the activation of Restriction (15), in the formulation V=-U and Equations (1) and (2), the multifractal-type conservation laws become a unique multifractal equation, a law of “diffusion” type
t ρ = λ ( d t ) [ 2 f ( α ) ] 1 l l ρ = σ l l ρ
As an immediate result, these evidentiated mechanisms manifests themselves as diffusion processes in a multifractal space but at diverse scale resolutions, which can be at large classified in the category of Fickian diffusion (when polymer relaxation time tr is much greater than the value td, characteristic solvent diffusion time) or, as non-Fickian-type diffusion (when trtd), with the same ease. If the medicament release is of the multifractal type and takes place in a perfect immersed state, the initial conditions and conditions at the boundary are accepted
t = 0 ,   α 2 < x < α 2 ,     ρ = ρ 0 t > 0 ,     x = ± α 2 ,     ρ = ρ 1
wherein ρ 0 is the primary drug density of the multifractal-type states, as part of the fractal-type “mechanism” used, and ρ 1   is the current medicament density according to polymer–fluid fractal interdependence. The equation solution in the given conditions may now be presented in the complex form as below [16]:
f = ρ t ρ = 2 ( σ t δ 2 ) 1 2 = { π 1 / 2 + n = 1 ( 1 ) n e r f c [ n δ 2 ( σ t ) 1 2 ] }
A precise expression was achieved for little values of time t
f = ρ t ρ = 2 ( σ t δ 2 ) 1 2 = c o n s t ( t ) 1 2
In these circumstances, the ratio ρ t ρ can be approved as an equivalent dissolved drug fraction, i.e., M t M ρ t ρ , where Mt is the drug quantity dissolved in time t and M is the drug quantity dissolved in the total time, considering that the pharmaceutical graduation shape has been already exhausted.
The function f depends on t, and σ plays the role of the parameter. In Figure 4, a 2D curve is obtained for a fixed value of the parameter. Thus, the black color represents the value σ = 1, the blue color represents σ = 2, the red color represents σ = 3, the green color represents σ = 4, and the olive color represents σ = 5. In Figure 4, a 3D function f = f (t, σ) is considered by three variables; one allows the variation of σ continuously with a parameter value between 0 and 10.
In Figure 5, a model confirmation of the 5-fluorouracil release out of the chitosan-founded matrix is depicted. The experimental results are fit tightly with the mathematical theoretical functions of multifractal conception replica. The obtained graphical figure demonstrated that the realized model is ready to successfully prophesy the drug-deliverance kinetics [15,16].
The legend of Figure 5, with the symbols for defining the experimental points and the colors for each type of fractal curve used in the case of the four tested substances, is presented above. As can be seen in Figure 5, from the graphic representation of cumulative drug release as a function of time, curves U1 and U4 are classic saturation curves, while curves U2 and U3 are increasing with time (tend to infinite values).

3. Evaluation of Polarized Optical Microscopy Pictures by Fractal Analysis

5-fluorouracil’s existence in chemical formulations was confirmed by means of polarized light microscopy (Figure 6). In the provided images, the obvious drug segregation into the hydrogels is clearly visible, with great density of crosslinking in compounds U1 and U2, in contrast with the hydrogel formulations of reduced crosslinking density (U4), about which several things are discussed. Thus, a birefringent comportment and a granular texture were declared; a specific feature of the crystal is its sub-micrometric dimensions, but it still falls below the detection limit of the analysis devices used [17].
The scale bar for the POM photographic images (from Figure 6: (a) U1; (b) U2; (c) U3) is 20 microns. The POM image in Figure 7 of experimentally produced formulations were processed according to fractal analysis standards, calculating the fractal dimension and lacunarity of each one. The processing method and the values obtained for fractal dimension and lacunarity, together with the voxel representation for each image separately are presented below, in paper continuation.
The fractal dimension and the lacunarity of the investigated POM images were calculated according to known fractal analysis procedures [18,19]. The same advanced mathematical utilities were successfully used in a previous work, but for images obtained by scanning electron microscopy (SEM), which will prove to be a complementary reading [20]. In order to obtain a high resolution, the calculation programs/software developed for the evaluation of some medical diseases with CT and MRI images were used. It is known that significant differences are noticed at the pixel level in order to give an early quantitative diagnosis (see [21,22,23]).
Following the numerical evaluations with the appropriate software of the selected image U1, the values of fractal dimension D = 1.7602, standard deviation s = ± σ 2 = ± 0.2026 (see [24]) and lacunarity Λ = 0.0215 were obtained [24,25,26], as in Table 1.
The graph in Figure 8 shows the results of the 2D box-count algorithm and the local fractal dimension of the investigated image U1 [25,27].
Figure 9 shows the verification of the selected U1 image area with the software Harmonic and Fractal Image Analyser Demo version 5.5.30 [28] of the fractal dimension for various ruler scales r.
Figure 10 represents a three-dimensional graph of the voxel representation for the image U1 from the modified area.
The POM image in Figure 11 of experimentally produced formulations were processed according to fractal analysis standards, calculating the fractal dimension and lacunarity of each one.
Following the numerical evaluations with the appropriate software of the selected image U2, the values of fractal dimension D = 1.7523, standard deviation s = ± σ 2 = ± 0.1949 (see [24]) and lacunarity Λ = 0.0363 were obtained [24,25,26], as in Table 2.
The graph in Figure 12 shows the results of the 2D box-count algorithm and the local fractal dimension of the investigated image U2 [25,27].
Figure 13 shows the verification of the selected U2 image area with the software Harmonic and Fractal Image Analyser Demo version 5.5.30 [28] of the fractal dimension for various ruler scales r.
Figure 14 represents the three-dimensional graph of voxel representation for the image U2 from the modified area.
Following the numerical evaluations with the appropriate software of the selected image U4, the values of fractal dimension D=1.7352, standard deviation s = ± σ 2 = ± 0.1831 (see [24]) and lacunarity Λ = 0.0385 were obtained [24,25,26], as in Table 3.
Figure 7a–d, Figure 11a–d and Figure 15a–d show the three transformations (b, c and d) of the original of Figure 7a, Figure 11a and Figure 15a. Thus, the following changes were obtained: (b) the grayscale version; (c) the gray scale version without luminance; (d) the binarized version. The thresholds of 14, 25, and 35 were used for binarization. In [Table 1, Table 2 and Table 3], the values of fractal dimension, standard deviation and lacunarity that were obtained following the evaluation of images U1, U2 and U3 from the original Figure 7a, Figure 11a and Figure 15a are shown.
The POM image in Figure 15 of experimentally produced formulations were processed according to fractal analysis standards, calculating the fractal dimension and lacunarity of each one.
The graph in Figure 16 shows the results of the 2D box-count algorithm and the local fractal dimension of the investigated image U4, [25,27].
Figure 17 shows the verification of the selected U4 image area with the software Harmonic and Fractal Image Analyser Demo version 5.5.30 [28] for the fractal dimension with various ruler scales r.
Figure 18 represents the three-dimensional graph of the voxel representation for the image U4 from the modified area.
Figure 8, Figure 9, Figure 12, Figure 13, Figure 16 and Figure 17 present the graphic representations of the 2D fractal dimension used for the linear regression slope calculation of the fractal dimension. Figure 10, Figure 14 and Figure 18 represent the three-dimensional graphs of voxel representation for the images U1, U2 and U4 from the modified area. The three coordinate axes are distributed as follows: the number of pixels on the ox axis, the number of pixels on the oy axis, and the gray level for each pixel on the oz axis. In 3D computer graphics, the voxel represents a numerical value associated with the regular grid of a three-dimensional coordinate space.

4. Conclusions

Suites of drug delivery arrangements were realized by 5-fluorouracil encapsulation in the hydrogel matrix configured by the reticulation of chitosan-based with 3,7-dimethyl-2,6-octadienal. The hydrogel demonstrated a special capability to anchor the medicament, ensuring its prolonged deliverance close to 24 h. The release speed was adjusted by fluctuating the crosslinking consistency, arriving at 90% for a great reticulation density and at 70 % for a reduced one. In the stated multi fractal paradigm, a mathematical model was evolved for the comprehension of the drug delivery dynamics, taking into account that these special comportments are closely related to the fact that the movement takes place on non-differentiable curves.
Worth mentioning is the fact that, for the one-dimensional case of the multifractal hydrodynamic formalism, the ratio between the differentiable and non-differentiable velocities of a distinct distance is contingent on time, in a homographic fashion. The simultaneous dynamic conditions suggest the synchronization of the drug deliverance techniques at both scale resolutions, formulated by the multifractal diffusion functions, since we have the diffusion process depending on the scale resolutions. The fitted curves of the theoretical model are in good agreement with the obtained experimental results.
The three images of the U1, U2 and U4 compounds, obtained by means of polarized light microscopy, were investigated according to the standards of fractal analysis. Thus, for image U1, we found the fractal dimension D = 1.7602 ± 0.2026 and lacunarity Λ = 0.0215 ; for image U2, we found the fractal dimension D = 1.7523±0.1949 and lacunarity Λ = 0.0363 ; and for Image U4, we found the fractal dimension D = 1.7352±0.1831 and lacunarity Λ = 0.0385 .

5. Materials and Methods

5.1. Materials

From the point of view of chemistry, the chitosan (CS) is a basic polysaccharide acquired by the deacetylation of chitin, one of the most plentiful native biopolymers in the world, with the exception of cellulose, evidently. In large part due to its high biocompatibility, biodegradability and reduced toxicity, chitosan has been utilized on a large scale in the pharmaceutical industry. In the last decade, chitosan utilization as a conveyor of continuous-release preparations to resolve the half-life period of chemotherapy medicaments has been announced [29]. In addition, chitosan-covered magnetic nanoparticles were utilized to extend the liberation of 5-fluorouracil (5-FU) [4]. CS (243 kDa, DA: 87%), 3,7-dimethyl-2,6-octadienal (95%), 5-FU, and phosphate tampon solution with a pH of 7.4 were acquisitioned from Aldrich chemical company and utilized as they were received.

5.2. Synthesis of the Hydrogel Formulation Preparations

The formulations used were fabricated through the hydrogelation in situ method of chitosan (CS) with 3,7-dimethyl-2,6-octadienal in the company of 5-fluorouracil (5-FU), which is a respected and classic protocol [3,4]. In a very short time, a 2% solution of 3,7-dimethyl-2,6-octadienal combined with 5-fluorouracil was gently dripped into a 3% chitosan-dissolved solution in aqueous acetic acid solution (1%). All substances were fabricated by Aldrich chemical company.
According to the observations made, it was accepted that 5-fluorouracil was fixed in the matricial structure as a crystalline component, and its size was assorted with the crosslinking density value. The in vitro drug release simulated in an environment that imitates the physiological habitat disclosed a 5-fluorouracil gradual deliverance, in close relation with the value of crosslinking density [23].
The drug and chitosan quantity were both retained steadily, while the 3,7-dimethyl-2,6-octadienal amount was diversified to touch various molar proportions of the amine/aldehyde groups, from the ratio of 1/1 to 4/1, in order to achieve hydrogels with different values of crosslinking consistency [7]. In the end, it was found that the hydrogelation period increased as the aldehyde quantity decreased. Thereby, the 1/1 molar ratio of amine/aldehyde functional classes immediately took place and continued sluggishly for 24 h for the 4/1 molar ratio of the same classes. Eventually, the acquired hydrogels were lyophilized and subjected to detailed analysis. Moreover, for the inventory, the formulations received the labels U1, U2, U3, U4, where the attached number was in agreement with the amino/aldehyde classes’ molar ratio.

5.3. Methods

The gelation period was visually established in the moment when the reaction blend was transformed out of its initial viscous state to a rubbery final state. We mention that the xerogels were realized by lyophilization of the adequate hydrogels, utilizing the Labconco FreeZone Freeze Dry System, (FreeZoner2.5 Liter Freeze Dry Systems) instrument, for a duration of 24 h at a temperature of −50 °C and pressure of 0.04 mbar. The formulations produced were described by FTIR spectroscopy using the FT-IR Bruker Vertex 70 Spectrophotometer(Art photonics GmbH, 12489 Berlin, Germany).
The drug-release kinetics of the detailed medicaments’ delivery systems were identified by recording the absorbance at 265 nm from the supernatant where the release was completed. The concentration of the solution was computed with the help of the Beer–Lambert law. The supernatant UV–vis spectra were recorded on a Horiba spectrophotometer, but the absorbance was accommodated on a previously drawn calibration curve. The 5-fluorouracil curve of calibration was drawn utilizing the absorption maximum spectrum at 265 nm.

Author Contributions

Conceptualization, V.-P.P. and M.-A.P.; methodology, V.-P.P.; software, V.-A.P.; validation, V.-P.P., M.-A.P. and V.-A.P.; formal analysis, V.-P.P., M.-A.P. and V.-A.P.; investigation, V.-A.P. and M.-A.P.; resources, V.-A.P. and M.-A.P.; data curation, V.-A.P.; writing—original draft preparation, V.-P.P.; writing—review and editing, M.-A.P. and V.-P.P.; visualization, V.-A.P.; supervision, V.-P.P.; project administration, V.-P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used to support the findings of this study cannot be accessed due to commercial confidentiality.

Acknowledgments

The co-authors M.A. Paun, V.A. Paun and V.P. Paun would like to thank Jenica Paun, for her continuous kind support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Prausnitz, M.R. Langer, R. Transdermal drug delivery. Nat. Biotechnol. 2008, 26, 1261–1268. [Google Scholar] [CrossRef] [PubMed]
  2. Patel, A.; Cholkar, K.; Agrahari, V.; Mitra, A.K. Ocular drug delivery systems: An overview. World J. Pharmacol. 2013, 2, 47–64. [Google Scholar] [CrossRef]
  3. Jacob, J.; Haponiuk, J.T.; Thomas, S.; Gopi, S. Biopolymer based nanomaterials in drug delivery systems: A review. Mater. Today Chem. 2018, 9, 43–55. [Google Scholar] [CrossRef]
  4. Xiong, S.; Marin, L.; Duan, L.; Cheng, X. Fluorescent chitosan hydrogel for highly and selectively sensing of p-nitrophenol and 2, 4, 6-trinitrophenol. Carbohydr. Polym. 2019, 225, 115253. [Google Scholar] [CrossRef] [PubMed]
  5. Pereira, L.M. Fractal Pharmacokinetics. Comput. Math. Methods Med. 2010, 11, 161–184. [Google Scholar]
  6. Kosmidis, K.; Argyrakis, P.; Macheras, P. Fractal kinetics in drug release from finite fractal matrices. J. Chem. Phys. 2003, 119, 6373–6377. [Google Scholar] [CrossRef] [Green Version]
  7. Le Mehaute, A.; Crepy, G. Introduction to transfer and motion in fractal media: The geometry of kinetics. Solid State Ionics 1983, 9-10, 17–30. [Google Scholar] [CrossRef]
  8. Barkai, E.; Metzler, R.; Klafter, J. From continuous time random walk to the fractional Fokker-Planck equation. Phys. Rev. E 2000, 61, 132–138. [Google Scholar] [CrossRef]
  9. Zhang, L.; He, G.; Yu, Y.; Zhang, Y.; Li, X.; Wang, S. Design of Biocompatible Chitosan/Polyaniline/Laponite Hydrogel with Photothermal Conversion Capability. Biomolecules 2022, 12, 1089. [Google Scholar] [CrossRef] [PubMed]
  10. Jiang, W.; Zhao, P.; Song, W.; Wang, M.; Yu, D.-G. Electrospun Zein/Polyoxyethylene Core-Sheath Ultrathin Fibers and Their Antibacterial Food Packaging Applications. Biomolecules 2022, 12, 1110. [Google Scholar] [CrossRef] [PubMed]
  11. Liu, H.; Wang, H.; Lu, X.; Murugadoss, V.; Huang, M.; Yang, H.; Wan, F.; Yu, D.-G.; Guo, Z. Electrospun structural nanohybrids combining three composites for fast helicide delivery. Adv. Compos. Hybrid Mater. 2022, 5, 1017–1029. [Google Scholar] [CrossRef]
  12. Higaki, K.; Yamashita, S.; Amidon, G.L. Time-dependent oral absorption models. J. Pharmacokinet. Pharmacodyn. 2001, 28, 109–128. [Google Scholar] [CrossRef] [PubMed]
  13. Karalis, V.; Tsantili-Kakoulidou, A.; Macheras, P. Quantitative structure–pharmacokinetic relationships for disposition parameters of cephalosporins. Eur. J. Pharm. Sci. 2003, 20, 115–123. [Google Scholar] [CrossRef]
  14. Chechetkin, V.R.; Lutovinov, V.S.; Turygin, A.Y. Multifractal structure of fully developed hydrodynamic turbulence. I. Kolmogorov’s third hypothesis revisited. J. Stat. Phys. 1990, 61, 573–588. [Google Scholar] [CrossRef]
  15. Basu, A.; Chakrabarti, B.K. Hydrodynamic descriptions for surface roughness in fracture front propagation. Philos Trans A Math Phys Eng Sci. 2019, 377, 20170387. [Google Scholar] [CrossRef] [Green Version]
  16. Kwak, K.; Yang, S. Developing a Multi-Dimensional Hydrodynamics Code with Astrochemical Reactions. IAU Gen. Assem. 2015, 29, 2250754. [Google Scholar]
  17. Iancu, R.; Irimiciuc, S.A.; Agop, M.; Frasila, M.; Paun, M.-A.; Paun, V.-A.; Paun, V.-P.; Stratulat, S. 5-fluorouracil release from chitosan-based matrix. Experimental and theoretical aspects. Mater. Plast. 2020, 57, 180–188. [Google Scholar] [CrossRef]
  18. O’Shaughnessy, B.; Procaccia, I. Analytical solutions for diffusion on fractal objects. Phys. Rev. Lett. 1985, 54, 455–458. [Google Scholar] [CrossRef]
  19. Richard, R. 1D Diffusion PDE. Introduction to Partial Differential Equations. 2019. Available online: http://personal.ph.surrey.ac.uk/~phs1rs/teaching/l3_pdes.pdf (accessed on 15 June 2022).
  20. Polowsky, P.J.; Tansman, G.; Kindstedt, P.S.; Hughes, J.M. Characterization and identification of surface crystals on smear-ripened cheese by polarized light microscopy. J. Dairy Sci. 2018, 101, 7714–7723. [Google Scholar] [CrossRef]
  21. Karperien, A.L.; Jelinek, H.F. Box-Counting Fractal Analysis: A Primer for the Clinician. In The Fractal Geometry of the Brain; Di Ieva, A., Ed.; Second Chapter; Springer Series in Computational Neuroscience; Springer Science + Business Media: New York, NY, USA, 2016. [Google Scholar]
  22. Karperien, A.; Jelinek, H.F.; Miloševic, N.T. Reviewing Lacunarity Analysis and Classification of Microglia in Neuroscience. In Proceedings of the 8th European Conference on Mathematical and Theoretical Biology, Kraków, Poland, 28 June–2 July 2011; pp. 1–6. [Google Scholar]
  23. Lungu, R.; Paun, M.-A.; Peptanariu, D.; Ailincai, D.; Marin, L.; Nichita, M.-V.; Paun, V.-A.; Paun, V.-P. Biocompatible Chitosan-Based Hydrogels for Bioabsorbable Wound Dressings. Gels 2022, 8, 107. [Google Scholar] [CrossRef]
  24. Nichita, M.V.; Paun, M.A.; Paun, V.A.; Paun, V.P. Fractal analysis of brain glial cells. Fractals dimension and lacunarity. Univ. Politeh. Buchar. Sci. Bull. Ser. A Appl. Math. Phys. 2019, 81, 273–284. [Google Scholar]
  25. Bordescu, D.; Paun, M.A.; Paun, V.A.; Paun, V.P. Fractal analysis of Neuroimagistic. Lacunarity degree, a precious indicator in the detection of Alzheimer’s disease. Univ. Politeh. Buchar. Sci. Bull. Ser. A Appl. Math. Phys. 2018, 80, 309–320. [Google Scholar]
  26. Postolache, P.; Borsos, Z.; Paun, V.A.; Paun, V.P. New Way in Fractal Analysis Of Pulmonary Medical Images. Univ. Politeh. Buchar. Sci. Bull. Ser. A Appl. Math. Phys. 2018, 80, 313–322. [Google Scholar]
  27. Scott, D.W. Statistics: A Concise Mathematical Introduction for Students, Scientists, and Engineers; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2020. [Google Scholar]
  28. Available online: http://imagesci.fch.vut.cz/includes/harfa_download.inc.php (accessed on 15 June 2022).
  29. Wang, X.; Ahmed, B.N.; Alvarez, S.G.; Tuttolomondo, M.V.; Hélary, C.; Desimone, M.F.; Coradin, T. Sol-gel encapsulation of biomolecules and cells for medicinal applications. Curr. Top. Med. Chem. 2015, 15, 223–244. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Three-dimensional plot representations of U ( ξ , η )   multifractal function for two distinct sets of variable values ξ and η: (a) ξ from −40 to 40; η from 0 to 10; (b) ξ from −5 to 25; η from 0 to 10.
Figure 1. Three-dimensional plot representations of U ( ξ , η )   multifractal function for two distinct sets of variable values ξ and η: (a) ξ from −40 to 40; η from 0 to 10; (b) ξ from −5 to 25; η from 0 to 10.
Gels 08 00661 g001
Figure 2. Three-dimensional plot representations of V ( ξ , η )   multifractal function for two distinct sets of variable values ξ and η: (a) ξ from −5 to 25; η from 0 to 10; (b) ξ from −40 to 40; η from 0 to 10.
Figure 2. Three-dimensional plot representations of V ( ξ , η )   multifractal function for two distinct sets of variable values ξ and η: (a) ξ from −5 to 25; η from 0 to 10; (b) ξ from −40 to 40; η from 0 to 10.
Gels 08 00661 g002
Figure 3. Three-dimensional plot representations of U/V multifractal function.
Figure 3. Three-dimensional plot representations of U/V multifractal function.
Gels 08 00661 g003
Figure 4. Two-dimensional contour (a) and three-dimensional plot (b) representations of f multifractal function.
Figure 4. Two-dimensional contour (a) and three-dimensional plot (b) representations of f multifractal function.
Gels 08 00661 g004
Figure 5. 5-fluorouracil release of experimental presentation [17]; the formulations are accommodated in multifractal theoretical replica (solid-colored lines).
Figure 5. 5-fluorouracil release of experimental presentation [17]; the formulations are accommodated in multifractal theoretical replica (solid-colored lines).
Gels 08 00661 g005
Figure 6. Representative POM images [17] of the formulations (a) U1; (b) U2; (c) U4.
Figure 6. Representative POM images [17] of the formulations (a) U1; (b) U2; (c) U4.
Gels 08 00661 g006
Figure 7. Primary processing of the selected image U1. (a)—original image (the entire portion); (b)—the grayscale version; (c)—the gray scale version without luminance; (d)—binarized version. A threshold of 14 was used for binarization.
Figure 7. Primary processing of the selected image U1. (a)—original image (the entire portion); (b)—the grayscale version; (c)—the gray scale version without luminance; (d)—binarized version. A threshold of 14 was used for binarization.
Gels 08 00661 g007
Figure 8. Box-count algorithm results: fractal dimension.
Figure 8. Box-count algorithm results: fractal dimension.
Gels 08 00661 g008
Figure 9. Graphic of fractal dimension for selected U1 image area.
Figure 9. Graphic of fractal dimension for selected U1 image area.
Gels 08 00661 g009
Figure 10. Voxel representation for the image U1.
Figure 10. Voxel representation for the image U1.
Gels 08 00661 g010
Figure 11. Primary processing of the selected image U2. (a)—original image (the entire portion); (b)—the grayscale version; (c)—the gray scale version without luminance; (d)—binarized version. A threshold of 25 was used for binarization.
Figure 11. Primary processing of the selected image U2. (a)—original image (the entire portion); (b)—the grayscale version; (c)—the gray scale version without luminance; (d)—binarized version. A threshold of 25 was used for binarization.
Gels 08 00661 g011
Figure 12. Box-count algorithm results: fractal dimension.
Figure 12. Box-count algorithm results: fractal dimension.
Gels 08 00661 g012
Figure 13. Graphic of fractal dimension for selected U2 image area.
Figure 13. Graphic of fractal dimension for selected U2 image area.
Gels 08 00661 g013
Figure 14. Voxel representation for the image U2.
Figure 14. Voxel representation for the image U2.
Gels 08 00661 g014
Figure 15. Primary processing of the selected image U4. (a)—original image (the entire portion); (b)—the grayscale version; (c)—the gray scale version without luminance; (d)—binarized version. A threshold of 35 was used for binarization.
Figure 15. Primary processing of the selected image U4. (a)—original image (the entire portion); (b)—the grayscale version; (c)—the gray scale version without luminance; (d)—binarized version. A threshold of 35 was used for binarization.
Gels 08 00661 g015
Figure 16. Box-count algorithm results: fractal dimension.
Figure 16. Box-count algorithm results: fractal dimension.
Gels 08 00661 g016
Figure 17. Graphic of fractal dimension for selected U4 image area.
Figure 17. Graphic of fractal dimension for selected U4 image area.
Gels 08 00661 g017
Figure 18. Voxel representation for the image U4.
Figure 18. Voxel representation for the image U4.
Gels 08 00661 g018
Table 1. Calculation of fractal parameters for picture U1.
Table 1. Calculation of fractal parameters for picture U1.
NameFractal DimensionStandard DeviationLacunarity
Image U11.7602±0.20260.0215
Table 2. Calculation of fractal parameters for picture U2.
Table 2. Calculation of fractal parameters for picture U2.
NameFractal DimensionStandard DeviationLacunarity
Image U21.7523±0.19490.0363
Table 3. Calculation of fractal parameters for picture U4.
Table 3. Calculation of fractal parameters for picture U4.
NameFractal DimensionStandard DeviationLacunarity
Image U41.7352±0.18310.0385
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Paun, M.-A.; Paun, V.-A.; Paun, V.-P. A Multifractal Vision of 5-Fluorouracil Release from Chitosan-Based Matrix. Gels 2022, 8, 661. https://doi.org/10.3390/gels8100661

AMA Style

Paun M-A, Paun V-A, Paun V-P. A Multifractal Vision of 5-Fluorouracil Release from Chitosan-Based Matrix. Gels. 2022; 8(10):661. https://doi.org/10.3390/gels8100661

Chicago/Turabian Style

Paun, Maria-Alexandra, Vladimir-Alexandru Paun, and Viorel-Puiu Paun. 2022. "A Multifractal Vision of 5-Fluorouracil Release from Chitosan-Based Matrix" Gels 8, no. 10: 661. https://doi.org/10.3390/gels8100661

APA Style

Paun, M. -A., Paun, V. -A., & Paun, V. -P. (2022). A Multifractal Vision of 5-Fluorouracil Release from Chitosan-Based Matrix. Gels, 8(10), 661. https://doi.org/10.3390/gels8100661

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop