Population pharmacokinetic/pharmacodynamic modeling of clopidogrel in Korean healthy volunteers and stroke patients

Population pharmacokinetic (PK) and pharmacodynamic (PD) modeling of clopidogrel was developed from pooled data from healthy volunteers (n = 44) and stroke patients (n = 35). The PK modeling used plasma concentrations of the clopidogrel metabolite (SR26334), and the PD modeling used platelet aggregation. The models were developed using NONMEM and evaluated via visual predictive check (VPC). Data were analyzed by 2-compartment modeling with Erlang's absorption and first-order elimination. There was no statistically significant covariate for each model parameter. The typical point estimates of PK were k(tr) (identical transfer rate constant) = 5.97 h(-1), k(e) (elimination rate constant) = 0.126 h(-1), k(d) (distribution rate constant) = 0.212 h(-1), V(2) (volume of central compartment) = 21.0 L, and V(3) (volume of peripheral compartment) = 38.8 L. The typical point estimates of PD were k(in) (input rate) = 27.9 h(-1), E(max) (maximum effect on input rate) = 0.292 h(-1), EC(5) (0) (median effective concentration) = 0.00629 ng/mL, and BASE (predose aggregation) = 66.7%. The final model was used to estimate individual parameters using patient data and showed good predictions using VPC.     

Simultaneous pharmacokinetic and pharmacodynamic modeling

The compartmental model approach previously proposed by Sheiner et al.for the simultaneous characterization of the pharmacokinetics and pharmacodynamics of a drug substance has been rigorously tested and extended in the present report. Pharmacokinetic and pharmacodynamic equations were derived for several commonly used models, as well as alternate models that relate effect to drug concentrations in a peripheral compartment. The versatility and flexibility of these models were tested using simulation and curve-fitting procedures. The utility as well as the limitations of this modeling procedure are discussed.     

Population pharmacokinetic/pharmacodynamic modeling of valproate in Chinese children with epilepsy

AIM： To establish a population pharmacokinetic/pharmacodynamic （PK/PD） model for valproate （VPA） in children with epilepsy in China, and promote reasonable use of antiepileptic drug in clinical practice. METHODS： Sparse data of VPA serum concentrations from 417 pediatric children were collected. These patients were divided into three groups： Population PK-Index group, n = 317; Population PK-Valid group, n = 100; 115 of the total 417 subjects were also included in the Population PD group. Population PK parameters of VPA were estimated based on the data from population PK-Index group. In the validation procedure, the serum concentrations of VPA from the population PK-Valid group were predicted by base and final models respectively. To assess the accuracy and precision of the predictions, mean prediction error （MPE）, mean squared prediction error （MSPE）, root mean squared prediction error （RMSPE）, weight-residues （WRES） and its 95 % confidence intervals （95 % CI） were all calculated, then compared between the two models. For population PD group, all of the ! 15 patients were VPA monotherapy. Efficacy of epilepsy treatment was divided into 5 grades according to the percentage of seizure frequency decreased （ PSFD% ）. The value of PSFD% 100%, 75% - 100%, 50% - 75%, 25% - 50%, or less than 25 % are corresponding to grade 1 to 5. For population PD group, the quantitative relationship between the VPA serum concentrations and the probability for its efficacy score was characterized by logistic regression.[第一段]     

Population pharmacokinetic and pharmacodynamic model of norvancomycin

去甲万古霉素群体药代动力学与药效学研究@张菁$Institute of Antibiotics, Huashan Hospital, Fudan University!Shanghai 200040, China @张婴元$Institute of Antibiotics, Huashan Hospital, Fudan University!Shanghai 200040, China @施耀国$Institute of Antibiotic     

Population pharmacokinetic and pharmacodynamic analysis of ribavirin in patients with chronic hepatitis C

The population pharmacokinetics of ribavirin were assessed in patients with chronic hepatitis C virus (HCV) infection treated with interferon alpha-2b and ribavirin in four clinical efficacy studies. The authors collected 3450 ribavirin serum concentrations from 1105 patients at different treatment weeks for inclusion in the analysis. Population factors included gender, age, body weight, serum creatinine, creatinine clearance, and previous interferon treatment history. Ribavirin apparent clearance (CL/F) was calculated from individual patients' daily doses divided by concentration values, and the influence of these factors was assessed by multiple regression. Body weight, gender, age, and serum creatinine affected CL/F. Population mean CL/F estimates were 17.9 L/h (female) and 21.5 L/h (male) assuming an age of 40 years and body weight of 70 kg. Ribavirin apparent clearance increased as a function of body weight and decreased at ages greater than 40 years. Serum creatinine had little influence on CL/F, which may reflect the relatively normal renal function of these patients. Total CL/F variability was approximately 28%. The four covariates in the model explained 27% of this variability, and were thus of limited clinical significance because of the substantial residual variability not accounted for by the model. In assessing the relationship between pharmacokinetics and pharmacodynamics, the week 4 hemoglobin nadir value was negatively associated with week 4 ribavirin concentrations. The percentage of reduction from baseline was positively associated with ribavirin concentrations, although these data were highly variable. Loss of HCV-RNA at 24 weeks after completion of treatment was considered a response to interferon and ribavirin treatment in a logistic regression analysis of clinical and pharmacokinetic variables and treatment response in the interferon-naive patients. Hepatitis C virus genotype, pretreatment HCV-RNA titer, duration of treatment period, week 4 ribavirin concentration, and patient age affected the likelihood of response. Higher ribavirin concentrations at treatment week 4 were associated with a higher response rate. Variables that have predictive value for treatment outcome in patients treated with interferon and ribavirin are similar to those previously reported for interferon monotherapy.     

Population pharmacokinetic/pharmacodynamic modeling of histamine response measured by histamine iontophoresis laser Doppler

The epicutaneous histamine (EH) test is the current gold standard method for the clinical evaluation of allergic conditions. However, the EH method is limited in providing an objective and qualitative assessment of histamine pharmacodynamic response. The histamine iontophoresis with laser Doppler (HILD) monitoring method, an alternative method, allows a fixed dose of histamine to be delivered and provides an objective, continuous, and dynamic measurement of histamine epicutaneous response in children and adults. However, due to the high sampling frequency (up to 40 Hz), the output files are usually too cumbersome to be directly used for further analysis. In this study, we developed an averaging algorithm that efficiently reduces the HILD data in size. The reduced data was further analyzed and a population linked effect pharmacokinetic/pharmacodynamic (PK/PD) model was developed to describe the local histamine response. The model consisted of a one-compartment PK model and a direct-response fractional maximum effect (Emax) model. The parameter estimates were obtained as follows: absorption rate constant (ka), 0.094/min; absorption lag time (Tlag), 2.72 min; partitioning clearance from local depot to systemic circulation (CLpar), 0.0006 L/min; baseline effect (E0), 13.1 flux unit; Emax, 13.4; concentration at half maximum effect (EC50) 31.1 mg/L. Covariate analysis indicated that age and race had significant influence on Tlag and EC50, respectively.     

Relevance of pharmacokinetic and pharmacodynamic modeling to clinical care of critically ill patients

Efficacious therapy is of utmost importance to save lives and prevent bacterial resistance in critically ill patients. This review summarizes pharmacokinetic (PK) and pharmacodynamic (PD) modeling methods to optimize clinical care of critically ill patients in empiric and individualized therapy. While these methods apply to all therapeutic areas, we focus on antibiotics to highlight important applications, as emergence of resistance is a significant problem. Nonparametric and parametric population PK modeling, multiple-model dosage design, Monte Carlo simulations, and Bayesian adaptive feedback control are the methods of choice to optimize therapy. Population PK can estimate between patient variability and account for potentially increased clearances and large volumes of distribution in critically ill patients. Once patient- specific PK data become available, target concentration intervention and adaptive feedback control algorithms can most precisely achieve target goals such as clinical cure of an infection or resistance prevention in stable and unstable patients with rapidly changing PK parameters. Many bacterial resistance mechanisms cause PK/PD targets for resistance prevention to be usually several-fold higher than targets for near-maximal killing. In vitro infection models such as the hollow fiber and one-compartment infection models allow one to study antibiotic-induced bacterial killing and emergence of resistance of mono- and combination therapies over clinically relevant treatment durations. Mechanism-based (and empirical) PK/PD modeling can incorporate effects of the immune system and allow one to design innovative dosage regimens and prospective validation studies. Mechanism-based modeling holds great promise to optimize mono- and combination therapy of anti-infectives and drugs from other therapeutic areas for critically ill patients.     

Population pharmacokinetic modeling of flurbiprofen

The paper is to report the establishment of a population pharmacokinetic model for flurbiprofen (FP), an active metabolite of flurbiprofen axetil (FA). 246 FP serum concentration and clinical data were perspectively collected from 23 general anaesthesia patients receiving FA intravenously before operation in Dentofacial Surgery and Otorhinolaryngology Department of the First Affiliated Hospital of Fujian Medical University. Population pharmacokinetic data analysis was performed using NONMEM software. The measure of Bootstrap was applied for internal validation, while Visual Predictive check was adopted for external validation. The data of FP correspond with two-compartment model. The body weight (WT) had conspicuous effect on clearance and volume of central compartment, while sex, age and daily dose of administration had no marked effect on pharmacokinetic parameter of FP. The basic model was described as follows: CL (L x h(-1)) = 1.28x EXP(ETA(1)), V1 (L) = 5.03x EXP(ETA(2)), Q (L x h(-1)) = 8.5 x EXP(ETA(3)), V2 (L) = 4.39 x EXP(ETA(4)). The final model was described as follows: CL (L x h(-1)) = 1.32 x (WT/60) x EXP(ETA(1)), V1 (L) = 5.23 x (WT/60) x EXP(ETA(2)), Q (L x h(-1)) = 8.45 x EXP(ETA(3)), V2 (L) = 4.37 x EXP(ETA(4)). The population typical value of CL, V1, Q and V2 were: 1.32 L x h(-1), 5.23 L, 8.45 L x h(-1) and 4.37 L, respectively. Bootstrap and visual predictive check show that the final model of FP is stable, effective and predictable. A novel population pharmacokinetic model is developed to estimate the individual pharmacokinetic parameter for patients intravenous injecting FA in terms of patients' characteristics and dosing history, and to design a prior dosage regimen.     

左旋多巴中国人群药动学模型的建立

目的:建立中国人群左旋多巴/苄丝肼复合制剂中左旋多巴的群体药动学模型。方法:前瞻性收集服用多巴丝肼片的帕金森病(PD)门诊患者稳态谷浓度97例102个血样和健康志愿者13例153个密集血样,高效液相色谱-电化学(HPLC-ECD)法测定左旋多巴(LD)血药浓度。应用NONMEM软件进行群体药动学数据分析,Bootstrap重复抽样用于模型的内部验证。另收集20例PD患者22个血样点作为验证组进行模型外部验证,计算最简模型和最终模型对验证组的平均预测误差(MPE)和平均绝对误差(MAE)对模型进行外部验证。结果:数据采用一房室模型拟合,年龄(AGE)对LD清除率有显著影响,性别(SEX)、体质量(WT)、给药剂量(TAMT)、合并用药不影响LD的药动学参数。LD的基础模型为:CL(CL/F)(L.h-1)=18.2×EXP[ETA(1)],V(V/F)(L)=48.4,ka(h-1)=2.13×EXP[ETA(2)];最终模型为:CL(CL/F)(L.h-1)=17.9×(55/AGE)0.59×(EXP[ETA(1)],V(V/F)(L)=47.5,ka(h-1)=2.14×EXP[ETA(2)]。CL、V、ka的群体典型值分别为17.9 L.h-1、47.5 L、2.14 h-1。Bootstrap重复抽样显示所建立的最终模型稳定、可靠,最终模型对验证组的MPE和MAE较最简模型有显著改善,显示模型有效,且有一定代表性。结论:根据患者的生理用药资料,结合上述模型,可估算个体药动学参数,为临床个体化给药提供参考。     

Population pharmacokinetic and pharmacodynamic modelling of artemisinin and mefloquine enantiomers in patients with falciparum malaria

OBJECTIVES: The aims of this study were to investigate whether artemisinin influences the pharmacokinetics of mefloquine enantiomers or vice versa and to model the antiparasitic effect of these drugs alone and in combination in Plasmodium falciparum malaria patients. METHODS: Forty-two male and female patients were randomised to treatment with either oral artemisinin 500 mg daily for 3 days followed by oral mefloquine 750 mg on day 4, oral artemisinin 500 mg daily for 3 days plus oral mefloquine 750 mg on day 1 or a single 750-mg oral dose of mefloquine. The data was modelled using NONMEM. RESULTS: All patients were successfully treated regardless of treatment. The fastest parasite clearance rates were observed in patients receiving artemisinin together with mefloquine on the first day of treatment. A pharmacodynamic model based on the life cycle of P. falciparum successfully described the efficacy of artemisinin, mefloquine and the combination. The time artemisinin concentration stays above a minimum inhibitory concentration was estimated to 2.97 h (relative standard error 4.7 h). The two mefloquine enantiomers exhibited different pharmacokinetics, with an oral clearance of 3.51 (7.9) l/h and 0.602 (6.9) l/h for RS-mefloquine and SR-mefloquine, respectively. In patients receiving only artemisinin the first 3 days, artemisinin oral clearance was 6.9-fold higher the last day of treatment compared with the first day. There was no difference in the pharmacokinetics of mefloquine enantiomers when mefloquine was given alone, in combination with artemisinin or after a 3-day regimen of artemisinin. There was a tendency towards, although non-significant, higher artemisinin concentrations when artemisinin was given together with mefloquine compared with when given alone. CONCLUSIONS: No significant pharmacokinetic interactions were observed after co-administration of artemisinin and mefloquine. The P. falciparum malaria pharmacodynamic model successfully described the antimalarial effect of artemisinin, mefloquine and a combination of the two drugs.     

Population pharmacokinetic and pharmacodynamic modeling of MNRP1685A in cynomolgus monkeys using two-target quasi-steady-state (QSS) model

MNRP1685A (anti-NRP1) is a fully human IgG1 monoclonal antibody against neuropilin-1 (NRP1), a protein necessary for blood vessel maturation. MNRP1685A binds to free membrane-bound NRP1 (mNRP1) and circulating NRP1 (cNRP1). Total cNRP1 increased in a dose-dependent manner following anti-NRP1 administration in mice, rats, and monkeys. The purpose of this study is to develop a mechanism-based model to simultaneously describe the kinetics of both unbound drug (MNRP1685A) and total cNRP1 in cynomolgus monkeys. Pharmacokinetic (PK) and pharmacodynamic (PD) profiles after single- or multiple-dose administrations were well described by the two-target quasi-steady-state (QSS) model. The estimated nonspecific clearance was 3.26 mL/day/kg and central compartment volume was 38.2 mL/kg. The maximum elimination rate for mNRP1-mediated disposition was 98.8 nM/day. The synthesis rate (3.8 nM/day), degradation rate constant (1.53 day⁻¹), and complex elimination rate constant (0.260 day⁻¹) for cNRP1 were also derived from the model. QSS constants were 6.94 nM for mNRP1 and 2.8 nM for cNRP1. The results suggest that cNRP1 has minimal effect on MNRP1685A PK while mNRP1 plays a major role in the target-mediated drug disposition. This finding is favorable as the desired pharmacological target is mNRP1, rather than cNRP1. The two-target QSS model provides mechanistic understanding of the nonlinear PK of MNRP1685A. Based on the model prediction, the free drug concentrations to maintain free mNRP1 and cNRP1 below 10% of baseline level are 10 and 20 μg/mL, respectively. This serves as a target concentration for clinical dose selection, assuming cynomolgus monkeys are predictive for humans.     

The place of pegvisomant in the management of acromegaly

Conventional treatments for acromegaly include surgery, radiotherapy, dopamine agonists and somatostatin (SMS) analogues, which effect disease control by lowering circulating growth hormone (GH). Due to variability in tumour characteristics, combinations of these treatment modalities leave a significant number of patients with sub-optimal serum GH and insulin-like growth factor-I (IGF-I) levels, which have been linked to increased morbidity and mortality. The GH receptor antagonist pegvisomant is a genetically engineered analogue of GH that prevents functional dimerisation of the growth hormone receptor (GHR); a process that is critical to GH action at the cellular level. A crucial amino acid substitution at Gly¹²⁰ to Arg¹²⁰ within the third alpha helix of the antagonist prevents functional GHR dimerisation. Pegvisomant represents a novel treatment for acromegaly as, unlike existing treatment modalities, the effectiveness of pegvisomant is independent of pituitary tumour characteristics. Initial clinical studies in patients with active acromegaly have demonstrated serum IGF-I normalisation in over 90% of patients receiving 20 mg per day, such that, in terms of serum IGF-I normalisation, pegvisomant now represents the most effective medical treatment for acromegaly. Although there are limited long-term data on the use of pegvisomant and questions regarding pituitary tumour growth and altered liver function remain, this therapy offers the prospect of serum IGF-I normalisation in the vast majority of patients with active acromegaly.     

Ocular pharmacokinetic/pharmacodynamic modeling for multiple anti-glaucoma drugs

We have constructed a new ocular pharmacokinetic pharmacodynamic (PK/PD) model for anti-glaucoma drugs to describe ocular hypotensive effects on intraocular pressure (IOP) after instillation of a combination of an alpha(1)-adrenergic antagonist, bunazosin, and a beta-adrenergic antagonist, timolol, into rabbits. This model was constructed by the combination of two ocular PK/PD models for bunazosin and timolol by including aqueous humor dynamics based on both action mechanisms. We also verified the reliability of this model by confirming the drug concentrations in aqueous humor and ocular hypotensive effects after instillation of the drug combination. The aqueous humor concentrations of timolol and bunazosin were determined by an HPLC, and ocular hypotensive effect-time profiles were measured using a telemetry system, which was able to record automatically detailed effects. The combined model could simulate the aqueous humor concentrations of both drugs and the additive IOP-lowering effect after instillation of the combination using the MULTI (RUNGE) program and PK/PD parameters which were obtained from ocular hypotensive effects after instillation of bunazosin alone or timolol alone. The theoretical concentration curves of both drugs in the aqueous humor and the theoretical ocular hypotensive effect curves almost agreed with both the observed concentrations and ocular hypotensive effects after instillation of the drug combination. These results indicate the reliability and usefulness of PK/PD modeling considering aqueous humor dynamics to predict IOP in multidrug therapy. This is the first study to develop a PK/PD model for multidrug therapy for the eye.     

丙戊酸清除率的群体药动学模型的建立

目的:建立中国癫痫病患者丙戊酸清除率的群体药动学模型.方法:前瞻性收集上海、北京两地4所医院服用丙戊酸的350名患者的稳态血药浓度(n=435).数据分析采用非线性混合效应模型.结果:最终模型为:CL(L·h-1)=0.0422·Dose·BSA0.269·1.41(如果合用卡马西平,否则为1)·1.37(如果合用苯妥因,否则为1)·(0.00735·PBS 0.807)(如果合用苯巴比妥,否则为1)·(Dose/950)(如果Dose大于950 mg·m-2·d-1,否则为1)·1.21(如果BSA大于1.7 m2,否则为1)·1.24(如果年龄小于6,否则为1).上式中Dose为日剂量(mg·m-2·d-1);BSA为体表面积(m2);PBS为苯巴比妥的日剂量(mg·m-2·d-1).结论:根据患者的生理用药资料,结合上述模型,可估算其清除率,为制定给药方案提供依据.     

Combined pharmacokinetic/pharmacodynamic (PK/PD) modeling

Several compartmental and noncompartmental approaches have been successfully applied to PK/PD modeling. Although great opportunities exist to expand on these techniques, application of these methods to relate concentrations of all drugs to their effects is possible.     

Population pharmacokinetic modeling of epoetin delta in pediatric patients with chronic kidney disease

This analysis quantifies the population pharmacokinetics of subcutaneous and intravenous epoetin delta, an epoetin produced in a human cell line, in pediatric patients with chronic kidney disease and estimates the effects of covariate factors on epoetin delta and epoetin alfa pharmacokinetic parameters. Erythropoietin serum concentration data, taken from a phase III study conducted in 60 patients aged 1 to 17 years, were best described by a 1-compartment model with first-order absorption and elimination. The typical point estimates were clearance (0.268 L/h), central volume of distribution (1.03 L), absorption rate constant (0.0554 h(-1)), and bioavailability (0.708) for a 35-kg male < or = 10 years who was predialysis and on subcutaneous epoetin delta treatment. Erythropoietin pharmacokinetic parameters were similar in pediatric patients as compared with adults when scaled by weight. The subcutaneous administration of epoetin alfa exhibited lower systemic bioavailability than subcutaneous administration of epoetin delta.     

Population pharmacokinetic and pharmacodynamic modelling of the antimalarial chemotherapy chlorproguanil/dapsone

AIMS: To determine the population pharmacokinetics of chlorproguanil, dapsone and the active metabolite of chlorproguanil, chlorcycloguanil; and to estimate the duration of parasitocidal activity for chlorpoguanil/dapsone against Plasmodium falciparum isolates of varying sensitivity. METHODS: Rich and sparse pharmacokinetic data were collected prospectively from: healthy volunteers (n=48) and adults (n=65) and children (n=68) suffering from P. falciparum malaria. All subjects received 2.0 mg kg-1 of chlorproguanil and 2.5 mg kg-1 of dapsone. RESULTS: The population pharmacokinetic parameter estimates for chlorproguanil were ka=00.09 h-1 (intersubject variability was 44%), CL/F=51.53 l h-1 (57%), CLD/F=54.67 l h-1, V1/F=234.40 l (50%) and V2/F=1612.75 l; for dapsone were ka=00.93 h-1, CL/F=1.99 l h-1 (72%) and V/F=76.96 l (48%); and for chlorcycloguanil were CLm/Fm=3.72 l h-1 kg-1 (67%) and Vm/Fm=12.76 l kg-1 (64%). For dapsone, CL/F and V/F were both significantly positively correlated with body weight. For a 10-kg child, the mean duration of parasitocidal activity for chlorproguanil/dapsone against the three most susceptible P. falciparum strains was 4.5 days [5th and 95th percentiles 2.4, 7.3] for W282; 5.9 days (3.6, 9.7) for ItG2F6; and 6.1 days (3.7, 10.1) for K39. For an isolate with the ile-164-leu mutation, V1/S, activity ranged from 0.8 days (0.0, 3.3) for a 10-kg child to 1.8 days (0.0, 4.0) for a 60-kg adult. CONCLUSIONS: Plasmodium falciparum malaria has no effect on the pharmacokinetic parameters for chlorproguanil, dapsone or chlorcycloguanil. Chlorproguanil/dapsone will probably prove to be ineffective against parasite strains with the mutation ile-164-leu, were these to become prevalent in Africa.