Population pharmacokinetics of perhexiline from very sparse, routine monitoring data

Using NONMEM, the population pharmacokinetics of perhexiline were studied in 88 patients (34 F, 54 M) who were being treated for refractory angina. Their mean +/- SD (range) age was 75 +/- 9.9 years (46-92), and the length of perhexiline treatment was 56 +/- 77 weeks (0.3-416). The sampling time after a dose was 14.1 +/- 21.4 hours (0.5-200), and the perhexiline plasma concentrations were 0.39 +/- 0.32 mg/L (0.03-1.56). A one-compartment model with first-order absorption was fitted to the data using the first-order (FO) approximation. The best model contained 2 subpopulations (obtained via the $MIXTURE subroutine) of 77 subjects (subgroup A) and 11 subjects (subgroup B) that had typical values for clearance (CL/F) of 21.8 L/h and 2.06 L/h, respectively. The volumes of distribution (V/F) were 1470 L and 260 L, respectively, which suggested a reduction in presystemic metabolism in subgroup B. The interindividual variability (CV%) was modeled logarithmically and for CL/F ranged from 69.1% (subgroup A) to 86.3% (subgroup B). The interindividual variability in V/F was 111%. The residual variability unexplained by the population model was 28.2%. These results confirm and extend the existing pharmacokinetic data on perhexiline, especially the bimodal distribution of CL/F manifested via an inherited deficiency in hepatic and extrahepatic CYP2D6 activity.

Population pharmacokinetics of magnesium in preeclampsia

OBJECTIVE: The purpose of this study was to determine the population pharmacokinetics of magnesium from sparse observational data in patients with preeclampsia. STUDY DESIGN: Serum magnesium concentrations (1-11 per patient) were obtained retrospectively from the records of 116 patients with preeclampsia who had a loading dose of magnesium sulfate (16 or 20 mmol), followed by a maintenance dose (1 mmol/h) over an average of 28 hours. Population clearance, volume of distribution, and the baseline magnesium concentration were estimated using the NONMEM program. RESULTS: The following population typical values, together with the interpatient variability (expressed as coefficient of variation) were obtained with the use of a 1-compartment model: systemic clearance, 4.28 L/h (37.3%); volume of distribution, 32.3 L (32.1%); baseline concentration, 0.811 mmol/L (18.5%). The average half-life was 5.2 hours. Clonus was not obtunded in 4 patients whose serum magnesium concentrations were similar to the average concentration of 1.7 mmol/L. The variability remaining unexplained after the population model was fitted to the data was 6.5% to 10.8%. CONCLUSION: This study extended knowledge of the pharmacokinetic disposition of magnesium in preeclampsia. The results are potentially useful for the calculation of loading and maintenance doses, particularly when the relationship between serum concentration and effect in preeclampsia is clarified.

Mechanisms of drought-induced population decline in an endangered wallaby

The bridled nailtail wallaby is restricted to one locality in central Queensland, Australia. The population declined severely during a major drought between 1991 and 1995. We investigated age-specific covariates of survival and proximate causes of mortality from 1994 to 1997, using mark-recapture and radio-tagging techniques at two study sites. Using a matrix population model, we also modelled the effect of drought on age-specific survival and the intrinsic rate of population increase,;,. The only significant covariate of survival for adults was a measure of health unrelated to drought. Rainfall, food, predator activity, year, sex and habitat were not associated with variation in adult survival. Juvenile survival was negatively affected by drought, and predation was the proximate cause of most juvenile deaths. The matrix projection model showed that the observed juvenile survivorship during the drought was low enough to have produced a population decline, although fecundity and survival of other age classes was high throughout the study. (C) 2001 Elsevier Science Ltd. All rights reserved.

Modelling the population dynamics of the Queensland fruit fly, Bactrocera (Dacus) tryoni: a cohort-based approach incorporating the effects of weather

Queensland fruit fly, Bactrocera (Dacus) tryoni (QFF) is arguably the most costly horticultural insect pest in Australia. Despite this, no model is available to describe its population dynamics and aid in its management. This paper describes a cohort-based model of the population dynamics of the Queensland fruit fly. The model is primarily driven by weather variables, and so can be used at any location where appropriate meteorological data are available. In the model, the life cycle is divided into a number of discreet stages to allow physiological processes to be defined as accurately as possible. Eggs develop and hatch into larvae, which develop into pupae, which emerge as either teneral females or males. Both females and males can enter reproductive and over-wintering life stages, and there is a trapped male life stage to allow model predictions to be compared with trap catch data. All development rates are temperature-dependent. Daily mortality rates are temperature-dependent, but may also be influenced by moisture, density of larvae in fruit, fruit suitability, and age. Eggs, larvae and pupae all have constant establishment mortalities, causing a defined proportion of individuals to die upon entering that life stage. Transfer from one immature stage to the next is based on physiological age. In the adult life stages, transfer between stages may require additional and/or alternative functions. Maximum fecundity is 1400 eggs per female per day, and maximum daily oviposition rate is 80 eggs/female per day. The actual number of eggs laid by a female on any given day is restricted by temperature, density of larva in fruit, suitability of fruit for oviposition, and female activity. Activity of reproductive females and males, which affects reproduction and trapping, decreases with rainfall. Trapping of reproductive males is determined by activity, temperature and the proportion of males in the active population. Limitations of the model are discussed. Despite these, the model provides a useful agreement with trap catch data, and allows key areas for future research to be identified. These critical gaps in the current state of knowledge exist despite over 50 years of research on this key pest. By explicitly attempting to model the population dynamics of this pest we have clearly identified the research areas that must be addressed before progress can be made in developing the model into an operational tool for the management of Queensland fruit fly. (C) 2003 Published by Elsevier B.V.

A D-optimal designed population pharmacokinetic study of itraconazole capsules and solution in adults with cystic fibrosis

Background: Oral itraconazole (ITRA) is used for the treatment of allergic bronchopulmonary aspergillosis in patients with cystic fibrosis (CF) because of its antifungal activity against Aspergillus species. ITRA has an active hydroxy-metabolite (OH-ITRA) which has similar antifungal activity. ITRA is a highly lipophilic drug which is available in two different oral formulations, a capsule and an oral solution. It is reported that the oral solution has a 60% higher relative bioavailability. The influence of altered gastric physiology associated with CF on the pharmacokinetics (PK) of ITRA and its metabolite has not been previously evaluated. Objectives: 1) To estimate the population (pop) PK parameters for ITRA and its active metabolite OH-ITRA including relative bioavailability of the parent after administration of the parent by both capsule and solution and 2) to assess the performance of the optimal design. Methods: The study was a cross-over design in which 30 patients received the capsule on the first occasion and 3 days later the solution formulation. The design was constrained to have a maximum of 4 blood samples per occasion for estimation of the popPK of both ITRA and OH-ITRA. The sampling times for the population model were optimized previously using POPT v.2.0.[1] POPT is a series of applications that run under MATLAB and provide an evaluation of the information matrix for a nonlinear mixed effects model given a particular design. In addition it can be used to optimize the design based on evaluation of the determinant of the information matrix. The model details for the design were based on prior information obtained from the literature, which suggested that ITRA may have either linear or non-linear elimination. The optimal sampling times were evaluated to provide information for both competing models for the parent and metabolite and for both capsule and solution simultaneously. Blood samples were assayed by validated HPLC.[2] PopPK modelling was performed using FOCE with interaction under NONMEM, version 5 (level 1.1; GloboMax LLC, Hanover, MD, USA). The PK of ITRA and OH‑ITRA was modelled simultaneously using ADVAN 5. Subsequently three methods were assessed for modelling concentrations less than the LOD (limit of detection). These methods (corresponding to methods 5, 6 & 4 from Beal[3], respectively) were (a) where all values less than LOD were assigned to half of LOD, (b) where the closest missing value that is less than LOD was assigned to half the LOD and all previous (if during absorption) or subsequent (if during elimination) missing samples were deleted, and (c) where the contribution of the expectation of each missing concentration to the likelihood is estimated. The LOD was 0.04 mg/L. The final model evaluation was performed via bootstrap with re-sampling and a visual predictive check. The optimal design and the sampling windows of the study were evaluated for execution errors and for agreement between the observed and predicted standard errors. Dosing regimens were simulated for the capsules and the oral solution to assess their ability to achieve ITRA target trough concentration (Cmin,ss of 0.5-2 mg/L) or a combined Cmin,ss for ITRA and OH-ITRA above 1.5mg/L. Results and Discussion: A total of 241 blood samples were collected and analysed, 94% of them were taken within the defined optimal sampling windows, of which 31% where taken within 5 min of the exact optimal times. Forty six per cent of the ITRA values and 28% of the OH-ITRA values were below LOD. The entire profile after administration of the capsule for five patients was below LOD and therefore the data from this occasion was omitted from estimation. A 2-compartment model with 1st order absorption and elimination best described ITRA PK, with 1st order metabolism of the parent to OH-ITRA. For ITRA the clearance (ClItra/F) was 31.5 L/h; apparent volumes of central and peripheral compartments were 56.7 L and 2090 L, respectively. Absorption rate constants for capsule (kacap) and solution (kasol) were 0.0315 h-1 and 0.125 h-1, respectively. Comparative bioavailability of the capsule was 0.82. There was no evidence of nonlinearity in the popPK of ITRA. No screened covariate significantly improved the fit to the data. The results of the parameter estimates from the final model were comparable between the different methods for accounting for missing data, (M4,5,6)[3] and provided similar parameter estimates. The prospective application of an optimal design was found to be successful. Due to the sampling windows, most of the samples could be collected within the daily hospital routine, but still at times that were near optimal for estimating the popPK parameters. The final model was one of the potential competing models considered in the original design. The asymptotic standard errors provided by NONMEM for the final model and empirical values from bootstrap were similar in magnitude to those predicted from the Fisher Information matrix associated with the D-optimal design. Simulations from the final model showed that the current dosing regimen of 200 mg twice daily (bd) would provide a target Cmin,ss (0.5-2 mg/L) for only 35% of patients when administered as the solution and 31% when administered as capsules. The optimal dosing schedule was 500mg bd for both formulations. The target success for this dosing regimen was 87% for the solution with an NNT=4 compared to capsules. This means, for every 4 patients treated with the solution one additional patient will achieve a target success compared to capsule but at an additional cost of AUD $220 per day. The therapeutic target however is still doubtful and potential risks of these dosing schedules need to be assessed on an individual basis. Conclusion: A model was developed which described the popPK of ITRA and its main active metabolite OH-ITRA in adult CF after administration of both capsule and solution. The relative bioavailability of ITRA from the capsule was 82% that of the solution, but considerably more variable. To incorporate missing data, using the simple Beal method 5 (using half LOD for all samples below LOD) provided comparable results to the more complex but theoretically better Beal method 4 (integration method). The optimal sparse design performed well for estimation of model parameters and provided a good fit to the data.

Paediatric population pharmacokinetics of itraconazole and its active metabolite hydroxy-itraconazole in cystic fibrosis and bone marrow transplant patients

Objectives: The aim of the study was to characterise the population pharmacokinetics (popPK) properties of itraconazole (ITRA) and its active metabolite hydroxy-ITRA in a representative paediatric population of cystic fibrosis (CF) and bone marrow transplant (BMT) patients. The goals were to determine the relative bioavailability between the two oral formulations, and to explore improved dosage regimens in these patients. Methods: All paediatric patients with CF taking oral ITRA for the treatment of allergic bronchopulmonary aspergillosis and patients undergoing BMT who were taking ITRA for prophylaxis of any fungal infection were eligible for the study. A minimum of two blood samples were drawn after the capsules and also after switching to oral solution, or vice versa. ITRA and hydroxy-ITRA plasma concentrations were measured by HPLC[1]. A nonlinear mixed-effect modelling approach (NONMEM 5.1.1) was used to describe the PK of ITRA and hydroxy-ITRA simultaneously. Simulations were used to assess dosing strategies in these patients. Results: Forty-nine patients (29CF, 20 BMT) were recruited to the study who provided 227 blood samples for the population analysis. A 1-compartment model with 1st order absorption and elimination best described ITRA kinetics, with 1st order conversion to hydroxy-ITRA. For ITRA, the apparent clearance (ClItra/F) and volume of distribution (Vitra/F) was 35.5L/h and 672L, respectively; the absorption rate constant for the capsule formulation was 0.0901 h-1 and for the oral solution formulation it was 0.959 h-1. The capsule comparative bioavailability (vs. solution) was 0.55. For hydroxy-ITRA, the apparent volume of distribution and clearance were 10.6 L and 5.28 L/h, respectively. Of several screened covariates only allometrically scaled total body weight significantly improved the fit to the data. No difference between the two populations was found. Conclusion: The developed popPK model adequately described the pharmacokinetics of ITRA and hydroxy-ITRA in paediatric patients with CF and patients undergoing BMT. High inter-patient variability confirmed previous data in CF[2], leukaemia and BMT[3] patients. From the population model, simulations showed the standard dose (5 mg/kg/day) needs to be doubled for the solution formulation and even 4 times more given of the capsules to achieve an adequate target therapeutic trough plasma concentration of 0.5 mg/L[4] in these patients.

Adopting turtle excluder devices in Australia and the United States: What are the differences in technology transfer, promotion, and acceptance?

Turtle excluder devices (TEDs) are being trialed on a voluntary basis in many Australian prawn (shrimp) trawl fisheries to reduce sea turtle captures. Analysis of TED introductions into shrimp trawl fisheries of the United States provided major insights into why conflicts occurred between shrimpers, conservationists, and government agencies. A conflict over the introduction and subsequent regulation of TEDs occurred because the problem and the solution were perceived differently by the various stakeholders. Attempts to negotiate and mediate the conflict broke down, resulting in litigation against the U.S. government by conservationists and shrimpers. Litigation was not an efficient resolution to the sea turtle-TED-trawl conflict but it appears that litigation was the only remaining path of resolution once the issue became polarized. We review two major Australian trawl fisheries to identify any significant differences in circumstances that may affect TED acceptance. Australian trawl fisheries are structured differently and good communication occurs between industry and researchers. TEDs are being introduced as mature technology. Furthermore, bycatch issues are of increasing concern to all stakeholders. These factors, combined with insights derived from previous conflicts concerning TEDs in the United Stares, increase the possibilities that TEDs will be introduced to Australian fishers with better acceptance.

Modelling dispersal behaviour on a fractal landscape

We use a spatially explicit population model to explore the population consequences of different habitat selection mechanisms on landscapes with fractal variation in habitat quality. We consider dispersal strategies ranging from random walks to perfect habitat selectors for two species of arboreal marsupial, the greater glider (Petauroides volans) and the mountain brushtail possum (Trichosurus caninus). In this model increasing habitat selection means individuals obtain higher quality territories, but experience increased mortality during dispersal. The net effect is that population sizes are smaller when individuals actively select habitat. We find positive relationships between habitat quality and population size can occur when individuals do not use information about the entire landscape when habitat quality is spatially autocorrelated. We also find that individual behaviour can mitigate the negative effects of spatial variation on population average survival and fecundity. (C) 1998 Elsevier Science Ltd. All rights reserved.

Swapping space for time and unfair tests of ecological models

Testing ecological models for management is an increasingly important part of the maturation of ecology as an applied science. Consequently, we need to work at applying fair tests of models with adequate data. We demonstrate that a recent test of a discrete time, stochastic model was biased towards falsifying the predictions. If the model was a perfect description of reality, the test falsified the predictions 84% of the time. We introduce an alternative testing procedure for stochastic models, and show that it falsifies the predictions only 5% of the time when the model is a perfect description of reality. The example is used as a point of departure to discuss some of the philosophical aspects of model testing.

Modelling life history strategies with capture-recapture data: Evolutionary demography of the water skink Eulamprus tympanum

Matrix population models, elasticity analysis and loop analysis can potentially provide powerful techniques for the analysis of life histories. Data from a capture-recapture study on a population of southern highland water skinks (Eulamprus tympanum) were used to construct a matrix population model. Errors in elasticities were calculated by using the parametric bootstrap technique. Elasticity and loop analyses were then conducted to identify the life history stages most important to fitness. The same techniques were used to investigate the relative importance of fast versus slow growth, and rapid versus delayed reproduction. Mature water skinks were long-lived, but there was high immature mortality. The most sensitive life history stage was the subadult stage. It is suggested that life history evolution in E. tympanum may be strongly affected by predation, particularly by birds. Because our population declined over the study, slow growth and delayed reproduction were the optimal life history strategies over this period. Although the techniques of evolutionary demography provide a powerful approach for the analysis of life histories, there are formidable logistical obstacles in gathering enough high-quality data for robust estimates of the critical parameters.

Why Scotland needs more than just a new migration policy

The Scottish Executive has adopted a policy to combat Scotland's declining population by encouraging inward migration. Using a multi-state population model this paper presents nine long-term population scenarios for Scotland using three net international migration levels and three fertility paths. The results show inward migration can slow population decline but makes little difference to population ageing. Without a higher fertility rate Scotland's population will become demographically unsustainable. Our simulations show that a higher fertility rate substantially reduces the future ageing.

A theory for optimal monitoring of marine reserves

Monitoring of marine reserves has traditionally focused on the task of rejecting the null hypothesis that marine reserves have no impact on the population and community structure of harvested populations. We consider the role of monitoring of marine reserves to gain information needed for management decisions. In particular we use a decision theoretic framework to answer the question: how long should we monitor the recovery of an over-fished stock to determine the fraction of that stock to reserve? This exposes a natural tension between the cost (in terms of time and money) of additional monitoring, and the benefit of more accurately parameterizing a population model for the stock, that in turn leads to a better decision about the optimal size for the reserve with respect to harvesting. We found that the optimal monitoring time frame is rarely more than 5 years. A higher economic discount rate decreased the optimal monitoring time frame, making the expected benefit of more certainty about parameters in the system negligible compared with the expected gain from earlier exploitation.

Should managed populations be monitored every year?

We often need to estimate the size of wild populations to determine the appropriate management action, for example, to set a harvest quota. Monitoring is usually planned under the assumption that it must be carried out at fixed intervals in time, typically annually, before the harvest quota is set. However, monitoring can be very expensive, and we should weigh the cost of monitoring against the improvement that it makes in decision making. A less costly alternative to monitoring annually is to predict the population size using a population model and information from previous surveys. In this paper, the problem of monitoring frequency is posed within a decision-theory framework. We discover that a monitoring regime that varies according to the state of the system call outperform fixed-interval monitoring This idea is illustrated using data for a red kangaroo (Macropits rufus) population in South Australia. Whether or not one should monitor in a given year is dependent on the estimated population density in the previous year, the uncertainty in that population estimate, and past rainfall. We discover that monitoring is-important when a model-based prediction of population density is very uncertain. This may occur if monitoring has not taken place for several years, or if rainfall has been above average. Monitoring is also important when prior information suggests that the population is near a critical threshold in population abundance. However, monitoring is less important when the optimal management action would not be altered by new information.

Dosing of gentamicin in patients with end-stage renal disease receiving hemodialysis

The aim of this study was to evaluate dosing schedules of gentamicin in patients with end-stage renal disease and receiving hemodialysis. Forty-six patients were recruited who received gentamicin while on hemodialysis. Each patient provided approximately 4 blood samples at various times before and after dialysis for analysis of plasma gentamicin concentrations. A population pharmacokinetic model was constructed using NONMEM (version 5). The clearance of gentamicin during dialysis was 4.69 L/h and between dialysis was 0.453 L/h. The clearance between dialysis was best described by residual creatinine clearance (as calculated using the Cockcroft and Gault equation), which probably reflects both lean mass and residual clearance mechanisms. Simulation from the final population model showed that predialysis dosing has a higher probability of achieving target maximum concentration (C-max) concentrations (> 8 mg/L) within acceptable exposure limits (area under the concentration-time curve [AUC] values > 70 and < 120 mg.h/L per 24 hours) than postdialysis dosing.

Population balance model of in vivo neutrophil formation following bone marrow rescue therapy

In this paper, we develop a simple four parameter population balance model of in vivo neutrophil formation following bone marrow rescue therapy. The model is used to predict the number and type of neutrophil progenitors required to abrogate the period of severe neutropenia that normally follows a bone marrow transplant. The estimated total number of 5 billion neutrophil progenitors is consistent with the value extrapolated from a human trial. The model provides a basis for designing ex vivo expansion protocols.