Evaluation of theophylline pharmacokinetics in a pediatric population using mixed effects models

A Retrospective Evaluation of Allometry, Population Pharmacokinetics, and Physiologically-Based Pharmacokinetics for Pediatric Dosing Using Clearance as a Surrogate

Physiologically-based pharmacokinetic models are increasingly applied for pediatric dose selection along with traditional methods such as allometry and population pharmacokinetic models. We report a retrospective evaluation of the three methods. Pediatric population pharmacokinetic models sourced from literature for a subset of eight compounds were used to predict clearances for children < 2 years when they were within the modeled age range (interpolation, N = 11) or including those outside the modeled age range (interpolation and extrapolation, N = 18). Pediatric/adult clearance ratios were evaluated with a strict performance criterion of 0.8-1.25 and with twofold criteria. For children > 2 years, 58-75% of the clinical studies (N = 10) met the strict criteria, and > 80% of the clinical studies were predicted within twofold by all three methods. For children < 2 years, physiologically-based pharmacokinetic, allometry with age-dependent exponents, and pediatric population pharmacokinetic models predict 54%, 82%, and 64% within twofold of the observed, respectively.

Population pharmacokinetics of theophylline in adult Chinese patients with asthma and chronic obstructive pulmonary disease

Background Theophylline has a narrow therapeutic range and large interindividual variability in blood levels, so a thorough understanding of its pharmacokinetic characteristics is essential. Population pharmacokinetic (PPK) approaches could achieve it and many PPK studies of theophylline have been reported in infants. However, none was conducted in Chinese adults and none has explored the effect of CYP1A2 genotypes on the PPK characteristics of theophylline in adults. Objective To evaluate the PPK characteristics of theophylline and to assess the possible influence of covariates, including CYP1A2 genotypes, on theophylline clearance in Chinese adult patients. Setting The study is conducted at the department of respiration in Zhujiang Hospital, Guangzhou, China. Methods Theophylline concentrations were obtained from eligible patients and were measured by high performance liquid chromatography. The polymorphisms of - 3860G > A, - 163C > A, C5347T (CYP1A2*1B) and G-3113A were genotyped using a direct sequencing method. Then, CYP1A2 genotypes, age, fat-free mass (FFM) and other covariates were used to develop a PPK model by NONMEM software. Bootstrap analysis was used to asses the accuracy and prediction of the PPK model. Main outcome measure The concentration and clearance of theophylline. Results A total of 134 theophylline concentrations from 95 patients were obtained. The final model was as follows: CL/F(L/h) = 4.530 × (FFM/56.1)^0.75 × 0.713^CYP1A2*1B, the inter-individual variability in clearance/bioavailability (CL/F) was 44.0%, and the residual variability was 9.8%. The final model was proved to be reliable by bootstrap analysis. Conclusion Theophylline clearance was significantly associated with FFM and CYP1A2*1B genotypes in Chinese adult patients.

Theophylline Population Pharmacokinetics and Dosing in Children Following Congenital Heart Surgery With Cardiopulmonary Bypass

Children undergoing cardiac surgery requiring cardiopulmonary bypass (CPB) frequently develop acute kidney injury due to renal ischemia. Theophylline, which improves renal perfusion via adenosine receptor inhibition, is a potential targeted therapy. However, children undergoing cardiac surgery and CPB commonly have alterations in drug pharmacokinetics. To help understand optimal aminophylline (salt formulation of theophylline) dosing strategies in this population, a population-based pharmacokinetic model was developed using nonlinear mixed-effects modeling (NONMEM) from 71 children (median age 5 months; 90% range 1 week to 10 years) who underwent cardiac surgery requiring CPB and received aminophylline as part of a previous randomized controlled trial. A 1-compartment model with linear elimination adequately described the pharmacokinetics of theophylline. Weight scaled via allometry was a significant predictor of clearance and volume. In addition, allometric scaled clearance increased with age implemented as a power maturation function. Compared to prior reports in noncardiac children, theophylline clearance was markedly reduced across age. In the final population pharmacokinetic model, optimized empiric dosing regimens were developed via Monte Carlo simulations. Doses 50% to 75% lower than those recommended in noncardiac children were needed to achieve target serum concentrations of 5 to 10 mg/L.

Pharmacokinetics and metabolism of natural methylxanthines in animal and man

Caffeine, theophylline, theobromine, and paraxanthine administered to animals and humans distribute in all body fluids and cross all biological membranes. They do not accumulate in organs or tissues and are extensively metabolized by the liver, with less than 2% of caffeine administered excreted unchanged in human urine. Dose-independent and dose-dependent pharmacokinetics of caffeine and other dimethylxanthines may be observed and explained by saturation of metabolic pathways and impaired elimination due to the immaturity of hepatic enzyme and liver diseases. While gender and menstrual cycle have little effect on their elimination, decreased clearance is seen in women using oral contraceptives and during pregnancy. Obesity, physical exercise, diseases, and particularly smoking and the interactions of drugs affect their elimination owing to either stimulation or inhibition of CYP1A2. Their metabolic pathways exhibit important quantitative and qualitative differences in animal species and man. Chronic ingestion or restriction of caffeine intake in man has a small effect on their disposition, but dietary constituents, including broccoli and herbal tea, as well as alcohol were shown to modify their plasma pharmacokinetics. Using molar ratios of metabolites in plasma and/or urine, phenotyping of various enzyme activities, such as cytochrome monooxygenases, N-acetylation, 8-hydroxylation, and xanthine oxidase, has become a valuable tool to identify polymorphisms and to understand individual variations and potential associations with health risks in epidemiological surveys.

Effect of gender on theophylline clearance in the asthmatic acute phase in Japanese pediatric patients

The effect of gender on theophylline clearance was investigated retrospectively in 96 Japanese pediatric patients (63 males and 33 females) ranging in age from 0.5 to 8 years and in weight from 6.3 to 36.8 kg. All patients received intravenous constant-rate infusion of aminophylline in the asthmatic acute phase. The theophylline clearances in males and females were 56.2+/-15.4 and 50.1+/-14.2 ml/h/kg for ages 0.5-<2 years, 58.7+/-18.8 and 48.3+/-6.5 ml/h/kg for ages 2-<4 years, and 65.7+/-12.0 and 52.1+/-16.8 ml/h/kg for ages 4-<9 years, respectively. At ages from 2 to 8 years, the theophylline clearance was 20% higher in males than in females (p<0.05). Our findings suggested that the initial dosage of theophylline should be adjusted according to the gender of pediatric patients and particularly in the case of infants.

Comparison of CYP1A2 and NAT2 phenotypes between black and white smokers

The lower incidence rate of transitional cell carcinoma of the urinary bladder in blacks than in whites may be due to racial differences in the catalytic activity of enzymes that metabolize carcinogenic arylamines in tobacco smoke. To examine this, we compared cytochrome P4501A2 (CYP1A2) and N-acetyltransferase-2 activities (NAT2) in black and white smokers using urinary caffeine metabolites as a probe for enzyme activity in a community-based study of 165 black and 183 white cigarette smokers. The paraxanthine (1,7-dimethylxanthine, 17X)/caffeine (trimethylxanthine, 137X) ratio or [17X+1,7-dimethyluric acid (17U)]/137X ratio was used as an indicator of CYP1A2 activity. The 5-acetyl-amino-6-formylamino-3-methyluracil (AFMU)/1-methylxanthine (1X) ratio indicated NAT2 activity. The odds ratio for the slow NAT2 phenotype associated with black race was 0.4; 95% confidence intervals 0.2-0.7. The putative combined low risk phenotype (slow CYP1A2/rapid NAT2) was more common in blacks than in whites (25% vs. 15%, P<0.02). There were no significant racial differences in slow and rapid CYP1A2 phenotypes, and in the combined slow NAT2/rapid CYP1A2 phenotype. Age, education, cigarette smoking amount, body mass index, GSTM1 and GSTM3 genotypes were unrelated to CYP1A2 and NAT2 activity. Intake of cruciferous vegetables (primarily broccoli), red meat, carrots, grapefruit and onions predicted CYP1A2 activity either for all subjects or in race-specific analyses. Carrot and grapefruit consumption was related to NAT2 activity. Collectively, these results indicated that phenotypic differences in NAT2 alone or in combination with CYP1A2 might help explain the higher incidence rates of transitional cell bladder cancer in whites.

Population pharmacokinetics of theophylline during paediatric extracorporeal membrane oxygenation

AIMS

To determine the population pharmacokinetics of theophylline during extracorporeal membrane oxygenation (ECMO) from routine monitoring data.

METHODS

Retrospective data were collected from 75 term neonates and children (age range 2 days to 17 years) receiving continuous infusions of aminophylline (mean rate 9.2 +/- 2.6 micro g kg-1 min-1) during ECMO. A total of 160 plasma concentrations (range 1-8 per patient), sampled at time intervals ranging from 10 h to 432 h, were included. Population PK analysis and model building were carried out using WinNonMix Professional (Version 2.0.1). Cross-validation was used to evaluate the validity and predictive accuracy of the model.

RESULTS

A one-compartment model with first order elimination combined with an additive error model was found to best describe the data. Of the covariables tested, bodyweight significantly influenced clearance and volume of distribution, whereas age was an important determinant of clearance, as adjudged by the differences in the -2 x log likelihood (P < 0.005) and the residual error value. The final model parameters were estimated as: clearance (l h-1) = 0.023 x bodyweight (kg) + 0.000057 x age (days) and volume of distribution (l) = 0.57 x bodyweight (kg). The interindividual variability in clearance and volume of distribution was 38% and 40%, respectively. The residual error corresponded to a standard deviation of 3.6 mg l-1. Cross-validation revealed a median (95% confidence interval) model bias of 9.4% (2.9, 16.5%) and precision of 29.5% (24.8, 36.0%).

CONCLUSIONS

The estimated clearance is significantly lower, and volume of distribution higher, than previously reported in non-ECMO patients of similar age. These differences are probably a result of the expanded circulating volume during ECMO and altered renal and hepatic physiology in this critically ill group. Large interindividual variability reflects the heterogeneous nature of patients treated on ECMO.

Drug-drug pharmacodynamic interaction detection by a nonparametric population approach. Influence of design and of interindividual variability

Population approaches are appealing methods for detecting then assessing drug-drug interactions mainly because they can cope with sparse data and quantify the interindividual pharmacokinetic (PK) and pharmacodynamic (PD) variability. Unfortunately these methods sometime fail to detect interactions expected on biochemical and/or pharmacological basis and the reasons of these false negatives are somewhat unclear. The aim of this paper is firstly to propose a strategy to detect and assess PD drug-drug interactions when performing the analysis with a nonparametric population approach, then to evaluate the influence of some design variates (i.e., number of subjects, individual measurements) and of the PD interindividual variability level on the performances of the suggested strategy. Two interacting drugs A and B are considered, the drug B being supposed to exhibit by itself a pharmacological action of no interest in this work but increasing the A effect. Concentrations of A and B after concomitant administration are simulated as well as the effect under various combinations of design variates and PD variability levels in the context of a controlled trial. Replications of simulated data are then analyzed by the NPML method, the concentration of the drug B being included as a covariate. In a first step, no model relating the latter to each PD parameter is specified and the NPML results are then proceeded graphically, and also by examining the expected reductions of variance and entropy of the estimated PD parameter distribution provided by the covariate. In a further step, a simple second stage model suggested by the graphic approach is introduced, the fixed effect and its associated variance are estimated and a statistical test is then performed to compare this fixed effect to a given value. The performances of our strategy are also compared to those of a non-population-based approach method commonly used for detecting interactions. Our results illustrate the relevance of our strategy in a case where the concentration of one of the two drugs can be included as a covariate and show that an existing interaction can be detected more often than with a usual approach. The prominent role of the interindividual PD variability level and of the two controlled factors is also shown.

The population pharmacokinetics of long-term methotrexate in rheumatoid arthritis

AIMS

Methotrexate is considered by many practitioners to be the agent of choice in the treatment of rheumatoid arthritis. The pharmacokinetics of methotrexate have been reported to exhibit significant intersubject variability. Therefore, this study was undertaken to evaluate the population pharmacokinetics of methotrexate during long-term administration in adults with rheumatoid arthritis.

METHODS

Methotrexate pharmacokinetics were evaluated in a 36 month study of 62 adults with rheumatoid arthritis. Patients received oral or intramuscular doses of methotrexate weekly with pharmacokinetic studies performed every 6 months. Data were analyzed with nonlinear mixed effects modeling.

RESULTS

Three thousand two hundred and sixty post oral or intramuscular dose serum methotrexate concentrations comprising 425 individual concentration vs time profiles were modeled using NONMEM. Covariates that significantly (P < 0.005) influenced the disposition of methotrexate were age (AGE, years), body weight (BW, kg), creatinine clearance (CL(CR), 1 h(-1)), gender (GEN; 0 = male, 1 = female), dose (DOSE, micromol), and fed vs fasted state (FED; 0 = fasted, 1 = fed). The final model describing the biexponential disposition of methotrexate was clearance(CL, 1 h(-1)) = (0.0810*BW + 0.257*CL(CR))*(1-(0.167*GEN); central volume (Vc, 1) = 0.311*BW; peripheral volume (Vp, 1) = 0.469*BW-0.169*AGE; intercompartmental clearance (Q, 1 h(-1)) = 4.27*(1-0.355*GEN); oral absorption rate constant (ka(p.o.), h(-1)) = 4.70-0.0439*DOSE*(1-0.507*FED); intramuscular absorption rate constant (ka(i.m.), h(-1)) = 0.122*DOSE; relative bioavailability (F) = 93.4%; and oral absorption lag time (LAG(p.o.), min) = 13.5. Pharmacokinetic parameters (%CV) for a typical fasted male subject in this study were CL, 7.341 h(-1) (27%); Vc, 23.51 (28%); Vp, 25.31 (31%); Q, 4.25 1 h(-1) (41%); ka(p.o.), 3.67 h(-1) (77%); and ka(i.m.), 3.09 h(-1) (44%).

CONCLUSIONS

The population pharmacokinetics of methotrexate in adults with rheumatoid arthritis were well described by this investigation. Substantial interpatient variability was explained by incorporating patient specific data into regression equations predicting pharmacokinetic parameters.

Clinical pharmacokinetics of vasodilators. Part II

Stimulating cardiac beta 1-adrenoceptors with oxyfedrine causes dilatation of coronary vessels and positive inotropic effects on the myocardium. beta 1-adrenergic agonists increase coronary blood flow in nonstenotic and stenotic vessels. The main indication for the use of the phosphodiesterase inhibitors pamrinone, mirinone, enoximone and piroximone is acute treatment of severe congestive heart failure. Theophylline is indicated for the treatment of asthma, chronic obstructive pulmonary disease, apnea in preterm infants ans sleep apnea syndrome. Severe arterial occlusive disease associated with atherosclerosis can be beneficially affected by elcosanoids. These drugs must be administered parenterally and have a half-life of only a few minutes. Sublingual or buccal preparations of nitrates are the only prompt method (within 1 or 2 min) of terminating anginal pain, except for biting nifedipine capsules. The short half-life (about 2.5 min) of nitroglycerin (glyceryl trinitrate) makes long term therapy impossible. Tolerance is a problem encountered with longer-acting nitric oxide donors. Knowledge of the pharmacokinetic properties of vasodilating drugs can prevent a too sudden and severe blood pressure decrease in patients with chronic hypertension. In considering the administration of a second dose, or another drug, the time necessary for the initially administered drug to reach maximal efficacy should be taken into account. In hypertensive emergencies urapidil, sodium nitroprusside, nitroglycerin, hydralazine and phentolamine are the drugs of choice, with the addition of beta-blockers during catecholamine crisis or dissecting aortic aneurysm. Childhood hypertension is most often treated with angiotensin-converting enzyme (ACE) inhibitors or calcium antagonists, primarily nifedipine. Because of the teratogenic risk involved with ACE inhibitors, extreme caution must be exercised when prescribing for adolescent females. The propagation of health benefits to breast-fed infants, combined with more women delaying pregnancy until their fourth decade, has entailed an increase in the need for hypertension management during lactation. Low dose hydrochlorothiazide, propranolol, nifedipine and enalapril or captopril do not pose enough of a risk of preclude breastfeeding in this group. The most frequently used antihypertensive agents during pregnancy are methyldopa, labetalol and calcium channel antagonists. Methyldopa and beta-blockers are the drugs of choice for treating mild to moderate hypertension. Prazosin and hydralazine are used to treat moderate to severe hypertension and hydralazine, urapidil or labetalol are used to treat hypertensive emergencies. The use of overly aggressive antihypertensive therapy during pregnancy should be avoided so that adequate uteroplacental blood flow is maintained. Methyldopa is the only drug accepted for use during the first trimester of pregnancy.

Aspects of theophylline clearance in children

Michaelis-Menten pharmacokinetic parameters for theophylline were estimated in a three-month infant following an accidental overdose of intravenous aminophylline. Fitting of time-concentration data was performed using nonlinear regression with MKMODEL. A mixed order elimination model was superior to a first order model. Parameter estimates were standardized to a 70 kg human using an allometric power model. Parameter estimates (SE) were: maximum rate of metabolism (Vmax) 71 (42) mg.h-1, Michaelis-Menten constant (Km) 32.3 (33.5) mg.l-1, volume of distribution (Vd) 46.9 (2.6)l. This Michaelis-Menten constant is lower than that reported for adults and consequently non-linear elimination will occur at lower plasma concentrations in infants than in adults. Theophylline clearance has traditionally been reported as directly proportional to body weight. This per kilogram model gives an erroneous impression that clearance is greatest in early childhood and then decreases with age until adult rates are reached in late adolescence. Age-related clearance values reported in the literature were reviewed using an allometric 3/4 power model. This size model demonstrates that clearance increases in infancy and reaches adult rates in the first one to two years of life.

Theophylline population pharmacokinetics from routine monitoring data in very premature infants with apnoea

1. Theophylline is commonly used in neonatology for the treatment and prophylaxis of apnoea of prematurity, and during ventilator weaning. 2. NONMEM was used to study the population pharmacokinetics of intravenous and oral theophylline from retrospective drug monitoring data in 82 premature neonates, weighing < 1500 g at birth, and < or = 32 weeks gestational age. 3. Clearance (CL), volume of distribution (V), and oral bioavailability (F1) from liquid preparations were modelled alone, and under the influence of demographic and clinical covariates, assuming a 1-compartment model with first-order elimination. 4. The final population models with influential co-variates were as follows: CL (1h-1) = 0.0000123 *body weight (g) + 0.000377 *postnatal age (days); V (1) = 0.000937 *body weight (g); F = 0.918. 5. The CL was lower and V was higher than previously reported for less premature neonates, term babies, and older children. 6. Predictive performance of the population models was evaluated by Bayesian forecasting in a similar, but independent cohort of 30 infants. There was statistically insignificant bias and imprecision between measured and predicted serum theophylline concentrations. 7. Based on the validated population models, recommended maintenance theophylline dosages are provided for infants aged between 2 and 50 days, and weighing 700 to 2000 g.

Gender differences in human pharmacokinetics and pharmacodynamics

Gender differences in pharmacokinetics and pharmacodynamics have long been recognized in animals. In humans, however, little attention has been paid to this field despite at least theoretical reasons to believe that gender may be an important variable in the processes of absorption, distribution, metabolism, and excretion. Gastric acid secretion, gastrointestinal blood flow, proportions of muscular and adipose tissue, amount of drug binding proteins, gender-specific cytochrome P450 isozymes, physiologic and hormonal changes during the menstrual cycle, and renal blood flow are several factors that may contribute to sex-related differences in pharmacokinetics. Clinical investigations have documented greater absorption and subsequent incorporation of iron into erythrocytes, and higher bioavailability of ethanol in females. Women have been shown to have a slower metabolism of mephobarbital and propranolol but an increased biotransformation of methylprednisolone, all three of which are metabolized by enzymes of the cytochrome P450 system. Lastly, the renal excretion of amantadine was inhibited significantly by quinidine and quinine in men but not in women. While gender-specific pharmacodynamic data are meager, evidence also supports the existence of sex-related differences. Women appear to be more prone to develop torsades de points from drugs such as quinidine and procainamide than men. A dimorphism in insulin sensitivity has been demonstrated with males having an enhanced response compared to females. Pharmacokinetic and pharmacodynamic sex-related differences exist and are complex. Future research efforts should be designed to provide more gender-specific information on drug disposition and clinical effect.

Pharmacokinetic variability in the adolescent: implications of body size and organ function for dosage regimen design

Although there are convincing clinical studies demonstrating important pharmacokinetic differences between the adult and pediatric patient, guidelines for adjusting the dosages of specific drugs are often based on empirical and limited data. Adult doses often provide the reference point for therapy in children with adjustment for body size. Both body weight and body surface area (BSA) are commonly used to adjust adult doses for pediatric patients but yield substantially different absolute doses. For the child age five years, the BSA-based dose will be more than 50 percent greater than the dose adjusted for body weight. The choice of weight or BSA is often arbitrary and may be an important confounding variable when evaluating drugs over wide ranges of age in pediatric patients or comparing two drugs for which doses are not adjusted in a similar manner. For example, comparative trials for HIV-infected pediatric patients are evaluating zidovudine dosed by BSA with zalcitabine dosed by body weight. Organ function changes with age in pediatric patients yet little information is available on the effects of maturation on hepatic or renal function and the consequences for drug disposition. The use of model substrates for organ function offers potential for elucidating organ function relative to maturation and, in particular, to those functional capacities associated with the onset of puberty and gender differentiation.

Some general estimation methods for nonlinear mixed-effects models

A nonlinear mixed-effects model suitable for characterizing repeated measurement data is described. The model allows dependence of random coefficients on covariate information and accommodates general specifications of a common intraindividual covariance structure, such as models for variance within individuals that depend on individual mean response and autocorrelation. Two classes of procedures for estimation in this model are described, which incorporate estimation of unknown parameters in the assumed intraindividual covariance structure. The procedures are straightforward to implement using standard statistical software. The techniques are illustrated by examples in growth analysis and assay development.

Digoxin pharmacokinetics in patients with high serum digoxin concentrations

Digitalis intoxication is a frequent iatrogenic effect in patients on treatment with digoxin. In the present study we evaluated the pharmacokinetic behaviour of digoxin and the factors responsible for intoxication by this drug in monitored patients exhibiting clinical signs of overdosing with serum levels > 2 ng/ml. A control group of patients was used as a reference whose population pharmacokinetic parameters obtained by a maximum likelihood method were: Vd = 542.92 +/- 274.53 (litre); Cl = 8.73 +/- 1.55 (litre/h) (mean +/- SD). Statistically significant differences (P < 0.001) were found between the mean Cl values in both groups of patients. The difference between the dose-level ratios established in both populations studied also proved to be significant (P < 0.001). Calculation of the optimum dose for each patient showed that the doses recommended in intoxicated patients should be three times lower than those used in the control population. A good correlation was found between the concentrations observed 24 h after administration and the mean concentrations observed at steady state predicted for both population groups. Multiple regression analysis showed that the variables with the greatest predictive value for clearance in intoxicated patients were age and renal function. The modifications observed in the pharmacokinetic behaviour and in the response to digoxin in this type of patient suggest systematic monitoring using pharmacokinetic and clinical criteria jointly.

Sparse data analysis

In recent years there has been a growing interest in techniques capable of analyzing sparse data, particularly gathered during Phase III clinical trials, and there is now pressure on manufacturers to obtain more kinetic and dynamic information from Phase III studies. Techniques for the analysis of sparse data are reviewed drawing on a number of examples taken from pharmacokinetic and pharmacodynamic experiments.

Smooth nonparametric maximum likelihood estimation for population pharmacokinetics, with application to quinidine

The seminonparametric (SNP) method, popular in the econometrics literature, is proposed for use in population pharmacokinetic analysis. For data that can be described by the nonlinear mixed effects model, the method produces smooth nonparametric estimates of the entire random effects density and simultaneous estimates of fixed effects by maximum likelihood. A graphical model-building strategy based on the SNP method is described. The methods are illustrated by a population analysis of plasma levels in 136 patients undergoing oral quinidine therapy.

Evaluation of hypothesis testing for comparing two populations using NONMEM analysis

In a simulation study of inference on population pharmacokinetic parameters, two methods of performing tests of hypotheses comparing two populations using NONMEM were evaluated. These two methods are the test based upon 95% confidence intervals and the likelihood ratio test. Data were simulated according to a monoexponential model and, in that context, power curves for each test were generated for (i) the ratio of mean clearance and (ii) the ratio of the population standard deviations of clearance. To generate the power curves, a range of these parameters was employed; other pharmacokinetic parameters were selected to reflect the variability typically present in a Phase II clinical trial. For tests comparing the means, the confidence interval tests had approximately the same power as the likelihood ratio tests and were consistently more faithful to the nominal level of significance. For comparison of the standard deviations, and when the volume of information available was relatively small, however, the likelihood ratio test was more able to detect differences between the two groups. These results were then compared to results on parameter estimation in order to gain insight into the question of power. As an example, the nonnormality of estimates of the ratio of standard deviations plays an important role in explaining the low power for the confidence interval tests. We conclude that, except for the situation of modeling standard deviations with only sparse information, NONMEM produces tests of significance that are effective at detecting clinically significant differences between two populations.

Bayesian parameter estimation and population pharmacokinetics

The widespread application of Bayesian parameter estimation in the area of therapeutic drug monitoring (TDM) has prompted the need for well conducted population studies to obtain relevant prior pharmacokinetic parameter estimates. In many cases the population has consisted of a relatively small number of subjects. This may be unavoidable for drugs used in cancer chemotherapy or in small, specific populations of patients. In contrast, information about drugs which are used extensively, such as the aminoglycosides, can be obtained by population studies which involve a large number of individuals. Indeed, this technique has proved particularly useful for determining parameter estimates which can be employed in neonatal TDM. Bayesian parameter estimation has been most frequently used for drugs with narrow therapeutic ranges such as the aminoglycosides, cyclosporin, digoxin, anticonvulsants (especially phenytoin), lithium and theophylline. However, the technique has now been extended to cytotoxic drugs, Factor VIII and warfarin. Bayesian methods have also been used to limit the number of samples required in more conventional pharmacokinetic studies with new drugs. Further advances in the use of these methods are likely to include measures of drug response and toxicity requiring population studies which also include relevant pharmacodynamic information.

Optimal sampling theory and population modelling: application to determination of the influence of the microgravity environment on drug distribution and elimination

Newer mathematical techniques such as optimal sampling theory and NON-linear Mixed Effects Modelling (NONMEM) allow the determination of pharmacokinetics and pharmacodynamics in populations of individuals previously believed to be "too ill" or "too difficult to study." Optimal sampling determines the most information-rich times to sample the system, allowing robust parameter estimates to be determined from the minimal number of samples. NONMEM, by taking the population as the unit of analysis, allows even fragmentary patient data sets to contribute to population parameter estimates. Obviously, the microgravity environment presents extreme logistical difficulties to the performance of traditional pharmacokinetic and pharmacodynamic studies. Examples of the validation of these techniques are presented, which indicates their likely utility in the important task of determining the influence of the microgravity environment on drug distribution and elimination, even accounting for the limitations of support which will be faced in this circumstance.

Steady-state population pharmacokinetics of sustained release theophylline in adult asthmatic patients

The steady-state population pharmacokinetics of theophylline were studied in 52 asthmatic adult patients who received sustained-release theophylline as armophylline or euphylline. A total of 92 steady-state plasma theophylline concentration-dosage pairs were analyzed using a nonlinear mixed effects model. The pharmacokinetic model used was a one-compartment open model with single path Michaelis-Menten elimination. Dosage was adjusted to body weight. The effects of age, gender, alcohol consumption, cigarette smoking, dosage form, concurrent treatment with beta-agonists or steroids, outpatient dosing, and plasma caffeine concentration on maximum elimination rate (Vm) and Michaelis constant for theophylline metabolism (Km) were investigated. Hypothesis testing produced a final model in which Km = 0.42 (mg/l), and Vm (mg/kg per day) was based on cigarette smoking and dosage form, with Vm = 7.54 + 2.01 (smoking) + 1.08 (euphylline). Estimated coefficients of variation for interindividual variability in Km and Vm were 162.6% and 48.1%, respectively. Residual variability in dosage rates was estimated as 0.90 mg/kg per day. The identification of factors influencing theophylline disposition should prove useful for the a priori design of theophylline dosage regimens and monitoring of drug levels during therapy.

A liquid chromatographic method for determination of theophylline in serum and capillary blood--a comparison

A simple, fast and reliable liquid chromatographic method for the determination of theophylline in serum and capillary blood after a solid phase extraction is described for therapeutic drug monitoring. The employment of capillary blood permits the determination of an individual drug profile and other pharmacokinetic studies in neonates and infants. There were no differences in venous- and capillary-blood levels but these values compared poorly with those in serum. An adjustment of the results by correction of the different volumes of serum and blood by haematocrit was unsuccessful. Differences in the binding of theophylline to erythrocytes could be an explanation for the differences in serum at blood levels of theophylline.