Population Pharmacodynamic Model Research Articles

Relationship of tafamidis binding site occupancy, transthyretin stabilization, and disease modification in tafamidis treated transthyretin amyloid cardiomyopathy patients

Abstract Introduction Tafamidis inhibits progression of transthyretin (TTR) amyloid cardiomyopathy (ATTR-CM) by binding TTR tetramer and inhibiting its dissociation to monomers which can denature and result in amyloid fibril formation and deposition in cardiac tissue. Purpose While the phase 3 ATTR-ACT clinical data clearly demonstrated efficacy, quantification of TTR stabilization in patients and the degree to which the approved dose captures the full potential of the mechanism has yet to be assessed. Methods Isothermal titration calorimetry and subunit fraction exchange were used to determine affinities to the two binding sites on TTR. These values were used to develop a model of tafamidis binding to TTR in plasma that was applied to individual patient data to calculate TTR binding site occupancy and the change in total TTR levels induced by TTR stabilization. Population pharmacodynamic (PD) models were developed for three measures of disease progression, plasma NT-proBNP levels, Kansas City Cardiomyopathy Questionnaire – Overall Score (KCCQ-OS), and six-minute walk test (6MWT) distance, to evaluate patient response with the degree of TTR occupancy. Results In vitro binding data of wild-type TTR confirmed tafamidis binds two sites of TTR with negative cooperativity and provided precise estimates of the binding affinity to TTR and albumin. Modeling individual patient data of tafamidis exposure and increased TTR plasma levels using the in vitro derived binding affinity values confirms single site binding is consistent with complete tetramer stabilization. Patients given 80 mg tafamidis meglumine, the clinically approved dose for ATTR-CM, had a 92% reduction in unbound, unstabilized TTR, which correlated with a 53% decrease in NT-proBNP elevation, a 56% decrease in KCCQ-OS worsening and a 49% reduced decline in the 6MWT. For 100% receptor occupancy and stabilization, the expected effects on these measures are 58%, 61%, and 54% for NT-proBNP, KCCQ-OS, and 6MWT, respectively. Conclusions These results demonstrate a quantitative relationship between TTR stabilization, the mechanism of action of tafamidis, and accepted laboratory and patient-based outcomes in ATTR-CM. These results also support the value of TTR stabilization as a clinically beneficial treatment option which maintains the protein in its physiologically active form within the body. Funding Acknowledgement Type of funding sources: Private company. Main funding source(s): This study was sponsored by Pfizer.

Optimising the Nadroparin Dose for Thromboprophylaxis During Hemodialysis by Developing a Population Pharmacodynamic Model Using Anti-Xa Levels.

The optimal nadroparin dose in patients undergoing hemodialysis is difficult to determine in clinical practice. Anti-Xa levels ≥0.4 IU/mL and <2.0 IU/mL are suggested to prevent thrombus formation within the extracorporeal circuit whilst minimizing bleeding risk. We aimed to characterize the variability in the association between dose and anti-Xa levels, identify patient and dialysis characteristics that explained this variability, and optimize nadroparin dosing based on the identified characteristics. Anti-Xa samples were collected in patients who received intravenous nadroparin as thromboprophylaxis during routine dialysis sessions. A population pharmacodynamic model was developed using non-linear mixed-effects modelling. The percentage of patients ≥0.4 IU/mL (efficacy) and <2.0 IU/mL (safety) was simulated for different doses, patient and dialysis characteristics. Patients (n =137) were predominantly receiving standard hemodialysis (84.7% vs. hemodiafiltration 15.3%) and had a mean bodyweight of 76.3 kg (±16.9). Lean body mass (LBM), mode of dialysis, and dialyzer partially explained between-subject variability in anti-Xa levels. Patients on hemodiafiltration and those receiving hemodialysis with a high LBM (≥80 kg) had a low probability (<29%) of anti-Xa levels ≥0.4 IU/mL during the entire dialysis session. All patients, except hemodialysis patients with a low LBM (<50 kg), had a high probability (>70%) of peak anti-Xa levels <2.0 IU/mL. Mainly patients receiving hemodiafiltration and those receiving hemodialysis with a high LBM can benefit from a higher nadroparin dose than currently used in clinical practice, while having anti-Xa levels <2.0 IU/mL.

Population Pharmacodynamic Analyses of Human Anti-Rabies Virus Monoclonal Antibody (Ormutivimab) in Healthy Adult Subjects.

Ormutivimab is the first recombinant human anti-rabies monoclonal antibody (rhRIG) approved for clinical application in China. In this study, a population pharmacodynamic (PPD) model was established to compare the neutralizing antibody activities of Ormutivimab and human rabies immunoglobulin (HRIG), alone or combined with human rabies vaccine (Vero), in a phase II clinical trial, and to recommend a target dose for the phase III trial. The model was verified to fit the PPD data well. The stability of the model was verified by the bootstrap method. The level of neutralizing antibodies in vivo increased rapidly after administration of Ormutivimab or HRIG. Neutralizing antibodies with a strong activity were produced at 7 days (Ormutivimab + vaccine) or 10 days (HRIG + vaccine) after induction by the vaccine in vivo. Compared to that induced by HRIG + vaccine, the level of the neutralizing antibodies induced by Ormutivimab + vaccine peaked higher and faster. The levels of neutralizing antibodies induced by Ormutivimab + vaccine and HRIG + vaccine were similar within 21 days after administration. According to these results and the safety data, 20 IU·kg−1 was recommended as the target dose in the confirmatory study of Ormutivimab. Registration: ClinicalTrials.gov #NCT02559921.

Strategies for developing Alzheimer’s disease treatments: application of population pharmacokinetic and pharmacodynamic models

Strategies for developing Alzheimer’s disease treatments: application of population pharmacokinetic and pharmacodynamic models

Rivaroxaban population pharmacokinetic and pharmacodynamic modeling in Iranian patients.

Although predictable pharmacokinetic and pharmacodynamic of rivaroxaban allow fixed dosing regimens without routine coagulation monitoring, there is still the necessity to monitor and predict the effects of rivaroxaban in specific conditions and different populations. The current study was designed and conducted to analyze the rivaroxaban population pharmacokinetics in Iranian patients and establish a pharmacokinetic/pharmacodynamic model to predict the relationship between rivaroxaban concentration and its anticoagulant activity. A sequential nonlinear mixed effect pharmacokinetic/pharmacodynamic modeling method was used to establish the relation between rivaroxaban concentration and anti-factor Xa activity, prothrombin time, and activated partial thromboplastin time (aPTT) as pharmacodynamic biomarkers in a population of sixty-nine Iranian patients under treatment with oral rivaroxaban. Rivaroxaban plasma concentration was quantified by a validated high-performance liquid chromatography-tandem mass spectrometry. The typical population values (inter-individual variability%) of the oral volume of distribution and clearance for a one-compartment model were 61.2 L (21%) and 3.68 L·h-1 (61%), respectively. Creatinine clearance and Child-Turcotte-Pugh score were found to affect the clearance. A direct link linear structural model best fitted the data for both prothrombin time and aPTT. The baseline estimates of aPTT and prothrombin time in the population were 35.0 (15%) and 12.6 (2%) seconds, respectively. The slope of the relationship between apTT, prothrombin time, and rivaroxaban concentration was 0.033 (28%) and 0.018 (54%) s·ml·ng-1 , respectively. The selected model for anti-factor Xa activity consisted of a direct link inhibitory Emax model with Hill coefficient. The maximum level of inhibition (Emax ) was 4 IU·ml-1 . The concentration of rivaroxaban producing 50% of the maximum inhibitory effect (EC50 ) was 180 (24%) ng·ml-1 , and Hill coefficient (γ) was 1.44 (108%). No covariates showed a statistically significant effect on PT and activated partial thromboplastin time prolonging properties and anti-factor Xa activity. Our results confirmed that pharmacokinetic/pharmacodynamic models similar to those of the other studies describe the relationship between the rivaroxaban concentration and its anticoagulant effect in Iranian patients. However, considerable differences were observed in the parameters of the pharmacodynamics-pharmacokinetic models with the results of other reports that can explain the unpredictable effects of rivaroxaban in some patients.

MO364: Population Pharmacodynamic Dose-Response Analysis of Serum Potassium Following Dosing With Sodium Zirconium Cyclosilicate

Abstract BACKGROUND AND AIMS Sodium zirconium cyclosilicate (SZC, LOKELMA®) is an orally administered suspension for treating hyperkalemia that permits selective, monovalent cation exchange and has a high exchange capacity for potassium cations in the gut leading to a net reduction in systemic potassium. The goals of the analysis were to describe the population pharmacodynamic response of serum potassium concentration (S-K) in hyperkalemic subjects treated with SZC; specifically, the longitudinal dose response of SZC for S-K in correction (0–48 h) and maintenance (&amp;gt;48 h) phases of treatment, and the differences in S-K response among populations identified by their intrinsic and extrinsic factors. METHOD A population pharmacodynamic model was developed using serum potassium, dose timing/amount, and patient demographic and medical history data pooled from 7 phase 2 and 3 clinical studies consisting of 2369 subjects with hyperkalemia (screening S K &amp;gt;5 mmol/L) with a total of 25 764 S-K observations. A virtual pharmacokinetic-pharmacodynamic model (K-PD) model was used to mimic a drug exposure without any measurements (none were available, drug was not absorbed). S-K over time was described by an indirect response model where potassium enters/leaves a compartment corresponding to serum potassium measurements. The impact of SZC exposure was modeled as an inhibition of the potassium production rate, consistent with the SZC mechanism of action to bind potassium and lower its absorption and appearance in circulation. RESULTS Diagnostics indicated good agreement between the observed data and the model, for mean S-K and change from baseline, and overall qualified the model as appropriate for predicting S-K changes following SZC dosing. The sigmoid exposure response on Kin was characterized by a high value for the EC50 (32.8 g) and a hill coefficient of 1.36. Not surprisingly, within the range of doses studied (0–15 g), the dose response appeared nearly linear. The placebo adjusted S-K change from baseline for QD maintenance treatment was 0.37 (95% CI –0.42 to –0.32) for 5 g, –0.75 (95% CI –0.82 to –0.65) for 10 g, and –1.06 (95% CI –1.17 to –0.89) for 15 g. Greater treatment response was associated with high S-K at baseline, more advanced age, lower body weight and lower eGFR. Furthermore, greater treatment response was also predicted for Black or African American and Asian patients, compared with Caucasian individuals, but they remained unworthy of dose adjustment. Lastly, dose-response simulations show no clinically meaningful difference between maintenance regimens of 5 g QOD and 2.5 g QD in the studied patient groups. CONCLUSION The population K-PD model of SZC adequately described the changes in serum potassium concentration during correction and maintenance phase dosing in the seven clinical studies. Greater treatment response was associated with high S-K at baseline, more advanced age, lower body weight, lower eGFR, Black or African American race and Asian race; however, the impact of each covariate was modest. The overall results of this analysis support the conclusion that no dose adjustment is needed for SZC based on any of the tested covariates.

A semimechanistic population pharmacokinetic and pharmacodynamic model incorporating autoinduction for the dose justification of TAS-114.

TAS‐114 is a dual deoxyuridine triphosphatase (dUTPase) and dihydropyrimidine dehydrogenase (DPD) inhibitor expected to widen the therapeutic index of capecitabine. Its maximum tolerated dose (MTD) was determined from a safety perspective in a combination study with capecitabine; however, its inhibitory effects on DPD activity were not assessed in the study. The dose justification to select its MTD as the recommended dose in terms of DPD inhibition has been required, but the autoinduction profile of TAS‐114 made it difficult. To this end, an approach using a population pharmacokinetic (PPK)/pharmacodynamic (PD) model incorporating autoinduction was planned; however, the utility of this approach in the dose justification has not been reported. Thus, the aim of this study was to demonstrate the utility of a PPK/PD model incorporating autoinduction in the dose justification via a case study of TAS‐114. Plasma concentrations of TAS‐114 from 185 subjects and those of the endogenous DPD substrate uracil from 24 subjects were used. A two‐compartment model with first‐order absorption with lag time and an enzyme turnover model were selected for the pharmacokinetic (PK) model. Moreover, an indirect response model was selected for the PD model to capture the changes in plasma uracil concentrations. Model‐based simulations provided the dose justification that DPD inhibition by TAS‐114 reached a plateau level at the MTD, whereas exposures of TAS‐114 increased dose dependently. Thus, the utility of a PPK/PD model incorporating autoinduction in the dose justification was demonstrated via this case study of TAS‐114.

Semimechanistic models to relate noxious stimulation, movement, and pupillary dilation responses in the presence of opioids.

Intraoperative targeting of the analgesic effect still lacks an optimal solution. Opioids are currently the main drug used to achieve antinociception, and although underdosing can lead to an increased stress response, overdose can also lead to undesirable adverse effects. To better understand how to achieve the optimal analgesic effect of opioids, we studied the influence of remifentanil on the pupillary reflex dilation (PRD) and its relationship with the reflex movement response to a standardized noxious stimulus. The main objective was to generate population pharmacodynamic models relating remifentanil predicted concentrations to movement and to pupillary dilation during general anesthesia. A total of 78 patients undergoing gynecological surgery under general anesthesia were recruited for the study. PRD and movement response to a tetanic stimulus were measured multiple times before and after surgery. We used nonlinear mixed effects modeling to generate a population pharmacodynamic model to describe both the time profiles of PRD and movement responses to noxious stimulation. Our model demonstrated that movement and PRD are equally depressed by remifentanil. Using the developed model, we changed the intensity of stimulation and simulated remifentanil predicted concentrations maximizing the probability of absence of movement response. An estimated effect site concentration of 2 ng/ml of remifentanil was found to inhibit movement to a tetanic stimulation with a probability of 81%.

A stochastic population pharmacodynamic model of QAP14 in the treatment of lung metastases of 4T1 breast cancer

A stochastic population pharmacodynamic model of QAP14 in the treatment of lung metastases of 4T1 breast cancer

Maxacalcitol Pharmacokinetic-Pharmacodynamic Modeling and Simulation for Secondary Hyperparathyroidism in Patients Receiving Maintenance Hemodialysis.

Maxacalcitol was approved in Taiwan in 2018 as the first active vitamin D3 injection for secondary hyperparathyroidism (SHPT) in patients on maintenance hemodialysis. However, no data from any clinical study with maxacalcitol in Taiwanese patients is available. This analysis aimed to evaluate the profiles of parathyroid hormone (PTH) and calcium (Ca) concentrations in Taiwanese SHPT patients on hemodialysis and maxacalcitol. We developed population pharmacokinetic (PK) and pharmacodynamic (PD) models using a modeling and simulation approach. The data for these analyses were obtained from two studies: a clinical pharmacology study in Japanese patients and an ethnic comparison study in healthy Japanese and -Taiwanese volunteers. We then conducted a simulation study with a PK-PD model comprising the PK and PD models developed here. Serum maxacalcitol concentration profile was modeled using a two-compartment model that took into consideration the distribution of concentrations below the lower limit of quantification. An ethnic difference in clearance was included in the PK model as a covariate. A PD model that used a PTH/Ca feedback loop best described the observed data. There were no significant differences in Ca or PTH concentrations between Taiwanese and Japanese based on the simulation results from our PK-PD model, even though maxacalcitol exposure was approximately 40% higher in Taiwanese than in Japanese. On the basis of these population PK and PD analyses and the clinical study conducted in Japan, there is no clinically relevant difference between Taiwanese and Japanese in terms of serum Ca or PTH levels.

Hospital Pharmacometrics for Optimal Individual Administration of Antimicrobial Agents for Anti-methicillin-resistant Staphylococcus aureus Infected Patients.

Therapeutic drug monitoring and target concentration intervention based on population pharmacokinetic and pharmacodynamic models has been strongly recommended for anti-methicillin-resistant Staphylococcus aureus (MRSA) agents in order to provide appropriate antimicrobial chemotherapy to each individual patient, and pharmacokinetic and pharmacodynamic analyses in hospitalized patients have been actively conducted, as evidenced with vancomycin. Teicoplanin, daptomycin, and linezolid have been the most studied antibiotics, using population pharmacokinetics of patients with MRSA. Infections caused by MRSA have higher severity and fatality rates than other antimicrobial-susceptible infections. Therefore, many medical facilities have been implementing infection control programs based on antimicrobial stewardship to prevent nosocomial infections and drug-resistant strains. Studies detailing pharmacometrics for these antibiotics have been reported to elucidate the pharmacokinetic and pharmacodynamic properties, to determine significant factors influencing variabilities between individuals, and to develop target concentration interventions and dosing regimens for adults, the elderly, patients with renal insufficiency including those on continuous renal replacement therapies, patients with low body weight, obese patients, and pediatric patients. This review presents the details of our recent research on the optimal dosing design of antimicrobial agents for the treatment of MRSA infection based on hospital pharmacometrics. In addition, the prospect of using modeling and simulation has shown major advantages in supporting dosing regimen selection.

ModVizPop: A shiny interface for empowering teams to perform interactive pharmacokinetic/pharmacodynamic simulations.

ModVizPop is an interactive and dynamic visualization tool developed for simulating differential equation–based population pharmacokinetic (PK) and pharmacodynamic (PD) models with variability. It has a built‐in PK/PD ordinary differential equations library of models to choose from alongside the ability to plug in a user‐defined model from a local or project directory. The user interface consists of several key inputs for performing the simulations as well as options to visualize the plots, perform simple noncompartmental analysis, and review inputs and model code. It also provides the ability to download the underlying model, plots, simulated data, or a comprehensive report consisting of all the key inputs and outputs of the simulations. The interface includes advanced features where users can overlay external data on a simulation, set a certain simulation scenario as a reference, or carry out sensitivity analysis–based simulations. This easy‐to‐use interface can serve as a valuable tool to project teams to evaluate potential scenarios facilitating collaborative decision making in the drug discovery and development paradigm.

Therapeutic Drug Monitoring of Antimicrobial Drugs in Neonates: An Opinion Article.

Background:Neonatal infections are associated with high morbidity and mortality rates. Optimal treatment of these infections requires knowledge of neonatal pharmacology and integration of neonatal developmental pharmacokinetics (PKs) of antimicrobial drugs in the design of dosing regimens for use with different gestational and postnatal ages. Population PK and pharmacodynamic models are used to personalize the use of these drugs in these fragile patients. The final step to further minimize variability in an individual patient is therapeutic drug monitoring (TDM), where the same population PK/pharmacodynamic models are used in concert with optimally drawn blood samples to further fine-tune therapy. The purpose of this article is to describe the present status and future role of model-based precision dosing and TDM of antimicrobial drugs in neonates.Methods:PubMed was searched for clinical trials or clinical studies of TDM in neonates.Results:A total of 447 articles were retrieved, of which 19 were concerned with antimicrobial drugs. Two articles (one aminoglycoside and one vancomycin) addressed the effects of TDM in neonates. We found that, in addition to aminoglycosides and vancomycin, TDM also plays a role in beta-lactam antibiotics and antifungal drugs.Conclusions:There is a growing awareness that, in addition to aminoglycosides and vancomycin, the use of beta-lactam antibiotics, such as amoxicillin and meropenem, and other classes of antimicrobial drugs, such as antifungal drugs, may benefit from TDM. However, the added value must be shown. New analytical techniques and software development may greatly support these novel developments.

Population pharmacodynamic modeling of intramuscular and oral dexamethasone and betamethasone effects on six biomarkers with circadian complexities in Indian women

Population pharmacokinetic/pharmacodynamic (PK/PD) analysis was performed for extensive data for differing dosage forms and routes for dexamethasone (DEX) and betamethasone (BET) in 48 healthy nonpregnant Indian women in a partial and complex cross-over design. Single doses of 6 mg dexamethasone phosphate (DEX-P), betamethasone phosphate (BET-P), or 1:1 mixture of betamethasone phosphate and acetate (BET-PA) were administered orally (PO) or intramuscularly (IM) where each woman enrolled in a two-period cross-over study. Plasma concentrations collected over 96 h were described with a two-compartment model with differing PO and IM first-order absorption inputs. Overall, BET exhibited slower clearance, similar volume of distribution, faster absorption, and longer persistence than DEX with BET acetate producing extremely slow absorption but full bioavailability of BET. Six biomarkers were assessed over a 24-h baseline period with four showing circadian rhythms with complex baselines. These baselines and the strong responses seen after drug dosing were fitted with various indirect response models using the Laplace estimation methods in NONMEM 7.4. Both the PK and six biomarker responses were well-described with modest variability likely due to the homogeneous ages, weights, and ethnicities of the women. The drugs either inhibited or stimulated the influx processes with some models requiring joint inclusion of drug effects on circadian cortisol suppression. The biomarkers and order of sensitivity (lowest IC50/SC50 to highest) were: cortisol, T-helper cells, basophils, glucose, neutrophils, and T-cytotoxic cells. DEX sensitivities were generally greater than BET with corresponding mean ratios for these biomarkers of 2.86, 1.27, 1.72, 1.27, 2.69, and 1.06. Overall, the longer PK (e.g. half-life) of BET, but lesser PD activity (e.g. higher IC50), produces single-dose response profiles that appear quite similar, except for the extended effects from BET-PA. This comprehensive population modeling effort provides the first detailed comparison of the PK profiles and six biomarker responses of five commonly used dosage forms of DEX and BET in healthy women.

The power of modelling pulsatile profiles

The quantitative description of individual observations in non-linear mixed effects models over time is complicated when the studied biomarker has a pulsatile release (e.g. insulin, growth hormone, luteinizing hormone). Unfortunately, standard non-linear mixed effects population pharmacodynamic models such as turnover and precursor response models (with or without a cosinor component) are unable to quantify these complex secretion profiles over time. In this study, the statistical power of standard statistical methodology such as 6 post-dose measurements or the area under the curve from 0 to 12 h post-dose on simulated dense concentration–time profiles of growth hormone was compared to a deconvolution-analysis-informed modelling approach in different simulated scenarios. The statistical power of the deconvolution-analysis-informed approach was determined with a Monte-Carlo Mapped Power analysis. Due to the high level of intra- and inter-individual variability in growth hormone concentrations over time, regardless of the simulated effect size, only the deconvolution-analysis informed approach reached a statistical power of more than 80% with a sample size of less than 200 subjects per cohort. Furthermore, the use of this deconvolution-analysis-informed modelling approach improved the description of the observations on an individual level and enabled the quantification of a drug effect to be used for subsequent clinical trial simulations.

Pharmacodynamic analysis of hypertension caused by lenvatinib using real-world postmarketing surveillance data.

Lenvatinib is a tyrosine kinase inhibitor of the vascular endothelial growth factor receptor used against nonoperative thyroid cancer; however, hypertension is a major dose‐limiting side effect. In this study, hypertension caused by lenvatinib was described through a novel population pharmacodynamic model using postmarketing surveillance data obtained in Japan. The model consists of two maximum effect model components based on the (1) concentration of lenvatinib in plasma and (2) cumulative area under the curve of lenvatinib. In addition, antihypertensive drug of either an angiotensin‐converting enzyme inhibitor/angiotensin receptor blocker or calcium channel blocker accounted for by lowering effect on diastolic blood pressure. Based on virtual simulations, the combination of antihypertensive drug and dose adjustment of lenvatinib showed a reduction in the probability of grade greater than or equal to 3 hypertension. The present model provides useful guidance in managing hypertension during treatment with lenvatinib in the real‐world setting.

Population pharmacokinetics of mycophenolic acid in paediatric patients.

Mycophenolic acid (MPA) is widely used in paediatric kidney transplant patients and sometimes prescribed for additional indications. Population pharmacokinetic or pharmacodynamic modelling has been frequently used to characterize the fixed, random and covariate effects of MPA in adult patients. However, MPA population pharmacokinetic data in the paediatric population have not been systematically summarized. The objective of this narrative review was to provide an up-to-date critique of currently available paediatric MPA population pharmacokinetic models, with emphases on modelling techniques, pharmacological findings and clinical relevance. PubMed and EMBASE were searched from inception of database to May 2020, where a total of 11 studies have been identified representing kidney transplant (n = 4), liver transplant (n = 1), haematopoietic stem cell transplant (n = 1), idiopathic nephrotic syndrome (n = 2), systemic lupus erythematosus (n = 2), and a combined population consisted of kidney, liver and haematopoietic stem cell transplant patients (n = 1). Critical analyses were provided in the context of MPA absorption, distribution, metabolism, excretion and bioavailability in this paediatric database. Comparisons to adult patients were also provided. With respect to clinical utility, Bayesian estimation models (n = 6) with acceptable accuracy and precision for MPA exposure determination have also been identified and systematically evaluated. Overall, our analyses have identified unique features of MPA clinical pharmacology in the paediatric population, while recognizing several gaps that still warrant further investigations. This review can be used by pharmacologists and clinicians for improving MPA pharmacokinetic-pharmacodynamic modelling and patient care.

Population pharmacokinetic and pharmacodynamic modeling of capecitabine and its metabolites in breast cancer patients.

The present study was performed to examine relationships between systemic exposure of capecitabine metabolites (5-FU, 5'-DFCR and 5'-DFUR) and toxicity or clinical response in patients with metastatic breast cancer. A population pharmacokinetic model for capecitabine and its three metabolites was built. Typical parameter values, characteristics of random distributions, associated with parameters, and covariates impact were estimated. Area under the curve (AUC) were computed for 5-FU and compared with grades of toxicity. Pharmacokinetic modeling was based on data collected on the first treatment cycle. Toxicity was assessed on the two first treatment cycles. The study was conducted in 43 patients. The population pharmacokinetic model (a one-compartment model per compound) was able to capture the very complex absorption process of capecitabine. Statistically significant covariates were cytidine deaminase, alkaline phosphatase and dihydrouracilemia (UH2)/uracilemia (U) ratio. UH2/U ratio was the most significant covariate on 5-FU elimination and CDA on the transformation of 5'-DFCR in 5'-DFUR. A trend was observed between 5-FU AUC and thrombopenia toxicity grades, but not with other toxicities. Best clinical response was not linked to systemic exposure of capecitabine metabolites. In our study, we propose a model able to describe, meanwhile, and its main metabolites, with a complex absorption process and inclusion of enzyme activity covariates such as CDA and UH2/U ratio. Trial registration Eudract 2008-004136-20, 2008/11/26.

A pharmacodynamic model of tidal volume and inspiratory sevoflurane concentration in children during spontaneous breathing.

High concentrations of sevoflurane causes respiratory depression, mainly due to the decrease in tidal volume (TV) during spontaneous ventilation. The purpose of this study was to identify clinical variables that affect the relationship between TV and sevoflurane concentration, and to establish a population pharmacodynamic modelling approach to TV and sevoflurane concentration in children. A prospective observational study involving 48 patients (≤ 6 years of age) scheduled to undergo general anesthesia using laryngeal mask airway was performed. When the inspiratory sevoflurane concentration reached 2 vol%, the vaporizer was increased to 4 vol% for 5 min, then sevoflurane was decreased to 2 vol% for 5 min. During the study period, TV, end-tidal carbon dioxide, and sevoflurane concentration were recorded every 30 s. Pharmacodynamic analysis using a sigmoid Emax model was performed to assess the TV-sevoflurane concentration relationship. To collapse hysteresis of the pharmacokinetic and pharmacodynamic relationship, the semicompartmental model was applied which does not require a structural model for equilibration delay causing the hysteresis. TV decreased with increasing inspiratory sevoflurane concentrations. Hysteresis between the TV and sevoflurane concentration was observed and was accounted for when the model was developed. Initial TV and maximal reduction in TV were related to body weight. The γ (a steepness of the concentration-response relation curve) was 8.78 and the keo, (a first-order rate constant determining the equilibrium between the end-tidal sevoflurane concentration and effect site sevoflurane concentration) was 2.27 min-1. Changes in TV were correlated with sevoflurane concentration with spontaneous breathing during sevoflurane anesthesia. The initial and maximal TV were related to body weight, in a pediatric population.

Population Pharmacodynamic Modeling Using the Sigmoid Emax Model: Influence of Inter-individual Variability on the Steepness of the Concentration\u2013Effect Relationship. a Simulation Study

The relationship between the concentration of a drug and its pharmacological effect is often described by empirical mathematical models. We investigated the relationship between the steepness of the concentration–effect relationship and inter-individual variability (IIV) of the parameters of the sigmoid Emax model, using the similarity between the sigmoid Emax model and the cumulative log-normal distribution. In addition, it is investigated whether IIV in the model parameters can be estimated accurately by population modeling. Multiple data sets, consisting of 40 individuals with 4 binary observations in each individual, were simulated with varying values for the model parameters and their IIV. The data sets were analyzed using Excel Solver and NONMEM. An empirical equation (Eq. (11)) was derived describing the steepness of the population-predicted concentration–effect profile (γ*) as a function of γ and IIV in C50 and γ, and was validated for both binary and continuous data. The tested study design is not suited to estimate the IIV in C50 and γ with reasonable precision. Using a naive pooling procedure, the population estimates γ* are significantly lower than the value of γ used for simulation. The steepness of the population-predicted concentration–effect relationship (γ*) is less than that of the individuals (γ). Using γ*, the population-predicted drug effect represents the drug effect, for binary data the probability of drug effect, at a given concentration for an arbitrary individual.