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.

The Role of Cannabinoids in the Neurobiology of Sensory Gating: A firing rate model study

Gating of sensory (e.g. auditory) information has been demonstrated as a reduction in the auditory-evoked potential responses recorded in the brain of both normal animals and human subjects. Auditory gating is perturbed in schizophrenic patients and pharmacologically by drugs such as amphetamine, phencyclidine or ketamine, which precipitate schizophrenic-like symptoms in normal subjects. The neurobiological basis underlying this sensory gating can be investigated using local field potential recordings from single electrodes. In this paper we use such technology to investigate the role of cannabinoids in sensory gating. Cannabinoids represent a fundamentally new class of retrograde messengers which are released postsynaptically and bind to presynaptic receptors. In this way they allow fine-tuning of neuronal response, and in particular can lead to so-called depolarization-induced suppression of inhibition (DSI). Our experimental results show that application of the exogenous cannabinoid WIN55, 212-2 can abolish sensory gating as measured by the amplitude of local field responses in rat hippocampal region CA3. Importantly we develop a simple firing rate population model of CA3 and show that gating is heavily dependent upon the presence of a slow inhibitory (GABAB) pathway. Moreover, a simple phenomenological model of cannabinoid dynamics underlying DSI is shown to abolish gating in a manner consistent with our experimental findings.

Role of climate on barn owls road-kill likelihood and the effect on population viability in future climate change

Roads represent a new source of mortality due to animal-vehicle risk of collision threatening log-term populations’ viability. Risk of road-kill depends on species sensitivity to roads and their specific life-history traits. The risk of road mortality for each species depends on the characteristics of roads and bioecological characteristics of the species. In this study we intend to know the importance of climatic parameters (temperature and precipitation) together with traffic and life history traits and understand the role of drought in barn owl population viability, also affected by road mortality in three scenarios: high mobility, high population density and the combination of previous scenarios (mixed) (Manuscript). For the first objective we correlated the several parameters (climate, traffic and life history traits). We used the most correlated variables to build a predictive mixed model (GLMM) the influence of the same. Using a population model we evaluated barn owl population viability in all three scenarios. Model revealed precipitation, traffic and dispersal have negative relationship with road-kills, although the relationship was not significant. Scenarios showed different results, high mobility scenario showed greater population depletion, more fluctuations over time and greater risk of extinction. High population density scenario showed a more stable population with lower risk of extinction and mixed scenario showed similar results as first scenario. Climate seems to play an indirect role on barn owl road-kills, it may influence prey availability which influences barn owl reproductive success and activity. Also, high mobility scenario showed a greater negative impact on viability of populations which may affect their ability and resilience to other stochastic events. Future research should take in account climate and how it may influence species life cycles and activity periods for a more complete approach of road-kills. Also it is important to make the best mitigation decisions which might include improving prey quality habitat.

Estimating Population Trends in American Woodcock (Scolopax Minor) Using Population Reconstruction Models

The American woodcock (Scolopax minor) population index in North America has declined 0.9% a year since 1968 prompting managers to identify priority information and management needs for the species (Sauer et al 2008). Managers identified a need for a population model that better informs on the status of American woodcock populations (Case et al. 2010). Population reconstruction techniques use long-term age-at-harvest data and harvest effort to estimate abundances with error estimates. Four new models were successfully developed using survey data (1999 to 2013). The optimal model estimates sex specific harvest probability for adult females at 0.148 (SE = 0.017) and all other age-sex cohorts at 0.082 (SE = 0.008) for the most current year 2013. The model estimated a yearly survival rate of 0.528 (SE = 0.008). Total abundance ranged from 5,206,000 woodcock in 2007 to 6,075,800 woodcock in 1999. This study represents the first population estimates of woodcock populations.

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.

Avaliação do stock de Abalones: uma pesca sustentável

Jornadas "Ciência nos Açores – que futuro? Tema Ciências Naturais e Ambiente", Ponta Delgada, 7-8 de Junho de 2013.

Evaluation of the exploitation potential of red swamp crayfish (Procambarus clarkii) in Portugal

Dissertação para obtenção do Grau de Mestre em Engenharia do Ambiente, perfil de Ordenamento do Território e Impactes Ambientais

Relationships between imatinib plasma levels and efficacy

Purpose: While imatinib has revolutionized the treatment of chronic myeloid leukaemia (CML) and gastrointestinal stromal tumors (GIST), its pharmacokinetic-pharmacodynamic relationships have been poorly studied. This study aimed to explore the issue in oncologic patients, and to evaluate the specific influence of the target genotype in a GIST subpopulation. Patients and methods: Data from 59 patients (321 plasma samples) were collected during a previous pharmacokinetic study. Based on a population model purposely developed, individual post-hoc Bayesian estimates of pharmacokinetic parameters were derived, and used to estimate drug exposure (AUC; area under curve). Free fraction parameters were deduced from a model incorporating plasma alpha1-acid glycoprotein levels. Associations between AUC (or clearance) and therapeutic response (coded on a 3-point scale), or tolerability (4-point scale), were explored by ordered logistic regression. Influence of KIT genotype on response was also assessed in GIST patients. Results: Total and free drug exposure correlated with the number of side effects (p < 0.005). A relationship with response was not evident in the whole patient set (with good-responders tending to receive lower doses and bad-responders higher doses). In GIST patients however, higher free drug exposure predicted better responses. A strong association was notably observed in patients harboring an exon 9 mutation or a wild type KIT, known to decrease tumor sensitivity towards imatinib (p < 0.005). Conclusions: Our results are arguments to further evaluate the potential benefit of a therapeutic monitoring program for imatinib. Our data also suggest that stratification by genotype will be important in future trials.

Semigroup analysis of structured populations with distributed states-at-birth arising in mathematical aquaculture

Motivated by the modelling of structured parasite populations in aquaculture we consider a class of physiologically structured population models, where individuals may be recruited into the population at different sizes in general. That is, we consider a size-structured population model with distributed states-at-birth. The mathematical model which describes the evolution of such a population is a first order nonlinear partial integro-differential equation of hyperbolic type. First, we use positive perturbation arguments and utilise results from the spectral theory of semigroups to establish conditions for the existence of a positive equilibrium solution of our model. Then we formulate conditions that guarantee that the linearised system is governed by a positive quasicontraction semigroup on the biologically relevant state space. We also show that the governing linear semigroup is eventually compact, hence growth properties of the semigroup are determined by the spectrum of its generator. In case of a separable fertility function we deduce a characteristic equation and investigate the stability of equilibrium solutions in the general case using positive perturbation arguments.

Relationship of imatinib-free plasma levels and target genotype with efficacy and tolerability

Imatinib has revolutionised the treatment of chronic myeloid leukaemia (CML) and gastrointestinal stromal tumours (GIST). Using a nonlinear mixed effects population model, individual estimates of pharmacokinetic parameters were derived and used to estimate imatinib exposure (area under the curve, AUC) in 58 patients. Plasma-free concentration was deduced from a model incorporating plasma levels of alpha(1)-acid glycoprotein. Associations between AUC (or clearance) and response or incidence of side effects were explored by logistic regression analysis. Influence of KIT genotype was also assessed in GIST patients. Both total (in GIST) and free drug exposure (in CML and GIST) correlated with the occurrence and number of side effects (e.g. odds ratio 2.7+/-0.6 for a two-fold free AUC increase in GIST; P<0.001). Higher free AUC also predicted a higher probability of therapeutic response in GIST (odds ratio 2.6+/-1.1; P=0.026) when taking into account tumour KIT genotype (strongest association in patients harbouring exon 9 mutation or wild-type KIT, known to decrease tumour sensitivity towards imatinib). In CML, no straightforward concentration-response relationships were obtained. Our findings represent additional arguments to further evaluate the usefulness of individualizing imatinib prescription based on a therapeutic drug monitoring programme, possibly associated with target genotype profiling of patients.

Developing methodologies to allow the robust prediction dynamic context dominated by perturbations

El principal objectiu del projecte era desenvolupar millores conceptuals i metodològiques que permetessin una millor predicció dels canvis en la distribució de les espècies (a una escala de paisatge) derivats de canvis ambientals en un context dominat per pertorbacions. En un primer estudi, vàrem comparar l'eficàcia de diferents models dinàmics per a predir la distribució de l'hortolà (Emberiza hortulana). Els nostres resultats indiquen que un model híbrid que combini canvis en la qualitat de l'hàbitat, derivats de canvis en el paisatge, amb un model poblacional espacialment explícit és una aproximació adequada per abordar canvis en la distribució d'espècies en contextos de dinàmica ambiental elevada i una capacitat de dispersió limitada de l'espècie objectiu. En un segon estudi abordarem la calibració mitjançant dades de seguiment de models de distribució dinàmics per a 12 espècies amb preferència per hàbitats oberts. Entre les conclusions extretes destaquem: (1) la necessitat de que les dades de seguiment abarquin aquelles àrees on es produeixen els canvis de qualitat; (2) el biaix que es produeix en la estimació dels paràmetres del model d'ocupació quan la hipòtesi de canvi de paisatge o el model de qualitat d'hàbitat són incorrectes. En el darrer treball estudiarem el possible impacte en 67 espècies d’ocells de diferents règims d’incendis, definits a partir de combinacions de nivells de canvi climàtic (portant a un augment esperat de la mida i freqüència d’incendis forestals), i eficiència d’extinció per part dels bombers. Segons els resultats dels nostres models, la combinació de factors antropogènics del regim d’incendis, tals com l’abandonament rural i l’extinció, poden ser més determinants per als canvis de distribució que els efectes derivats del canvi climàtic. Els productes generats inclouen tres publicacions científiques, una pàgina web amb resultats del projecte i una llibreria per a l'entorn estadístic R.