Population pharmacokinetic analysis of mycophenolic acid in renal transplant recipients following oral administration of mycophenolate mofetil

Aim To develop a population pharmacokinetic model for mycophenolic acid in adult kidney transplant recipients, quantifying average population pharmacokinetic parameter values, and between- and within-subject variability and to evaluate the influence of covariates on the pharmacokinetic variability. Methods Pharmacokinetic data for mycophenolic acid and covariate information were previously available from 22 patients who underwent kidney transplantation at the Princess Alexandra Hospital. All patients received mycophenolate mofetil 1 g orally twice daily. A total of 557 concentration-time points were available. Data were analysed using the first-order method in NONMEM (version 5 level 1.1) using the G77 FORTRAN compiler. Results The best base model was a two-compartment model with a lag time (apparent oral clearance was 271 h(-1), and apparent volume of the central compartment 981). There was visual evidence of complex absorption and time-dependent clearance processes, but they could not be successfully modelled in this study. Weight was investigated as a covariate, but no significant relationship was determined. Conclusions The complexity in determining the pharmacokinetics of mycophenolic acid is currently underestimated. More complex pharmacokinetic models, though not supported by the limited data collected for this study, may prove useful in the future. The large between-subject and between-occasion variability and the possibility of nonlinear processes associated with the pharmacokinetics of mycophenolic acid raise questions about the value of the use of therapeutic monitoring and limited sampling strategies.

Toward better outcomes with tacrolimus therapy: Population pharmacokinetics and individualized dosage prediction in adult liver transplantation

Patient outcomes in transplantation would improve if dosing of immunosuppressive agents was individualized. The aim of this study is to develop a population pharmacokinetic model of tacrolimus in adult liver transplant recipients and test this model in individualizing therapy. Population analysis was performed on data from 68 patients. Estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F) using the nonlinear mixed effects model program (NONMEM). Factors screened for influence on these parameters were weight, age, sex, transplant type, biliary reconstructive procedure, postoperative day, days of therapy, liver function test results, creatinine clearance, hematocrit, corticosteroid dose, and interacting drugs. The predictive performance of the developed model was evaluated through Bayesian forecasting in an independent cohort of 36 patients. No linear correlation existed between tacrolimus dosage and trough concentration (r(2) = 0.005). Mean individual Bayesian estimates for CL/F and V/F were 26.5 8.2 (SD) L/hr and 399 +/- 185 L, respectively. CL/F was greater in patients with normal liver function. V/F increased with patient weight. CL/F decreased with increasing hematocrit. Based on the derived model, a 70-kg patient with an aspartate aminotransferase (AST) level less than 70 U/L would require a tacrolimus dose of 4.7 mg twice daily to achieve a steady-state trough concentration of 10 ng/mL. A 50-kg patient with an AST level greater than 70 U/L would require a dose of 2.6 mg. Marked interindividual variability (43% to 93%) and residual random error (3.3 ng/mL) were observed. Predictions made using the final model were reasonably nonbiased (0.56 ng/mL), but imprecise (4.8 ng/mL). Pharmacokinetic information obtained will assist in tacrolimus dosing; however, further investigation into reasons for the pharmacokinetic variability of tacrolimus is required.

Prospective evaluation of a D-optimal designed population pharmacokinetic study

Recently, methods for computing D-optimal designs for population pharmacokinetic studies have become available. However there are few publications that have prospectively evaluated the benefits of D-optimality in population or single-subject settings. This study compared a population optimal design with an empirical design for estimating the base pharmacokinetic model for enoxaparin in a stratified randomized setting. The population pharmacokinetic D-optimal design for enoxaparin was estimated using the PFIM function (MATLAB version 6.0.0.88). The optimal design was based on a one-compartment model with lognormal between subject variability and proportional residual variability and consisted of a single design with three sampling windows (0-30 min, 1.5-5 hr and 11 - 12 hr post-dose) for all patients. The empirical design consisted of three sample time windows per patient from a total of nine windows that collectively represented the entire dose interval. Each patient was assigned to have one blood sample taken from three different windows. Windows for blood sampling times were also provided for the optimal design. Ninety six patients were recruited into the study who were currently receiving enoxaparin therapy. Patients were randomly assigned to either the optimal or empirical sampling design, stratified for body mass index. The exact times of blood samples and doses were recorded. Analysis was undertaken using NONMEM (version 5). The empirical design supported a one compartment linear model with additive residual error, while the optimal design supported a two compartment linear model with additive residual error as did the model derived from the full data set. A posterior predictive check was performed where the models arising from the empirical and optimal designs were used to predict into the full data set. This revealed the optimal'' design derived model was superior to the empirical design model in terms of precision and was similar to the model developed from the full dataset. This study suggests optimal design techniques may be useful, even when the optimized design was based on a model that was misspecified in terms of the structural and statistical models and when the implementation of the optimal designed study deviated from the nominal design.

Genetic diversity and population structure of Tasmanian devils, the largest marsupial carnivore

Genetic diversity and population structure were investigated across the core range of Tasmanian devils (Sarcophilus laniarius; Dasyuridae), a wide-ranging marsupial carnivore restricted to the island of Tasmania. Heterozygosity (0.386-0.467) and allelic diversity (2.7-3.3) were low in all subpopulations and allelic size ranges were small and almost continuous, consistent with a founder effect. Island effects and repeated periods of low population density may also have contributed to the low variation. Within continuous habitat, gene flow appears extensive up to 50 km (high assignment rates to source or close neighbour populations; nonsignificant values of pairwise F-ST), in agreement with movement data. At larger scales (150-250 km), gene flow is reduced (significant pairwise F-ST) but there is no evidence for isolation by distance. The most substantial genetic structuring was observed for comparisons spanning unsuitable habitat, implying limited dispersal of devils between the well-connected, eastern populations and a smaller northwestern population. The genetic distinctiveness of the northwestern population was reflected in all analyses: unique alleles; multivariate analyses of gene frequency (multidimensional scaling, minimum spanning tree, nearest neighbour); high self-assignment (95%); two distinct populations for Tasmania were detected in isolation by distance and in Bayesian model-based clustering analyses. Marsupial carnivores appear to have stronger population subdivisions than their placental counterparts.

Using the effect size to model change in preference values from descriptive health status

Objectives: This pilot study describes a modelling approach to translate group-level changes in health status into changes in preference values, by using the effect size (ES) to summarize group-level improvement. Methods: ESs are the standardized mean difference between treatment groups in standard deviation (SD) units. Vignettes depicting varying severity in SD decrements on the SF-12 mental health summary scale, with corresponding symptom severity profiles, were valued by a convenience sample of general practitioners (n = 42) using the rating scale (RS) and time trade-off methods. Translation factors between ES differences and change in preference value were developed for five mental disorders, such that ES from published meta-analyses could be transformed into predicted changes in preference values. Results: An ES difference in health status was associated with an average 0.171-0.204 difference in preference value using the RS, and 0.104-0.158 using the time trade off. Conclusions: This observed relationship may be particular to the specific versions of the measures employed in the present study. With further development using different raters and preference measures, this approach may expand the evidence base available for modelling preference change for economic analyses from existing data.

Quantitative justification for target concentration intervention - parameter variability and predictive performance using population pharmacokinetic models for aminoglycosides

Aims [1] To quantify the random and predictable components of variability for aminoglycoside clearance and volume of distribution [2] To investigate models for predicting aminoglycoside clearance in patients with low serum creatinine concentrations [3] To evaluate the predictive performance of initial dosing strategies for achieving an aminoglycoside target concentration. Methods Aminoglycoside demographic, dosing and concentration data were collected from 697 adult patients (>=20 years old) as part of standard clinical care using a target concentration intervention approach for dose individualization. It was assumed that aminoglycoside clearance had a renal and a nonrenal component, with the renal component being linearly related to predicted creatinine clearance. Results A two compartment pharmacokinetic model best described the aminoglycoside data. The addition of weight, age, sex and serum creatinine as covariates reduced the random component of between subject variability (BSVR) in clearance (CL) from 94% to 36% of population parameter variability (PPV). The final pharmacokinetic parameter estimates for the model with the best predictive performance were: CL, 4.7 l h(-1) 70 kg(-1); intercompartmental clearance (CLic), 1 l h(-1) 70 kg(-1); volume of central compartment (V-1), 19.5 l 70 kg(-1); volume of peripheral compartment (V-2) 11.2 l 70 kg(-1). Conclusions Using a fixed dose of aminoglycoside will achieve 35% of typical patients within 80-125% of a required dose. Covariate guided predictions increase this up to 61%. However, because we have shown that random within subject variability (WSVR) in clearance is less than safe and effective variability (SEV), target concentration intervention can potentially achieve safe and effective doses in 90% of patients.

Validation of a food-frequency questionnaire assessment of carotenoid and vitamin E intake using weighed food records and plasma biomarkers: The method of triads model

Background: Reliability or validity studies are important for the evaluation of measurement error in dietary assessment methods. An approach to validation known as the method of triads uses triangulation techniques to calculate the validity coefficient of a food-frequency questionnaire (FFQ). Objective: To assess the validity of an FFQ estimates of carotenoid and vitamin E intake against serum biomarker measurements and weighed food records (WFRs), by applying the method of triads. Design: The study population was a sub-sample of adult participants in a randomised controlled trial of beta-carotene and sunscreen in the prevention of skin cancer. Dietary intake was assessed by a self-administered FFQ and a WFR. Nonfasting blood samples were collected and plasma analysed for five carotenoids (alpha-carotene, beta-carotene, beta-cryptoxanthin, lutein, lycopene) and vitamin E. Correlation coefficients were calculated between each of the dietary methods and the validity coefficient was calculated using the method of triads. The 95% confidence intervals for the validity coefficients were estimated using bootstrap sampling. Results: The validity coefficients of the FFQ were highest for alpha-carotene (0.85) and lycopene (0.62), followed by beta- carotene (0.55) and total carotenoids (0.55), while the lowest validity coefficient was for lutein (0.19). The method of triads could not be used for b- cryptoxanthin and vitamin E, as one of the three underlying correlations was negative. Conclusions: Results were similar to other studies of validity using biomarkers and the method of triads. For many dietary factors, the upper limit of the validity coefficients was less than 0.5 and therefore only strong relationships between dietary exposure and disease will be detected.

Canine model for investigating the impact of oral enrofloxacin on commensal coliforms and colonization with multidrug-resistant Escherichia coli

A model was developed in dogs to determine the impact of oral enrofloxacin administration on the indigenous coliform population in the gastrointestinal tract and subsequent disposition to colonization by a strain of multidrug-resistant Escherichia coli (MDREC). Dogs given a daily oral dose of 5 mg enrofloxacin kg(-1) for 21 consecutive days showed a significant decline in faecal coliforms to levels below detectable limits by 72 In of administration. Subsequently, faecal coliforms remained suppressed throughout the period of enrofloxacin dosing. Upon termination of antibiotic administration, the number of excreted faecal coliforms slowly returned over an 8-day period, to levels comparable to those seen prior to antibiotic treatment. Enrofloxacin-treated dogs were more effectively colonized by MDREC, evidenced by a significantly increased count of MDREC in the faeces (7.1 +/- 1.5 log(10) g(-1)) compared with non-antibiotic-treated dogs (5.2 +/- 1.2; P = 0.003). Furthermore, antibiotic treatment also sustained a significantly longer period of MDREC excretion in the faeces (26.8 +/- 10.5 days) compared with animals not treated with enrofloxacin (8.5 +/- 5.4 days; P = 0.0215). These results confirm the importance of sustained delivery of an antimicrobial agent to maintain and expand the colonization potential of drug-resistant bacteria in vivo, achieved in part by reducing the competing commensal coliforms in the gastrointestinal tract to below detectable levels in the faeces. Without in vivo antimicrobial selection pressure, commensal coliforms dominated the gastrointestinal tract at the expense of the MDREC population. Conceivably, the model developed could be used to test the efficacy of novel non-antibiotic strategies aimed at monitoring and controlling gastrointestinal colonization by multidrug-resistant members of the Enterobacteriaceae that cause nosocomial infections.

Monitoring green turtles (Chelonia mydas) at a coastal foraging area in Baja California, Mexico: multiple indices to describe population status

From June 1995 to August 2002 we assessed green turtle (Chelonia mydas) population structure and survival, and identified human impact, at Bahia de los Angeles, a large bay that was once the site of the greatest sea turtle harvest rates in the Gulf of California, Mexico. Turtles were captured live with entanglement nets and mortality was quantified through stranding surveys and flipper tag recoveries. A total of 14,820 netting hours (617.5 d) resulted in 255 captures of 200 green turtles. Straight-carapace length and mass ranged from 46.0-100.0 cm (mean = 74.3 +/- 0.7 cm) and 14.5-145.0 kg (mean = 61.5 +/- 1.7 kg), respectively. The size-frequency distribution remained stable during all years and among all capture locations. Anthropogenic-derived injuries ranging from missing flippers to boat propeller scars were present in 4% of captured turtles. Remains of 18 turtles were found at dumpsites, nine stranded turtles were encountered in the study area, and flipper tags from seven turtles were recovered. Survival was estimated at 0.58 for juveniles and 0.97 for adults using a joint live-recapture and dead-recovery model (Burnham model). Low survival among juveniles, declining annual catch per unit effort, and the presence of butchered carcasses indicated human activities continue to impact green turtles at this foraging area.

Of sheep and rain: large-scale population dynamics of the red kangaroo

1. We analysed time-series data from populations of red kangaroos (Macropus rufus, Desmarest) inhabiting four areas in the pastoral zone of South Australia. We formulated a set of a priori models to disentangle the relative effects of the covariates: rainfall, harvesting, intraspecific competition, and domestic herbivores, on kangaroo population-growth rate. 2. The statistical framework allowed for spatial variation in the growth-rate parameters, response to covariates, and environmental variability, as well as spatially correlated error terms due to shared environment. 3. The most parsimonious model included all covariates but no area-specific parameter values, suggesting that kangaroo densities respond in the same way to the covariates across the areas. 4. The temporal dynamics were spatially correlated, even after taking into account the potentially synchronizing effect of rainfall, harvesting and domestic herbivores. 5. Counter-intuitively, we found a positive rather than negative effect of domestic herbivore density on the population-growth rate of kangaroos. We hypothesize that this effect is caused by sheep and cattle acting as a surrogate for resource availability beyond rainfall. 6. Even though our system is well studied, we must conclude that approximating resources by surrogates such as rainfall is more difficult than previously thought. This is an important message for studies of consumer-resource systems and highlights the need to be explicit about population processes when analysing population patterns.

Approximating persistence in a general class of population processes

We provide a general framework for estimating persistence in populations which may be affected by catastrophic events, and which are either unbounded or have very large ceilings. We model the population using a birth-death process modified to allow for downward jumps of arbitrary size. For such processes, it is typically necessary to truncate the process in order to make the evaluation of expected extinction times (and higher-order moments) computationally feasible. Hence, we give particular attention to the selection of a cut-off point at which to truncate the process, and we present a simple method for obtaining quantitative indicators of the suitability of a chosen cut-off. (c) 2005 Elsevier Inc. All rights reserved.

Increase of a Caribbean leatherback turtle Dermochelys coriacea nesting population linked to long-term nest protection

The leatherback turtle Dermochelys coriacea is considered to be at serious risk of global extinction, despite ongoing conservation efforts. Intensive long-term monitoring of a leatherback nesting population on Sandy Point (St. Croix, US Virgin Islands) offers a unique opportunity to quantify basic population parameters and evaluate effectiveness of nesting beach conservation practices. We report a significant increase in the number of females nesting annually from ca. 18-30 in the 1980s to 186 in 2001, with a corresponding increase in annual hatchling production from ca. 2000 to over 49,000. We then analyzed resighting data from 1991 to 2001 with an open robust-design capture-mark-recapture model to estimate annual nester survival and adult abundance for this population. The expected annual survival probability was estimated at ca. 0.893 (95% CL 0.87-0.92) and the population was estimated to be increasing ca. 13% pa since the early 1990s. Taken together with DNA fingerprinting that identify mother-daughter relations, our findings suggest that the increase in the size of the nesting population since 1991 was probably due to an aggressive program of beach protection and egg relocation initiated more than 20 years ago. Beach protection and egg relocation provide a simple and effective conservation strategy for this Northern Caribbean nesting population as long as adult survival at sea remains relatively high. (c) 2005 Elsevier Ltd. All rights reserved.

Population pharmacokinetics and association between A77 1726 plasma concentrations and disease activity measures following administration of leflunomide to people with rheumatoid arthritis

Aims To investigate the concentration-effect relationship and pharmacokinetics of leflunomide in patients with rheumatoid arthritis (RA). Methods Data were collected from 23 RA patients on leflunomide therapy (as sole disease modifying antirheumatic drug (DMARD)) for at least 3 months. Main measures were A77 1726 (active metabolite of leflunomide) plasma concentrations and disease activity measures including pain, duration/intensity of morning stiffness, and SF-36 survey. A population estimate was sought for apparent clearance (CL/F ) and volume of distribution was fixed (0.155 l kg(-1)). Factors screened for influence on CL/F were weight, age, gender and estimated creatinine clearance. Results Significantly higher A77 1726 concentrations were seen in patients with less swollen joints and with higher SF-36 mental summary scores than in those with measures indicating more active disease (P < 0.05); concentration-effect trends were seen with five other disease activity measures. Statistical analysis of all disease activity measures showed that mean A77 1726 concentrations in groups with greater control of disease activity were significantly higher than those in whom the measures indicated less desirable control (P < 0.05). There was large between subject variability in the dose-concentration relationship. A steady-state infusion model best described the pharmacokinetic data. Inclusion of age as a covariate decreased interindividual variability (P < 0.01), but this would not be clinically important in terms of dosage changes. Final parameter estimate (% CV interindividual variability) for CL/F was 0.0184 l h(-1) (50%) (95% CI 0.0146, 0.0222). Residual (unexplained) variability (% CV) was 8.5%. Conclusions This study of leflunomide in patients using the drug clinically indicated a concentration-effect relationship. From our data, a plasma A77 1726 concentration of 50 mg l(-1) is more likely to indicate someone with less active disease than is a concentration around 30 mg l(-1). The marked variability in pharmacokinetics suggests a place for individualized dosing of leflunomide in RA therapy.

A simulation model of cereal-legume intercropping systems for semi-arid regions II. Model application

Smallholder farmers in Africa practice traditional cropping techniques such as intercropping. Intercropping is thought to offer higher productivity and resource milisation than sole cropping. In this study, risk associated with maize-bean intercropping was evaluated by quantifying long-term yield in both intercropping and sole cropping in a semi-arid region of South Africa (Bloemfontein, Free State) with reference to rainfall variability. The crop simulation model was run with different cultural practices (planting date and plant density) for 52 summer crop growing seasons (1950/1951-2001/2002). Eighty-one scenarios, consisted of three levels of initial soil water, planting date, maize population, and bean population, were simulated. From the simulation outputs, the total land equivalent ratio (LER) was greater than one. The intercrop (equivalent to sole maize) had greater energy value (EV) than sole beans, and the intercrop (equivalent to sole beans) had greater monetary value (MV) than sole maize. From these results, it can be concluded that maize-bean intercropping is advantageous for this semi-arid region. Soil water at planting was the most important factor of all scenario factors, followed by planting date. Irrigation application at planting, November/December planting and high plant density of maize for EV and beans for MV can be one of the most effective cultural practices in the study region. With regard to rainfall variability, seasonal (October-April) rainfall positively affected EV and MV, but not LER. There was more intercrop production in La Nina years than in El Nino years. Thus, better cultural practices may be selected to maximize maize-bean intercrop yields for specific seasons in the semi-arid region based on the global seasonal outlook. (c) 2004 Elsevier B.V. All rights reserved.

Mentoring for population health in general practice divisions

This paper describes the implementation and evaluation of a three-way model of service development mentoring. This population health mentoring program was funded by the Commonwealth Department of Health and Ageing to enable staff from eight Divisions of General Practice in South Australia to gain a sound understanding of population health concepts relevant to their workplace. The distinguishing features of service development mentoring were that the learning was grounded within an individual's work setting and experience; there was an identified population health problem or issue confronting the Division of General Practice; and there was an expectation of enhanced organisational performance. A formal evaluation found a consensus among all learners that mentoring was a positive and worthwhile experience, where they had achieved what they had set out to do. Mentors found the model of learning agreeable and effective. Division executive officers recognised enhanced skills among their "learner" colleagues, and commented positively on the benefits to their organisations through the development of well researched and relevant projects, with the potential to improve the efficiency of their population health activities.