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.

Searching for convergence in phylogenetic Markov chain Monte Carlo

Markov chain Monte Carlo (MCMC) is a methodology that is gaining widespread use in the phylogenetics community and is central to phylogenetic software packages such as MrBayes. An important issue for users of MCMC methods is how to select appropriate values for adjustable parameters such as the length of the Markov chain or chains, the sampling density, the proposal mechanism, and, if Metropolis-coupled MCMC is being used, the number of heated chains and their temperatures. Although some parameter settings have been examined in detail in the literature, others are frequently chosen with more regard to computational time or personal experience with other data sets. Such choices may lead to inadequate sampling of tree space or an inefficient use of computational resources. We performed a detailed study of convergence and mixing for 70 randomly selected, putatively orthologous protein sets with different sizes and taxonomic compositions. Replicated runs from multiple random starting points permit a more rigorous assessment of convergence, and we developed two novel statistics, delta and epsilon, for this purpose. Although likelihood values invariably stabilized quickly, adequate sampling of the posterior distribution of tree topologies took considerably longer. Our results suggest that multimodality is common for data sets with 30 or more taxa and that this results in slow convergence and mixing. However, we also found that the pragmatic approach of combining data from several short, replicated runs into a metachain to estimate bipartition posterior probabilities provided good approximations, and that such estimates were no worse in approximating a reference posterior distribution than those obtained using a single long run of the same length as the metachain. Precision appears to be best when heated Markov chains have low temperatures, whereas chains with high temperatures appear to sample trees with high posterior probabilities only rarely. [Bayesian phylogenetic inference; heating parameter; Markov chain Monte Carlo; replicated chains.]

Bayesian population pharmacokinetic analysis of sirolimus

Objective: It is usual that data collected from routine clinical care is sparse and unable to support the more complex pharmacokinetic (PK) models that may have been reported in previous rich data studies. Informative priors may be a pre-requisite for model development. The aim of this study was to estimate the population PK parameters of sirolimus using a fully Bayesian approach with informative priors. Methods: Informative priors including prior mean and precision of the prior mean were elicited from previous published studies using a meta-analytic technique. Precision of between-subject variability was determined by simulations from a Wishart distribution using MATLAB (version 6.5). Concentration-time data of sirolimus retrospectively collected from kidney transplant patients were analysed using WinBUGS (version 1.3). The candidate models were either one- or two-compartment with first order absorption and first order elimination. Model discrimination was based on computation of the posterior odds supporting the model. Results: A total of 315 concentration-time points were obtained from 25 patients. Most data were clustered at trough concentrations with range of 1.6 to 77 hours post-dose. Using informative priors, either a one- or two-compartment model could be used to describe the data. When a one-compartment model was applied, information was gained from the data for the value of apparent clearance (CL/F = 18.5 L/h), and apparent volume of distribution (V/F = 1406 L) but no information was gained about the absorption rate constant (ka). When a two-compartment model was fitted to the data, the data were informative about CL/F, apparent inter-compartmental clearance, and apparent volume of distribution of the peripheral compartment (13.2 L/h, 20.8 L/h, and 579 L, respectively). The posterior distribution of the volume distribution of central compartment and ka were the same as priors. The posterior odds for the two-compartment model was 8.1, indicating the data supported the two-compartment model. Conclusion: The use of informative priors supported the choice of a more complex and informative model that would otherwise have not been supported by the sparse data.

Distribution of an asymmetrical copepod, Hatschekia plectropomi, on the gills of Plectropomus leopardus

Hatschekia plectropomi, an ectoparasitic copepod found on the gills, infected Plectropomus leopardus from Heron Island Reef with 100% prevalence (n = 32) and a mean +/- S.E. infection intensity of 131.9 +/- 22.1. The distribution of 4222 adult female parasites across 32 individual host fish was investigated at several organizational levels ranging from the level of holobranch pairs to that of individual filaments. Parasites demonstrated a site preference for the two central holobranchs (2 and 3). Along the lengths of hemibranchs, filaments near the dorsal and ventral ends and those in the proximity of the bend region were rarely occupied. The probability of coming into contact with a suitable attachment site and the ability to withstand ventilation forces at that site were proposed as the major factors affecting distribution. Two H. plectropomi morphotypes were identified based on the direction of body curvature. Regardless of morphotype, 99.9% of individuals were attached such that the convex side of the body was oriented towards the oncoming ventilating water currents. Further, 93.3% of individuals attached to the posterior faces of filaments, leading to a predictable pattern of attachment for this species. It is suggested that the direction of body curvature develops in response to the direction of the ventilating water currents. (c) 2006 The Fisheries Society of the British Isles.

Distribution kinetics of solutes in the isolated in-situ perfused rat head using the multiple indicator dilution technique and a physiological two-barrier model

The purpose of this study was to determine the pharmacokinetics of [C-14]diclofenac, [C-14]salicylate and [H-3]clonidine using a single pass rat head perfusion preparation. The head was perfused with 3-[N-morpholino] propane-sulfonic acid-buffered Ringer's solution. Tc-99m-red blood cells and a drug were injected in a bolus into the internal carotid artery and collected from the posterior facial vein over 28 min. A two-barrier stochastic organ model was used to estimate the statistical moments of the solutes. Plasma, interstitial and cellular distribution volumes for the solutes ranged from 1.0 mL (diclofenac) to 1.6 mL (salicylate), 2.0 mL (diclofenac) to 4.2 mL (water) and 3.9 mL (salicylate) to 20.9 mL (diclofenac), respectively. A comparison of these volumes to water indicated some exclusion of the drugs from the interstitial space and salicylate from the cellular space. Permeability-surface area (PS) products calculated from plasma to interstitial fluid permeation clearances (CLPI) (range 0.02-0.40 mL s(-1)) and fractions of solute unbound in the perfusate were in the order: diclofenac>salicylate >clonidine>sucrose (from 41.8 to 0.10 mL s(-1)). The slow efflux of diclofenac, compared with clonidine and salicylate, may be related to its low average unbound fraction in the cells. This work accounts for the tail of disposition curves in describing pharmacokinetics in the head.

The distribution and morphological characteristics of cholinergic cells in the brain of monotremes as revealed by ChAT immunohistochemistry

The present study employs choline acetyltransferase (ChAT) immunohistochemistry to identify the cholinergic neuronal population in the central nervous system of the monotremes. Two of the three extant species of monotreme were studied: the platypus (Omithorhynchus anatinus) and the short-beaked echidna (Tachyglossus aculeatus). The distribution of cholinergic cells in the brain of these two species was virtually identical. Distinct groups of cholinergic cells were observed in the striatum, basal forebrain, habenula, pontomesencephalon, cranial nerve motor nuclei, and spinal cord. In contrast to other tetrapods studied with this technique, we failed to find evidence for cholinergic cells in the hypothalamus, the parabigeminal nucleus (or nucleus isthmus), or the cerebral cortex. The lack of hypothalamic cholinergic neurons creates a hiatus in the continuous antero-posterior aggregation of cholinergic neurons seen in other tetrapods. This hiatus might be functionally related to the phenomenology of monotreme sleep and to the ontogeny of sleep in mammals, as juvenile placental mammals exhibit a similar combination of sleep elements to that found in adult monotremes. Copyright (C) 2002 S. Karger AG, Basel.

Liquid distribution as a means to describing the granule growth mechanism

Copper concentrate (chalcopyrite) was granulated in a rotating drum with a diameter of 0.3 m and a length of 0.2 m. Water was used as the binder and it was sprayed onto the powder bed with a nozzle. This material exhibited induction type behaviour, which was defined by Iveson and Litster [AIChE J. 44 (1998) 1510]. Induction type behaviour is characterized by the occurrence of an induction stage, during which the granules are gradually being compacted and little or no growth occurs. At the end of this induction stage, binder liquid is squeezed from the interior of the granules onto the granule surface and the granules are then surface-wet. This results in a rapid growth rate of the granules. Different types of experiments were conducted. The influence of the nozzle pressure and the distance from the nozzle to the powder bed on the growth behaviour of the granules as well as on the binder distribution was examined. The results of these experiments led to the postulation of a modified mechanism for induction type behaviour: it was found that after the binder was delivered, there were large granules containing a high amount of binder and small granules containing less binder. During the induction stage, the granules are compacted and binder liquid continuously appears at the surface of the large granules. These wet spots that are continuously being formed pick up the dry and small granules. When all the small granules have been picked up, further expulsion of binder liquid onto the granules' surface results in granules that remain surface-wet. This phenomenon marks the end of the induction stage and it coincides with the disappearance of the small granules. The hypothesis was tested by selectively removing the smaller granules during an experiment. As expected, this resulted in a shorter induction time.

Distribution of major, minor, and trace metals in lake environments of Antarctica

The concentrations of major, minor and trace metals were measured in water samples collected from five shallow Antarctic lakes (Carezza, Edmonson Point (No 14 and 15a), Inexpressible Island and Tarn Flat) found in Terra Nova Bay (northern Victoria Land, Antarctica) during the Italian Expeditions of 1993-2001. The total concentrations of a large suite of elements (Al, As, Ba, Ca, Cd, Ce, Co, Cr, Cs, Cu, Fe, Ga, Gd, K, La, Li, Mg, Mn, Mo, Na, Nd, Ni, Pb, Pr, Rb, Sc, Si, Sr, Ta, Ti, U, V, Y, W, Zn and Zr) were determined using spectroscopic techniques (ICP-AES, GF-AAS and ICP-MS). The results are similar to those obtained for the freshwater lakes of the Larsemann Hills, East Antarctica, and for the McMurdo Dry Valleys. Principal Component Analysis (PCA) and Cluster Analysis (CA) were performed to identify groups of samples with similar characteristics and to find correlations between the variables. The variability observed within the water samples is closely connected to the sea spray input; hence, it is primarily a consequence of geographical and meteorological factors, such as distance from the ocean and time of year. The trace element levels, in particular those of heavy metals, are very low, suggesting an origin from natural sources rather than from anthropogenic contamination.

Distribution, expression, and motif variability of ankyrin domain genes in Wolbachia pipientis

The endosymbiotic bacterium Wolbachia pipientis infects a wide range of arthropods, in which it induces a variety of reproductive phenotypes, including cytoplasmic incompatibility (CI), parthenogenesis, male killing, and reversal of genetic sex determination. The recent sequencing and annotation of the first Wolbachia genome revealed an unusually high number of genes encoding ankyrin domain (ANK) repeats. These ANK genes are likely to be important in mediating the Wolbachia-host interaction. In this work we determined the distribution and expression of the different ANK genes found in the sequenced Wolbachia wMel genome in nine Wolbachia strains that induce different phenotypic effects in their hosts. A comparison of the ANK genes of wMel and the non-CI-inducing wAu Wolbachia strain revealed significant differences between the strains. This was reflected in sequence variability in shared genes that could result in alterations in the encoded proteins, such as motif deletions, amino acid insertions, and in some cases disruptions due to insertion of transposable elements and premature stops. In addition, one wMel ANK gene, which is part of an operon, was absent in the wAu genome. These variations are likely to affect the affinity, function, and cellular location of the predicted proteins encoded by these genes.

Cytoplasmic incompatibility in Drosophila populations: influence of assortative mating on symbiont distribution

Cytoplasmic incompatibility is known to occur between strains of both Drosophila simulans and D. melanogaster. Incompatibility is associated with the infection of Drosophila with microorganismal endosymbionts. This paper reports survey work conducted on strains of D. simulans and D. melanogaster from diverse geographical locations finding that infected populations are relatively rare and scattered in their distribution. The distribution of infected populations of D. simulans appears to be at odds with deterministic models predicting the rapid spread of the infection through uninfected populations. Examination of isofemale lines from four localities in California where populations appear to be polymorphic for the infection failed to find evidence for consistent assortative mating preferences between infected and uninfected populations that may explain the basis for the observed polymorphism.

Effects of memory load and distraction on performance and event-related slow potentials in a visuospatial working memory task

Brain electrical activity related to working memory was recorded at 15 scalp electrodes during a visuospatial delayed response task. Participants (N = 18) touched the remembered position of a target on a computer screen after either a 1 or 8 sec delay. These memory trials were compared to sensory trials in which the target remained present throughout the delay and response periods. Distracter stimuli identical to the target were briefly presented during the delay on 30% of trials. Responses were less accurate in memory than sensory trials, especially after the long delay. During the delay slow potentials developed that were significantly more negative in memory than sensory trials. The difference between memory and sensory trials was greater at anterior than posterior electrodes. On trials with distracters, the slow potentials generated by memory trials showed further enhancement of negativity whereas there were minimal effects on accuracy of performance. The results provide evidence that engagement of visuospatial working memory generates slow wave negativity with a timing and distribution consistent with frontal activation. Enhanced brain activity associated with working memory is required to maintain performance in the presence of distraction. © 1997 by the Massachusetts Institute of Technology

Sequence and tissue distribution of chicken cellular nucleic acid binding protein cDNA

We sequenced cDNAs coding for chicken cellular nucleic acid binding protein (CNBP). Two slightly different variations of the open reading frame were found, each of which translates into a protein with seven zinc finger domains. The longest transcript contains an in-frame insert of 3 bp. The sequence conservation between chick CNBP cDNAs with human, rat and mouse CNBP cDNAs is extreme, especially in the coding region, where the deduced amino acid sequence identity with human, rat and mouse CNBP is 99%. CNBP-like transcripts were also found in various tissues from insect, shrimp, fish and lizard. Regions with remarkable nucleotide conservation were also found in the 3' untranslated region, indicating important functions for these regions. Quantitative reverse transcription polymerase chain reaction (RT-PCR) indicated that in the chick, CNBP is present in all tissues examined in approximately equal ratios to total RNA. RT-PCR of total RNA isolated from different phyla indicate CNBP-like proteins art widespread throughout the animal kingdom. The extraordinary level of conservation suggests an important physiological role for CNBP. (C) 1997 Elsevier Science Inc.