Pharmacokinetic Delivery and Metabolizing Rate of Nicardipine Incorporated in Hydrophilic and Hydrophobic Cyclodextrins Using Two-Compartment Mathematical Model

The dispersion routes of cyclodextrin complexes with nicardipine (NC), such as hydrophilic hydroxypropyl-β-cyclodextrin (NC/HPβCD) and hydrophobic triacetyl-β-cyclodextrin (NC/TAβCD), through the body for controlled drug delivery and sustained release have been examined. The two-compartment pharmacokinetic model described the mechanisms of how the human body handles with ingestion of NC-cyclodextrin complexes in gastrointestinal tract (GI), distribution in plasma, and their metabolism in the liver. The model showed that drug bioavailability was significantly improved after oral administration of cyclodextrin complexes. The mathematical significance of this study to predict nicardipine delivery using pharmacokinetic two-compartment mathematical model with linear ordinary differential equations (ODE) approach represents a valuable tool to emphasize its effectiveness and metabolizing rate and diminish the side effects.


Introduction
Formulation of drugs with different drug carrier materials to control and improve drug release is a rapidly growing area in the pharmaceutical sciences. One example of this application is the usage of cyclodextrins (CDs) as excipients for a sustained release of various chemical compounds after their oral administration [1,2]. Various mathematical models were proposed to study the CD release profiles and match them to the experimentally-determined kinetic data [3,4]. Among those methods were Weibull distribution along with zero-order, Higuchi, and Korsmeyer-Peppas that best describe drug release phenomena with major applicability to the versatile experimental data [5,6].
In this context, the experimental and theoretical approaches are increasing importance of the investigation of complexation mechanisms and pharmacokinetic properties of CD-drug conjugates, as well as their delivery mechanisms in the different compartments of the human body [7][8][9].
Nicardipine (NC) is a calcium-channel blocker used in the treatment of different cardiovascular and cerebral diseases; these diseases are mainly hypertension, subarachnoid hemorrhage, and acute strokes [10,11]. This drug was used as a model to simulate the limited bioavailability effects (up to 40%) and rapid biotransformation in the liver with a very short half-life time [12].
The study of Fernandes and coauthors had reported that nicardipine release rate might be controlled by the formulation of this drug with hydroxypropyl--cyclodextrin (NC/HP CD) and triacetyl--cyclodextrin (NC/TA CD) derivatives [1]. The nicardipine delivery can be modeled by using two-compartment pharmacokinetic model developed previously in the works of Spitznagel and Shonkwiler and Herod [13,14]. This approach was already implemented successfully to investigate the transfer and elimination of lead compound from blood and bone tissues and assess the influence of bone-lead depots [15].
The current study attempts to predict the nicardipine release rate based on its formulation with hydrophilic and hydrophobic cyclodextrins in the gastrointestinal tract (GI) and plasma using pharmacokinetic two-compartment mathematical model with linear ordinary differential equations 2 The Scientific World Journal (ODE) approach. Additionally, the metabolizing rate analysis was conducted, and the results were compared to the experimentally-acquired data using the referred model.

Computational Methods
The experimental pharmacokinetic parameters, such as maximum plasma concentration ( max ), the time to reach max ( max ), and the elimination half-life ( 1/2 ) together with dissolution profiles, were taken from in vivo study of Fernandes and coauthors [1]. The authors analyzed the drug release rate in rabbits ( = 6) using high performance liquid chromatography with UV detection.
To simulate the nicardipine dose-dependent release and its metabolic degradation, we used computer linear algebra algorithms implemented in the Maple 15 software package (Maplesoft GmbH, Aachen, Germany). The pharmacokinetic two-compartment model included two different compartments: the first compartment ({ ( )} or / ) describes the drug concentration in the GI, and the second compartment ({ ( )} or / ) describes its concentration in plasma in dosage-dependent manner (long-range) by following differential equation: where ( ) is the drug dosage in g/mL, DR is the dissolution rate, and 1/2 value is the half-life parameter either for the GI or plasma. The dissolution rate was calculated from experimentally-determined profiles in the acidic medium (pH 1.2) for nicardipine and its complexes as where complex and drug are the dissolution time in minutes either for nicardipine or its complexes in the acidic environment. The Runge-Kutta method was used to solve numerically the ODE by using a trial step at the midpoint of an interval to cancel out lower-order error terms [16,17]. This method was already proven for a variety of pharmacokinetic simulations as a very robust technique [18][19][20].
The ODE approach was also applied to solve the firstorder irreversible series of drug translocations and modifications to obtain symbolically the exponential matrix as a function of time [21]; the independent variable uses the following scheme: where GI and Plasma are the estimated rate constants either in the GI or plasma. Therefore, in the Maple software package for a linear system of simultaneous coupled first-order differential equations, is × 1 matrix, is a × matrix, and exp( ) is a × matrix, and so on [21]. The governing equation for this reaction can be represented as where and are the concentrations for analyzing species either in the GI or plasma (short-range). The value was set to 1.0 ng/mL as the initial concentration. The concentration of metabolite species at any given time ( ( )) is The estimated rate constants (Table 1) were calculated from experimentally-determined 1/2 values for nicardipine and its complexes as For a compartment model, we used the diagonal terms of the matrix and the column sums, which were negative or zero. Therefore, the solution has a tendency to be driven by drug dosage ( ( )).

Results and Discussion
According to our model, nicardipine and its CD complexes were taken orally at periodic time intervals to provide a pulse of dosage delivered to the gastrointestinal tract. From there, the substances entered into the bloodstream in concentration-dependent and dissolution-dependent manner. The drug was assumed to be distributed into the bloodstream compartment after its partial or complete elimination from the GI implying the bioexponential drug disposition and distribution.
To establish the absorption rate as the 1/2 value of NC, NC/HP CD, and NC/TA CD substances, we specified dissolution-dependent drug dose taking into account that the NC dissolution was observed within 180 min due to high partitioning coefficient (log = 3.94) [22], after 15  Hazai et al. [23] ensured a significant enhancement of the low NC dissolution and solubility because of surfactantlike properties of the carrier diminishing drug-medium interfacial tension [24,25]. The starting ( ) value was calculated for each substance as ( )/12 for NC and ( )/32 for NC/TA CD, following the logic that these substances have 12-and 32-fold decrease in dissolution in comparison to hydrophilic cyclodextrin formulation. However, the final dosage value for NC/HP CD was actually devised from precalculated concentration (0.49 g/mL) using initial ( ). Further, this value was diminished to ( )/12 to fit its predicted max parameter to the experimentally-determined results (39.65 ± 1.11 g/mL).
Since the maximal plasma levels of NC were already observed after 30 min of oral administration, we used a 30 min half-life in the GI tract for nicardipine such as a 1/2 of 2 ln(2). For the NC/HP CD and NC/TA CD complexes, we used empirically determined 1/2 values of ln(2)/4 and ln(2)/2.2 since these compounds have much longer half-lives in the GI tract in comparison to NC [24]. The plasma 1/2 parameters were ln(2)/1.6 for NC, ln(2)/5.49 for NC/HP CD, and ln(2)/7.38 for NC/TA CD according to their experimentally-determined half-life rates [1]. We set (0) = (0) = 0 for initial conditions to adjust the starting point where no substance is present either in the GI tract or circulatory system. The Runge-Kutta method was implemented to produce NC concentration. Interestingly, the NC concentration in the bloodstream was predicted to be oscillating every 6 hours with a time to achieve the maximum plasma level (54.6 ng/mL) and the GI level (60.1 ng/mL) at half-life times of 30 min in the GI and 1.6 hours in the plasma, respectively. On the other hand, the NC concentration in the GI was superimposed on a gradually increasing level of the plasma concentration starting from a 6-hour time point.
The NC/HP CD concentration in the bloodstream was predicted to be oscillating every 24 hours to achieve the maximum GI level (81.5 ng/mL) and plasma level (41.5 ng/mL) at half-life times of 4 hours in the GI and 5.49 hours in the plasma (exp. value in plasma: 5.49 ± 0.17 hours). The NC/TA CD concentration in the bloodstream was predicted to be oscillating also every 24 hours to reach the maximum plasma level (21.6 ng/mL) and the GI level (28.9 ng/mL) at half-life times of 2.2 hours in the GI and 7.38 hours in the plasma (exp. value in plasma: 7.38 ± 1.61 hours). On the other hand, the complexes did not show sharp peak plasma concentrations similar to those observed by other researchers [1,26]. All predicted maximum plasma concentrations were in agreement with experimentally-derived data for analyzed substances: 69.64 ± 8.57 ng/mL for NC, 39.65 ± 1.11 ng/mL for NC/HP CD, and 18.25 ± 0.74 ng/mL for NC/TA CD, respectively. The GI and plasma concentration levels for NC and the GI concentrations for NC/HP CD and NC/TA CD after oral administration showed sharp peaks and abrupt decrease of concentration, especially for NC due to a short elimination half-life (1.6 hours). The NC drug kinetics in plasma were already previously described by a similar model revealing a very short 1/2 of distribution phase [27].
The plasma concentration levels for both formulations were significantly prolonged to 24 hours in retention time because of the improvement of NC dissolution profile upon complexation with cyclodextrins in simulated intestinal or gastric fluids with higher stability constants [24].
In the drug designing process, it is important to keep the concentration as a uniform function to reach the limit cycle as fast as possible. Therefore, the half-life parameters in the GI or plasma can be adjusted to prolong and optimize the "time release" mechanism. The plasma 1/2 value in our model represents characteristic properties of the analyzed compound. Because the { ( )} value is ( )and ( )dependent as the level of the drug in the circulatory system, it should take longer time period for a drug to be effective.
Subsequently, all the studied substances were metabolized in the liver as a result of first-order series of parallel irreversible reactions represented by a linear system of coupled ODE [28].
The dissolution profile of NC/TA CD was determined as 4.38-fold less calculated from the GI(NC) / GI(NC/TA CD) ratio than this profile for the NC drug alone. Therefore, initial concentration parameter ( ) was set to /4.38. The parameter was set to 2. 5 for the same compound. The process of NC, NC/HP CD, and NC/TA CD conversion into their metabolites might be referred to as an initial value problem because the initial conditions of the dependent variables must have a known which determines how the dependent variable changes with time as shown in Algorithm 2.
The short-range pharmacokinetic profiles of plasma concentrations for analyzed compounds were precalculated to determine maximum plasma concentration ( pre max ) and the time to reach it ( pre max ), such as 0.58 ng/mL and 0.07 min for NC, 0.43 ng/mL and 0.4 min for NC/HP CD, and 0.18 ng/mL and 0.44 min for NC/TA CD using experimentally-determined rate constants ( Table 1). The rate constants were optimized to produce the max values close to the experimentally-determined ones for the analyzed substances, such as GI(Plasma) /7.14 for NC, GI(Plasma) /4.88 for NC/HP CD, and GI(Plasma) /5.38 for NC/TA CD, in comparison to those from the experimental study [1]. Further, these refined rate constants ( GI * and Plasma * ) were implemented in the model to fit both predicted and experimental parameters ( Table 2).
Since nicardipine undergoes extensive bio-transformation of the N-benzyl side-chain on position 3 of the molecule under liver control and oxidation to the analogous pyridine metabolite [29], it is important to assess and compare the NC metabolizing rates with its cyclodextrin formulated complexes. Under these conditions, the stability of the NCcyclodextrin complex plays an important role in extending a protection against the oxidation due to the fact that the 6 The Scientific World Journal > restart; Use linear algebra package "(linalg)" in Maple. > with(linalg); > with(plots); Enter the governing equations (5) and (6) as follow:  Tables 1 and 2 and the independent variable . Use the "pars" command to illustrate the results: > pars:={k1,k2}; Next, substitute these values for the kinetic constants into the concentrations store as , , and : > Ca:=subs(pars,ca); > Cb:=subs(pars,cb); > Cc:=subs(pars,cc); Plot the concentration profiles for these values of the kinetic constants: > plot(Ca,t=0..10); > plot(Cb,t=0..10); > plot(Cc,t=0..10); Algorithm 2 drug hydrolysis in CD complex is slower [30]. Therefore, to adapt the metabolic body conditions and reach the equilibrium, the absorption of NC in the GI was very rapid resulting in a decrease of its GI concentration level after 2 min, 10 min for NC/HP CD, and 25 min for NC/TA CD formulations (Figure 3(a)). The NC plasma concentration reached its maximum level (0.59 ng/mL) after 0.52 min. After this threshold, the NC plasma concentration was dropped down significantly reflecting drug rapid disappearance in plasma after 5 min. On the contrary, the NC/HP CD and NC/TA CD plasma concentrations were reached with their maximum levels (0.43 ng/mL and 0.18 ng/mL) only after 1.98 and 2.34 min, respectively (Figure 3(b)). Subsequently, they were also dropped down to undetectable levels after 15 and 20 min, indicating the sustained-release behavior The Scientific World Journal of these complexes with significant 1/2 parameters. This observation was experimentally confirmed as the plasma concentrations of the NC formulations were maintained at a relatively constant level for a long time (up to 24 hours) after their oral administration [1]. The plateau obtained in the NC metabolite concentration rate was more evident even after 5.5 min time in comparison to 15 min for NC/HP CD and 26 min for NC/TA CD (Figure 3(c)). It is interesting to note that the plateau was reached after significant delay time for the NC formulations, which suggested a time increase in their metabolizing rates corresponding to a greater extent of their plasma levels compared to the drug alone. Similarly, the NC plasma peak concentration did not appear after the administration of its complexes, indicating that the side effects may also be reduced.

Conclusion
The NC release rate is the half-life-dependent process that can be adjusted to prolong and optimize the "time release" mechanism by controlling the formulation ratio of its hydrophilic (NC/HP CD) or hydrophobic (NC/TA CD) complexes. Because { ( )} is ( )and ( )-dependent as the level of the drug in the circulatory system, it should take a longer period of time for a drug to be effective. With subsequent 8 The Scientific World Journal equal dose levels for specified chemical compounds and when taken orally, the periodic plasma concentrations were significantly prolonged for the NC formulations, especially for its hydrophobic form, and maintained thereafter. The plasma 1/2 value in our model represents characteristic properties of the analyzed compound, which affects not just plasma concentration levels of the drug but its metabolizing rates.