Navigating the landscape of parameter identifiability methods: A workflow recommendation for model development.

In pharmacometric modeling, it is often important to know whether the data is sufficiently rich to identify the parameters of a proposed model. While it may be possible to assess this based on the results of a model fit, it is often difficult to disentangle identifiability issues from other model fitting and numerical problems. Furthermore, it can be of value to ascertain identifiability beforehand. This paper compares four methods for parameter identifiability, namely Differential Algebra for Identifiability of SYstems (DAISY), the sensitivity matrix method (SMM), Aliasing, and the Fisher information matrix method (FIMM). We discuss the characteristics of the methods and apply them to a set of applications, consisting of frequently used PK model structures, with suitable dosing regimens and sampling times. These applications were selected to validate the methods and demonstrate their usefulness. While traditional identifiability analysis provides a categorical result [PLoS One, 6, 2011, e27755; CPT Pharmacometrics Syst Pharmacol, 8, 2019, 259; Bioinformatics, 30, 2014, 1440], we argue that in practice a continuous scale better reflects the limitations on the data and is more informative. The methods were generally consistent in their evaluation of the applications. The Fisher information matrix method seemed to provide the most consistent answers. All methods provided information on the parameters that were unidentifiable. Some of the results were unexpected, indicating identifiability issues where none were foreseen, but could be explained upon further analysis. This illustrated the usefulness of identifiability assessment.

A population K-PD model analysis of long-term testosterone inhibition in prostate cancer patients undergoing intermittent androgen deprivation therapy.

Intermittent androgen deprivation therapy with gonadotropin-releasing-hormone (GnRH) agonists can prevent or delay disease progression and development of castration resistant prostate cancer for subpopulations of prostate cancer patients. It may also reduce risk and severity of side effects associated with chemical castration in prostate cancer (PCa) patients. One of the earliest comprehensively documented clinical trials on this was reported in a Canadian patient population treated with leuprorelin preceded by a lead-in with cyproterone acetate. A systems-based mixed effect analysis of testosterone response in active and recovery phases allows inference of new information from this patient population. Efficacy of androgen deprivation therapy is presumed to depend on a treshold value for testosterone at the nadir, below which no additional beneficial effects on PSA reponse can be expected, and occurance of testosterone breakthroughs during active therapy. The present analysis results in a mixed effect model, incorporating GnRH receptor activation, testosterone turnover and feedback mechanisms, describing and predicting testosterone inhibition under intermittent androgen deprivation therapy on the individual and population level, during multiple years of therapy. Testosterone levels in these patients decline over time with an estimated first order rate constant of 0.083 year-1(T1/2 = 8.4 y), with a substantial distribution among this patient population, compared to the general population. PCa patients leaving the trial due to unmanageble PSA relapse appear to have slightly higher testosterone levels at the nadir than sustained responders. These findings are expected to contribute to an increased understanding of the role of testosterone in long term disease progression of prostate cancer.

Pharmacokinetic and pharmacodynamic modeling of the metastin/kisspeptin analog, TAK-448, for its anti-tumor efficacy in a rat xenograft model.

TAK-448 is the investigational metastin/kisspeptin analog, which is known to have an anti-tumor effect through suppression of androgen hormones (luteinizing hormone and testosterone) levels. This study developed pharmacokinetic-pharmacodynamic (PK/PD) models of TAK-448 and leuprorelin acetate (TAP-144) in a rat vertebral-cancer of the prostate (VCaP) androgen-sensitive prostate cancer xenograft model to quantitatively assess and compare the anti-tumor effects of both drugs. A potential contribution of the hormone-independent direct effects of TAK-448 to the tumor growth inhibition was also investigated in the in vivo rat xenograft model, because our in vitro experiments revealed that TAK-448 may also directly suppress VCaP cellular proliferation. The PK/PD model successfully described the time course of tumor growth inhibition after drug treatment as well as the development of resistance to the inhibition of androgen hormones, following drug treatment or castration. The EC50 of the hormone-dependent inhibitory effect of TAK-448 was much lower than that of TAP-144, and TAK-448 also has a faster onset of anti-tumor effect than TAP-144, demonstrating that TAK-448 has a stronger overall anti-tumor effect than TAP-144. In addition, model inference, by incorporating a hormone-independent inhibition pathway of TAK-448 into the PK-PD model, suggested that such a direct inhibition pathway for TAK-448 cannot be excluded, as also indicated by in vitro studies, but its EC50 would be approximately three orders of magnitude higher than that of the hormone-dependent pathway. This study helps to understand the potential and mechanism of TAK-448 as a prostate cancer treatment.

Influence of body composition on disposition of the highly fat distributed compound as analysed by physiologically based pharmacokinetic (PBPK) modeling and simulation.

A recent study suggested that the pharmacokinetics (PK) of highly fat distributed compounds can be affected by acute changes in the volume of adipose tissue. The present study investigates possible influences of body composition on the disposition of the highly lipophilic compound TAK-357 in two rat strains. Physiologically based PK (PBPK) modeling and simulation was applied on single and multiple dose PK data of TAK-357 in obese Wistar fatty rats and Wistar lean rats having approximately 45% and 13% body fat, respectively. The observed effects of an elevated fat mass in Wistar fatty rats on the plasma concentrations appeared to be partly compensated for by other differences between the two rat strains. A decrease in the tissue to blood partition coefficients under high body fat conditions was identified as another factor contributing to the difference in PK. A higher lipid content in the plasma in high body fat animals may result in relatively lower tissue to blood partition coefficients. PBPK-based simulations indicate that the plasma concentrations of lipophilic compounds in high body fat conditions can differ by up to two-times at steady-state. This confirms that there is only a small impact of body composition change on the plasma concentration of highly lipophilic drugs and that the need for therapeutic dose adjustments may be limited.

Impact of acute fat mobilisation on the pharmacokinetics of the highly fat distributed compound TAK-357, investigated by physiologically based pharmacokinetic (PBPK) modeling and simulation.

In a dog toxicokinetic study, an unusual plasma concentration increase of the highly lipophilic compound TAK-357 was observed 2 weeks after termination of a 2-week repeated dosing in one dog with acute body weight loss. The present study investigates the cause of this increase. A physiologically based pharmacokinetic (PBPK) model was constructed using the rat and dog pharmacokinetic data. Using the constructed model, the TAK-357 concentration profile in the case of body weight change was simulated. The PBPK model-derived simulation suggested that redistribution from adipose tissues to plasma due to a loss of body fat caused the observed concentration increase of TAK-357 in dog plasma. The analysis demonstrates that the disposition of a highly lipophilic and fat-distributed compound can be affected by acute changes in adipose tissue mass. PBPK modeling and simulation proved to be efficient tools for the quantitative hypothesis testing of apparently atypical PK phenomena resulting from acute physiological changes.

Evaluation of the long-term durability and glycemic control of fasting plasma glucose and glycosylated hemoglobin for pioglitazone in Japanese patients with type 2 diabetes.

BACKGROUND: This study applied a pharmacodynamic model-based approach to evaluate the long-term durability and glycemic control of pioglitazone in comparison with other oral glucose-lowering drugs in Japanese type 2 diabetes mellitus (T2DM) patients. SUBJECTS AND METHODS: Japanese T2DM patients were enrolled in a prospective, randomized, open-label, blinded-end point study and received pioglitazone with or without other oral glucose-lowering drugs (excluding another thiazolidinedione [TZD]) (n=293) or oral glucose-lowering drugs excluding TZD (n=294). Treatment was adjusted to achieve glycosylated hemoglobin (HbA1c) <6.9%, and samples for fasting plasma glucose (FPG) and HbA1c were collected over 2.5-4 years. A simultaneous cascading indirect response model structure was applied to describe the time course of FPG and HbA1c. HbA1c levels were described using both an FPG-dependent and an FPG-independent function. To account for titration, drug effects for both treatment groups were implemented using a time-dependent Emax model. RESULTS: Pioglitazone was superior in both time to maximum effect and the magnitude of reduction achieved in FPG and HbA1c. A greater reduction in median FPG (-21 mg/dL vs. -9 mg/dL) was observed with pioglitazone (P<0.05). Maximum drug effect for FPG was predicted to occur earlier (11 months) for pioglitazone than for the control group (14 months). The simulated additional reduction in FPG and HbA1c achieved with pioglitazone was predicted to be maintained beyond the currently observed study duration. CONCLUSIONS: Pioglitazone was found to result in improved glycemic control and durability compared with control treatment. This model-based approach enabled the quantification of differences in FPG and HbA1c for both treatment groups and simulation to evaluate longer-term durability on FPG and HbA1c.

Pharmacokinetics of recombinant human antithrombin in delivery and surgery patients with hereditary antithrombin deficiency.

Population pharmacokinetic (PK) analyses were conducted to refine dosing recommendations for recombinant human anti-thrombin therapy in surgery and delivery patients with hereditary antithrombin deficiency (HD). Single-dose PK data from patients with HD and nonlinear mixed-effects modeling were used to devise a dosing regimen to target antithrombin (AT) activity levels between 80% and 120% of normal. External validation with data from a phase 3 trial confirmed the correctness of a covariate-free model for surgery patients, but dosing adjustment was necessary for delivery patients. After different covariates were tested, the model was updated to incorporate the influential covariate, delivery. Simulations were used to develop a therapeutic drug-monitoring scenario that results in steady state AT activity levels within the target range as quickly as practically feasible. Data from a second clinical trial provided additional external validation and confirmed the accuracy of the dosing model for both groups of patients.

A model-based approach to analyze the influence of UGT2B15 polymorphism driven pharmacokinetic differences on the pharmacodynamic response of the PPAR agonist sipoglitazar.

The pharmacokinetics of sipoglitazar, a peroxisome proliferator activated receptor agonist, are subject to high inter-individual variability resulting from a polymorphism of the UGT2B15 genotype. The aim of the current analysis was to apply a PK-PD model-based approach to evaluate the influence of UGT2B15 driven pharmacokinetic differences on the clinical response. Efficacy and safety of sipoglitazar compared to placebo were assessed in Type 2 Diabetes Mellitus patients in two Phase II randomized, double-blind studies (sipoglitazar once daily: 8, 16, 32, or 64 mg; sipoglitazar twice daily: 16 or 32 mg; rosiglitazone 8 mg once daily and placebo for 13 weeks; n = 780). Changes in fasting plasma glucose (FPG) and glycosylated hemoglobin (HbA1c) levels over time were described as a function of individual drug exposure using a simultaneous, cascading indirect response model structure. The effects on FPG and HbA1c could successfully be described for placebo, rosiglitazone, and sipoglitazar treated groups in all three UGT2B15 genotypes. Differences in drug effects between genotypes were fully explained by differences in drug exposure. The current PK-PD analysis confirms that UGT2B15 genotype is a major determinant for differences in FPG and HbA1c response to sipoglitazar treatment between Type 2 Diabetes mellitus patients, due to related differences in drug exposure.

Evaluation of the impact of UGT polymorphism on the pharmacokinetics and pharmacodynamics of the novel PPAR agonist sipoglitazar.

Sipoglitazar is a peroxisome proliferator-activated receptor α, δ, and γ agonist. During phase I, a wide distribution of clearance between individuals was observed. Hypothesized to result from a polymorphism in the uridine 5'-diphospate-glucuronosyltransferase (UGT)2B15 enzyme, pharmacogenetic samples were collected from each individual for genotyping UGT2B15 in a subsequent phase I trial in healthy subjects (n = 524) and in 2 phase II trials in type 2 diabetes subjects (n = 627), total genotype frequency was as follows: *1/*1 (22%), *1/*2 (51%), and *2/*2 (27%). The impact of genotype on exposure was assessed using a pharmacokinetic modeling approach; the influence of genotype on efficacy was evaluated using 12-week HbA1c change from baseline. Model analysis demonstrated UGT2B15 genotype accounted significantly for the variability in sipoglitazar clearance; however, a small fraction of subjects had a clearance that could not be explained entirely by genotype. HbA1c drop increased with daily drug dose. When stratified by both dose and genotype, HbA1c drop was larger in the UGT2B15*2/*2 compared with UGT2B15*1/*1 and UGT2B15*1/*2 genotypes (P < .05). In summary, UGT2B15 genotype is a strong predictor for sipoglitazar clearance; a greater clinical response observed in the UGT2B15*2/*2 genotype appears to confirm this. However, overlap in individual rates of clearance across genotypes remains after accounting for genotype.

Predictive Performance of a Recently Developed Population Pharmacokinetic Model for Morphine and its Metabolites in New Datasets of (Preterm) Neonates, Infants and Children.

BACKGROUND AND OBJECTIVE: Model validation procedures are crucial when models are to be used to develop new dosing algorithms. In this study, the predictive performance of a previously published paediatric population pharmacokinetic model for morphine and its metabolites in children younger than 3 years (original model) is studied in new datasets that were not used to develop the original model. METHODS: Six external datasets including neonates and infants up to 1 year were obtained from four different research centres. These datasets contained postoperative patients, ventilated patients and patients on extracorporeal membrane oxygenation (ECMO) treatment. Basic observed versus predicted plots, normalized prediction distribution error analysis, model refitting, bootstrap analysis, subpopulation analysis and a literature comparison of clearance predictions were performed with the new datasets to evaluate the predictive performance of the original morphine pharmacokinetic model. RESULTS: The original model was found to be stable and the parameter estimates were found to be precise. The concentrations predicted by the original model were in good agreement with the observed concentrations in the four datasets from postoperative and ventilated patients, and the model-predicted clearances in these datasets were in agreement with literature values. In the datasets from patients on ECMO treatment with continuous venovenous haemofiltration (CVVH) the predictive performance of the model was good as well, whereas underprediction occurred, particularly for the metabolites, in patients on ECMO treatment without CVVH. CONCLUSION: The predictive value of the original morphine pharmacokinetic model is demonstrated in new datasets by the use of six different validation and evaluation tools. It is herewith justified to undertake a proof-of-principle approach in the development of rational dosing recommendations - namely, performing a prospective clinical trial in which the model-based dosing algorithm is clinically evaluated.

Predictive performance of a recently developed population pharmacokinetic model for morphine and its metabolites in new datasets of (preterm) neonates, infants and children.

BACKGROUND AND OBJECTIVE: Model validation procedures are crucial when models are to be used to develop new dosing algorithms. In this study, the predictive performance of a previously published paediatric population pharmacokinetic model for morphine and its metabolites in children younger than 3 years (original model) is studied in new datasets that were not used to develop the original model. METHODS: Six external datasets including neonates and infants up to 1 year were obtained from four different research centres. These datasets contained postoperative patients, ventilated patients and patients on extracorporeal membrane oxygenation (ECMO) treatment. Basic observed versus predicted plots, normalized prediction distribution error analysis, model refitting, bootstrap analysis, subpopulation analysis and a literature comparison of clearance predictions were performed with the new datasets to evaluate the predictive performance of the original morphine pharmacokinetic model. RESULTS: The original model was found to be stable and the parameter estimates were found to be precise. The concentrations predicted by the original model were in good agreement with the observed concentrations in the four datasets from postoperative and ventilated patients, and the model-predicted clearances in these datasets were in agreement with literature values. In the datasets from patients on ECMO treatment with continuous venovenous haemofiltration (CVVH) the predictive performance of the model was good as well, whereas underprediction occurred, particularly for the metabolites, in patients on ECMO treatment without CVVH. CONCLUSION: The predictive value of the original morphine pharmacokinetic model is demonstrated in new datasets by the use of six different validation and evaluation tools. It is herewith justified to undertake a proof-of-principle approach in the development of rational dosing recommendations - namely, performing a prospective clinical trial in which the model-based dosing algorithm is clinically evaluated.

Pharmacokinetic-pharmacodynamic analysis of the static allodynia response to pregabalin and sildenafil in a rat model of neuropathic pain.

The objective of this study was to develop a pharmacokinetic-pharmacodynamic (PK-PD) model of the static allodynia response to pregabalin with and without sildenafil in a chronic constriction injury model of neuropathic pain. Six treatment groups were evaluated every 30 min for 6 h. Rats were treated with either 1) a saline infusion; 2) a 2-h pregabalin infusion at 4 mgxkg(-1)xh(-1); 3) a 2-h pregabalin infusion at 10 mgxkg(-1)xh(-1); 4) a 2.2-mg loading dose + 12 mgxkg(-1)xmin(-1) infusion of sildenafil; 5) a 2-h pregabalin infusion at 1.6 mgxkg(-1)xh(-1) with sildenafil; and 6) a 2-h infusion of pregabalin at 4 mgxkg(-1)xh with sildenafil. The static allodynia endpoint was modeled by using three population PD approaches: 1) the behavior of the injured paw using a three-category ordinal logistic regression model; 2) paw withdrawal threshold (PWT) (g) between the injured and uninjured paw using the Hill equation with a baseline function; and 3) the baseline normalized difference in PWT between the injured and uninjured paw. The categorical model showed a significant shift in the concentration-response relationship of pregabalin to lower concentrations with concomitant sildenafil. Likewise, the continuous PK-PD models demonstrated a reduction in the EC(50) of pregabalin necessary for PD response in the presence of sildenafil. The difference-transformed PD model resulted in a 54.4% (42.3-66.9%) decrease in EC(50), whereas the percentage-transformed PD model demonstrated a 53.5% (42.7-64.3%) shift. It is concluded from these studies that there is a synergistic PD interaction between pregabalin and sildenafil.

Pharmacokinetics of clindamycin in pregnant women in the peripartum period.

The study presented here was performed to determine the pharmacokinetics of intravenously administered clindamycin in pregnant women. Seven pregnant women treated with clindamycin were recruited. Maternal blood and arterial and venous umbilical cord blood samples were obtained. Maternal clindamycin concentrations were analyzed by nonlinear mixed-effects modeling with the NONMEM program. The data were best described by a linear three-compartment model. The clearance and the volume of distribution at steady state were 10.0 liters/h and 6.32 x 10(3) liters, respectively. Monte Carlo simulations were performed to determine the area under the concentration curve (AUC) for the free (unbound) drug (f) in maternal serum for 24 h divided by the MIC (fAUC(0-24)/MIC). At a MIC of 0.5 mg/liter, which is the EUCAST breakpoint, the attainment at the lower 95% confidence interval (CI) was 24.6 if the level of protein binding was 65%, and this value concurred well with the target value of 27. However, for higher degrees of protein binding, as has been described in the literature, the attainment was lower, down to 10.2 for a protein binding level of 85% (lower 95% CI). The concentrations in umbilical cord blood were lower than those in maternal blood. The concentration-time profiles in maternal serum indicate that the level of exposure to clindamycin may be too low in these patients. Together with the lower concentrations in umbilical cord blood, this finding suggests that the current dosing regimen may not be adequate to protect all neonates from group B streptococcal disease.

Population pharmacokinetic model of the pregabalin-sildenafil interaction in rats: application of simulation to preclinical PK-PD study design.

PURPOSE: Preliminary evidence has suggested a synergistic interaction between pregabalin and sildenafil for the treatment of neuropathic pain. The focus of this study was to determine the influence of sildenafil on the pharmacokinetics (PK) of pregabalin with the objective of informing the design of a quantitative pharmacodynamic (PD) study. METHODS: The pharmacokinetics were determined in rats following 2-hr intravenous infusions of pregabalin at doses of 4 mg/kg/hr and 10 mg/kg/hr with and without a sildenafil bolus (2.2 mg) and steady state infusion (12 mg/kg/hr for 6 h). This PK model was utilized in a preclinical trial simulation with the aim of selecting the optimal sampling strategy to characterize the PK-PD profile in a future study. Eight logistically feasible PK sampling strategies were simulated in NONMEM and examined through trial simulation techniques. RESULTS: A two-compartment population PK model best described pregabalin pharmacokinetics. Significant model covariates included either a binary effect of sildenafil administration (30.2% decrease in clearance) or a concentration-dependent effect due to sildenafil's active metabolite. CONCLUSIONS: Analysis of simulations indicated that three post-PD samples had the best cost/benefit ratio by providing a significant increase in the precision (and minor improvement in bias) of both PK and PD parameters compared with no PK sampling.

Morphine glucuronidation in preterm neonates, infants and children younger than 3 years.

BACKGROUND AND OBJECTIVE: A considerable amount of drug use in children is still unlicensed or off-label. In order to derive rational dosing schemes, the influence of aging on glucuronidation capacity in newborns, including preterms, infants and children under the age of 3 years was studied using morphine and its major metabolites as a model drug. METHODS: A population pharmacokinetic model was developed with the nonlinear mixed-effects modelling software NONMEM V, on the basis of 2159 concentrations of morphine and its glucuronides from 248 infants receiving intravenous morphine ranging in bodyweight from 500 g to 18 kg (median 2.8 kg). The model was internally validated using normalized prediction distribution errors. RESULTS: Formation clearances of morphine to its glucuronides and elimination clearances of the glucuronides were found to be primarily influenced by bodyweight, which was parameterized using an allometric equation with an estimated exponential scaling factor of 1.44. Additionally, a postnatal age of less than 10 days was identified as a covariate for formation clearance to the glucuronides, independent of birthweight or postmenstrual age. Distribution volumes scaled linearly with bodyweight. CONCLUSIONS: Model-based simulations show that in newborns, including preterms, infants and children under the age of 3 years, a loading dose in microg/kg and a maintenance dose expressed in microg/kg1.5/h, with a 50% reduction of the maintenance dose in newborns younger than 10 days, results in a narrow range of morphine and metabolite serum concentrations throughout the studied age range. Future pharmacodynamic investigations are needed to reveal target concentrations in this population, after which final dosing recommendations can be made.

Pharmacokinetics of amoxicillin in maternal, umbilical cord, and neonatal sera.

The pharmacokinetics of amoxicillin were studied in umbilical cord and neonatal sera relative to maternal concentrations in prevention of neonatal group B streptococcus infection. The subjects were 44 pregnant women receiving amoxicillin as 1 or 2 g as an intravenous infusion. To measure the concentrations, blood samples were obtained from the mother, the arterial and venous umbilical cord, and the neonate. The pharmacokinetics were characterized by a five-compartment model by using nonlinear mixed-effects (population) modeling. The population estimates for the clearance, central volume of distribution, and the two peripheral maternal volumes of distribution were 19.7 +/- 0.99 liters/h, 6.40 +/- 0.61 liters, and 5.88 +/- 0.83 liters (mean +/- standard error), respectively. The volume of distribution of the venous umbilical cord and the neonatal volume of distribution were 3.40 liters and 11.9 liters, respectively. The pharmacokinetic parameter estimates were used to simulate the concentration-time profiles in maternal, venous umbilical cord, and neonatal sera. The peak concentration in the venous umbilical cord serum was 18% of the maternal peak concentration. It was reached 3.3 min after the maternal peak concentration. The concentration-time profile in neonatal serum was determined by the profile in venous umbilical cord serum, which in turn depended on the profile in maternal serum. Furthermore, the simulated concentrations in maternal, venous umbilical cord, and neonatal sera exceeded the MIC for group B streptococcus for more than 90% of the 4-h dosing interval. In a first approximation, the 2-g infusion to the mother appears to be adequate for the prevention of group B streptococcal disease. However, to investigate the efficacy of the prophylaxis, further studies of the interindividual variability in pharmacokinetics are indicated.

Amoxicillin pharmacokinetics in pregnant women with preterm premature rupture of the membranes.

OBJECTIVE: This study was undertaken to study the pharmacokinetics of intravenously administered amoxicillin in pregnant women with preterm premature rupture of the membranes (PPROM). STUDY DESIGN: Healthy women with PPROM were recruited and treated with amoxicillin (2 g initially and 1 g subsequently). Blood samples were obtained from the opposite arm and concentrations determined with the use of high-pressure liquid chromatography. Nonlinear mixed-effects modeling was performed in nonlinear mixed effect (population) modeling. RESULTS: The pharmacokinetics of 17 patients was described by a 3-compartment model. Clearance and volume of distribution at steady state were 22.8 L/h and 21.4 L/h, respectively, similar to values in nonpregnant individuals. There was little variability between patients. No relationship was observed between values of individual pharmacokinetic parameters and various covariates. CONCLUSION: The pharmacokinetics of amoxicillin in pregnant patients with PPROM similar to nonpregnant individuals. Given the small interindividual variability in pharmacokinetics, no dose adjustments are required to account for differences between subjects under normal circumstances.

Pharmacokinetics of penicillin G in infants with a gestational age of less than 32 weeks.

The pharmacokinetics of penicillin G were studied in 20 preterm neonates with a gestational age of less than 32 weeks on day 3 of life by using a population approach performed with the nonlinear mixed effects modeling program NONMEM. The derived population estimates and the correlation matrix of these estimates were used to perform Monte Carlo simulations and obtain the probability of target attainment (PTA). The pharmacokinetics of penicillin G were best described by a two-compartment pharmacokinetic model. The population estimates of the central volume of distribution, the peripheral volume of distribution, the intercompartmental clearance, and the total body clearance were 0.359 +/- 0.06 liter, 0.152 +/- 0.03 liter, 0.774 +/- 0.28 liter/h, and 0.103 +/- 0.01 liter/h (mean +/- standard error), respectively. The terminal half-life was 3.9 h. Clearance increased significantly with increasing birth weight. Assuming the percentage of time that the concentration of unbound drug remained above the MIC of 50% for preterm neonates, the susceptibility breakpoint based on a 100% PTA was < or =4 mg/liter, simulating the current dosing regimen of 50,000 U/kg every 12 h. This regimen is therefore adequate for the treatment of common infections in neonates on the third day of life.

Application of the convection-dispersion equation to modelling oral drug absorption.

Models of systemic drug absorption after oral administration are frequently based on a direct or a delayed first-order rate process. In practice, the use of the first-order approach to predict drug concentrations in blood plasma frequently yields a considerable mismatch between predicted and measured concentration profiles. This is particularly true for the upswing of the plasma concentration after oral administration. The current investigation explores an alternative model to describe the absorption rate based on the convection-dispersion equation describing the transport of chemicals through the GI tract. This equation is governed by two parameters, transport velocity and dispersion coefficient. One solution of this equation for a specific set of initial and boundary conditions was used to model absorption of paracetamol in a 22-year-old man after oral administration. The GI-tract passage rate in this subject was influenced by co-administration of drugs that stimulate or delay gastric emptying. The transport-limited absorption function is more accurate in describing the plasma concentration versus time curve after oral administration than the first-order model. Additionally, it provides a mechanistic explanation for the observed curve through the differences in GI-tract passage rate.

Pharmacokinetics and pharmacodynamics of midazolam and metabolites in nonventilated infants after craniofacial surgery.

BACKGROUND: Because information on the optimal dose of midazolam for sedation of nonventilated infants after major surgery is scant, a population pharmacokinetic and pharmacodynamic model is developed for this specific group. METHODS: Twenty-four of the 53 evaluated infants (aged 3-24 months) admitted to the Pediatric Surgery Intensive Care Unit, who required sedation judged necessary on the basis of the COMFORT-Behavior score and were randomly assigned to receive midazolam, were included in the analysis. Bispectral Index values were recorded concordantly. Population pharmacokinetic and pharmacodynamic modeling was performed using NONMEM V (GloboMax LLC, Hanover, MD). RESULTS: For midazolam, total clearance was 0.157 l/min, central volume was 3.8 l, peripheral volume was 30.2 l, and intercompartmental clearance was 0.30 l/min. Assuming 60% conversion of midazolam to 1-OH-midazolam, the volume of distribution for 1-OH-midazolam and 1-OH-midazolamglucuronide was 6.7 and 1.7 l, and clearance was 0.21 and 0.047 l/min, respectively. Depth of sedation using COMFORT-Behavior could adequately be described by a baseline, postanesthesia effect (Emax model) and midazolam effect (Emax model).The midazolam concentration at half maximum effect was 0.58 mum with a high interindividual variability of 89%. Using the Bispectral Index, in 57% of the infants the effect of midazolam could not be characterized. CONCLUSION: In nonventilated infants after major surgery, midazolam clearance is two to five times higher than in ventilated children. From the model presented, the recommended initial dosage is a loading dose of 1 mg followed by a continuous infusion of 0.5 mg/h during the night for a COMFORT-Behavior of 12-14 in infants aged 1 yr. Large interindividual variability warrants individual titration of midazolam in these children.