Three-bottle exponential decay infusion : A novel intravenous infusion strategy for drugs which conform to three-compartment model with observations on propofol kinetics

The aim of this study is to develop a novel three-bottle exponential decay infusion method to achieve the steady state concentration rapidly for drugs which conform to three-compartment model. The mathematical model and dosing regimen of three-bottle exponential decay infusion were developed. The input concentration of this method is a multi-exponential decay variable. With this method, blood concentration increased very rapidly to plateau value in the early phase of infusion and the target plateau concentration can be achieved much more quicker (approximately 10 to 35 times) than conventional constant rate infusion. This procedure was verified in rabbits and dogs and compared with conventional constant rate infusion with observations on propofol kinetics. This new quantitative methodology may provide a novel intravenous infusion strategy for drugs which conform to three-compartment model and can be used for the program design of infusion pump. Article Info Received: 21 September 2012 Accepted: 22 September 2012 Available Online: 12 October 2012 DOI: 10.3329/bjp.v7i4.11908 Cite this article: Huang XH, Xie HT, Huang JH, Sun RY. Three-bottle exponential decay infusion: A novel intravenous infusion strategy for drugs which conform to three-compartment model with observations on propofol kinetics. Bangladesh J Pharmacol. 2012; 7: 223-35. This work is licensed under a Creative Commons Attribution 3.0 License. You are free to copy, distribute and perform the work. You must attribute the work in the manner specified by the author or licensor. In our previous study, we successfully developed a two -bottle exponential decay infusion method for drugs which conform to two-compartment model and verified it in dogs and healthy volunteers (Song et al., 1983). In the present study, a quantitative methodology on the basis of the three-compartment open-system model is described by which steady state concentration can be achieved much more rapidly than conventional CRI method by exponentially changing the infusion concentration during the infusion process. This quantitative method may provide a novel intravenous infusion strategy for drugs which conform three-compartment model and can be used for the program design of infusion pump. Materials and Methods Mathematical model and dosing regimen Figures 1 show the mathematical model of the threebottle exponential decay infusion. Related equations of dosing regimen and the concentration-time courses are described in Table I. Mathematical proof of these equations was attached in appendix. For the purposes of this discussion, it is postulated that the body behaves as a three-compartment open-system model with components shown in Figure 1. The model where elimination occurs from a central compartment which is reversibly connected to two peripheral compartments. The rates of transferring between compartments and the rate of drug elimination from the central compartment are considered to follow first-order kinetics. When carrying out this method, three infusion bottles (A, B and P) are connected in series. The drug concentration (Ca) of bottle A were set equal to that of the conventional CRI method. The intial concentration of bottle B (Cb0) is set higher than Ca and the intial concentration of bottle P (Cp0) is set higher than Cb0. The concentrations, and drug volumes and amounts of these three infusion bottles can easily calculated according to the equations in Table I. And this dosing regimen can be easily written in a computer program. From equation 1, we can see that if the infusion rates of these three bottle remain constant, the input concentration of TBEDI is a multi-exponential decay variable. Initially, the concentration in bottle P (the highest concentration), was directly infused to the central compartment. The concentrations in bottle B and P are diluted in an exponential way all the time. As the infusion time tends to infinite or lasts long enough, the concentrations in the three bottles can be approximately equal to that of Ca (the lowest concentration). At the very beginning of TBEDI, the blood concentration can increased very rapidly to plateau concentration and maintained at this level as long as desired. Experimental procedure Rabbits (♀ and ♂, n = 6, weighing 2.0 ± 0.2 kg) and dogs (♀ and ♂, n = 3, weighing 5.0 ± 0.2 kg) were obtained from the Animal Breeding Center of Wannan Medical College (Anhui Province, China). The study protocol was approved by the Ethical Committee for Animal Experimentation of Wannan Medical College. Propofol emulsion for intravenous injection (10 mg/1 mL, 200 mg per vial) was purchased from Zeneca company. 5% Glucose was obtained from the Pharmaceutical Department of Yijishan Hosptical Affiliated to Wannan Medical College. 224 Bangladesh J Pharmacol 2012; 7: 223-235 Figure 1: Diagram of three-bottle exponential decay infusion. Ro, kb and kp are the infusion rate constants of bottles A, B and P, respectively; d is the infusion rate; Xa, Xb0 and Xp0 are the drug amounts of bottles A, B and P, respectively; Vb and Vp are the drug volumes of bottles B and P, respectively. X1, X2 and X3 are the drug amounts of central compartment and two peripheral compartments, respectively. K12 and K13 are distribution rate constants from central compartment to the peripheral compartments. K21 and K31 are distribution rate constants from the peripheral compartments to the central compartment; K10 is elimination rate constant for loss of drug from the central compartment


Introduction
Constant rate infusion (CRI) is very commonly used in clinical practice.For drugs with a long elimination halflife, the time at which the steady state is achieved might be very long.To overcome this drawback, a loading dose is frequently used to achieve the steady state as early as possible.It is easy to calculate the loading dose and infusion rate for one-compartment model drugs, but such calculations are rather complicated for two-or three-compartment model drugs (Gibaldi et al., 1982).Some scholars have proposed different dosing regimens to solve such problem for two-compartment drugs (Kruger-thiemer et al., 1968;Wagner et al., 1974;Boyes et al., 1971;Mitento et al., 1972;Zhu et al., 1980;Sun et al., 1983;Sun et al., 2004;Xie et al., 2011).Each method has its own advantages and disadvantages (Xie et al., 2011).However, the pharmacokinetic features of some drugs, such as propofol (Puri et al., 2012;Ye et al., 2012), raltitrexed (Blair et al., 2004;Clarke et al., 2000), thiopental (Russo et al., 1998), midazolam (Hardmeier et al., 2012), and digoxin (Krishna et al., 2011) etc. should be effectively described with the use of the three -compartment open-system model.For such drugs, the target plateau concentrations often can not be achieved rapidly enough by conventional CRI method, which may delays the therapeutic effects (van den Nieuwenhuyzen et al., 2000;Xie et al., 2011).According to the literature search using MEDLINE, there is hardly satisfactory quantitative methodology to address this problem.
In our previous study, we successfully developed a two -bottle exponential decay infusion method for drugs which conform to two-compartment model and verified it in dogs and healthy volunteers (Song et al., 1983).In the present study, a quantitative methodology on the basis of the three-compartment open-system model is described by which steady state concentration can be achieved much more rapidly than conventional CRI method by exponentially changing the infusion concentration during the infusion process.This quantitative method may provide a novel intravenous infusion strategy for drugs which conform three-compartment model and can be used for the program design of infusion pump.

Mathematical model and dosing regimen
Figures 1 show the mathematical model of the threebottle exponential decay infusion.Related equations of dosing regimen and the concentration-time courses are described in Table I.Mathematical proof of these equations was attached in appendix.For the purposes of this discussion, it is postulated that the body behaves as a three-compartment open-system model with components shown in Figure 1.The model where elimination occurs from a central compartment which is reversibly connected to two peripheral compartments.The rates of transferring between compartments and the rate of drug elimination from the central compartment are considered to follow first-order kinetics.
When carrying out this method, three infusion bottles (A, B and P) are connected in series.The drug concen-tration (Ca) of bottle A were set equal to that of the conventional CRI method.The intial concentration of bottle B (Cb0) is set higher than Ca and the intial concentration of bottle P (Cp0) is set higher than Cb0.The concentrations, and drug volumes and amounts of these three infusion bottles can easily calculated according to the equations in Table I.And this dosing regimen can be easily written in a computer program.
From equation 1, we can see that if the infusion rates of these three bottle remain constant, the input concentration of TBEDI is a multi-exponential decay variable.Initially, the concentration in bottle P (the highest concentration), was directly infused to the central compartment.The concentrations in bottle B and P are diluted in an exponential way all the time.As the infusion time tends to infinite or lasts long enough, the concentrations in the three bottles can be approximately equal to that of Ca (the lowest concentration).At the very beginning of TBEDI, the blood concentration can increased very rapidly to plateau concentration and maintained at this level as long as desired.

Experimental procedure
Rabbits (♀ and ♂, n = 6, weighing 2.0 ± 0.2 kg) and dogs (♀ and ♂, n = 3, weighing 5.0 ± 0.2 kg) were obtained from the Animal Breeding Center of Wannan Medical College (Anhui Province, China).The study protocol was approved by the Ethical Committee for Animal Experimentation of Wannan Medical College.Propofol emulsion for intravenous injection (10 mg/1 mL, 200 mg per vial) was purchased from Zeneca company.5% Glucose was obtained from the Pharmaceutical Department of Yijishan Hosptical Affiliated to Wannan Medical College.Rabbits and dogs were administered by three kinds of drug administration methods (i.e., intravenous injecttion, TBEDI and CRI) with a one-week wash-out period for each method.The experimental procedures were as follows: First, propofol in a dose of 10 mg/kg was injected in rabbits and dogs to study the kinetic behavior of this drug and obtain the individual pharmacokinetic parameters.The information obtained was used to calculate the dosing regimen of TBEDI and CRI.Blood samples were collected before and at 1,2,3,4,5,6,8,10,15,20,25,30,35,40,45,50,60 min following the bolus intravenous injection of propofol.Then two infusion methods, TBEDI and CRI, were carried out respectively and compared.When using TBEDI method, The different concentrations and drug volumes for the bottle A, B and P were calculated according to the equations in Table I based on individual pharmacokinetic parameters.When using CRI method, only one bottle (bottle A) was used.Blood samples were drawn according to the exact time schedules (0,2,4,6,8,10,12,15,18,20,25,30 min for TBEDI,and 0,2,6,10,15,20,25,30 min for CRI.).When carrying out and comparing these two infusion methods, the desired steady state plasma drug concentration was set at 3 μg/mL.The infusion device was connected to a constant flow pump and the flow rate was adjusted at 1 mL/min.During the above procedures, the ampule of propofol was shaked for at least 1 min before use and the appropriate dose was aspirated accurately according to the body weight of the rabbit or dog.The rabbits or dogs were fixed properly on the operating table and all the animals had venous canula inserted for sampling.Venous blood propofol concentrations were determined by HPLC.DAS 3.0 program was used to analyze the pharmacokinetic data.DAS is the most widely used

Equations of dosing regimen for three-bottle exponential decay infusion
The time course of plasma concentration in bottle P : The time course of plasma concentration in central compartment : When , Initial drug concentration in bottle B: Initial drug concentration in bottle P: Drug volume of bottle P: Initial drug amount in P bottle: Xp0 = Cp0×Vp ( 10 ) pharmacokinetic software in mainland China.
Blood samples were collected in heparinized centrifugation tubes and stored at 4°C until assay.After adding 0.1 mL blood sample and 0.1 mL acetonitrile in 1 mL plastic centrifugation tube, the mixture then was blended on an inversion mixer for 10 sec and centrifugated at 15,000 g for 5 min.A 100 uL volume of the supernatant was taken for HPLC analysis.HPLC system (Waters Assoc.)consisted of a Waters 515 MS system controller, a Waters 717 plus automatic sample injector and a 474 fluorimetric detector.Chromatographic conditions: The HPLC mobile phase consisted of acetonitrile and distilled water (65:35, v/v).pH was adjusted to 4.0 using acetic acid.Chromatographic separation was carried out by a 200 x 4.6 mm column packed with silane particle (C18) which was used at ambient temperature of 30°C.The excitation and emission wavelenths were set at 216 and 370 nm respectively.The linearity of the curve in the range of 5 to 5,000 ng/mL and the relative standard deviations in this concentration range were less than 10%.The limit of quantification proved to be 0.4 mg/L.Plasma drug concentrations were calculated according to the peak area of the samples.

Results
The plasma concentration-time profiles of propofol following the bolus intravenous injection are shown in Figure 2. The pharmacokinetics of propofol were best fitted to a three-compartment open model.The main parameters are summarized in Table II.
Table III shows the dosing regimen of the present study.The different concentrations and drug volumes for the bottle A, B and P were calculated and dispensed according to the equations in Table I and individual pharmacokinetic parameters in Table II.The results are presented graphically for each animal.Figure 3 and 4 show the predicted and measured propofol concentration-time profiles of TBEDI and conventional CRI methods in rabbits and dogs.In general, the observed plasma values agree well with the predicted concentrations.The results demonstrate that the time of obtain 90 per cent plateau concentration is approximately 4 to 6 min by TBEDI for all the subjects.On the other hand, for conventional mono-bottle infusion, it would require up to 30 min to reach 50 per cent and estimated approximately 59 to 143 min to reach 90 per cent of plateau concentration.The target concentration was maintained for 30 min in our study and could be as long as desired.The result also demonstrates that the optimal anesthetic target, or adequate depth and satisfactory quality of anesthesia were achieved with maintain concentration of between 2.5 and 4.7 μg/mL of propofol in rabbits and dogs.With this method, the speed of achieving the 90 per cent plateau concentration was more than 10 to 35 times faster than that of the conventional constant rate infusion in the present study.

Discussion
Many therapeutic agents can be applied more predictably and effectively if a constant plasma concentration is maintained.Many clinical situations often demand that the steady-state concentration be reached as rapidly as possible.It is often desirable in such cases to administer an intravenous loading dose just before starting the intravenous infusion.For two-compartment drugs, the calculation of loading dose and fusion rate is relatively complicated.Some theoretical and experimental studies have been published attempting to provide possible solutions.Each method has its own advantages and disadvantages.For example, Boyes's method is relatively safe, but this approach has the drawback of a minimum concentration which may be well below that required, followed by a slow return to plateau values (Boyes et al., 1971;Sun et al., 2004).With the method developed by Mitenko (Mitenko et al., 1972), the concentration during the early infusion phase is very high, sometimes even the 1.5 to 5 times to target concentration, which may often exceed the safe limits of  drugs with small therapeutic indexes.Some researches proposed regimens similar to these two methods, but still could not solve the problem of avoiding the high concentrations during early infusion phase (Wagner et al., 1974;Chau, 1977;Zhu, 1980).In our previous study (Song et al., 1983), we successfully developed a doublebottle exponential decay infusion method with observations on gentamicin kinetics.With this method, the rate of infusion is exponential with time and steady-state concentration can be obtained very rapidly for one-or two-compartment models.Actually, two-bottle exponential decay infusion can be seen as the particular case of the TBEDI method described in this study.This procedure was verified in 12 dogs and 8 normal male volunteers and was more safe compared with Boyes', Mitenko' and Zhu' methods.
However, there are a number of effective but highly toxic drugs that do manifested multi-exponential characters which are too simplistic to explain using one-or two-compartment model.In particular, it has been pointed out that the pharmacokinetic characteristics of many anesthetic drugs can be effectively described with the use of three-compartment open-system model (Nieuwenhuyzen et al., 2000).In such cases, there is considerable potential for inadequate or excessively high drug concentration in the blood in clinical practice.
On the other hand, with conventional CRI infusion method, target plateau value can be achieved slowly only after 4-5 elimination half-lives have passed which is clearly unsuitable for anaesthetic drugs in many clinical situations.According to the literature search using MEDLINE, there is hardly satisfactory quantitative methodology to address this problem.
A possible solution is offered in this paper.This method can be useful for drugs which conform to threecompartment open model, especially those with long half-lives and small therapeutic index.Only a programmable hand calculator is needed to calculate the Take rabbit R1 as an example: Va dosing regimen.With this method, a specified steadystate plasma concentration can be rapidly obtained by combining three bottles in series with exponential decay of drug concentration.The time needed to achieve the steady state can be more than 10 to 35 times faster than that of the conventional constant rate infusion.This method is especially suitable for anesthetic agents, such as propofol and alfentanil.Optimal therapy with such drugs requires appropriate dosing regimen to safely obtain the desired clinical effects In clinical application, to ensure the liquid medicines in the bottles B and P be blend completely, a electromagnetic stirring device could be used.And constant rate infusion pump could be used to control the infu-sion rate to be remained constant and ensure accurate delivery.However, according to our experi-ences in previous study (Song et al., 1983), even if only the clamps are used to control the infusion rate in clinical practice, the satisfactory result could also be obtained.
Actually, even if the liquid medicines in bottles B and P can not be diluted completely, the concentrations during the early infusion phase of this method are still much lower than those of Mitento' and Zhu' method (Song and Sun, 1983;Xie et al., 2011).It should be also noted that, when applying this method, in the early phase of infusion the drug amount into the body is larger than that of CRI, and drug concentrations rise relatively quickly.So, for intense irritant drugs and those with large toxic effects on cardiovascular or central nervous system, this method should be used with caution.In addition, The three-bottle exponential decay infusion scheme may also be applied to maintain the concentration of drug in a particular compartment of a multi-compartment model (e.g., the compartment representing the site of drug effect), but the mathematic models may be more complicated.

Conclusion
The three-bottle exponential decay infusion quantitative method, with which drug concentration plateaus can be achieved much more rapidly compared with conventional constant rate infusion, may provide a novel intravenous infusion strategy for drugs which conform to three-compartment model and can be used for the program design of infusion pump.

Figure 1 :
Figure 1: Diagram of three-bottle exponential decay infusion.Ro, kb and kp are the infusion rate constants of bottles A, B and P, respectively; d is the infusion rate; Xa, Xb0 and Xp0 are the drug amounts of bottles A, B and P, respectively; Vb and Vp are the drug volumes of bottles B and P, respectively.X1, X2 and X3 are the drug amounts of central compartment and two peripheral compartments, respectively.K12 and K13 are distribution rate constants from central compartment to the peripheral compartments.K21 and K31 are distribution rate constants from the peripheral compartments to the central compartment; K10 is elimination rate constant for loss of drug from the central compartment

Figure 2 :
Figure 2: Measured plasma propofol concentrations in dogs following the bolus intravenous injection at 10 mg/kg.(A: six rabbits; B: three dogs)

Figure 4 :
Figure 4: The comparison of predicted and measured propofol concentration-time profiles of TBEDI and MBI methods in each experimental animal.R1 to R6 and D1 to D3 denote six rabbits and three dogs, respectively.The symbols are defined in Figure 3

Figure 4 :
Figure 4: The comparison of predicted and measured propofol concentration-time profiles of TBEDI and MBI methods in each experimental animal.R1 to R6 and D1 to D3 denote six rabbits and three dogs, respectively.The symbols are defined in Figure 3 (Cont.)