`
`BIOPHARMACEUTICS & DRUG DISPOSITION, VOL. 11, 93—l05 (1990)
`
`ABSORPTION OF DOXYCYCLINE FROM A
`CONTROLLED RELEASE PELLET
`FORMULATION:' THE INFLUENCE OF FOOD ON
`
`BIOAVAILABILITY
`
`DESMOND B. WILLIAMS*, WILLIAM J. O‘REILLYT, GARTII I3OEHM* AND MICHAEL J. STORY*
`
`*F.H. Faulding and Co. Ltd, Adelaide, South Australia and
`1‘School of Pharmacy, South Australian Institute of Technology
`
`ABSTRACT
`
`A three—way crossover study was performed to compare the bioavailability of a new
`pelletised doxycycline product administered either with food or without food and a
`reference product taken without food.
`Four different methods were used to calculate pharmacokinetic parameters from the
`data. The sums of squares, Akaike’s Information Criterion (AIC), and the ranges for the
`parameters obtained were used for comparison. Good fits to the data were obtained when
`all four methods were used, each with a lag time. The two compartment open model was
`the most efficient method for describing the data. The one compartment open model was
`the least efficient, particularly with respect to predicting the peak concentration of doxy-
`cycline in plasma. All the models gave similar rank order results with respect to bioavail—
`ability differences between the three treatments.
`Analysis of the data by different methods suggests that pelletised doxycycline is
`bioequivalent to the reference product when taken in the absence of food. A standardized
`feeding regimen affected the rate, but not extent of absorption of doxycycline from the
`pelletised formulation.
`
`KEY WORDS Doxycycline Absorption Bioavailability Food
`
`INTRODUCTION
`
`This study had three objectives. The first was to examine the bioavailability of
`two doxycycline products in fasted subjects and the second was to compare the
`effects of food on the bioavailability of one product. The third objective was to
`use the plasma doxycycline data to compare different pharmacokinetic methods
`for evaluation of bioavailability.
`
`*Reprint requests to: Dr Desmond B. Williams, F. H. Faulding and Co. Ltd, P.O. Box 746,
`Salisbury, South Australia, 5108, Australia.
`
`0142~2782/90/020093~13$06.50
`© 1990 by John Wiley & Sons, Ltd.
`
`Received2 October 1987
`Revised 7 April 1989
`
`Amneal 1068
`
`Amneal v. Supernus
`IPRZO13-00371
`
`
`
`L;._
`
`Doxycycline is a lipid soluble antibiotic and is generally considered to be well
`absorbed (greater
`than 90%).” When compared with other
`tetracycline
`derivatives, the absorption of doxycycline from the gastrointestinal tract appears
`to be less affected by food. Welling et al.3 found that ingestion with test meals
`reduced the absorption of tetracycline and doxycycline by about 50 and 20 per
`cent, respectively.
`Other comparative bioavailability studies have been performed with doxycy-
`cline.“ Insufficient time (up to 24 h) was allowed for collection of blood
`following administration of doxycycline for pharmacokinetic analysis and no
`attempt was made to fit a pharmacokinetic modgl to the data in these studies.
`When pharmacokinetic analysis has been performed, a one compartment open
`model has usually been assumed.3'7’8 However, a two compartment open model
`was assumed when doxycycline was administered by intravenous infusion.9
`Numerous methods have been proposed to compare the bioavailability of
`various formulations of the same drug.” The measurement of relative bioavail-
`ability is normally determined by two variables, the rate of absorption, and the
`extent of absorption of drug from the product.
`to the
`These are often determined by fitting a compartmental model
`concentration of drug in blood versus time data. Some of the model parameters
`(e. g. absorption rate constant) are used as the basis for comparing the rate aspect
`of drug absorption. The extent of absorption is measured by a comparison of
`areas enclosed by the concentration of drug in blood versus time profiles (AUC).
`More recently, techniques for evaluation of bioavailability known as either
`non—compartmental or model-independent (i.e. independent of compartment
`models) methods have been developed.” Differences in the extent of absorption
`are evaluated from area measurements and the rates of absorption are compared
`by determining the time of peak concentration of drug in the blood. Useful
`parameters, such as the mean residence time (MRT), have been developed which
`allow some combination of rate and extent of absorption into a combined
`concept.“
`In this paper, a number of these methods were applied to the extensive data
`sets generated in the study. It was hoped to show that a variety of methods, if
`properly applied, will provide useful and equivalent information about the
`relative bioavailability of different pharmaceutical products.
`
`EXPERIMENTAL
`
`Subject protocol
`
`the study was performed by Biodecision
`The experimental part of
`Laboratories, Pittsburg, Pennsylvania. Twenty-four healthy male subjects
`entered the study. Four subjects (subjects 1, 12, 18, and 22) withdrew for reasons
`
`
`
`DOXYCYCLINE ABSORPTION
`
`95
`
`unrelated to the study. The twenty subjects who completed the study were from
`19 to 35 years (255 i 5-1 years, mean :2 SD) and weighed between 609 and
`864 kg (72-5 i 8-0 kg). Each subject was given a physical examination, a medical
`history was taken, and informed consent was obtained. Before each dose, the
`subjects were restricted to at least a 10 h fast, followed by a high protein, low fat
`FDA recommended diet for 48 h after each dose. The subjects were not
`permitted to take either antibiotics for at least 15 days or other drug products for
`at least 7 days, prior to the start of the study.
`One hundred milligrams of doxycycline was administered as doxycycline
`hyclate to each subject during each of three phases. Water (180 ml) was taken
`with each dose. Doxycycline was administered as the pelletised product (Doryx
`Capsules®, doxycycline hyclate delayed-release capsules, USP, F.H. Faulding &
`Co. Ltd, Batch 3EDl57) with a standard breakfast for the first treatment.
`The pelletised product and the reference product (Vibramycin Capsules®,
`Pfizer Laboratories, Batch 23091) were administered‘ 2h before a standard
`breakfast for the second and third treatments, respectively. All three treatments
`were administered during each phase according to a Latin square design and
`each phase was separated by 1 week.
`Blood samples were collected (Vacutainer®, Beckton Dickinson and Co.) by
`venepuncture immediately prior to drug administration and at 0-5, 1, 1-5, 2, 3, 4,
`6, 8, 12, 16, 24, 36, and 48 h following drug administration. The volume of each
`collection was
`10 ml, except
`for
`the first collection, which was 20 ml.
`Immediately after blood collection, the samples were stored in an ice bath and
`protected from light. Following centrifugation, the plasma was collected and
`stored at —20° until analysis. Urine was collected at -1 to 0, 0 to 1, 1 to 2, 2 to 4, 4
`to 6, 6 to 8, 8 to 12, 12 to 24, and 24 to 48h following drug administration.
`During the collection period, all urine samples were stored at 4°. After
`measurement of the urine volume and pH, a 15 ml aliquot was taken from each
`sample and frozen until required for analysis.
`
`Determination of doxycycline in biologicalfluids
`
`Doxycycline concentrations in plasma and urine were measured by a micro-
`biological agar diffusion assay using Bacillus cereus ATCC 11778 as the test
`organism.”
`‘T
`Standards were prepared from USP doxycycline hyclate reference material
`using 0-1 M potassium phosphate buffer, pH 4-5, as diluent. Plasma standards
`were prepared by spiking buffer standards with either 1:5 or 1:10 plasma which
`was free from doxycycline. Low concentration and high concentration controls
`were prepared by spiking plasma and urine with doxycycline. The controls were
`divided into aliquots and stored with the subjects’ samples.
`A set of standards and controls were processed on each analysis day. The
`samples were diluted with 0-1 M potassium phosphate buffer, pH 4-5, to fall
`
`
`
`within the linear portion of the standard curve. The concentrations of
`doxycycline in plasma samples were determined by interpolation on curves
`obtained from standards containing the same proportion of plasma as the
`diluted sample. Urine samples were interpolated on a curve obtained from buffer
`standards. The measured concentrations were corrected for dilution.
`The lowest concentration of doxycycline which could be quantified with an
`acceptable degree of precision was 0-0125 mg 1-1. Any sample below this
`concentration was reported as zero. Day—to-day reproducibility of the assay was
`determined by analysing control samples. The coefficient of variation for the
`plasma assay was 5-5 per cent at a concentration of 0-60 mg 1”‘ and 9-9 per cent at
`a concentration of 1.19 mg 1-1. In urine, the coefffcient of variation was 11-3 per
`cent at a concentration of 0-51 mg 1'1 and 13-7 per cent at a concentration of
`1.13 mg 1”‘.
`Recovery was determined by comparing the measured concentration and the
`expected cdhcentration of each standard. Recovery on a typical analysis day in
`plasma ranged from 86 to 112 per cent and averaged 100 per cent. In urine,
`recovery ranged from 86 to 113 per cent and averaged 100 per cent.
`The assay was linear over the concentration range 000 to 0-70 mg l_1.
`
`‘Data analysis
`
`The data were subjected to curve peeling techniques by standard graphical
`methods and the use of ESTRIP.” Two and three term polyexponential
`equations were fitted to the data. The ESTRIP progam generated estimates of
`the coefficients and indices of the equations and also estimates of the lag time
`(nag) before absorption occurs.
`_
`Areas under the curves of concentration of drug in plasma versus time
`(AUC38) up to the last sampling time were estimated by the trapezoidal method.
`The extrapolated areas to infinity (AUC.3°) were calculated by dividing the last
`measured concentration by the terminal phase rate constant. The MRT was
`calculated from
`
`MRT = AUMC/AUC°6’
`
`(1)
`
`where AUMC is the area under the first moment curve.”’14 The AUMC was
`calculated by the method of Riegelman and Collier.”
`When the values of the coefficients, indices, and nag had been determined for
`each subject, the equation describing absorption and biexponential elimination
`was used to estimate the maximum concentration and the time at which the
`maximum concentration occurred. This was achieved by entering the equation
`into a computer program (available on request) which calculated the
`concentration at increasing increments of time near the peak.
`The renal clearance of doxycycline was calculated from
`
`
`
`DOXYCYCLINE ABSORPTION
`
`A48
`CIR = ——_
`
`AUC. BW
`
`,
`
`t
`
`97
`
`(2)
`
`where A43 is the amount of docycycline excreted in urine in 48 h and BW is the
`body weight.
`A number of compartmental models were fitted to the data as described
`below.
`
`One compartment model (model I). The simple one compartment model with
`first order input and output") and including a lag time for absorption was fitted
`to the data. In the fitting procedure, the absorption rate constant, kg, elimination
`rate constant, kc, nag, and FDo/ Vb were entered as parameters, where F is the
`fraction of dose absorbed, Do is the dose administered and Vb is the apparent
`Volume of distribution of the drug. Initial Values of these parameters were
`calculated using ESTRIP. Final Values of these parameters were obtained from
`the best fit
`to the data as described under ‘Model fitting and statistical
`procedures’.
`The maximum concentration and the time at which the maximum
`
`concentration occurred for each subject were calculated from standard
`equations. 10
`
`Two compartment open model with first order input (model 2). A standard
`triexponential equationlo with a tlag incorporated was fitted to the data. The
`Values obtained for slope constants, A1 and A2, were used as initial estimates of oz
`and B, the complex constants of the two compartment model. The smallest
`index Value obtained in the fitting procedure was assumed to be ,8. The Value of
`oz was equated with the index Value closest to the known values of oz for
`doxycycline given by the intravenous route (mean = O-80h"1, reference 8). The
`remaining index value (in most cases the largest)‘ was assigned as ka. In the
`with—food treatment,
`the or Value used as a starting estimate in the fitting
`procedure was an average of the at Values obtained with the same subject after
`fitting the other two treatments. In the fitting procedure, best fits were obtained
`for k.., oz, ,8, FD.,/ Va, /C21, and”1“iag, where k21 is a distribution rate constant and Vc
`is the Volume of the central compartment.
`The constants kl; and km were calculated from oz, B, and km using standard
`methods.”
`The Values obtained for the maximum concentration and the time at which
`the maximum concentration occurred were estimated as described for the non-
`
`compartmental parameters.
`
`Two compartment open model with zero-order input (model 3). This model
`assumes that absorption occurs at a constant rate and that the absorption
`
`
`
`process terminates at a definite point followed by an elimination phase“)
`
`(a) During the absorption phase equation 246 in Gibaldi and Perrier,” was used
`with tag incorporated. The zero-order absorption rate constant, kc, was
`calculated as the dose absorbed (Do) divided by the time of absorption:
`
`k =.
`0
`
`° _
`T‘ [lag
`
`I
`
`(3)
`
`i
`
`where Tis the time at which absorption stops. This is an approximate estimate
`of koz, as the absorption of doxycycline is generally slightly less than 100 per
`cent.
`
`(b) During the p‘ost—absorptive phase, equation 257 of Gibaldi and Perrier” was
`used. In fitting these equations to the data, 0:, kzi, Vc/ F and mg were entered into
`the equations as parameters which were varied until the best fit was obtained.
`The dose of drug administered (Do), and the time for the end of absorption (D
`were entered as constants. The values of the starting parameters were those
`obtained for each subject from the fitting of the two compartment equation
`assuming first order absorption. T was determined as the time of peak
`concentration of doxycycline in plasma since this was the point at which
`absorption ceased and a decline in the concentration of doxycycline in plasma
`began.
`
`Modelfitting and statistical procedures
`
`The equations described above were fitted to the unweighted data by a Monte
`Carlo parameter estimation technique known as REVOL.” The appropriate
`equations were written into the program and parameter estimates as described
`above were used to start the fitting procedure. In each case, the procedure was
`continued until a satisfactory fit was obtained. A number of criteria were used to
`determine the best fit. First, a fitted curve was required to be as free as possible
`from large regions of systematic error. Second, the smallest Values obtainable for
`the error sum (ES) as defined in by”
`
`1': 1 Mn
`
`(4)
`
`where n is the number of data points, yr is the ith observation and ym is the
`calculated value for that observation.
`
`The Akaike Information Criterion (AIC) is
`determining the goodness of fit:”
`
`the third guide useful for
`
`AIC=nlnSS+2p
`
`(5)
`
`
`
`
`
`DOXYCYCLINE ABSORPTION
`
`99
`
`where n is the number of data points, SS is the sum of the squared deviation
`between observed and calculated observations and p is the number of parameters
`fitted in the model. The AIC may also be used to compare the relative suitability
`of different models.” The lowest value of AIC generally indicates the model
`which best fits the data.
`
`Analysis of variance (ANOVA) used for complete crossover design was used
`on the various parameters from the individual treatments. Differences between
`individual treatments were determined using Tukey’s ‘honestly significant
`
`Table 1. Model independent parameters for doxycycline
`
`Parameters*
`
`Treatment means
`Pellet product
`Pellet
`product with without food
`food
`(A)
`
`(B)
`
`Reference
`product
`without food
`(C)
`
`ANOVAT
`
`27
`(4)
`
`31
`(6)
`
`39
`(10)
`
`0-33
`(0-08)
`
`0-55
`(0-24)
`
`1-8
`(0-5)
`
`1-2
`(0-5)
`
`1-7
`(0-4)
`
`22
`(4)
`
`27
`(4)
`
`31
`(5)
`
`38
`(7)
`
`0-33
`(0-09)
`
`0-27
`(0-17)
`
`-
`
`1-4
`(0-6)
`
`1-1
`(0-5)
`
`1-8
`(0-4)
`
`22
`(4)
`
`nsi
`
`ns
`
`ns
`
`ns
`
`p <0-001
`(A>B>C)"
`
`-- p <0-001
`(A>B=C)
`
`0-002 < p < 0-005
`A>B=C
`
`0-002 <p <0-005
`(A <C,A=B,B=C)
`
`ns
`
`AUC38 (mg h 1“)
`
`AUCB’ (mg h 1*)
`
`Au‘t,“(mg)
`
`cm (m1min"l kg”)
`
`rm 0“)
`
`tmax (h)
`
`[max — flag (h)
`
`0,, max (mg 1“)
`
`MRT (11)
`
`26
`(3)§
`
`29
`(4)
`
`38
`(8)
`
`0-34
`(0-08)
`
`1-1
`(0-37)
`
`3-1
`(0-9)
`
`1-9
`(0-9)
`
`1-4
`(0-6)
`
`22
`(3)
`
`*See text for definition of symbols.
`TAnalysis of variance.
`;tN0t significant.
`§Standard deviation.
`||Tukey’s test.
`
`
`
`
`
`
`
`PlasmaDoxycycline(mg/I)
`
`o
`
`10
`
`20
`
`30 S
`
`40
`
`so
`
`Time (h)
`
`Figure 1. Plasma concentrations of doxycycline following administration of the pellet product with a
`standard meal to 20 subjects. The data are shown as means and standard deviations and the solid line
`was generated by assuming a two compartment open model with first order absorption
`
`difference’ test.” All computing procedures were run on either a RSTS system or
`a Sirius 1 microcomputer.
`
`RESULTS AND DISCUSSION
`
`Visual examination of the plasma doxycycline concentration versus time plots
`(Figures 1, 2, and 3), area analysis of individual curves (AUC4(f and AUC, Table
`1) and the total drug recovered in urine (Alf, Table 1) all support the View that
`the extent of absorption of each product is equivalent. The MRT values (Table
`1) are also consistent with this result.
`As usual in the evaluation of bioavailability by area analysis, it is assumed that
`averaged clearance (and hence Volume of distribution) is unchanged between
`treatments.” The lack of variation in renal clearance between treatments (Table
`
`‘(J
`
`i
`
`is
`
`
`
`DOXYCYCLINE ABSORPTION
`
`101
`
`
`
`
`
`PlasmaDoxycycline(mg/I)
`
`Time (h)
`
`Figure 2. Plasma concentrations of doxycycline following administration of the pellet product 2h
`before a standard breakfast. The data and fitted curve were as described‘ in Figure l
`
`1) would support the assumption of constancy of total clearance in this study.
`Area analysis is related to the total amount of drug passing through the system
`but gives no information about the rate of absorption of drug from the gut. The
`rate of absorption is related to the time of peak drug concentration in blood
`(tmax) which is estimated from the data plots. The Values of tmax indicate that food
`significantly slows down themabsorption of the pellet product while there is no
`marked difference between the absorption rates of the products when fasting.
`The same conclusion is obtained when the tmax figures are corrected for an initial
`delay in drug release (tlag) calculated by the polynomial curve fitting program.
`The corrected Values (tmax — tlag, Tables 1 and 2) indicate the same differences in
`absorption rate as tmax. The tlag Values (Tables 1 and 2) indicate a delay in the
`initiation of absorption by food for the pellet product. In fasting subjects, the tiag
`for the pellet product was longer than the reference product. There was no
`correlation between either the urinary pH or urinary Volume and doxycycline
`elimination.
`
`
`
`Compartment models and bioavailability
`
`Three compartment models were fitted to the data. The one compartment
`model (model 1) and the two compartment model with zero order input (model
`3) showed a considerably poorer fit
`to the data compared with the two
`compartment model with first order input (model 2). Hence the latter is used in
`this discussion to evaluate bioavailability. Similar bioavailability determinations
`were made with the former models.
`
`Goodness of fit for the three models was determined by three methods. (a)
`Systematic error was least evident with model 2, particularly in the region of
`peak drug concentration and during the post peak distribution phase. (b) The
`error sum estimates (ES) were smallest for model 2. (c) The AIC Values were
`smallest for model 2.
`
`
`
`
`
`PlasmaDoxycycline(mg/I)
`
`0
`
`10
`
`20
`
`30
`
`40
`
`50
`
`Time (m
`
`Figure 3. Plasma concentrations of doxycycline following administration of the reference product 2 h
`before a standard breakfast. The data and fitted curve were as described in Figure 1
`
`
`
`DOXYCYCLINE ABSORPTION
`
`103
`
`Table 2. Parameters obtained assuming a two compartment open model with first order
`absorption for docycycline disposition
`
`Parameters*
`
`ka 01"‘)
`
`01 (11-1)
`
`13 (h“1)
`
`FD / V9 (mg 1"‘)
`
`[C10 01”‘)
`
`1621 (h")
`
`/C12 01”‘)
`
`tlag (11)
`
`tmax (h)
`
`tmax ~ tlag (h)
`
`Treatment means
`Pellet product
`Pellet
`product with without food
`food
`(A)
`
`(B)
`
`_
`Reference
`product
`without food
`(C)
`
`ANOVAT
`
`1-4
`(0-9)::
`
`0-60
`(0-20)
`
`0-046
`(0-007)
`
`2-2
`(0-5)
`
`0-075
`(0-021)
`
`0-37
`(0-12)
`
`0-22
`(0-13)
`
`1-13
`(0-35)
`
`3-1
`(0-8)
`
`2-00
`(0-75)
`
`2-6
`(1-3)
`
`0-79
`(0-32)
`
`0-047
`(O-008)
`
`2-6
`(1-0)
`
`0-084
`(0-037)
`
`0-44
`(0-11)
`
`0-31
`(0-26)
`
`0-53
`(0-23)
`
`1-6
`(0-4)
`
`-
`
`1-08
`(0-33)
`
`2-4
`(1-1)
`
`0-84
`(0-27)
`
`0-002 <p <0-005
`(A <B=C)§
`
`0-01 <p <0-05
`(A <B <c)
`
`0-046
`(0-006)
`
`ns 1'
`
`2-7
`(0-6)
`
`0-087
`(0-020)
`
`0-45
`(0-16)
`
`0-35
`(0-16)
`
`0-30
`(0-14)
`
`1-5
`(0-7)
`
`1-13
`(0-71)
`
`ns
`
`ns
`
`ns
`
`ns
`
`p < <0-001
`(A> 13> C)
`
`p < <0-001
`A > B=C
`
`p <.,_0-001
`A > B=C
`
`0-002 <p <0-005
`1-9
`1-8
`1-5
`cm (mg 1*‘)
`A <13=c turn
`(0-4)
`(0-4)
`(0-2)
`
`*See text for definition of symbols.
`‘(Analysis of variance.
`1Standard deviation.
`fiTukey’s test.
`Not significant.
`
`The parameters obtained using model 2 (Table 2) which were relevant to the
`absorption rate of docycycline (i.e., ka, tlag, tmax, tmax_t]ag, and Cm“) all showed the
`same relative differences between treatments as observed with the model
`
`
`
`independent analysis (Table 1). Those parameters concerned with distribution
`and elimination processes (i.e.,
`,8, km, km, and kzi) showed no significant
`difference between treatments. The only exception was on which was significantly
`different between treatments in the same order as ka varied. This discrepency
`probably reflects the difficulty of separating oz and kn, particularly where oz and
`ka approach each other in magnitude. The pharmacokinetic models used in this
`study assumed a continuous elimination phase, however, some subjects showed
`discontinuous jumps in concentration during the elimination phase. It has been
`reported that this observation is ‘consistent with a discontinuous absorption
`mechanism occurring,
`such as
`from enterohepatic cycling through the
`gall-bladde1‘.22
`t
`We
`
`CONCLUSIONS
`
`Each of the pharmacokinetic models tested suggests that doxycycline is
`absorbed at the same rate and extent from the pelletised product and the non-
`pelletised product when they are taken in the absence of food. Doxycycline is
`absorbed to an equivalent extent but at a slower rate when taken in the form of
`the pelletised product with food.
`'
`The results from this study suggest that different pharmacokinetic models may
`be used to draw similar conclusions about the bioavailability of doxycycline
`following oral administration from different dosage forms. Comparison of the
`systematic errors, error sum estimates, and AIC values suggested that the two
`compartment open model which included a lag time and a first order absorption
`process best fitted the plasma doxycycline data.
`
`ACKNOWLEDGEMENTS
`
`The authors are indebted to Dr R. Don Brown for the ESTRIP program, and to
`the staff at Biodecisions Laboratories, Pittsburg, Pennsylvania,
`for their
`cooperation in performing this study.
`
`REFERENCES
`
`l. M. Schach von Wittenau, Chemotherapy (supp1.), 13, 41 (1968).
`2. B. A. Cunha, C. M. Sibley and A. M. Ristuccia, Ther. Drug Manit., 4, 115 (1982).,
`3. P. G. Welling, P. A. Koch, C. C. Lau and W. A. Craig, Antimicrob. Agents Chemother., 11,462
`(1977).
`-
`4. J. D. Arcilla, J. L. Fiore, O. Resnick, J. W. Nadelmann, J. L. Huth and W. M. Traetel, Curr.
`Ther’. Res., 16, 1126 (1974).
`5. R. Kitzes Cohen, F. Grauer and E. Weisenberg, Israel J. Med. Sell, 16, 545 (1980).
`
`
`
`DOXYCYCLINE ABSORPTION
`
`I05
`
`»—-o\ooo\1_c~«
`
`._.._.
`
`A.—S. Malmborg, Chemotherapy, 30, 76 (1984).
`. M. Gibaldi, Chemotherapia, 12, 265 (1967).
`. G. Ceccarelli, R. Rossoni, F. Ronita and A. Naddeo, Chemotherapy, 16, 1 (1971).
`. T. C. Raghuram and K. Krishnaswamy, Br J Clin Pharmacol, 14, 785 (1982).
`. M. Gibaldi and D. Perrier, Pharmacokinetics, 2nd edn. Marcel Dekker, New York, 1982.
`. K. Yamaoka, T. Nakagawa and T. Uno, J. Pharmacokinet. Bz'opharm., 6, 547 (1978).
`Code of Federal Regulations, U.S. Government Printing Office, Washington, 1984, p. 244.
`R. D. Brown and J. E. Manno, J. Pharm. Scz'., 67, 1687 (1978).
`. S. Riegelman and P. Collier, J. Pharmacokinet. Biopharm, 8, 509 (1980).
`. M. Gibaldi, D. Perrier, Pharmacokinetics, lst edn. Marcel Dekker, New York, 1975.
`P. Koeppe and C. Hamann, Compuz. Prog. Biomed., 12, 121 (1980).
`K. Yamaoka, T. Nakagawa and T. Uno, J. Pharmacokinet. Bz'opharm., 6, 165 (1978).
`R. P. Sokal and F. J. Rohlf, Biometry, 2nd edn. W. H. Freeman, San Francisco, 1981.
`. J. M. Jaffe, J. L. Colaizzi, R. I. Poust and R. H. McDonald, J. Pharmacokinet. Bz'opharm., 1
`267 (1973).
`. J. M. Jaffe, R. I. Poust, S. L. Feld and J. L. Colaizzi, J. Pharm. Sci., 63, 1256 (1974).
`. J. G. Wagner and C. M. Metzler, J. Pharm. Sci., 56, 658 (1967).
`. P. Veng Pedersen and R. Miller, J. Pharm. Sci., 69, 204 (I980).
`
`‘€358BESEGEBE
`
`3
`
`haw