<<

Br J clin Pharmac 1994; 37: 545-551

Mixed linear and non-linear disposition of lazabemide, a reversible and selective inhibitor of B

T. W. GUENTERT', N. H. G. HOLFORD2, J. P. PFEFEN' & J. DINGEMANSE1 'Department of Clinical Pharmacology, F. Hoffmann-La Roche Ltd, CH-4002 Basel, Switzerland and 2Department of Pharmacology and Clinical Pharmacology, University of Auckland, New Zealand

1 Single oral doses (100-300 mg) and multiple oral doses (100-350 mg 12 hourly for 7 days) of lazabemide were administered to 35 young and 40 elderly healthy subjects. Plasma concentrations of unchanged drug were determined to study the dose-concentration relationship. 2 The elimination phase time course of lazabemide concentrations indicated concen- tration-dependent elimination after both single and multiple dosing. Nevertheless, maximum concentrations and areas under concentration-time curves increased almost proportionally with dose and accumulation after chronic dosing was less than a factor of 2; steady-state concentrations were achieved by the third day of dosing. The apparent half-life determining accumulation was approximately 8-9 h. 3 Drug absorption commenced rapidly after a dose; two components to the absorp- tion process were dectable in young subjects possibly due to simultaneous admin- istration of multiple tablets at the higher doses. 4 Observations after single and multiple dosing were described with a compartmen- tal model allowing for parallel saturable (population mean ± s.d.: maximum elimi- nation rate VmaxIF: 2.8 ± 1.4 mg h-'; concentration at half-maximum elimination K_: 36 + 19 pg 1-1) and first-order (CL/F 16 ± 3.8 1 h-') elimination pathways. No important difference between the young and the elderly subjects was noted in absorption or disposition parameters of lazabemide.

Keywords lazabemide Ro 19-6327 MAO-inhibition population pharmacokinetics modelling

Introduction Monoamine oxidases (MAO) are mitochondrial carboxamide; Ro 19-6327) is an investigational, selec- enzymes found in most mammalian tissues and are tive and reversible inhibitor of MAO-B [3]. During classified into type A and type B according to sub- single- and mulitple-dose tolerability studies in young strate specificity [1]. Owing to their involvement and elderly volunteers, plasma concentrations and in the degradation of neurotransmitter amines inhibition of platelet MAO-B activity were deter- (e.g. MAO-A: , , noradrenaline, mined to obtain initial information on the drug in ; MAO-B: dopamine) inhibitors of these humans and to enable studies on the relationship enzymes are clinically used or are being evaluated in between dose, drug concentrations and drug effects. the treatment of a number of diseases of the central The results from two different studies are reported. nervous system. While MAO-A inhibitors are pri- The first included healthy young subjects receiving marily used in the treatment of depression and certain single doses of 100-300 mg. In the second study phobic-anxiety states, MAO-B inhibitors are thought elderly subjects received active drug for 7 days at to be useful in Parkinson's disease as an adjunct to doses of 100 to 350 mg every 12 h. This report levodopa therapy or to slow down the progression of describes the search for a model to describe the dose- the disease [2]. concentration relationship and the population pharma- Lazabemide (N-(2-aminoethyl)-5-chloro-2-pyridine- cokinetic parameters of the most satisfactory model.

Correspondence: Professor Dr T. W. Guentert, Department of Clinical Pharmacology, F. Hoffmann-La Roche Ltd, CH-4002 Basel, Switzerland 545 546 T. W. Guentert et al.

Methods morning of day 1 and then twice a day (08.00 h, 20.00 h) for 7 days starting 48 h after the initial dose Clinical studies (days 3-9). The last dose was given in the morning of study day 10. Two different (same as single dose The clinical studies were carried out according to the study) tablet formulations (50, 100 mg) were used. principles laid down in the Declaration of Helsinki Serial blood samples of 10 ml were collected after following approval of the Institutional Ethics Com- the morning doses on days 1 and 10 (pre-dose; 0.5, 1, mittee and written informed consent from the sub- 1.5, 2, 4, 6, 8, 12, 24, 36, 48 h post-dose) and addi- jects. All volunteers understood that they were free to tional samples were obtained before the morning withdraw from the study at any time. doses on days 5, 7 and 9. Subjects were admitted dur- Both clinical studies reported were designed as ing the entire study period until 72 h after the last fixed sequence ascending dose, double-blind, parallel dose. group investigations with placebo control because pharmacological effects were also recorded. Prior to Blood collection procedure enrolment in the study, subjects underwent a physical examination and laboratory tests. Volunteers were Blood samples were collected via an indwelling excluded if these tests were abnormal, if they smoked venous catheter or by venipuncture into polypropy- more than 10 cigarettes/day, if they had a significant lene tubes containing EDTA as anticoagulant. food or drug allergy or if they had any history of drug Samples were centrifuged (3000 g for 10 min at 40 C) abuse or drug intolerance. within 30 min of collection and the plasma was Subjects were admitted to the clinical research cen- stored at -20° C until analyzed for unchanged drug. tre in the evening before the first drug administration. After an overnight fast they received drug or match- Chemical analyses ing placebo at 08.00 h together with 150-200 ml water. After the first (and, in the case of the multiple For determination of lazabemide 1 ml plasma was dose study, also the last) administration of the test combined in a glass tube with 50 pl of a solution con- the volunteers fasted for at least a further taining the internal standard (0.4 pg bromo-analogue 4 h. They did not lie down until 4 h after dosing of lazabemide ml-') in 0.4 N NaOH. After extraction except for the evaluation of vital signs and ECG (head over head 300 rev min-1; 15 min) with 10 ml t- recordings unless adverse events required a supine butyl-methyl ether/1-butanol (8:2, v:v), the organic position. phase (approximately 9 ml) was transferred into another glass tube containing 350 p1 of diluted Single-dose study phosphoric acid (1 ml o-phosphoric acid 85% w/v in 500 ml water). Back-extraction and centrifugation Thirty-six male volunteers between 19 and 36 years was then followed by a quantitative transfer of the old and within ± 15% of ideal body weight (mean 78 acidic phase into a new glass tube avoiding carry- kg, range 59-99 kg) were assigned to one of six over of the organic phase. After addition of 30 p1 of a active dose levels (100, 150, 200, 220, 250, 300 mg). sodium phosphate buffer (225 ml 1 N NaOH + 275 ml Three different tablet formulations (20, 50, 100 mg) of an aqueous solution containing per litre: 20 g were combined to achieve the designated active drug Na2HPO4 and 2.9 g NaH2PO4, pH 7.4), 200 p1 of dose. Blood (10 ml) was collected before and after fluorescamine solution (0.5 mg fluorescamine/10 ml drug administration (0.25, 0.5 1, 1.5, 2, 4, 6, 8, 10, acetonitrile) were added slowly under constant mix- 12, 16, 24, 32, 48 h). Subjects remained in the clini- ing. After evaporation of the acetonitrile in vacuo at cal research centre under medical observation until 400 C, 270 pl of the remaining solution was injected 32 h after drug administration and returned for post- into an h.p.l.c. system. study screening and the last blood sample at 48 h. Chromatographic separation of lazabemide from internal standard and interfering peaks was achieved Multiple-dose study with a Lichrocart Superspher 100 RP-18e cartridge (250 x 4 mm i.d., Merck) kept at a constant tempera- Twenty-one male and 19 female volunteers between ture of 600 C. The mobile phase (flow rate: 1.3 ml the ages of 60 and 78 years and within -15 to +25% min-') was a mixture of 744 g aqueous buffer solu- of ideal body weight (mean 73 kg, range 52-99 kg) tion (18.4 g NaH2PO4 1-l water adjusted to pH 5.45 participated in the study. Intake of medication com- with 1 N NaOH), 204 g acetonitrile, 38.4 g tetra- monly used by an elderly population was allowed hydrofurane. The column effluent was monitored by a during the study. Dose alterations were avoided Waters 470 fluorescence detector (2 excitation 370 during the study period as far as possible. However, nm, emission 485 nm). subjects requiring medication with a pronounced Along with every batch of study samples, calibra- effect on the central nervous system or with an tion samples (0.5-1000 pg 1-1') prepared daily by expected influence on the pharmacokinetics of other spiking drug-free plasma and aliquots of quality con- drugs (e.g. enzyme inhibitors/inducers, high doses of trol samples (1, 4, 40, 400, 1000 pg 1-1) prepared in antacids) were excluded from the study. Groups of sufficient volume before the entire analytical period eight subjects were administered one of five active were analyzed. drug doses (100, 150, 200, 250, 350 mg) on the Concentrations in quality control and study Pharmacokinetics of lazabemide 547 samples were calculated from the calibration curve Parameter estimation using peak height ratios of active drug/internal stan- dard. Concentrations in study samples were accepted Models were defined and parameters estimated using only if results from a minimum of 75% of quality MKMODEL [4]. Differential equations were solved control samples deviated less than 15% from their using a variable 4th/Sth order Runge-Kutta method known true value. Samples with drug concentrations [5] with a relative local error tolerance of 10-6 to above 1000 pg 1-F were diluted with drug-free plasma 10-12. The variance of the ith measured lazabemide and reanalyzed. concentration was predicted using a power variance The minimum quantifiable concentration of lazabe- model: mide was 1 pg 1-1 with an inter-day coefficient Vari = SD * (Yi + VO) of variation ranging from 4% (1000 pg I-') to 12% A (1 ,g 1-1). where SD is a scale factor for the error, Yi is the ith model prediction and VO the background variance. The power coefficient PWR was set to a constant Pharmacokinetic models with value 2. Absorption was modelled by either a single input or Model discrimination two sequential inputs with a lag-time prior to the start of each input. Each input process was tested as either Models of different structure and complexity were a first-order or zero-order process. Disposition was assessed by inspection of graphs of predictions and described by distribution into one or two compart- observations. With models that produced similar ments with elimination from the central compartment goodness of fit by visual criteria, the Schwarz crite- by a concentration-dependent (Michaelis-Menten, rion [6] was used to determine the most parsimonious mixed order) process either alone or in combination model. with a parallel first-order elimination pathway. The possible model parameters included a lag-time (tlag) Non-compartmental parameters and duration of zero-order input (TKO) or first-order absorption rate constant for each input; the fraction of Maximum concentration (Cmax) and the time of its the dose absorbed by the initial process (fl); inter- occurrence (tmax) were determined directly from the compartmental clearance (CLic) and volume terms observed time-concentration profile. Owing to the (V1, V2) for distribution; maximum elimination rate non-linear appearance of the concentration profile (Vmax), concentration at half-maximal elimination (Figures 2, 3) terminal half-lives were not calculated. (Km) and first-order clearance (CL) for elimination. The area under the time-concentration curve from 0 Vss was calculated from the sum of V1 and V2. The to 12 or 24 h (AUC(0,12), AUC(0,24)) was calcu- general model is depicted in Figure 1. lated from the observed concentrations by the linear Multiple dose accumulation was predicted from the trapezoidal rule. Concentrations could not be mea- results of the single dose study (Table 2). Concentra- sured after 24 h following 100 mg single doses; with tions after 15 doses given at 12 h intervals were higher doses 24 h concentrations were close to the obtained from a simulation using the mixed-order quantitation limit and samples collected 32 h after elimination model (single zero-order input into a drug administration were unmeasurable. Only after single compartment without lag-time). From these 300 mg single and 200 mg multiple doses could con- concentrations the AUC(0,12) was calculated and centrations be regularly measured up to 32-36 h. The compared with the AUC(0,12) obtained from accumulation behaviour of the drug after mul- observed concentrations after the first dose. tiple dosing was determined by comparing Cmax and AUC(0,12) after the first dose to the respective values after the last dose during the multiple dosing fli tlagl, study. Using the resulting accumulation factors (Rac) TKO1 or kai the drug's apparent half-life determining accumula- tion (tl, app) was calculated in every individual by rearranging the equation Rac = 1/(1 -e-(ln2/tAaPP)T), where X is the dosing interval (12 h).

Statistical procedures Figure 1 Pharmacokinetic model tested for description of the absorption and disposition of lazabemide after single Dose-normalized AUC within the single dose and and multiple dosing (tiag: lag-time; TKO: duration of zero- multiple dose studies and across studies at the same order input; ka first-order absorption rate constant; f1 = dose was compared using the Peritz F-test procedure fraction of the dose absorbed by process 1; V1/V2: volumes which includes adjustment for multiple comparisons of central/peripheral compartment; CLj,: inter- compartmental clearance; Vmax: maximum elimination rate; [7]. Kin: concentration at Vmax/2; CL: first-order elimination Using the expectation maximization (EM) algo- clearance). rithm proposed by Racine-Poon [8] the population 548 T. W. Guentert et al.

mean and standard deviation of the pharmacokinetic noted in any parameter after the first dose given to model parameter estimates of greatest interest young and elderly volunteers. (Vmax/F, Ki, CL/F) were obtained from the individ- Plasma drug concentrations in the single dose trial ual subject parameter estimates after generating the were described initially by a model allowing for covariance matrix with the full one- and two-com- only concentration-dependent elimination. Zero-order partment models. inputs were preferred over first-order processes to accommodate the rapid rise in concentrations after a dose. Irregularities around the peak could be mod- elled using two drug inputs, the second usually start- Results ing after a short lag time (Table 2). Although these models were well suited to describe drug concentra- The drug was rapidly absorbed with no detectable tions after a single dose, they predicted a much delay for the initial absorption process (tlag = 0). Two higher accumulation of drug than was actually seen components to the absorption process were detectable after multiple dosing. At a dosing level of 250 mg after the first dose in one-third of the younger sub- twice daily the AUC was expected to be four times jects in the single dose study (Figure 2) but in only higher after 7 days of dosing than after the first dose. one subject in the multiple dose study. They occurred The observed accumulation ratio was in most cases more often at high than at low doses. The model sug- less than 2-fold. This discrepancy between prediction gested that in the 12 subjects where double peaks and observation was clearly due to model were seen the first drug input started immediately misspecification. Therefore, for the description of after dosing and explained 78% (range 50-96%) of multiple dose data, the model was expanded to allow the total dose administered. Maximum concentrations for a second, concentration-independent elimination normalized to 100 mg ranged after single doses from clearance. Figures 3 and 4 show observed and pre- 365 to 544 jg 11- (Table 1) and after the multiple dicted plasma drug concentrations for two subjects. dosing regimen from 585 to 789 ig l-1. Model parameters Vmax/F, Ki,, CL/F and VSs/F (Table The terminal portions of the plasma concentration- 3) did not depend on the dose. The error model para- time profiles of the MAO-B inhibitor were character- meter VO was fixed at a negligibly small value (10-1o) istic of concentration-dependent elimination after and the PWR parameter at 2; the standard deviation single (Figure 2) and multiple dosing (Figure 3). In scale factor was determined as 20 ± 61 jig 1-1. A one- spite of this non-linearity in individual concentration- compartment model was sufficient to describe the sin- time profiles, Cmax and AUC increased almost propor- gle dose data in young volunteers but for many of the tionally with the dose (Table 1). There was no elderly subjects (28 out of 37) a two-compartment significant difference in dose normalized Cmax either in the single dose or in the multiple dose study (first and last doses). For the normalized AUC(0,12) and - 10 000 AUC(0,24) there was a slight trend to increase with -*a) 1000 increasing doses (Table 1). However, the only a) significant differences were detected in the multiple E 100 a) dose study (AUC(0,12) day 1: 100 mg vs 350 mg; .0 day 10: 100 mg vs 250 and 350 mg). Steady-state co 10 plasma drug concentrations were achieved already on C 1 the third E day of chronic dosing. The overall mean cu 0.1 accumulation ratio last/first dose for Cmax was 1.47 ± L 100 200 300 0.36 (range 0.86-2.4) and for AUC(0,12) it was 1.61 Time (h) ± 0.31 (range 0.97-2.2). No relevant difference was Figure 3 Observed (0) and predicted ( ; single input/dual elimination pathway model) plasma concentrations of lazabemide following multiple doses of 200 mg twice daily to an elderly subject.

m 1000 -1200 -3500- -n 100 C 900

a) 10 600 2500

N -o 300 E 1 'u-CD 1500 < 300 10 20 E , 30 N w w # X CD cn 0 1 L L- 8 TimefI..(h) - Cu 0 10 20 30 500 ~4= Time (h) E CD 0 L 0 1 0 20 30 40 Figure 2 Observed (0) and predicted ( ; dual input/single elimination pathway model) plasma Time (h) concentrations of lazabemide following a single dose of Figure 4 Comparison of plasma concentrations of 250 mg to a young subject. The figure with log-lin scale lazabemide in a healthy elderly volunteer after first (day 1; illustrates the concentration-dependent elimination; on the -0-) and last (day 10; -O-) dose during a multiple dosing lin-lin scale the dual input process is readily seen. regimen of 350 mg twice daily. Pharmacokinetics of lazabemide 549

00 - 0 r- t ON. Cl4 0I 000o ON 00 "t v:) 00 ON 00 o Cl Cl4 ON "z en \C _ tf r - C _)-n 00

o Cl 1- o o 00 Cli 0) In -.4 I- In r- r- -o _l _ 0t 111 CZ 1.) I-I

0 In ~00 - ', Cl4 Ce) 0 ~, --400 - I-- 0) C-0 Cl L. 0 -4 In r- c- rN t IC 00 0 ci In

e- ~ -. M "0 Ce) V) - Cef) 'Itt r2 ~ 00C CZ0) lO S on 0 E 00 In Cl .0C)co 0 o00 -o

.) >C Ien CO II4 N 0) XQI - _ Cs o O -O 00 ON Cl E o U:/C ":t

E 0 ci/ -0 00 - -o InUU: 00 -4 0Z 0Z -oo 0-_ -4 CO 0u

0)) CQ

"0 -0 o~~~~~~~~~~~~ I- _ I_0 O) 0 0) 0) C)00) V .: VD CZ CZ on r o _ 00 "0"0oA "0 ,~0"0°n Q C,) C) C) 000 :C-l CO C'-) a) C) t C W -t -0 '-0o ,:C O ( (0 1- _o m~~~ (0E .Wr_ 0~~ ~Z E m'-A E 2O 2T2) 2 550 T. W. Guentert et al. model was required to provide an adequate fit to the The failure of the single elimination pathway multiple dose data. model to predict the behaviour after multiple doses The combined mixed- and first-order elimination pointed to the existence of an additional elimination model predicted the observations after the first dose process. Concentration-dependent elimination is well and the concentration profile after the last dose known in the commonly acceptable range of equally well in 37 out of 40 multiple dose data sets. concentrations but the time course of elimination However, trough concentrations during the accumula- appears to be first-order at concentrations 10 times tion phase tended in many cases to be lower than pre- greater than that (reviewed by Holford [9]). The addi- dicted. In three subjects concentrations after the last tion of a first-order elimination path in parallel with dose were systematically lower than expected the concentration-dependent, mixed-order elimination (approximately 25%), based on the predictions after process provided a good description of the time the first dose. course of lazabemide not only after the first dose but also during accumulation to steady state and elimination of the last dose. First-order elimination Discussion dominates drug accumulation while the concentra- tion-dependent process determines the drug concen- A concentration-dependent elimination pathway for tration profile in the range of concentrations lazabemide was suggested by inspection of plasma associated with major changes in MAO-B inhibition drug concentration-time curves and supported by a [10]. single saturable elimination pathway model describ- Renal clearance accounts for approximately 60% ing data following a single dose. The maximum elim- of clearance of unchanged drug in the dose range of ination capacity estimated from the single pathway 50-150 mg (unpublished results); metabolic clearance model (365 mg day-l) predicted that oral doses leads mainly to a single acidic metabolite formed by greater than 360 mg day- would accumulate as long oxidative deamination of the side chain. The exact as drug was administered and did not reach steady mechanism for the saturable elimination of the drug state. However, plasma drug concentration-time is not confirmed at present. However, the observa- curves obtained during the multiple dose study tions are compatible with deamination of the mol- showed that accumulation was essentially complete ecule by the MAO-B enzyme which is itself inhibited and steady state was reached after less than 4 days by lazabemide. After sufficiently large doses lazabe- even with daily doses up to 700 mg day-' (Figure 3). mide's kinetics are principally determined by con-

Table 2 Single dose study: single elimination pathway parameter estimates I Mean ± s.d. Parameter Two stage EM Minimum Maximum VmaXIF (mg h-1) 15.2 ± 5.9 14.9 5.7 7.1 27.6 Km (gg 1-) 199 ± 135 187 124 10 493 Vss/F (1) 232 ± 39 155 332 TKO1 (h) 0.62 ± 0.42 0.21 2.2 f12 0.93 ± 0.13 - 0.5 1 tlag2 (h) 1.5 ± 1.5 0.25 5.7 TKo2(h) 1.6±1.1 0.3 3.6 'Estimates based on 35 data sets; 12 data sets required dual input model. 2fi in 12 data sets requiring dual input model: mean 78%, range 50-96%.

Table 3 Multiple dose study: dual elimination pathway parameter estimates1 Mean ± s.d. Parameter Two stage EM Minimum Maximum Vmax/F (mg h-') 3.72 ± 2.24 2.83 ± 1.35 0.75 9.16 Km (lg I-') 65 ± 58 36 ± 19 2.6 232 CL/F (1 h-1) 15.3 ± 4.13 16.1 ± 3.8 6.9 23.6 Vss/F (1) 229 ± 36 166 297 V2/F (1) 92.8 ± 42 38.0 203 CLic (lh-') 170 ± 172 TKO (h) 0.90 ± 0.41 - 0.22 1.9 'Estimates based on 37 data sets; 28 data sets required two- compartment disposition model. Pharmacokinetics of lazabemide 551 centration-independent disposition processes. The expected [11]. It is realized that several assumptions observation of overall dose linearity in the key are made in these calculations; based on present evi- kinetic parameters Cmax and AUC in the present dence these are not unreasonable. studies suggests that the non-linear elimination path- The rate of absorption of lazabemide was very way is only a small component of the overall disposi- rapid without any detectable delay. However, in some tion of the drug at doses of 200 to 400 mg day-'. The cases not all of the drug was released at the same ability to detect dose-dependent changes in dose time from the formulations. In 12 subjects in the normalized AUCs depends on the dose range studied. single dose trial double peaks were clearly seen. The Predicted departures from linearity at the doses fact that they were only seen when multiple tablets studied are of no clinical relevance. However, they were given suggests that this reflects separation of are expected to become more important if doses much tablets within the g.i. tract. The pharmacokinetic lower than 100 mg are studied. model developed needed to accommodate the irregu- Lazabemide was administered orally in both larities in absorption profiles in order to provide a studies. Therefore, several of the model parameters good description of observed concentrations at all determined are confounded by the extent of bio- times. This was achieved by allowing the dose to be availability. Although absolute bioavailability of the absorbed by two sequential processes. administered oral formulations is not known, almost Based on the observed accumulation ratio of Cmax complete extent of bioavailability of the drug is and AUC(0,12) the apparent half-life of lazabemide expected. After oral administration of radiolabelled determining accumulation is 7.5 ± 3.2 h and 8.8 ± 2.5 compound over 95% of the dose was found in the h, respectively. As expected this is somewhat shorter urine (unpublished results), confirming almost com- than the half-life (11 ± 3.9 h) predicted from the first- plete uptake of drug from the gastrointestinal tract order clearance and Vss- into the protal blood. Because 55% of an oral dose of The single and multiple dose studies reported were 100 mg is eliminated unchanged in the urine, renal performed in two different populations of young and clearance of unchanged drug has been estimated to be elderly healthy volunteers, both using very similar 14.4 1 h-1 (unpublished results in healthy Japanese study protocols with respect to single dosing. When subjects). The numerical similarity of the previously results across studies are compared the influence of found renal clearance and the first-order clearance age on the pharmacokinetics of lazabemide can be found in our study (16.1 1 h-1) may be fortuitous, but evaluated. From such a comparison it appears that no it might also be suggestive of a mechanistic explana- important changes occur in absorption or disposition tion. The saturable pathway would then be respon- parameters of lazabemide with progressing age in sible for any first-pass after oral dosing. healthy people. The possibility exists that in the With a blood/plasma partition ratio of 1.2 (unpub- elderly more of a distribution phase is seen in the lished results), the blood Km can be calculated to be concentration-time profiles compared with young around 43 ig 1-1. Assuming an hepatic blood flow of subjects. 90 1 h-', average portal venous blood concentrations during absorption of 100 mg over 1 h would be The SAS code for the expectation maximization algorithm around 1100 ,ug 1-1. The average intrinsic clearance was kindly supplied by Dr A. Racine-Poon, Ciba-Geigy, by the mixed-order pathway during absorption would Basel, Switzerland. Dr Ludger Banken, from F. Hoffmann- therefore be 2.5 1 h-'. From this value for total blood La Roche, gave valuable assistance in applying this code to CL/F an hepatic extraction ratio of less than 3% is our results.

References

1 McDaniel KD. Clinical pharmacology of monoamine 7 Marcus R, Peritz E, Gabriel KR. On closed testing pro- oxidase inhibitors. Clin Neuropharmacol 1986; 9: cedures with special reference to ordered analysis of 207-234. variance. Biometrika 1976; 63: 655-660. 2 Fowler CJ. Selective inhibitors of monoamine oxidase 8 Racine-Poon A. A Bayesian approach to nonlinear ran- types A and B and their clinical usefulness. Drugs of dom effects models. Biometrics 1985; 41: 1015-1023. the Future 1982; 7: 501-517. 9 Holford NHG. Clinical pharmacokinetics of ethanol. 3 Da Prada M, Kettler R Keller HH, Kyburz E, Burkard Clin Pharmacokin 1987; 13: 273-292. WP. Ro 19-6327: A novel highly selective and 10 Holford NHG, Guentert TW, Dingemanse J, Kettler R. reversible MAO-B inhibitor. Acta Pharmac Tox 1987; Pharmacodynamics of lazabemide, a reversible and 60 (Suppl): 10. selective inhibitor of . Br J clin 4 Holford, NHG. MKMODEL. Biosoft, Cambridge, Eng- Pharmac 1994; 37: 553-557. land, 1990. 11 Wilkinson GR, Shand DG. A physiological approach to 5 Fehlberg E. Low-order classical Runge-Kutta formulas hepatic drug clearance. Clin Pharmac Ther 1975; 18: including stepsize control and their applications to 377-390. some heat transfer problems. NASA Technical Report 1969; R-315. (Received 14 October 1993, 6 Schwarz G. Estimating the dimension of a model. Ann accepted 31 January 1994) Stat 1978; 6: 461-464.