Some considerations on the design of population pharmacokinetic studies

The goal of this manuscript is to introduce a framework for consideration of designs for population pharmacokinetic orpharmacokinetic-pharmacodynamic studies. A standard one compartment pharmacokinetic model with first-order input and elimination is considered. A series of theoretical designs are considered that explore the influence of optimizing the allocation of sampling times, allocating patients to elementary designs, consideration of sparse sampling and unbalanced designs and also the influence of single vs. multiple dose designs. It was found that what appears to be relatively sparse sampling (less blood samples per patient than the number of fixed effects parameters to estimate) can also be highly informative. Overall, it is evident that exploring the population design space can yield many parsimonious designs that are efficient for parameter estimation and that may not otherwise have been considered without the aid of optimal design theory.

Estimation of pharmacokinetic parameters from non-compartmental variables using Microsoft Excel((R))

This study was conducted to develop a method, termed 'back analysis (BA)', for converting non-compartmental variables to compartment model dependent pharmacokinetic parameters for both one- and two-compartment models. A Microsoft Excel((R)) spreadsheet was implemented with the use of Solver((R)) and visual basic functions. The performance of the BA method in estimating pharmacokinetic parameter values was evaluated by comparing the parameter values obtained to a standard modelling software program, NONMEM, using simulated data. The results show that the BA method was reasonably precise and provided low bias in estimating fixed and random effect parameters for both one- and two-compartment models. The pharmacokinetic parameters estimated from the BA method were similar to those of NONMEM estimation.

Analysis of population pharmacokinetic data using NONMEM and WinBUGS

The aim of this report is to describe the use of WinBUGS for two datasets that arise from typical population pharmacokinetic studies. The first dataset relates to gentamicin concentration-time data that arose as part of routine clinical care of 55 neonates. The second dataset incorporated data from 96 patients receiving enoxaparin. Both datasets were originally analyzed by using NONMEM. In the first instance, although NONMEM provided reasonable estimates of the fixed effects parameters it was unable to provide satisfactory estimates of the between-subject variance. In the second instance, the use of NONMEM resulted in the development of a successful model, albeit with limited available information on the between-subject variability of the pharmacokinetic parameters. WinBUGS was used to develop a model for both of these datasets. Model comparison for the enoxaparin dataset was performed by using the posterior distribution of the log-likelihood and a posterior predictive check. The use of WinBUGS supported the same structural models tried in NONMEM. For the gentamicin dataset a one-compartment model with intravenous infusion was developed, and the population parameters including the full between-subject variance-covariance matrix were available. Analysis of the enoxaparin dataset supported a two compartment model as superior to the one-compartment model, based on the posterior predictive check. Again, the full between-subject variance-covariance matrix parameters were available. Fully Bayesian approaches using MCMC methods, via WinBUGS, can offer added value for analysis of population pharmacokinetic data.

Population pharmacokinetics of metformin in late pregnancy

The pharmacokinetic disposition of metformin in late pregnancy was studied together with the level of fetal exposure at birth. Blood samples were obtained in the third trimester of pregnancy from women with gestational diabetes or type 2 diabetes, 5 had a previous diagnosis of polycystic ovary syndrome. A cord blood sample also was obtained at the delivery of some of these women, and also at delivery of others who had been taking metformin during pregnancy but from whom no blood had been taken. Plasma metformin concentrations were assayed by a new, validated, reverse-phase HPLC method, A 2-compartment, extravascular maternal model with transplacental partitioning of drug to a fetal compartment was fitted to the data. Nonlinear mixed-effects modeling was performed in'NONMEM using FOCE with INTERACTION. Variability was estimated using logarithmic interindividual and additive residual variance models; the covariance between clearance and volume was modeled simultaneously. Mean (range) metformin concentrations in cord plasma and in maternal plasma were 0.81 (range, 0.1-2.6) mg/L and 1.2 (range, 0. 1-2.9) mg/L, respectively. Typical population values (interindividual variability, CV%) for allometrically scaled maternal clearance and volume of distribution were 28 L/h/70 kg (17.1%) and 190 L/70 ka (46.3%), giving a derived population-wide half-life of 5.1 hours. The placental partition coefficient for metformin was 1.07 (36.3%). Neither maternal age nor weight significantly influenced the pharmacokinetics. The variability (SD) of observed concentrations about model-predicted concentrations was 0.32 mg/L. The pharmacokinetics were similar to those in nonpregnant patients and, therefore, no dosage adjustment is warranted. Metformin readily crosses the placenta, exposing the fetus to concentrations approaching those in the maternal circulation. The sequelae to such exposure, ea, effects on neonatal obesity and insulin resistance, remain unknown.

Pharmacokinetic-pharmacodynamic modelling of QT interval prolongation following citalopram overdoses

Aims To develop a pharmacokinetic-pharmacodynamic model describing the time-course of QT interval prolongation after citalopram overdose and to evaluate the effect of charcoal on the relative risk of developing abnormal QT and heart-rate combinations. Methods Plasma concentrations and electrocardiograph (ECG) data from 52 patients after 62 citalopram overdose events were analysed in WinBUGS using a Bayesian approach. The reported doses ranged from 20 to 1700 mg and on 17 of the events a single dose of activated charcoal was administered. The developed pharmacokinetic-pharmacodynamic model was used for predicting the probability of having abnormal combinations of QT-RR, which was assumed to be related to an increased risk for torsade de pointes (TdP). Results The absolute QT interval was related to the observed heart rate with an estimated individual heart-rate correction factor [alpha = 0.36, between-subject coefficient of variation (CV) = 29%]. The heart-rate corrected QT interval was linearly dependent on the predicted citalopram concentration (slope = 40 ms l mg(-1), between-subject CV = 70%) in a hypothetical effect-compartment (half-life of effect-delay = 1.4 h). The heart-rate corrected QT was predicted to be higher in women than in men and to increase with age. Administration of activated charcoal resulted in a pronounced reduction of the QT prolongation and was shown to reduce the risk of having abnormal combinations of QT-RR by approximately 60% for citalopram doses above 600 mg. Conclusion Citalopram caused a delayed lengthening of the QT interval. Administration of activated charcoal was shown to reduce the risk that the QT interval exceeds a previously defined threshold and therefore is expected to reduce the risk of TdP.

Prediction of subcellular localization using sequence-biased recurrent networks

Motivation: Targeting peptides direct nascent proteins to their specific subcellular compartment. Knowledge of targeting signals enables informed drug design and reliable annotation of gene products. However, due to the low similarity of such sequences and the dynamical nature of the sorting process, the computational prediction of subcellular localization of proteins is challenging. Results: We contrast the use of feed forward models as employed by the popular TargetP/SignalP predictors with a sequence-biased recurrent network model. The models are evaluated in terms of performance at the residue level and at the sequence level, and demonstrate that recurrent networks improve the overall prediction performance. Compared to the original results reported for TargetP, an ensemble of the tested models increases the accuracy by 6 and 5% on non-plant and plant data, respectively.

SARS coronavirus nucleocapsid immunodominant T-cell epitope cluster is common to both exogenous recombinant and endogenous DNA-encoded immunogens

Correspondence between the T-cell epitope responses of vaccine immunogens and those of pathogen antigens is critical to vaccine efficacy. In the present study, we analyzed the spectrum of immune responses of mice to three different forms of the SARS coronavirus nucleocapsid (N): (1) exogenous recombinant protein (N-GST) with Freund's adjuvant; (2) DNA encoding unmodified N as an endogenous cytoplasmic protein (pN); and (3) DNA encoding N as a LAMP-I chimera targeted to the lysosomal MHC II compartment (p-LAMP-N). Lysosomal trafficking of the LAMP/N chimera in transfected cells was documented by both confocal and immunoelectron microscopy. The responses of the immunized mice differed markedly. The strongest T-cell IFN-gamma and CTL responses were to the LAMP-N chimera followed by the pN immunogen. In contrast, N-GST elicited strong T cell IL-4 but minimal IFN-gamma responses and a much greater antibody response. Despite these differences, however, the immunodominant T-cell ELISpot responses to each of the three immunogens were elicited by the same N peptides, with the greatest responses being generated by a cluster of five overlapping peptides, N76-114, each of which contained nonameric H2(d) binding domains with high binding scores for both class I and, except for N76-93, class II alleles. These results demonstrate that processing and presentation of N, whether exogenously or endogenously derived, resulted in common immunodominant epitopes, supporting the usefulness of modified antigen delivery and trafficking forms and, in particular, LAMP chimeras as vaccine candidates. Nevertheless, the profiles of T-cell responses were distinctly different. The pronounced Th-2 and humoral response to N protein plus adjuvant are in contrast to the balanced IFN-gamma and IL-4 responses and strong memory CTL responses to the LAMP-N chimera. (C) 2005 Elsevier Inc. All rights reserved.

Arf6-independent GPI-anchored protein-enriched early endosomal compartments fuse with sorting endosomes via a Rab5/phosphatidylinositol-3 '-kinase-dependent machinery

In the process of internalization of molecules from the extracellular milieu, a cell uses multiple endocytic pathways, consequently generating different endocytic vesicles. These primary endocytic vesicles are targeted to specific destinations inside the cell. Here, we show that GPI-anchored proteins are internalized by an Arf6-independent mechanism into GPI-anchored protein-enriched early endosomal compartments (GEECs). Internalized GPI-anchored proteins and the fluid phase are first visualized in GEECs that are acidic, primary endocytic structures, negative for early endosomal markers, Rab4, Rab5, and early endosome antigen (EEA)1. They subsequently acquire Rab5 and EEA1 before homotypic fusion with other GEECs, and heterotypic fusion with endosomes containing cargo from the clathrin-dependent endocytic pathway. Although, the formation of GEECs is unaffected by inhibition of Rab5 GTPase and phosphatidylinositol-3'-kinase (PI3K) activity, their fusion with sorting endosomes is dependent on both activities. Overexpression of Rab5 reverts PI3K inhibition of fusion, providing evidence that Rab5 effectors play important roles in heterotypic fusion between the dynamin-independent GEECs and clathrin- and dynamin-dependent sorting endosomes.

Development of a dosage strategy in patients receiving enoxaparin by continuous intravenous infusion using modelling and simulation

Aim To develop an appropriate dosing strategy for continuous intravenous infusions (CII) of enoxaparin by minimizing the percentage of steady-state anti-Xa concentration (C-ss) outside the therapeutic range of 0.5-1.2 IU ml(-1). Methods A nonlinear mixed effects model was developed with NONMEM (R) for 48 adult patients who received CII of enoxaparin with infusion durations that ranged from 8 to 894 h at rates between 100 and 1600 IU h(-1). Three hundred and sixty-three anti-Xa concentration measurements were available from patients who received CII. These were combined with 309 anti-Xa concentrations from 35 patients who received subcutaneous enoxaparin. The effects of age, body size, height, sex, creatinine clearance (CrCL) and patient location [intensive care unit (ICU) or general medical unit] on pharmacokinetic (PK) parameters were evaluated. Monte Carlo simulations were used to (i) evaluate covariate effects on C-ss and (ii) compare the impact of different infusion rates on predicted C-ss. The best dose was selected based on the highest probability that the C-ss achieved would lie within the therapeutic range. Results A two-compartment linear model with additive and proportional residual error for general medical unit patients and only a proportional error for patients in ICU provided the best description of the data. Both CrCL and weight were found to affect significantly clearance and volume of distribution of the central compartment, respectively. Simulations suggested that the best doses for patients in the ICU setting were 50 IU kg(-1) per 12 h (4.2 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). The best doses for patients in the general medical unit were 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 100 IU kg(-1) per 12 h (8.3 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). These best doses were selected based on providing the lowest equal probability of either being above or below the therapeutic range and the highest probability that the C-ss achieved would lie within the therapeutic range. Conclusion The dose of enoxaparin should be individualized to the patients' renal function and weight. There is some evidence to support slightly lower doses of CII enoxaparin in patients in the ICU setting.

Interferon-induced, antiviral human MxA protein localizes to a distinct subcompartment of the smooth endoplasmic reticulum

Human MxA protein belongs to the superfamily of dynamin-like large GTPases that are involved in intracellular membrane trafficking. MxA is induced by interferons-alpha/beta (IFN-alpha/beta) and is a key component of the antiviral response against RNA viruses. Here, we show that MxA localizes to membranes that are positive for specific markers of the smooth endoplasmic reticulum, such as Syntaxin17, but is excluded from other membrane compartments. Overexpression of MxA leads to a characteristic reorganization of the associated membranes. Interestingly, Hook3, mannose-6-phosphate receptor, and Lamp-1, which normally accumulate in cis-Golgi, endosomes, and lysosomes, respectively, also colocalized with MxA, indicating that these markers were redistributed to the MxA-positive compartment. Functional assays, however, did not show any effect of MxA on endocytosis or the secretory pathway. The present results demonstrate that MxA is an IFN-induced antiviral effector protein that resembles the constitutively expressed large GTPase family members in its capacity to localize to and reorganize intracellular membranes.

Postnatal developmental expression of the PDZ scaffolds Na+-H+ exchanger regulatory factors 1 and 2 in the rat cochlea

Sensory transduction in the mammalian cochlea requires the maintenance of specialized fluid compartments with distinct ionic compositions. This is achieved by the concerted action of diverse ion channels and transporters, some of which can interact with the PDZ scaffolds, Na+-H+ exchanger regulatory factors 1 and 2 (NHERF-1, NHERF-2). Here, we report that NHERF-1 and NHERF-2 are widely expressed in the rat cochlea, and that their expression is developmentally regulated. Reverse transcription/polymerase chain reaction (RT-PCR) and Western blotting initially confirmed the RNA and protein expression of NHERFs. We then performed immunohistochemistry on cochlea during various stages of postnatal development. Prior to the onset of hearing (P8), NHERF-1 immunolabeling was prominently polarized to the apical membrane of cells lining the endolymphatic compartment, including the stereocilia and cuticular plates of the inner and outer hair cells, marginal cells of the stria vascularis, Reissner's epithelia, and tectorial membrane. With maturation (P21, P70), NHERF-1 immunolabeling was reduced in the above structures, whereas labeling increased in the apical membrane of the interdental cells of the spiral limbus and the inner and outer sulcus cells, Hensen's cells, the inner and outer pillar cells, Deiters cells, the inner border cells, spiral ligament fibrocytes, and spiral ganglion neurons (particularly type II). NHERF-1 expression in strial basal and intermediate cells was persistent. NHERF-2 immunolabeling was similar to that for NHERF-1 during postnatal development, with the exception of expression in the synaptic regions beneath the outer hair cells. NHERF-1 and NHERF-2 co-localized with glial fibrillary acidic protein and vimentin in glia. The cochlear localization of NHERF scaffolds suggests that they play important roles in the developmental regulation of ion transport, homeostasis, and auditory neurotransmission.

The β domain is required for Vps4p oligomerization into a functionally active ATPase

Endocytic and biosynthetic trafficking pathways to the lysosome/vacuole converge at the prevacuolar endosomal compartment. During transport through this compartment, integral membrane proteins that are destined for delivery to the lysosome/vacuole lumen undergo multivesicular body (MVB) sorting into internal vesicles formed by invagination of the endosomal limiting membrane. Vps4 is an AAA family ATPase which plays a key role in MVB sorting and facilitates transport through endosomes. It possesses an N-terminal microtubule interacting and trafficking domain required for recruitment to endosomes and an AAA domain with an ATPase catalytic site. The recently solved 3D structure revealed a P domain, which protrudes from the AAA domain, and a final C-terminal alpha-helix. However, the in vivo roles of these domains are not known. In this study, we have identified motifs in these domains that are highly conserved between yeast and human Vps4. We have mutated these motifs and studied the effect on yeast Vps4p function in vivo and in vitro. We show that the P domain of the budding yeast Vps4p is not required for recruitment to endosomes, but is essential for all Vps4p endocytic functions in vivo. We also show that the P domain is required for Vps4p homotypic interaction and for full ATPase activity. In addition, it is required for interaction with Vta1p, which works in concert with Vps4p in vivo. Our studies suggest that assembly of a Vps4p oligomeric complex with full ATPase activity that interacts with Vta1p is essential for normal endosome function.

The key role of CD40 ligand in overcoming tumor-induced dendritic cell dysfunction

Overcoming dendritic cell (DC) dysfunction is a prerequisite for successful active immunotherapy against breast cancer. CD40 ligand (CD40L), a key molecule in the interface between T-lymphocytes and DCs, seems to be instrumental in achieving that goal. Commenting on our data that CD40L protects circulating DCs from apoptosis induced by breast tumor products, Lenahan and Avigan highlighted the potential of CD40L for immunotherapy. We expand on that argument by pointing to additional findings that CD40L not only rescues genuine DCs but also functionally improves populations of immature antigen-presenting cells that fill the DC compartment in patients with breast cancer.

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

Objective: The objective of the study was to characterise the population pharmacokinetic properties of itraconazole and its active metabolite hydroxyitraconazole in a representative paediatric population of cystic fibrosis and bone marrow transplant (BMT) patients and to identify patient characteristics influencing the pharmacokinetics of itraconazole. The ultimate goals were to determine the relative bioavailability between the two oral formulations (capsules vs oral solution) and to optimise dosing regimens in these patients. Methods: All paediatric patients with cystic fibrosis or patients undergoing BMT at The Royal Children's Hospital, Brisbane, QLD, Australia, who were prescribed oral itraconazole for the treatment of allergic bronchopulmonary aspergillosis (cystic fibrosis patients) or for prophylaxis of any fungal infection (BMT patients) were eligible for the study. Blood samples were taken from the recruited patients as per an empirical sampling design either during hospitalisation or during outpatient clinic visits. ltraconazole and hydroxy-itraconazole plasma concentrations were determined by a validated high-performance liquid chromatography assay with fluorometric detection. A nonlinear mixed-effect modelling approach using the NONMEM software to simultaneously describe the pharmacokinetics of itraconazole and its metabolite. Results: A one-compartment model with first-order absorption described the itraconazole data, and the metabolism of the parent drug to hydroxy-itraconazole was described by a first-order rate constant. The metabolite data also showed one-compartment characteristics with linear elimination. For itraconazole the apparent clearance (CLitraconazole) was 35.5 L/hour, the apparent volume of distribution (V-d(itraconazole)) was 672L, the absorption rate constant for the capsule formulation was 0.0901 h(-1) and for the oral solution formulation was 0.96 h-1. The lag time was estimated to be 19.1 minutes and the relative bioavailability between capsules and oral solution (F-rel) was 0.55. For the metabolite, volume of distribution, V-m/(F (.) f(m)), and clearance, CL/(F (.) fm), were 10.6L and 5.28 L/h, respectively. The influence of total bodyweight was significant, added as a covariate on CLitraconazoie/F and V-d(itraconazole)/F (standardised to a 70kg person) using allometric three-quarter power scaling on CLitraconazole/F, which therefore reflected adult values. The unexplained between-subject variability (coefficient of variation %) was 68.7%, 75.8%, 73.4% and 61.1% for CLitraconazoie/F, Vd(itraconazole)/F, CLm/(F (.) fm) and F-rel, respectively. The correlation between random effects of CLitraconazole and Vd((itraconazole)) was 0.69. Conclusion: The developed population pharmacokinetic model adequately described the pharmacokinetics of itraconazole and its active metabolite, hydroxy-itraconazole, in paediatric patients with either cystic fibrosis or undergoing BMT. More appropriate dosing schedules have been developed for the oral solution and the capsules to secure a minimum therapeutic trough plasma concentration of 0.5 mg/L for these patients.