Abstract

The main aim of the paper was to simulate the drug release by a multifractal theoretical model, as a valuable method to assess the drug release mechanism. To do this, drug delivery films were prepared by mixing poly(vinyl alcohol boric acid) (PVAB) and diclofenac (DCF) sodium salt drug in different mass ratios from 90/10 to 70/30, in order to obtain drug delivery systems with different releasing rates. The different drug content of the three systems was confirmed by energy-dispersive spectroscopy (EDAX) analysis, and the encapsulation particularities were investigated by scanning electron microscopy (SEM), atomic force microscopy (AFM), and polarized optical microscopy (POM) techniques. The ability of the PVAB matrix to anchor the DCF was assessed by Fourier transform infrared (FTIR) spectroscopy. The in vitro release of the diclofenac sodium salt from the formulations was investigated in biomimetic conditions ( and 37°C) by UV-Vis spectroscopy, measuring the absorbance of the drug at 275 nm and fitting the results on a previously drawn calibration curve. An estimation of the drug release kinetics was performed by fitting three traditional mathematical models on experimental release data. Further, the drug delivery was simulated by the fractal theory of motion, in which the release dynamics of the polymer-drug complex system is described through various Riccati-type “regimes.” To explain such dynamics involved multifractal self-modulation in the form of period doubling, quasiperiodicity, intermittency, etc., as well as multifractal self-modulation of network type. Standard release dynamics were explained by multifractal behaviors of temporary kink type. The good correlation between the traditional mathematical models and the new proposed theoretical model demonstrated the validity of the multifractal model for the investigation of the drug release.

1. Introduction

Drug delivery is one of the most important fields related to medicine and healthcare, in which materials innovation plays a crucial role [1, 2]. Comprising different approaches and different systems, drug delivery is aimed at transporting pharmaceutical active compounds within a living organism, creating the frame for its therapeutic effect. The versatility and multitude of systems used in drug delivery comes from the wide range of both drugs and drugs carriers [38], considering that the same active molecule, encapsulated in a different carrier, gets a new life, due to the intrinsic properties of the carriers but also due to the interactions established between the matrix and the drug. Moreover, the existence of plenty of administration routes such as oral, nasal, ocular, and intravenous makes the field of drug delivery a complex one in which continuous development and innovation are required [9, 10].

Diclofenac (DCF) is a nonsteroid anti-inflammatory drug (NSAID), frequently used for the amelioration of acute or chronic painful conditions. It is efficient for the treatment of acute postoperative pain [11], osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, extra-articular rheumatism [12], pancreatitis [13], urological disses [14], and so on. Moreover, recent evidences showed that DCF is also able to inhibit tumor angiogenesis of some cancer cell lines [15]. It appears that the common root of its anti-inflammatory and antitumor efficiency results from its ability to inhibit the cyclooxygenase enzyme [14]. DCF, which is a potent inhibitor of COX-2 and prostaglandin E2 synthesis, displays a range of effects on the immune system, angiogenic cascade, chemo- and radio-sensitivity, and tumor metabolism. Nevertheless, the inhibition of this enzyme was also associated with side effects, such as gastrointestinal disorders, liver hypersensitivity reaction, immune hemolytic anemia, and immune thrombocytopenia. Systematic investigations revealed that these side effects are caused by interplay of metabolic and immunologic factors, and their incidence is higher for the patients who take DCF at a dose of more than 150 mg/day [16]. Moreover, the low bioavailability of DCF reached by oral administration is a real concern, as its presence in environment raises health problems for aquatic organisms, plants, and mammals and led to the development of drug-resistant strains [17]. In this context, there is a constant concern to find new pathways for an efficient DCF administration, which should minimize the side effects. To this end, the design of formulations for DCF prolonged release keeps the promise of a potential solution. Many systems were designed and investigated, based on the encapsulation of DCF into a large realm of polymeric matrixes, based on carboxymethyl cellulose [18], silk fibroin membranes [19], chitosan [6], polyvinyl alcohol [2023], poly(D,L-lactic acid-co-glycolic acid) [24], poly(D,L-lactide-co-glycolide) [2528], poly(epsilon-capro/D,L-lactide) [29], poly(maleic anhydride-alt-2-methoxyethyl vinyl ether) [30, 31], pectin [32], and so on. These data showed that a prolonged release of diclofenac can be reached from polymeric matrices able to develop intermolecular forces with diclofenac molecules [6, 33]. With all these in mind, we designed novel formulations based on diclofenac and poly(vinyl alcohol boric acid) [4]. The rational design was raised from the idea of strong intermolecular forces between the electron-deficient boron atom of PVAB and esteric groups of diclofenac, which should assure a good anchoring of the drug into the matrix [3437] and thus a slow release. The in vitro release has been investigated in an environment mimicking the body fluids, and the drug release kinetics was assessed by fitting the data on three traditional mathematical models. The next step in the investigation of such novel systems was the mathematical modeling of the diclofenac release by the fractal theory of motion, a new approach which offers various advantages when compared with the other well-known empirical and semiempirical models.

2. Materials and Methods

Diclofenac (DCF) sodium salt, polyvinyl alcohol boric acid (PVAB, ), and ethanol (>99.8%) were purchased from Sigma Aldrich and used without any previous purification, while double-distilled water was obtained in the laboratory.

2.1. Preparation of the Drug Delivery Systems

The drug delivery systems were obtained by solvent-induced phase separation method, by three consecutive cycles of 30 min stirring/5 min vortex applied to a mixture of PVAB and DCF solutions in different mass ratios [38, 39]. More precisely, to 10 mL PVAB solution 7.5% in water, were added 1.666, 3.75, and 6.42 mL solution of DCF in ethanol (5%), leading to three different formulations, D1, D2, and D3, containing different amounts of drugs. After the third cycle, the solutions were casted in Petri dishes and dried at room temperature, leading to self-standing flexible films (Figure 1).

2.2. The Evaluation of the DCF Release Kinetics

The ability of the obtained formulations to release the encapsulated drug was monitored in vitro, according to a previously described procedure [4]. Pieces of films from the D1-D3 formulations, containing the same amount of drug (14 mg), were used in the in vitro release determinations. The samples had sizes between ~1 cm2 and ~3 cm2 depending on drug content. They were immersed in 10 mL PBS, . In order to mimic the in vivo conditions, the temperature was kept constant at 37°C all over the experiment. From time to time, 2 mL from the supernatant was removed and replaced with fresh PBS. The absorbance of the DCF drug in the supernatant was evaluated by UV-Vis spectroscopy and correlated with concentration through a previously drawn calibration curve, on the DCF absorption maximum from 275 nm [6]. All the measurements were done in triplicate.

The absorbance of the drug was measured on a PerkinElmer Lambda 35 UV-Vis spectrophotometer.

The kinetic data were fitted on three mathematical models, further described as follows: (a)Zero-Order Model. , where is the amount of drug dissolved in the time and is the zero-order release constant(b)Higuchi Model. , where is the amount of drug released in the time and is the Higuchi dissolution constant(c)Korsmeyer-Peppas Model. , where is the fraction of drug released in the time , is the release rate constant, and is the release exponent

2.3. Mathematical Model
2.3.1. Release Dynamics through Multifractal Functions

Classical, commonly used models are usually founded on the otherwise unjustified supposition that variables describing the dynamics of any polymer-drug complex system are differentiable [40, 41] (for details, see the models of zero-order, Higuchi, Korsmeyer-Peppas, etc. [6]). Thus, the success of the abovementioned models should be understood as sequential, on domains in which differentiability and integrability are still valid. The differentiable and integrable mathematical procedures are otherwise inadequate when the dynamics of any polymer-drug complex system involves both nonlinearity and chaoticity. However, in order to describe such dynamics, but still employing differential mathematical procedures, it is necessary to introduce the scale resolution explicitly into the expression of physical variables and implicitly into the expression of fundamental equations that govern these dynamics. This means that any variable dependent on space and time coordinates, in a classical sense, will depend both on the space and time coordinates and on the scale resolutions in the new mathematical sense (that of nondifferentiability and nonintegrability). In other words, instead of operating with a variable described through a nondifferentiable function, we will work with approximations of this mathematical function, obtained by its mediation at various scale resolutions. As a result, any variable designed to describe the dynamics of any polymer-drug complex system will function as the limit of a family of mathematical functions, this being nondifferentiable for null scale resolution and differentiable for nonzero scale resolutions [42, 43].

This method of describing the dynamics of any polymer-drug complex system clearly implies the development of new geometrical structures and also of new mathematical models for which the motion laws, invariant to spatial and temporal transformations, are integrated with scale laws, invariant to spatial and temporal scale transformations. In our opinion, such a geometrical structure can be based on the concept of a “multifractal,” and the corresponding mathematical model can be based on the fractal theory of motion in an arbitrary and constant fractal dimension. Applications of the model in the analysis of the dynamics of complex systems are given in [4453]. For complex systems with release, the dynamics analysis is given in [5458].

The fundamental assumption of our model is the one that the dynamics of structural units of any polymer-drug complex system will be described by continuous but nondifferentiable motion curves (multifractal motion curves). These multifractal motion curves exhibit the property of self-similarity in every point, which can be translated into a property of holography (every part reflects the hole). Basically, we are discussing about “holographic implementations of structural units’ dynamics of any polymer-drug complex system” through multifractal “regimes” of Riccati-type equations (describing the dynamics of structural units of any polymer-drug complex system by using Riccati-type equations at various scale resolutions).

Obviously, logistic-type equations (which are a particular case of Riccati-type equations) for describing the dynamics of the structural units of any complex polymer-drug system have been used [59]. However, their applicability is limited considering the absence of the scale resolution from the description of such dynamics. A first result within the meaning of the above was given by us in [60]. In what follows, we will analyze the dynamics in complex polymer-drug systems using Riccati-type equations at various scale resolutions (multifractal Riccati-type equations).

Let us consider now the multifractal Riccati-type equation: where is the multifractal mass of drug released at time , is the nonmultifractal time, with the role of affine parameter of the release curves, and , , and are the multifractal parameters characteristic of the release dynamics.

These parameters depend on the intrinsic structure of the polymeric matrix and of the drug, on the polymer-drug matrix interactions, on the complex release system-biostructure interactions (which imply the immunity system of the biostructure), and on the scale resolution.

2.3.2. Multifractal Transient Phenomena in Drug Delivery

With regard to the solution of the multifractal Riccati-type equation (1), we must first notice that the roots of the polynomial can be written as

Performing the homographic transformation it results through direct calculus that is a solution of the linear and homogenous first-order equation which allows the solution

Therefore, if the initial condition is conveniently expressed, the general solution of (1) can be found by writing (4) as where is a real multifractal constant which characterizes the solution. By using (3), we can write this solution in real terms, as which highlights a self-modulation of the pulsation-type characteristic known as the Stoler-type transformation [61, 62] implying a complex form for this parameter. In Figure 2, we present this self-modulation phenomenon through time dependences, for various values of and .

It results that the release dynamics are explained by multifractal self-modulation (at different scale resolutions given by the maximum values of ) in the form of period doubling, quasiperiodicity, damped oscillations, intermittency, etc. (specificities dictated by the Fourier transforms of the analyzed quantities).

From a physical point of view, the fractal representation of the drug release mechanism showcased different release scenarios. These scenarios describe behaviors at a nondifferentiable (mesoscopic) scale from a local perspective encompassing all types of interactions between drug-polymer-medium. Our model predicts the drug release scenarios in the transition regime before the interactions reach an equilibrium drug release regime. In Figures 2(a) and 2(b), it can be seen a periodic release of specific quantities of the drug over a controllable period of time. The period doubling significance is important as each high concentration of drug is succeeded by a lower dose. In Figures 2(c) and 2(d), it can be seen that alternatives make these dynamics into more complex behaviors. It was observed that a useful damped periodical scenario do control the high doses of the specific drug. The intermittence release of drugs was based on a controllable modulating frequency (chosen as a control parameter). Finally, we achieved a modulated response of the drug release systems which can be controlled in both concentration and frequency of release. In the fractal paradigm, we can obtain various scenarios for the drug release mechanism, which help us to understand the behavior of drugs that may benefit from such complex scenarios.

The 3D and 2D (contour plot) dependences of on and , for a constant value of , are shown in Figure 3. In such a situation, the release dynamics involve multifractal self-modulation (at different nondifferentiable scale resolutions given by the maximum values of ) of the network dynamics. The overall evolution of our system is clearly presented in Figure 3. There we can see the damped drug release scenario seen for various concentrations of the particular drug (controlled in our model by ). For various nondifferentiable scale resolutions, we obtained different scenarios such as modulated release seen in Figure 3(b) while an intermittent release can be seen in Figure 3(c). The latter scenarios care complex release scenario and depend strongly on the physical properties of the polymer matrix and the biological release conditions. Finally, our model predicts the presence of unwanted regimes with a quasichaotic release. It is worth noticing that the complete chaotic dynamic is never reached, instead the adjustments made through the control parameter will force the system in a period doubling state. Therefore, even when the dynamic can seemingly be chaotic it can be easily rectified towards a more controllable state.

3. Results and Discussions

Three different formulations were obtained by combining a polymeric matrix-poly(vinyl alcohol boric acid) and diclofenac sodium salt anti-inflammatory drug, in different mass ratios from 70/30% to 90/10%. By casting the solutions into Petri dishes and slowly evaporating the water, self-standing flexible films were obtained, which were further characterized in terms of morphology, anchoring forces, and drug release rates.

3.1. Morphologic Characterization

The EDAX analysis confirmed the obtaining of three drug delivery systems with different amounts of drug, by the increase of the percentage of sodium and chlorine starting from D1 to D3 (Figures 4(a)4(c)). The corresponding SEM images indicated that DCF was homogenously encapsulated into the polymeric matrix as micrometric crystals, whose density depends on the mass ratio between PVAB polymer and DCF drug [39]. Therefore, in the case of the D1, which contains only 10% DCF, fewer micrometric geometric shapes were observed in comparison with the case of D3, in which the density was much higher (Figures 4(d)4(f)).

In order to have an insight about the morphology at nanoscale, the samples were further analyzed by AFM, which revealed the presence of nanograins in the case of all samples (Figure 5).

To confirm the crystalline nature of the encapsulated drug, the D1-D3 samples were analyzed by polarized light microscopy (Figure 6). As can be observed in Figure 7, the sample presented birefringent rectangular shapes characteristic to low molecular weight crystals [63], attributed to drug crystals, while the polymeric matrix did not present any birefringence, due to its amorphous nature. It can be observed that as the drug content increased, the size of the crystals decreased and their density increased, an aspect also seen for other systems based on PVAB matrix [3436]. The explanation is related to the increased viscosity of the system containing a higher drug amount, due to the higher density of physical forces developed between electron-deficient PVAB matrix and DCF containing rich electron heteroatoms. The higher density hindered the diffusion of the DCF molecules and thus the growth of crystals.

3.2. Anchoring of the DCF into the PVAB Matrix

The PVAB polymer has been chosen as polymer matrix for DCF, anticipating strong interactions which should anchor the drug, prompting thus a prolonged release. In order to investigate them, the FTIR spectra were recorded for the formulations and for their components (Figure 6). In the spectra of both DCF and PVAB, multiple bands appeared in accordance with their structures as follows: both the DCF drug and PVAB polymer presented between 3750 and 2750 cm-1 two bands at 3386 and 3262 cm-1 in the case of DCF and a band with two maxima at 3277 and 3177 cm-1 in the case of the polymer, which correspond to NH and OH stretching vibrations, respectively [64]. In the fingerprint region, in the DCF spectrum, bands were observed corresponding to the stretching vibrations of the C-Cl bond at 743 cm-1 and vibrations of the carboxylic units at 1572 cm-1. The PVAB polymer presented also multiple peaks, the more significant ones being at 1330 and 1375 cm−1, due to the B-O bond stretching vibration and at ~1188 cm−1, due to B-OH deformation. In the spectra of the formulations, important changes appeared in both peak position and number. Therefore, in the 3750-2750 cm-1 spectral region, the two maxima from the polymer spectrum shifted to higher wavenumbers, indicating significant changes in the hydrogen bonds, indicating new interactions with the DCF molecules. The bands corresponding to B-OH deformation and C-Cl vibrations disappeared from the formulation spectrum indicating their implication in new strong interactions, the most probably coordinative bonds of the electron-deficient boron atom with electron-rich atoms of diclofenac (e.g. nitrogen, chloride). These spectral modifications clearly indicate that due to the structural peculiarities of both PVAB and DCF drug, strong interactions between them occurred in the formulations, the DCF molecules being “anchored” in the PVAB matrix by strong physical interactions, similar to liquid crystals previously reported by us [3436].

3.3. The Release Kinetics

The release kinetics of the DCF drug from the formulations was investigated in biomimetic conditions, in PBS, at a pH of 7.4 and 37°C, according to a procedure previously described in the literature (Figure 8) [6].

As it could be observed, the samples released the encapsulated drug with different rates, which depend on the amount of DCF used in their preparation. Therefore, the D1 sample released the drug with a slower rate due to the fact that it contains only 10% drug in 90% PVAB polymer. This mass ratio between the components very probably led to a better dispersion and anchoring of the drug molecules in the polymeric matrix, influencing in a positive manner the drug release rate. When the amount of the drug increased, up to 30%, the releasing rate increased, leading to an increased burst effect. In all three cases, regardless the amount of the encapsulated drug, the release profile proceeded in two stages: in the first 4 hours, a burst effect was observed, while in the next hours (8-144 hours), the drug release rate decreased, leading to a sustained release (Figure 8, Table 1).

The kinetic data obtained from the release experiments were fitted with three mathematical models as follows: zero-order, Higuchi, and Korsmeyer-Peppas, on both release stages (Figure 9). The high values obtained for the correlation coefficient (>0.98) in the first stage for the three mathematical models revealed that the DCF release from the drug delivery systems was a process controlled mainly by two factors—drug dissolution and its diffusion through the polymeric matrix. In the second stage, the kinetic models fitted pretty well only for the D1 sample. This indicates that the drug release was mainly governed by the matrix erosion in the second-stage process which was also confirmed with the naked eye.

To assess the performances of the PVAB matrix for drug encapsulation, compared to other polymers, we realized a literature survey on other drug delivery systems based on biopolymers and DCF. As can be seen in Table 2, PVAB has the ability to encapsulate bigger amounts of drug, which were released in a more prolonged manner. This confirmed that the PVAB ability to strongly anchor DCF is beneficial to encapsulate a high amount of drugs and further to favor its prolonged release.

In our mathematical model, the drug release mechanism was developed in the framework of the fractal theory of motion in which dynamics were described by means of fractal functions. More precisely, the various interactions of the polymer-drug-medium as a complex system are described through various Riccati-type “regimes” (multifractal Riccati-type equations). We obtained several release kinetic scenarios at a nondifferentiable (mesoscopic) scale. These scenarios mimic the usual chaos transition scenarios: period doubling, quasiperiodicity, damped oscillations, intermittency, and quasichaoticity. Moreover, we can find specific oscillation modes of the complex polymer-drug networks, where the kink-type solution corresponds to a steady macroscopic behavior. The kink solution will be used to validate the above-presented experimental data.

Indeed, the multifractal Riccati-type equation (1) has the bounded solution:

If we set , then a relation that can describe the standard dynamics of release at various differentiable scale resolutions, through a convenient choice of variables and parameters (Figure 10). It results that standard release dynamics are explained by temporary kink-type multifractal behaviors (for details on nonlinear kink solutions, see [6769]). We notice also that such a particular solution of our model fits well the experimental data. These results underline the generality of the fractal approach and its wide range in terms of drug release scenarios, strongly depending on the physical properties of the polymer matrix and medium properties.

4. Conclusions

A series of formulations based on poly(vinyl alcohol boric acid) and diclofenac sodium salt were prepared, and the drug release kinetics was investigated by three traditional mathematical models and a new theoretical multifractal model proposed by us. The drug was encapsulated into the polymeric matrix as micrometric and submicrometric crystals, anchored by physical forces. This assured a prolonged release of the DCF in two stages, reaching a period of 5 days in the case of the sample which contained the 10% drug. The fitting of the three traditional mathematical models on experimental release data suggested that drug release was governed by the DCF dissolution and diffusion in the first stage (four hours) and by the matrix erosion in the second one (five days). The modeling of the drug release by the fractal theory of motion evidenced several transient scenarios for a nondifferentiable (mesoscopic) scale in the form of period doubling, quasiperiodicity, damped oscillations, intermittency, and quasichaoticity. Usually, such behaviors were assimilated to the chaos transition scenarios of the complex system dynamics. These scenarios are specific oscillation modes of complex polymer-drug networks, where the kink-type solution corresponds to a steady macroscopic behavior. In such a situation, through a convenient choice of variables and release parameters, the theoretical model was validated on the basis of experimental data.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The research leading to these results has received funding from the Romanian National Authority for Scientific Research, MEN-UEFISCDI grant, project number PN-III-P1-1.2-PCCDI-2017-0569, contract no. 10PCCDI/2018. Access to research infrastructure developed in the “Petru Poni” Institute of Macromolecular Chemistry through the European Social Fund for Regional Development, Competitiveness Operational Programme Axis 1, Project InoMatPol (ID P_36_570, Contract 142/10.10.2016, cod MySMIS: 107464) is gratefully acknowledged.