Document Type : Original Research Article

**Authors**

Chemical engineering department, Engineering faculty, University of Tehran, Tehran, Iran

**Abstract**

**Objective(s)****:** The quantitative calculation of release data is more facil when mathematics come to help. mathematically modeling could aid optimizing and amending the delivery systems design. Aim of this study is to find out the isoniazid release kinetic. **Methods:** In this work degradable temperature sensitive dextran-hydroxy ethyl methacrylate- poly-N-isopropyl acryl amide (Dex-HEMA-PNIPAAm) nanogels which were synthesized by UV polymerization were loaded by Isoniazid. The Isoniazid release amounts taken from in vitro studies at two different temperatures, below and upper lower criticalsolution temperature (LCST) were mathematically modeled to investigate the kinetic of drug release. Mathematically inquiry of release phenomenon of Isoniazid makes it easy to predict and recognize the influence of delivery device laying out parameters on release kinetic formulation. The modeling was performed using model dependent methods, such as zero order, first order, Higuchi, Korsmeyer- Pepas, Hixon and Crowel. **Results: **The best fitted model showing the highest determination coefficient (R2) was Korsmeyer-Pepas which means predominant release mechanism is controlled by diffusion. **Conclusions:** The Isoniazid release pattern of most samples was combination of swelling, diffusion and degradation.

**Keywords**

**INTRODUCTION**

Intelligent drug delivery systems, showing brilliant features as pre-determined time delivery, pre-planned rate and self controlling, are widely used in various areas of pharmaceutics world. Such delivery systems, mostly employing polymers to design new systems by help of nanotechnology, attracted large amount of attention in different medical applications. Size controlling, surface specification and release are among the most important requirements of nano-sized delivery systems(1–5).

Drug release includes steps of therapeutic agent freedom from drug which is then exposed to variety of mechanisms like absorption, distribution, metabolisms and exertion. Drug may be released through diffusion, dissolution, erosion, leaching and mixed of them(6–9).

Drug release occurs in varied way. Sometimes release is immediate in which drug releases immediately. Delayed and extended are two types of modified release that drug will release after a while by delay and in an extended period, respectively. Ultimately, pulsatile release a type of controlled release makes drug free at specific time intervals.

In drug release studies kinetic is a key factor because it can be related to the drug concentration in plasma which is the available facile parameter to determine. Applying kinetic model is valuable in clarifying release mechanism which is helpful in the drug release control and design(8–15).

There are vast numbers of kinetic models illustrating drug release which can be categorised in three distinct groups:

- Statistical approaches

ü Exploratory data analysis method

ü Repeated measurement design

ü Multivariate approach (ANOVA based model, analysis of variance)

- Model dependent approaches

ü Diffusion model (Fickʼs law)

ü Zero order model

ü First order model

ü Higuchi model

ü Korsmeyer – Peppas model

ü Hixon – Crowell model

ü Weibull model

ü Baker – Lonsdale model

ü Hopfenberg model

ü Gompertz model

ü Sequential layer model

- Model independent approaches

ü Ratio factors

ü Fit factors (similarity factor, difference factor, resign index)

The other way of categorisation is empirical and semi-empirical methods like Peppas and Sahlin; and mechanistic and empirical methods(8,10,12,13,15–26).

In this present study we will mathematically model the release kinetic of Isoniazid release which is an antibiotic using in curing diseases like tuberculosis from five series of Dex-HEMA-PNIPAAm nano/micro gels. The comprehensive synthesis method, characterization, drug loading and other worth full relative data were narrated in our previous publication (27,28). The model dependent category is chosen to evaluate and predict the effect of cross linking agent amount and polymers kinds on the type of release. At the moment obtained release data will fit by the most important models which have been introduced above and the results will relate to the structure of the samples.

** **

**METHODS**

Investigation of Isoniazid release kinetic from Dex-HEMA-PNIPAAm nanogels has been performed by applying nine kinetic models on release data obtained from five series of sample which were synthesized, characterized and loaded by Isoniazid previously in two different temperature, 25 ºC and 37 ºC. The determination is accomplished by comparison of determination coefficient; the higher the determination coefficient the most suitable model is achieved. Unexpurgated information of initial materials quantities and ratios, also Isoniazid release amounts were narrated before, here just a brief glance to the necessary data for comparison and interpretation of the samples behaviour and their structure is presented in (figures 1 and 2). As it is obvious a portion of release occurred as burst release at the beginning. Burst release may happen due to the high specific surface of Isoniazid carrier gels which is exposed to the media, instability of nanogels, exec amount of Isoniazid or just under the surface freedom. Release kinetic data were fitted to the following models which are explained completely and in accordance with regression coefficient of determination (R2) models were compared and studied.

*Zero order kinetic model*

Zero order release reveals an ideal delivery of drugs due to constant remaining of drug in blood plasma during delivery process. This release model points out the constant release regime which will be applicable in the case of heart or blood pressure maintenance, pain control and etc. This model can be represented by the following equation:

Where Q_{t }is the amount of drug release at time t, Q_{0 IS }initial amount of drug in solution which is often equal to zero, K_{0} is the zero order constant and finally t is time. To investigate the release kinetic, experimental in vitro cumulative amounts were plotted versus time.

*First order kinetic model*

First order kinetic model was first employed by Gibaldi and Feldman (1967) and Wagner (1969). Mentioned model also has been used for absorption and elimination of some drugs (1982). However first order kinetic model cannot clarify the mechanisms due to theories. In this case drug release will express by the following equation:

Where Q_{t }is the amount of drug release at time t, Q_{0 IS }initial amount of drug in solution, K is the first order constant and t is time. To investigate the release kinetic, experimental in vitro log cumulative percentage of drug remaining were plotted versus time. The result will be a straight line having the slope of -K/2.303.

*Higuchi model*

Among all delivery systems hydrophilic ones are mostly used. Such systems can swell and make a layer on the drug surface and as a consequence will prevent more water entrance and drug release. In this case swelling and erosion/degradation both are responsible of drug release rate control. Hence, often time dependent behaviour especially drug release reduction as the time pass, are observed in hydrophilic delivery systems. In 1961 and 1963 Higuchi proposed his models for water soluble and low soluble drugs in solid and semi-solid matrix. Simplified Higuchi model can be expressed by the following equation:

Where Q_{t }is the amount of drug release at time t, K_{H} is the Higuchi release rate constant and t is time. To investigate the release kinetic, experimental in vitro log cumulative percentage of drug remaining were plotted versus square root of time.

*Korsmeyer – Peppas model*

Korsmeyer and Peppas (1984) and Ritger and Peppas (1987) presented their equation analysing both Fickian and non-Fickian drug release from swelling and non-swelling polymeric matrix. Fickian release used in steady cases in which concentration does not change by time and it is just dependent to place. In the Fickian release in un- steady situations the concentration gradient in a special place is changing due to the time, in the other cases release follows non-Fickian regime. By fitting the first 60% of drug release to the model release mechanism can be fined out. The equation is represented as:

Where is the fraction of drug released at time t, K is the release rate constant, n is the release exponent that is shown in (table 1, 2) and t is time. In this model some factors are considered such as water diffusion in to the delivery system, swelling, gel formation, drug diffusion out, dissolution/erosion of polymeric matrix. Plot may be catches by log cumulative percentage drug release vs. Log time to investigate the in vitro release kinetic. When burst release is possible b parameter is added to the equation.

*Weibull model*

This model can be adapted to dissolution release and can be shown by the following equation:

Where Q is the amount of drug released at time t, is the total amount of drug released, T is the lag time and t is time.

*Gompertz model*

Proportional equation showing this model is as follow:

Where Q_{t} is the amount of drug released at time t, is the maximum amount of drug released, t is time, a is unreleased amount of drug at t=1 and b is release rate per unit of time.

*Baker- Lonsdale model*

Baker and Lonsdale model was derived from Higuchi model in 1974 for spherical matrix and is represented as follow:

Where Q_{t} is the amount of drug released at time t, is the amount of drug released at infinite time, K is the release rate and t is time. Plot respected to root of time inverse to study release kinetic of Isoniazid from nanogels.

* *

*Hopfenberg model *

Hopfenberg model is used for the cases of release from eroding polymer surface in which surface area keep constant during the process. For spherical geometry following equation express the mentioned model:

Where Q_{t} is the amount of drug released at time t, is the total amount of drug released, is the initial amount of drug in matrix, K is rate constant, t is time and a is system half thickness (radius).

*Hixson – Crowell model*

Hixson – Crowell model explain the release from systems in which surface area and diameter change. He suggested (1931) that the changes in diameter and area are relative to cube root of the volume and is described by the following equation:

Where Q_{t }is the remaining amount of drug in delivery particles at time t, Q_{0 IS }initial amount of drug in delivery particle, K is the rate constant and finally t is time. This model hypothesis is that the release rate is not controlled by the drug diffusion, but is controlled by particles erosion, dissolution and degradation. To compare the release kinetic of Isoniazid by this model cube root of drug percentage remaining in matrix plotted vs. Time(8–17,19–24,29).

** **

**RESULTS AND DISCUSSION**

Peer group comparisons of samples in each condition, temperature 25 ºC and 37ºC (considering and eliminating burst release), were performed and outcomes are revealed in (Tables 3 to 5).

(Table 3) shows the release data fitted to the chosen models at 25ºC. This temperature is below LCST temperature of the nanogels in which gels are in their swollen forms. The three most proper models fitted to data based on having closest R2to 1 were Korsmeyer-Pepas, Baker-Lonsdale, and Hopfenberg respectively. This order represent that the most release controlling mechanism in the temperatures below phase transition temperature is diffusion; as the water penetrate in gels structure make them swell and in this situation the drug particles diffuse out of swollen nanogels. Then after the degradation is the other mechanism of Isoniazid release from nanogels. The values of n (n < 0.5) reveal that diffusion pattern is of kind of Fickian diffusion.

At 37ºC nanogels are in the temperature upper LCST, hence there are in their deswelling form. The burst release is much sharper here than 25 ºC. As the gels were located in 37ºC they transited from swelling form to deswelling form and as a sequence large amount of Isoniazid came out suddenly. Noting the determination coefficient of models, only Korsmeyer-Pepas model fitting results gave R2> 0.9, which confirm that considering burst results the only release mechanism in nanogels are diffusion and hither also Fickian diffusion is dominant. (Table 4) shows the Determination coefficient of models considering burst release at 37 º C.

Just as mentioned above a reason for sharp burst was the phase transition. Here neglecting the first hour release that mean omitting burst release, again data were fitted to the selected models. This time more models exist having R2> 0.9. From the comparison of average amount of R2it is find out that the total first four best models are Baker-Lonsdale, Higuchi, Korsmeyer-Pepas and Hopfenberg, respectively. Now the first controlling mechanism is combination of swelling and degradation of spherical particles. Then after diffusion take parts in Isoniazid release. This means that here drug release does not obey one rule and both degradation and diffusion are responsible of release. For each sample individually the order is a bit unlike. Korsmeyer-Pepas, Baker-Lonsdale, Weibull, Korsmeyer-Pepas, Korsmeyer-Pepas proportionally belongs to sample 1 to 5 in an orderly manner. Diffusion is the main mechanism of Isoniazid from sample 1, in which Dex-HEMA/PNIPAAm ratio is 0.5, so because of small amount of hydrolysable ester group degradation is not dominant. More after cross linking agent value of sample 1 is equal to 1 that let the sample swell easily and drug molecules could travel out of nanogels. In sample 2 Dex-HEMA/PNIPAAm ratio and cross linking agent value are 1 and 2 respectively. In this case quantities of degradable parts are more in the structure and swelling characteristic of gel is fewer due to high cross links and as results proved Baker-Lonsdale best fitted for this sample that emphasis on swelling and degradation mechanism of Isoniazid release. (Table 5) shows the determination coefficient of models not considering burst release at 37 º C.

In samples 3 and 4 Dex-HEMA/PNIPAAm ratio is 2 but cross linking agent value in sample 3 is 1 and in sample 4 in 3. Both samples contain equal and the maximum of hydrolysable ester groups’ amounts but they release kinetic follow different regimes. In sample 3 Isoniazid release agree with Weibull model and sample 4 confirm the Korsmeyer-Pepas. As it is obvious the amount of cross links and as a consequence swelling ratio of nanogels can influence the release pattern. The more swelling happen the more degradation and the less diffusion observe. Finally sample 5 has less hydrolysable group in its structure and simultaneously the most cross links and so the less swelling which could lead to the diffusion. The data given in table convince the Fickian diffusion as a best linear fitted release kinetic model.

** **

**CONCLUSIONS**

Reviewing the various kinetic pattern open ups the correlation among drug release and structure of delivery systems mathematically. In this attempt, focuses were on Isoniazid release kinetic from Dex-HEMA-PNIPAAm nanogels. In vitro drug release investigation revealed controlled Isoniazid release containing a burst release at the beginning for all samples and all conditions. Considering that, these bursts were sharper upper phase transition temperature than ambient temperature. The Isoniazid release pattern of most samples was combination of swelling, diffusion and degradation. Results best fitted with Korsmeyer-Pepas model with Fickian diffusion.

** **

**ACKNOWLEDGMENTS**

The authors would like to acknowledge the financial support of University of Tehran Science and Technology Park for this research. This study was part of a PhD thesis under Tehran University.

** **

**CONFLICTS OF INTEREST**

The authors declare that there are no conflicts of interest.