Mixed Effects Modelling (mixed + effects_modelling)

Distribution by Scientific Domains

Kinds of Mixed Effects Modelling

  • nonlinear mixed effects modelling


  • Selected Abstracts


    Effect of obesity on serum amiodarone concentration in Japanese patients: population pharmacokinetic investigation by multiple trough screen analysis

    JOURNAL OF CLINICAL PHARMACY & THERAPEUTICS, Issue 3 2009
    H. Fukuchi MS
    Summary Objective:, To evaluate the influence of obesity on pharmacokinetics of amiodarone (AMD) using Non-Linear Mixed Effects Modelling (nonmem) in Japanese patients treated with oral therapy. Method:, Serum concentrations of AMD were determined by high performance liquid chromatography. One hundred and fifty-one trough concentrations from 23 patients receiving repetitive oral AMD were collected. Body mass index (BMI) and body fat percentage were measured. Results:, Estimates generated using nonmem indicated that the clearance of AMD was influenced by BMI, age and daily dosage of AMD. The final pharmacokinetic model was CL (L/h) = 0·16 · TBW · 0·53AGE , 65 · 0·78BMI , 25 · DD0·51, Vd (L) = 10·2 · TBW, where CL is total body clearance, TBW is total body weight (kg), DD (mg/kg/day) is daily dosage of AMD, AGE (years) ,65 = 1 for patient was 65 years old or over and 0 otherwise, BMI (kg/m2) ,25 = 1 for patient was 25 kg/m2 or over and 0 otherwise and Vd is apparent volume of distribution. The clearance of AMD decreased significantly by 22·3% with a BMI higher than 25 kg/m2. The clearance of AMD also decreased significantly by 46·9% when patient age was more than 65 years. Conclusion:, Population pharmacokinetic analysis confirms that obesity affects the pharmacokinetics of AMD. [source]


    Population modelling of the effect of inogatran, at thrombin inhibitor, on ex vivo coagulation time (APTT) in healthy subjects and patients with coronary artery disease

    BRITISH JOURNAL OF CLINICAL PHARMACOLOGY, Issue 1 2001
    Marie Cullberg
    Aims, The purpose of this study was to characterize the relationship between the degree of anticoagulation, assessed by APTT, and the plasma concentration of inogatran in healthy subjects and in patients with coronary artery disease. Methods, Data from five phase I studies in 78 healthy males and two phase II multicentre studies in 948 patients of both sexes with unstable angina pectoris or non-Q-wave myocardial infarction were evaluated. A total of 3296 pairs of concentration-APTT samples were obtained before, during, and after intravenous infusions of inogatran. Mixed effects modelling was used for population pharmacodynamic analysis of the drug effect and for describing the variability in baseline APTT. Results, The population mean baseline APTT was 29 s, but large variations between individuals (s.d. 3.6 s) were observed. The variability between studies (1.3 s) and centres (1.8 s) were of less importance, though statistically significant. APTT increased in a nonlinear manner with increasing inogatran concentration and the relationship was well described by a combined linear and Emax model. A significant part of the overall variability could be ascribed to the APTT reagent and equipment used at the different study centres. These method-dependent differences were compensated for by including the lower limit of the normal reference range as a covariate, affecting both baseline and Emax, in the model. For the typical healthy subject and patient, the method-corrected population mean parameters were: APTTbaseline 35 and 31 s, slope 8.0 and 5.8 s l µmol,1, Emax 36 and 34 s, and EC50 0.54 and 0.72 µmol l,1, respectively. The model predicted plasma concentration needed to double the APTT from the baseline value was 1.25 and 1.45 µmol l,1 in the healthy volunteer and patient, respectively. Conclusions, The nonlinear relationship between APTT and inogatran concentration in plasma was well described by a combined linear and Emax model. Pooling of data was made possible by incorporating a centre-specific characteristic of the assay method in the model. Patients had lower baseline APTT and appeared to have less pronounced effect of inogatran than young healthy subjects. [source]


    Population pharmacokinetic investigation of disopyramide by mixed effect modelling using routine clinical pharmacokinetic data in Japanese patients

    JOURNAL OF CLINICAL PHARMACY & THERAPEUTICS, Issue 4 2005
    E. Yukawa PhD
    Summary Objective:, To estimate the population pharmacokinetic parameters of disopyramide using non-linear mixed effects modelling. Method:, A total of 148 serum levels from 109 patients (61 males and 48 females) receiving disopyramide were collected. Results:, The final pharmacokinetic model was Cl (L/h) = 3·75·TBW0·567·AGE,0·374·Conc,0·719·1·48DOSE , 5, Vd (L/kg) = 4·13 and ka (h,1) = 0·363, where Cl is total body clearance, Vd is apparent volume of distribution, ka is absorption rate constant, TBW is total bodyweight (kg), AGE is age (years), Conc is the concentration of disopyramide (,g/mL), and DOSE , 5 = 1 for patient received 5 mg/kg/day of disopyramide dosage or over and 0 otherwise. Conclusion:, Application of the findings in this study to patient care may permit selection of an appropriate initial maintenance dosage to achieve target disopyramide concentrations and the desired therapeutic effect. [source]


    Clonidine disposition in children: a population analysis

    PEDIATRIC ANESTHESIA, Issue 6 2007
    AL Potts
    Background:, There are few data describing clonidine population pharmacokinetics in children (0,15 years) despite common use. Current paediatric data, described in terms of elimination half-life or Cmax and Tmax, poorly explain variability in drug responses among individuals representative of those in whom the drug will be used clinically. Methods:, Published data from four studies investigating clonidine PK after intravenous, rectal and epidural administration (n = 42) were combined with an open-label study undertaken to examine the pharmacokinetics of IV clonidine 1,2 ,g·kg,1 bolus in children after cardiac surgery (n = 30, EC approval granted). A population pharmacokinetic analysis of clonidine time-concentration profiles (380 observations) was undertaken using nonlinear mixed effects modelling. Estimates were standardised to a 70 kg adult using allometric size models. Results:, Children had a mean age of four (SD 3.6 years, range 1 week,14 years) year and weight 17.8 (SD 12.6, range 2.8,60 kg). A two compartment disposition model with first order elimination was superior to a one compartment model. Population parameter estimates (between subject variability) were clearance (CL) 14 (CV 28.3%) 1 h,1·70 kg,1, central volume of distribution (V1) 56.7 (67.5%) l·70 kg,1, inter-compartment clearance (Q) 143 (19.1%) l h,1 70 kg,1 and peripheral volume of distribution (V2) 123 (72.8%) l.70kg,1. Clearance at birth was 4.7 l·h,1·70kg,1 and matured with a half-time of 25.5 weeks to reach 85% adult rate by 1 year of age. The volumes of distribution, but not clearance, were increased after cardiac surgery (V1 180%, V2 117%). There was a lag time of 2.6 (CV 64%) min before absorption began in the rectum. The absorption half-life from the epidural space was slower than that from the rectum 1.04 (CV31%) h vs 0.28 (CV24%) h. The relative bioavailability of epidural and rectal clonidine was unity (F = 1). Conclusions:, Clearance in neonates is approximately one third that described in adults, consistent with immature clearance pathways. Maintenance dosing, which is a function of clearance, should be reduced in neonates and infants when using a target concentration approach. [source]


    Decadal dynamics of tree cover in an Australian tropical savanna

    AUSTRAL ECOLOGY, Issue 6 2009
    CAROLINE E. R. LEHMANN
    Abstract Spatio-temporal variation in tropical savanna tree cover remains poorly understood. We aimed to quantify the drivers of tree cover in tropical mesic savannas in Kakadu National Park by relating changes in tree cover over 40 years to: mean annual rainfall, fire activity, initial tree cover and prior changes in tree cover. Aerial photography, acquired in 1964, 1984 and 2004, was obtained for fifty sites in Kakadu that spanned a rainfall gradient from approximately 1200 to 1600 mm. The remotely sensed estimates of tree cover were validated via field survey. Linear mixed effects modelling and multi-model inference were used to assess the strength and form of the relationships between tree cover and predictor variables. Over the 40 years, tree cover across these savannas increased on average by 4.94 ± 0.88%, but was spatio-temporally variable. Tree cover showed a positive albeit weak trend across the rainfall gradient. The strength of this positive relationship varied over the three measurement times, and this suggests that other factors are important in controlling tree cover. Tree cover was positively related to prior tree cover, and negatively correlated with fire activity. Over 20 years tree cover was more likely to increase if (i) tree cover was initially low or (ii) had decreased in the previous 20-year interval or (iii) there had been fewer fires. Across the examined rainfall gradient, the greater variability in fire activity and inherently higher average tree cover at the wetter latitudes resulted in greater dynamism of tree cover compared with the drier latitudes. This is consistent with savanna tree cover being determined by interactions between mean annual rainfall, tree competition and frequent fire in these tropical mesic savannas. [source]


    Characterisation of the human liver in vitro metabolic pattern of artemisinin and auto-induction in the rat by use of nonlinear mixed effects modelling

    BIOPHARMACEUTICS AND DRUG DISPOSITION, Issue 2 2003
    Ulrika S.H. Svensson
    Abstract Aims: The aims of the study were to characterise the metabolic pattern of artemisinin in human and rat liver microsomes and to assess the magnitude of auto-induction in the rat. Methods: 14C-artemisinin was incubated with human liver microsomes and with liver microsomes from rats pretreated with oral artemisinin or placebo. The metabolic fate of 14C-artemisinin in microsomes from human B-lymphoblastoid cell lines transformed with CYP2A6, CYP2B6 and CYP3A4 was also investigated. The human liver microsome data and the rat liver microsomes data were analysed by nonlinear mixed effects modelling and naïve pooling using NONMEM, respectively. Results: Four metabolites were radiometrically detected in experiments with rat liver microsomes. The model that best described the data involved three primary metabolites of which one metabolite was further metabolised to a secondary metabolite. The formation of the four metabolites was induced 2.8, 7.2, 4.8 and 2.5-fold, respectively, in liver microsomes from rats pre-treated with artemisinin. Three metabolites were formed in human liver microsomes; having the same retention times as three of the metabolites formed in the rat. The final model consisted of two primary metabolites and a secondary metabolite with CYP2B6 and CYP2A6 influencing the formation rates of the major and minor primary metabolites, respectively. Conclusions: CYP2B6 and CYP2A6 activities described variability in the formation of the major and minor primary metabolites, respectively, in human liver microsomes. All artemisinin metabolic pathways in rat liver microsomes were induced in artemisinin pretreated animals. We suggest modelling as a method for the discrimination and detection of more complex metabolic patterns from in vitro metabolism rate data. Copyright © 2002 John Wiley & Sons, Ltd. [source]


    Population pharmacokinetics of intravenously and orally administered docetaxel with or without co-administration of ritonavir in patients with advanced cancer

    BRITISH JOURNAL OF CLINICAL PHARMACOLOGY, Issue 5 2010
    Stijn L. W. Koolen
    WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT , Docetaxel is an approved drug for the treatment of cancer of various primary origins. , An oral docetaxel regimen is warranted because of patient convenience and the opportunity to investigate more schedule intensive treatment regimens. , Co-administration of ritonavir significantly enhanced the apparent oral bioavailability of docetaxel. WHAT THIS STUDY ADDS , This study demonstrates that ritonavir increased the absorption of docetaxel after oral administration. , Furthermore, we showed that the clearance of docetaxel was inhibited in a concentration dependent manner. , The developed model will be used for further development of an oral docetaxel regimen. AIM Docetaxel has a low oral bioavailability due to affinity for P-glycoprotein and cytochrome P450 (CYP) 3A4 enzymes. Inhibition of the CYP3A4 enzymes by ritonavir resulted in increased oral bioavailability. The aim of this study was to develop a population pharmacokinetic (PK) model and to evaluate and quantify the influence of ritonavir on the PK of docetaxel. METHODS Data from two clinical trials were included in the data analysis, in which docetaxel (75 mg m,2 or 100 mg) had been administered intravenously or orally (10 mg or 100 mg) with or without co-administration of oral ritonavir (100 mg). Population modelling was performed using non-linear mixed effects modelling. A three-compartment model was used to describe the i.v. data. PK data after oral administration, with or without co-administration of ritonavir, were incorporated into the model. RESULTS Gut bioavailability of docetaxel increased approximately two-fold from 19 to 39% (CV 13%) with ritonavir co-administration. The hepatic extraction ratio and the elimination rate of docetaxel were best described by estimating the intrinsic clearance. Ritonavir was found to inhibit in a concentration dependent manner the intrinsic clearance of docetaxel, which was described by an inhibition constant of 0.028 µg ml,1 (CV 36%). A maximum inhibition of docetaxel clearance of more then 90% was reached. CONCLUSIONS A PK model describing both the PK of orally and intravenously administered docetaxel in combination with ritonavir, was successfully developed. Co-administration of ritonavir lead to increased oral absorption and reduced elimination rate of docetaxel. [source]


    Population pharmacokinetic analysis of cilostazol in healthy subjects with genetic polymorphisms of CYP3A5, CYP2C19 and ABCB1

    BRITISH JOURNAL OF CLINICAL PHARMACOLOGY, Issue 1 2010
    Hee-Doo Yoo
    WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT , The interindividual variability of the pharmacokinetic parameters of cilostazol is relatively large. , Cilostazol undergoes extensive hepatic metabolism via the P450 enzymes, primarily CYP3A and, to a lesser extent, CYP2C19. , Indeed, <1% of the administered dose of cilostazol is excreted unchanged in the urine. WHAT THIS STUDY ADDS , A population pharmacokinetic analysis of cilostazol was conducted to evaluate the impact of CYP3A, CYP2C19 and ABCB1 polymorphisms on cilostazol disposition in vivo. , Genetic polymorphisms of CYP3A5 and CYP2C19 explain the substantial interindividual variability in the pharmacokinetics of cilostazol. , ABCB1 genotypes do not to appear to be associated with the disposition of cilostazol. AIMS To investigate the influence of genetic polymorphisms in the CYP3A5, CYP2C19 and ABCB1 genes on the population pharmacokinetics of cilostazol in healthy subjects. METHODS Subjects who participated in four separate cilostazol bioequivalence studies with the same protocols were included in this retrospective analysis. One hundred and four healthy Korean volunteers were orally administered a single 50- or 100-mg dose of cilostazol. We estimated the population pharmacokinetics of cilostazol using a nonlinear mixed effects modelling (nonmem) method and explored the possible influence of genetic polymorphisms in CYP3A (CYP3A5*3), CYP2C19 (CYP2C19*2 and CYP2C19*3) and ABCB1 (C1236T, G2677T/A and C3435T) on the population pharmacokinetics of cilostazol. RESULTS A two-compartment model with a first-order absorption and lag time described the cilostazol serum concentrations well. The apparent oral clearance (CL/F) was estimated to be 12.8 l h,1. The volumes of the central and the peripheral compartment were characterized as 20.5 l and 73.1 l, respectively. Intercompartmental clearance was estimated at 5.6 l h,1. Absorption rate constant was estimated at 0.24 h,1 and lag time was predicted at 0.57 h. The genetic polymorphisms of CYP3A5 had a significant (P < 0.001) influence on the CL/F of cilostazol. When CYP2C19 was evaluated, a significant difference (P < 0.01) was observed among the three genotypes (extensive metabolizers, intermediate metabolizers and poor metabolizers) for the CL/F. In addition, a combination of CYP3A5 and CYP2C19 genotypes was found to be associated with a significant difference (P < 0.005) in the CL/F. When including these genotypes, the interindividual variability of the CL/F was reduced from 34.1% in the base model to 27.3% in the final model. However, no significant differences between the ABCB1 genotypes and cilostazol pharmacokinetic parameters were observed. CONCLUSIONS The results of the present study indicate that CYP3A5 and CYP2C19 polymorphisms explain the substantial interindividual variability that occurs in the metabolism of cilostazol. [source]


    Population pharmacokinetic analysis of varenicline in adult smokers

    BRITISH JOURNAL OF CLINICAL PHARMACOLOGY, Issue 5 2009
    Patanjali Ravva
    WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT? , Several clinical pharmacology studies have characterized the pharmacokinetics of varenicline in young adult and elderly smokers and subjects with impaired renal function. , Varenicline pharmacokinetics is linear over the recommended dose range. , Varenicline total clearance is linearly related to its renal clearance. , Both are progressively reduced as renal function declines, which results in a progressive increase in varenicline systemic exposure and prolonged half-life. WHAT THIS STUDY ADDS? , This work provides an integrated model-based analysis of varenicline pharmacokinetics across multiple studies in the target patient population. , The model describes the impact of patient-specific covariates, such as renal function, and provides a rationale for dose adjustment. , The resulting model also provides a means to predict individual-specific drug exposures to clinical responses in subsequent analyses. AIMS To characterize the population pharmacokinetics of varenicline and identify factors leading to its exposure variability in adult smokers. METHODS Data were pooled from nine clinical studies consisting of 1878 subjects. Models were developed to describe concentration,time profiles across individuals. Covariates were assessed using a full model approach; parameters and bootstrap 95% confidence intervals (CI) were estimated using nonlinear mixed effects modelling. RESULTS A two-compartment model with first-order absorption and elimination best described varenicline pharmacokinetics. The final population parameter estimates (95% CI) were: CL/F, 10.4 l h,1 (10.2, 10.6); V2/F, 337 l (309, 364); V3/F, 78.1 l (61.9, 98.9); Q/F, 2.08 l h,1 (1.39, 3.79); Ka, 1.69 h,1 (1.27, 2.00); and Alag, 0.43 h (0.37, 0.46). Random interindividual variances were estimated for Ka[70% coefficient of variation (CV)], CL/F (25% CV), and V2/F (50% CV) using a block covariance matrix. Fixed effect parameters were precisely estimated [most with % relative standard error < 10 and all with % relative standard error < 25], and a visual predictive check indicated adequate model performance. CL/F decreased from 10.4 l h,1 for a typical subject with normal renal function (CLcr = 100 ml min,1) to 4.4 l h,1 for a typical subject with severe renal impairment (CLcr = 20 ml min,1), which corresponds to a 2.4-fold increase in daily steady-state exposure. Bodyweight was the primary predictor of variability in volume of distribution. After accounting for renal function, there was no apparent effect of age, gender or race on varenicline pharmacokinetics. CONCLUSIONS Renal function is the clinically important factor leading to interindividual variability in varenicline exposure. A dose reduction to 1 mg day,1, which is half the recommended dose, is indicated for subjects with severe renal impairment. [source]