Population total and unbound pharmacokinetics and pharmacodynamics of ciprofol and M4 in subjects with various renal functions.

AIMS: The aim of this study was to develop a population pharmacokinetic (PK) model to simultaneously describe both total and unbound concentrations of ciprofol and its major glucuronide metabolite, M4, and to link it to the population pharmacodynamics (PD) model in subjects with various renal functions. METHODS: A total of 401 and 459 pairs of total and unbound plasma concentrations of ciprofol and M4, respectively, as well as 2190 bispectral index (BIS) data from 24 Chinese subjects with various renal functions were available. Covariates that may potentially contribute to the PK and PD variability of ciprofol were screened using a stepwise procedure. The optimal ciprofol induction dosing regimen was determined by model-based simulations. RESULTS: The PK of unbound ciprofol could best be described by a three-compartment model, while a two-compartment model could adequately describe unbound M4 PK. The concentrations of total and unbound ciprofol and M4 were linked using a linear protein binding model. The relationship between plasma concentrations of ciprofol and BIS data was best described by an inhibitory sigmoidal Emax model with a two-compartment biophase distribution compartment. Hemoglobin was the identified covariate determining the central compartment clearance of ciprofol; uric acid was a covariate affecting the central compartment clearance of M4 and protein binding rate, kB . The included covariates had no effect on the PD of ciprofol. Simulation results indicated that the label-recommended dose regimen was adequate for anaesthesia induction. CONCLUSIONS: The developed model fully characterized the population PK and PD profiles of ciprofol. No dose adjustment is required in patients with mild and moderate renal impairment.

Population pharmacokinetics and pharmacodynamics of Tranexamic acid in women undergoing caesarean delivery.

AIMS: The population pharmacokinetics (PK) and pharmacodynamics (PD) of tranexamic acid (TXA) have not been studied to prevent postpartum haemorrhage (PPH) in pregnant women. It is unclear which TXA dose assures sufficient PPH prevention. This study investigated population PK/PD of TXA in pregnant women who underwent caesarean delivery to determine the optimal prophylactic doses of TXA for future studies. METHODS: We analysed concentration (PK) and maximum lysis (PD) data from 30 pregnant women scheduled for caesarean delivery who received 5, 10 or 15 mg/kg of TXA intravenously using population approach. RESULTS: TXA PK was best described by a two-compartment model with first-order elimination and the following parameters: clearance (between-subject variability) of 9.4 L/h (27.7%), central volume of 10.1 L (47.4%), intercompartmental clearance of 22.4 L/h (66.7%), peripheral volume of 14.0 L (13.1%) and additive error of 1.4 mg/L. The relationship between TXA concentration and maximum lysis was characterized by a sigmoid Emax model with baseline lysis of 97%, maximum inhibition of 89%, IC50 of 6.0 mg/L (65.3%), hill factor of 8.5 (86.3%) and additive error of 7.3%. Simulations demonstrated that 500 and 650 mg of TXA maintained therapeutic targets for 30 minutes and 1 hour, respectively, in 90% of patients. CONCLUSION: This is the first population PK and PD study of TXA in pregnant women undergoing caesarean delivery. Our analysis suggests that a 650 mg dose provides adequate PPH prophylaxis up to 1 hour, which is less than the currently used 1000 mg of TXA in pregnant women.

Efficacy of DPP-4 inhibitors, GLP-1 analogues, and SGLT2 inhibitors as add-ons to metformin monotherapy in T2DM patients: a model-based meta-analysis.

AIMS: The aim of the present study was to quantitate the hypoglycaemic effects of dipeptidyl peptidase-4 inhibitors (DPP-4i), glucagon-like peptide-1 receptor agonists (GLP-1r) and sodium glucose cotransporter 2 inhibitors (SGLT2i) as add-on treatments to metformin monotherapy in patients with type 2 diabetes mellitus (T2DM) using a model-based meta-analysis (MBMA). METHODS: A systematic literature search of public databases was conducted to develop models that describe the time courses of the fasting plasma glucose (FPG)- and haemoglobin A1c (HbA1c)-lowering effects of three antidiabetic classes using NONMEM 7.3.0. RESULTS: Seventy-six publications were eligible for this study, and 873 FPG and 1086 HbA1c values were collected. We developed a physiological indirect response model that described the time courses of FPG and HbA1c and simulated reductions in these values 90 days after the initiation of add-on treatments. FPG and HbA1c reductions with once weekly exenatide, liraglutide and dulaglutide were greater than those with other drugs. Mean changes from baseline FPG and HbA1c with these drugs were as follows: exenatide (-22.5 and -16.6%), liraglutide (-22.1 and -16.3%), and dulaglutide (-19.3 and -14.3%). The hypoglycaemic effects of DPP-4i and SGLT2i were similar. CONCLUSIONS: Once weekly exenatide, liraglutide and dulaglutide provided better hypoglycaemic effects among the antidiabetic drugs analysed. Long-acting GLP-1r appears to be more useful for T2DM patients inadequately controlled with metformin monotherapy.

Pharmacokinetic-pharmacodynamic (dipeptidyl peptidase-4 inhibition) model to support dose rationale in diabetes patients, including those with renal impairment, for once-weekly administered omarigliptin.

AIMS: To characterize the population pharmacokinetics (PK) and pharmacodynamics (PD) of the once-weekly dipeptidyl peptidase-4 (DPP-4) inhibitor omarigliptin in healthy subjects and patients with type 2 diabetes mellitus, and use these models to support the dosing recommendation for patient labelling including patients with renal impairment. METHODS: PK and PD were assessed from a total of 9827 omarigliptin concentrations collected from 1387 healthy subjects and patients participating in Phase 1, 2 and 3 studies examining single- or multiple-dose weekly administration of omarigliptin at doses ranging from 0.25 to 400 mg. Population PK and PD analyses were performed using nonlinear mixed effect modelling. RESULTS: A semi-mechanistic 2-compartment model with linear unbound clearance and concentration-dependent binding of omarigliptin to the DPP-4 enzyme in both the central and peripheral compartments adequately described omarigliptin PK. Key covariates on omarigliptin PK included reduced unbound clearance with renal impairment. A direct effect sigmoid maximum inhibitory efficacy model adequately described the relationship between omarigliptin plasma concentrations and DPP-4 inhibition. These models supported the current Japan label instructions that the approved omarigliptin 25-mg once-weekly dose be halved in patients with severe renal impairment and in those with end-stage renal disease. Also, if patients missed a dose, the next dose of omarigliptin should be taken as soon as remembered up to and including the day before the next scheduled dose. No other clinically important covariates were identified. CONCLUSION: The models in the present analysis adequately described PK and PD characteristics of omarigliptin and supported the dosing and administration section of the omarigliptin label.

Concentration-response modeling of ECG data from early-phase clinical studies to assess QT prolongation risk of contezolid (MRX-I), an oxazolidinone antibacterial agent.

The effects of contezolid (MRX-I, an oxazolidinone antibacterial agent) on cardiac repolarization were evaluated retrospectively using a population modeling approach in a Phase I study incorporating single ascending dose, multiple ascending dose, and food effect assessments. Linear mixed effect models were used to assess the relationships between MRX-I plasma concentrations and QT/QTc/∆QTc (baseline-adjusted), in which different correction methods for heart rate have been included. The upper bound of the one-sided 95% confidence interval (CI) for predicted ∆∆QTc was < 10 ms (ms) at therapeutic doses of MRX-I. Model performance/suitability was determined using diagnostic evaluations, which indicated rationality of one-stage concentration-QT model, as well as C-QT model suggested by Garnett et al. The finding demonstrated that MRX-I may have no clinical effects on the QT interval. Concentration-QT model may be an alternative to conventional thorough QT studies.

A model-based approach to assess the exposure-response relationship of Lorenzo's oil in adrenoleukodystrophy.

AIMS: X-linked adrenoleukodystrophy (X-ALD) is a peroxisomal disorder, most commonly affecting boys, associated with increased very long chain fatty acids (C26:0) in all tissues, causing cerebral demyelination and adrenocortical insufficiency. Certain monounsaturated long chain fatty acids including oleic and erucic acids, known as Lorenzo's oil (LO), lower plasma C26:0 levels. The aims of this study were to characterize the effect of LO administration on plasma C26:0 concentrations and to determine whether there is an association between plasma concentrations of erucic acid or C26:0 and the likelihood of developing brain MRI abnormalities in asymptomatic boys. METHODS: Non-linear mixed effects modelling was performed on 2384 samples collected during an open label single arm trial. The subjects (n = 104) were administered LO daily at ~2-3 mg kg(-1) with a mean follow-up of 4.88 ± 2.76 years. The effect of erucic acid exposure on plasma C26:0 concentrations was characterized by an inhibitory fractional Emax model. A Weibull model was used to characterize the time-to-developing MRI abnormality. RESULTS: The population estimate for the fractional maximum reduction of C26:0 plasma concentrations was 0.76 (bootstrap 95% CI 0.73, 0.793). Our time-to-event analyses showed that every mg l(-1) increase in time-weighted average of erucic acid and C26:0 plasma concentrations was, respectively, associated with a 3.7% reduction and a 753% increase in the hazard of developing MRI abnormality. However, the results were not significant (P = 0.5344, 0.1509, respectively). CONCLUSIONS: LO administration significantly reduces the abnormally high plasma C26:0 concentrations in X-ALD patients. Further studies to evaluate the effect of LO on the likelihood of developing brain MRI abnormality are warranted.

Target attainment analysis and optimal sampling designs for population pharmacokinetic study on piperacillin/tazobactam in neonates and young infants.

OBJECTIVES: Population pharmacokinetic (popPK) analyses for piperacillin/tazobactam in neonates and infants of less than 2 months of age have been performed by our group previously. The results indicate that a dose of 44.44/5.56 mg/kg piperacillin/tazobactam every 8 or 12 h may not be enough for controlling infection in this population. In order to determine the appropriate dosing regimen and to provide a rationale for the development of dosing guidelines suitable for this population, further popPK studies of piperacillin/tazobactam would need to be conducted. The aim of the present study was to determine the appropriate dosing regimen and optimal sampling schedules in neonates and infants of less than 2 months of age. METHODS: Pharmacodynamic profiling of piperacillin using Monte Carlo simulation was performed to explore the target attainment probability of different dosing regimens for infections caused by different isolated pathogens. D-optimal designs for piperacillin and tazobactam were conducted separately, and the times that overlapped were chosen as the final sampling scheme for future popPK studies in neonates and young infants of less than 2 months of age. RESULTS: Our findings revealed that compared to the current empirical piperacillin/tazobactam dose regimen (50 mg/kg every 12 h by 5-min infusion in our hospital), the clinical outcome could be improved by increasing doses, increasing administration frequency, and prolonging intravenous infusion in neonates and infants of less than 2 months of age. The following optimal sampling windows were chosen as the final sampling scheme: 0.1-0.11, 0.26-0.29, 0.97-2.62, and 7.95-11.9 h administered every 12 h with 5-min infusion; 0.1-0.12, 0.39-0.56, 2.86-4.95, and 8.91-11.8 h administered every 12 h with 3-h infusion; 0.1-0.11, 0.22-0.29, 0.91-1.96, and 5.56-7.93 h administered every 8 h with 5-min infusion; 0.1-0.11, 0.38-0.48, 2.54-3.82, and 6.86-7.93 h administered every 8 h with 3-h infusion; 0.1-0.11, 0.25-0.28, 0.84-1.69, and 4.55-5.94 h administered every 6 h with 5-min infusion; and 0.1-0.11, 0.37-0.54, 3.13-3.72, and 5.57-5.99 h administered every 6 h with 3-h infusion. CONCLUSIONS: The dosing regimen and sampling schedules proposed in this study should be evaluated in future popPK studies of piperacillin/tazobactam in neonates and infants. To the best of our knowledge, this is the first study that combined optimal sampling design with Monte Carlo simulation for designing popPK studies of piperacillin/tazobactam.

Predicting the parasite killing effect of artemisinin combination therapy in a murine malaria model.

OBJECTIVES: To develop a mechanism-based model that describes the time course of the malaria parasite in infected mice receiving a combination therapy regimen of dihydroartemisinin and piperaquine. METHODS: Total parasite density-time data from Swiss mice inoculated with Plasmodium berghei were used for the development of population models in S-ADAPT. The mice were administered a single intraperitoneal dose of 30 mg/kg dihydroartemisinin, 10 mg/kg piperaquine phosphate or a combination of both antimalarials at 64 h post-inoculation. In a separate study, mice received multiple dihydroartemisinin doses (5 × 10 mg/kg or 30 mg/kg dihydroartemisinin followed by two 10 mg/kg doses). Parasite recrudescence after treatment was defined using a model that incorporated each erythrocytic stage of the P. berghei life cycle. RESULTS: The disposition of dihydroartemisinin and piperaquine was described by a one-compartment and two-compartment model, respectively. The estimated clearance was 1.95 L/h for dihydroartemisinin and 0.109 L/h for piperaquine. A turnover model described the parasite killing curve after single-agent dosing, with an estimated mean IC50 of 0.747 µg/L for dihydroartemisinin and 16.8 µg/L for piperaquine. In addition, the rate of parasite killing by dihydroartemisinin was almost 50-fold faster than for piperaquine. Parameters from the monotherapy models adequately described the parasite density-time curve following dihydroartemisinin/piperaquine combination therapy or multiple-dose regimens of dihydroartemisinin. CONCLUSIONS: This study has developed mechanistic models that describe the parasite-time curve after single, multiple or combination dosing of antimalarials to mice. These structural models have potential application for pre-clinical investigations to design and refine artemisinin-based combination therapy dosage regimens.

Population pharmacodynamic analysis of LDL-cholesterol lowering effects by statins and co-medications based on electronic medical records.

AIMS: HMG-CoA reductase inhibitors are available for use in low density lipoprotein-cholesterol (LDL-C) lowering therapy. The purposes of this study were to develop a population pharmacodynamic (PPD) model to describe the time course for the LDL-C lowering effects of statins and assess the efficacy of combination therapy based on electronic medical records. METHODS: Patient backgrounds, laboratory tests and prescribed drugs were collected retrospectively from electronic medical records. Patients who received atorvastatin, pitavastatin or rosuvastatin were enrolled. A physiological indirect response model was used to describe the changes observed in LDL-C concentrations. The PPD analysis was performed using nonmem 7.2.0 with the first order conditional estimation method with interaction (FOCE-INTER). RESULTS: An indirect response Imax model, based on the 2863 LDL-C concentrations of 378 patients, successfully and quantitatively described the time course for the LDL-C lowering effects of three statins. The combination of ezetimibe, a cholesterol absorption inhibitor, decreased the LDL synthesis rate (Kin ) by 10.9%. A simulation indicated that the combined treatment of ezetimibe with rosuvastatin (2.5 mg day(-1) ) led to superior clinical responses than those with high doses of rosuvastatin (5.0 mg day(-1) ) monotherapy, even in patients with higher baseline LDL-C concentrations prior to the treatment. CONCLUSIONS: A newly constructed PPD model supported previous evidence for the beneficial effects of ezetimibe combined with rosuvastatin. In addition, the established framework is expected to be applicable to other drugs without pharmacokinetic data in clinical practice.

Population Pharmacodynamic Models of Risperidone on PANSS Total Scores and Prolactin Levels in Schizophrenia.

Currently, research predominantly focuses on evaluating clinical effects at specific time points while neglecting underlying patterns within the treatment process. This study aims to analyze the dynamic alterations in PANSS total scores and prolactin levels in patients with schizophrenia treated with risperidone, along with the influencing covariates. Using data from an 8-week randomized, double-blind, multicenter clinical trial, a population pharmacodynamic model was established for the PANSS total scores of and prolactin levels in patients treated with risperidone. The base model employed was the Emax model. Covariate selection was conducted using a stepwise forward inclusion and backward elimination approach. A total of 144 patients were included in this analysis, with 807 PANSS total scores and 531 prolactin concentration values. The PANSS total scores of the patients treated with risperidone decreased over time, fitting a proportionally parameterized sigmoid Emax model with covariates including baseline score, course of the disease, gender, plasma calcium ions, and lactate dehydrogenase levels. The increase in prolactin levels conformed to the ordinary Emax model, with covariates encompassing course of the disease, gender, weight, red blood cell count, and triglyceride levels. The impacts of the baseline scores and the course of the disease on the reduction of the PANSS scores, as well as the influence of gender on the elevation of prolactin levels, each exceeded 20%. This study provides valuable quantitative data regarding PANSS total scores and prolactin levels among patients undergoing risperidone treatment across various physiological conditions.

Pharmacodynamic Age Structured Population Model For Cell Trafficking.

OBJECTIVES: Cell trafficking encompasses movement of the immune system cells (e.g., granulocytes, lymphocytes) between the blood and the extravascular tissues (e.g., lymph nodes). Corticosteroids are known to suppress cell trafficking. The age-structured cell population models introduce the transit time as a structure that allows one to quantify the distribution of times the immune cells spend in the blood and the extravascular tissues. The objective of this work is to develop an age-structured cell population model describing drug effects on cell trafficking and to implement the model in pharmacometric software to enable parameter estimation and simulations. METHODS: We adopted the well-known McKendrick age-structured population model to describe the age distributions in two cell populations: blood cells and cells in the extravascular space. The hazard of cell recirculation from the extravascular tissues was age dependent and described by the Weibull function with the shape ν and scale ß parameters. The drug effect on cell trafficking was modeled as the product of the Emax function of the drug plasma concentration and the Weibull hazard. The model was implemented in NONMEM 7.5.1. The model was applied to the basophil data in 34 healthy subjects who received a single intramuscular or oral dose of 6 mg dexamethasone (DEX). A recently published pharmacokinetic model was applied to describe DEX plasma concentration. Typical values of parameter estimates were further used to simulate the DEX effect of the basophil mean transit time in the extravascular tissues. RESULTS: Simulations of basophil time courses for varying ν demonstrated that the rebound in the blood count data following drug administration is only possible for ν >1. The estimates of model parameters were ν = 3.02, ß = 0.00863 1/h, and IC50 = 7.47 ng/mL. The calculated baseline mean transit times of basophils in the blood 7.2 h and extravascular tissues 104.9 h agree with the values reported in the literature. CONCLUSIONS: We introduced an age-structured population model to describe cell trafficking between the blood and extravascular tissues. The model was adopted to account for the inhibitory drug effect on the cell recirculation. We showed that the age structure is essential to explain the rebound observed in the blood count response to a single dose drug administration. The model was validated using the basophil responses to DEX treatment in healthy subjects.

PKPD Modeling and Simulations to Support Biosimilar Development of Biphasic Insulin Aspart 30.

This paper presents an analysis of data from a comparative study of biosimilarity in terms of pharmacokinetics and pharmacodynamics in healthy volunteers using a hyperinsulinemic euglycemic clamp for reference and test biphasic insulin aspart 30 (BIAsp 30). As a result of the study, one of the secondary pharmacodynamic (PD) endpoints did not satisfy the classical criterion of 80%-125% (the lower limit for PD parameter area under the glucose infusion rate-time curve [ AUC GIR 0 - t ${\rm{AUC}}_{{\rm{GIR}}_{0 - {\rm{t}}}}$ ] turned out to be 79.5%). The main hypothesis explaining this result is that the sample size is insufficient to conduct a PD test with 90% statistical power, since the sample size has been calculated based on the coefficient of variation (CV) of pharmacokinetic (PK) parameters. To test this hypothesis, population PKPD (popPKPD) modeling and subsequent simulations of the required number of PD profiles were used. Two popPKPD models were constructed (a one-compartment double simultaneous absorption model for PK and an effect compartment Emax model for PD) to describe the PKPD data of reference and test insulins. As a result, using real data along with model-based simulation data, a biosimilarity test for PD was performed, and the lower limit for AUC GIR 0 - t ${\rm{AUC}}_{{\rm{GIR}}_{0 - {\rm{t}}}}$ became 82.6%, while the CV decreased from 31.7% to 24.1%. Thus, popPKPD modeling and simulations have been shown to be effective in interpreting and supporting the results of clinical biosimilarity trials.

Population pharmacokinetic/pharmacodynamic analysis of AK111, an IL-17A monoclonal antibody, in subjects with moderate-to-severe plaque psoriasis.

AK111 is an innovative IL-17A antibody, presenting high affinity to IL-17A and showing similar pharmacokinetic (PK) characteristics to those of typical immunoglobulin (Ig) G1 antibodies. To optimize the dosage regimen for phase 2/3 clinical trials, PK and pharmacodynamics (PD) of AK111 were first characterized in Chinese moderate-to-severe plaque psoriasis patients in a phase 1b study. AK111 PK serum sample and Psoriasis Area and Severity Index (PASI) score data were collected from 48 moderate-to-severe psoriasis patients in this study. Non-linear mixed-effects modeling was used for the population PK/PD analysis. A one-compartment model with a first-order absorption and a first-order elimination best described the PK behavior of AK111. The apparent systemic clearance was 0.182 L/day, and the central volume was 6.65 L. The exposure-response relationship was characterized using an indirect response model. The pharmacological effect of AK111 was described in the form of inhibiting the formation of psoriatic plaque, whereas placebo was quantified in the form of promoting the degradation of psoriatic skin lesions. The maximum effect of drug effect (Imax) and placebo effect (PLBmax) was 1 and 0.429, respectively. The rate constant for psoriatic plaque production (Kin) was 0.474 PASI/day and psoriatic plaque loss (Kout) was 0.024 day-1. The body surface area (BSA) affected by psoriasis was identified as a significant covariate on K o u t . The simulation results confirmed that all of the predicted PASI90 response rates at week 12 were higher than 60% at 150 and 300 mg dose levels with different regimens and could reach higher than 80% at week 24. We hope this first PK/PD study of AK111 in Chinese moderate-to-severe plaque psoriasis patients will be of help in the further clinical development of AK111 and provide a reference to the dosage optimization for similar antibodies with a long half-life.

Recent lessons learned from population pharmacokinetic studies of mycophenolic acid: physiological, genomic, and drug interactions leading to the prediction of drug effects.

INTRODUCTION: Mycophenolic acid (MPA) is a widely used immunosuppressant in transplantation and autoimmune disease. Highly variable pharmacokinetics have been observed with MPA, but the exact mechanisms remain largely unknown. AREAS COVERED: The current review provided a critical, comprehensive update of recently published population pharmacokinetic/dynamic models of MPA (n = 16 papers identified from PubMed and Embase, inclusive from January 2017 to August 2021), with specific emphases on the intrinsic and extrinsic factors influencing the pharmacology of MPA. The significance of the identified covariates, potential mechanisms, and comparisons to historical literature have been provided. EXPERT OPINION: While select covariates affecting the population pharmacokinetics of MPA are consistently observed and mechanistically supported (e.g. cyclosporine and post-transplant time on MPA clearance), some variables have not been regularly reported and/or lacked mechanistic explanation (e.g. diarrhea and several genetic polymorphisms). Very few pharmacodynamic models were available, pointing to the need to extrapolate pharmacokinetic findings. Ideal models of MPA should consist of: i) utilizing optimal sampling points to allow the characterizations of absorption, re-absorption, and elimination phases; ii) characterizing unbound/total MPA, MPA metabolites, plasma/urinary concentrations, and genetic polymorphisms to facilitate mechanistic interpretations; and iii) incorporating actual outcomes (e.g. rejection, leukopenia, infections) and pharmacodynamic data (e.g. inosine-5'-monophosphate dehydrogenase activities) to establish clinical relevance. We anticipate the field will continue to expand in the next 5 to 10 years.

Classification Tree Analysis Based On Machine Learning for Predicting Linezolid-Induced Thrombocytopenia.

Linezolid-induced thrombocytopenia is related to linezolid exposure, baseline platelet count and patient background. Although the relationship usually reflects the time of onset of thrombocytopenia, if the platelet maturation process is taken into account, the platelet decrease can be considered to have started at the beginning of treatment. To predict linezolid-induced thrombocytopenia, classification and regression tree (CART) analysis based on machine learning has been applied to identify predictive factors and cutoff values. We examined 74 patient data with or without linezolid-induced thrombocytopenia. Linezolid concentration and platelet count change, baseline platelet count, age, body weight and creatinine clearance estimate were evaluated as predictive factors for linezolid-induced thrombocytopenia. CART analysis selected the final tree containing two cutoff values: a platelet count reduction to less than 2.3% from baseline at 96 h after the initial dose and a linezolid concentration greater than or equal to 13.5 mg/L at 96 h after the initial dose. The targets for therapeutic intervention were concluded to be the linezolid concentration and the platelet change from baseline at 96 h after the initial dose. These cutoff values occur prior to the onset of thrombocytopenia and should be monitored to avoid linezolid-induced thrombocytopenia.

Corrigendum: Population Pharmacokinetic/Pharmacodynamic Model-Guided Dosing Optimization of a Novel Sedative HR7056 in Chinese Healthy Subjects.

[This corrects the article DOI: 10.3389/fphar.2018.01316.].

Antiplatelet Effect of a Pulaimab [Anti-GPIIb/IIIa F(ab)2 Injection] Evaluated by a Population Pharmacokinetic-pharmacodynamic Model.

BACKGROUND: Cardiovascular disease has one of the highest mortality rates among all the diseases. Platelets play an important role in the pathogenesis of cardiovascular diseases. Platelet membrane glycoprotein GPIIb/IIIa antagonists are the most effective antiplatelet drugs, and pulaimab is one of these. The study aims to promote individual medication of pulaimab [anti-GPIIb/IIIa F(ab)2 injection] by discovering the pharmacological relationship among the dose, concentration, and effects. The goal of this study is to establish a population pharmacokineticpharmacodynamic model to evaluate the antiplatelet effect of intravenous pulaimab injection. METHODS: Data were collected from 59 healthy subjects who participated in a Phase-I clinical trial. Plasma concentration was used as the pharmacokinetic index, and platelet aggregation inhibition rate was used as the pharmacodynamic index. The basic pharmacokinetics model was a two-compartment model, whereas the basic pharmacodynamics model was a sigmoid-EMAX model with a direct effect. The covariable model was established by a stepwise method. The final model was verified by a goodness-of-fit method, and predictive performance was assessed by a Bootstrap (BS) method. RESULTS: In the final model, typical population values of the parameters were as follows: central distribution Volume (V1), 183 L; peripheral distribution Volume (V2), 349 L; Central Clearance (CL), 31 L/h; peripheral clearance(Q), 204 L/h; effect compartment concentration reaching half of the maximum effect (EC50), 0.252 mg/L; maximum effect value (EMAX), 54.0%; and shape factor (γ), 0.42. In the covariable model, thrombin time had significant effects on CL and EMAX. Verification by the goodness-of-fit and BS methods showed that the final model was stable and reliable. CONCLUSION: A model was successfully established to evaluate the antiplatelet effect of intravenous pulaimab injection that could provide support for the clinical therapeutic regimen.

Population Pharmacokinetic/Pharmacodynamic Model-Guided Dosing Optimization of a Novel Sedative HR7056 in Chinese Healthy Subjects.

HR7056 is a new benzodiazepine, showing more faster acting onset and recovery than currently available short-acting sedatives. To avoid inadequate anesthesia and predict return of cognition, allowing for immediate neurological evaluation, HR7056 pharmacokinetics and pharmacodynamics were characterized in Chinese healthy subjects. We report on modeling of the data and simulations of dosage regimens for future study. Up to 63 subjects were evaluated, using Bispectral Index (BIS) and Modified Observer's Assessment of Alertness/Sedation (MOAA/S) as pharmacodynamics endpoints. A three-compartment model best described HR7056 pharmacokinetics. Total clearance was 1.49 L min-1, central volume was 2.1 L, inter-compartmental clearances were 0.96 and 0.27 L min-1, respectively. The population mean pharmacodynamic parameters were as follows: BIS, E0: 95.3; IC50: 503 ng mL-1; γ: 1.5; ke0: 0.0855 min-1; Imax: 47.9 and MOAA/S, IC50: 436 ng mL-1; γ: 1.5; ke0: 0.05 min-1; Imax: 27.9. The model simulation will enable maintenance doses to be given more accurately for future study. Clinical Trial Registration: identifier: NCT01970072.

Characterization of changes in HbA1c in patients with and without secondary failure after metformin treatments by a population pharmacodynamic analysis using mixture models.

The objective of the present study was to develop a population pharmacodynamic (PPD) model to describe the glycated hemoglobin (HbA1c)-lowering effects of metformin in type 2 diabetes mellitus patients with and without secondary failure and to characterize changes in HbA1c levels in the two subpopulations using a mixture model. Information on patients was collected retrospectively from electronic medical records. In this study, the mixture model was used to characterize the bimodal effects of metformin. A PPD analysis was performed using NONMEM 7.3.0. A physiological indirect response model, based on 829 HbA1c levels of 69 patients, described the time course for the HbA1c-lowering effects of metformin. Evidence for the different effectiveness of metformin subpopulations was provided using the mixture model. In the final PPD model, the inhibition effect was constant over a study duration in a patient subpopulation without secondary failure. In contrast, the inhibition effect decreased as a function of time after start of metformin treatment in a subpopulation with secondary failure. These results indicated that HbA1c improvements appeared to deteriorate over time in patients with secondary failure. In a PPD analysis of metformin, it was possible to assign patients with secondary failure using the mixture model.

Prediction of the human in vivo antiplatelet effect of S- and R-indobufen using population pharmacodynamic modeling and simulation based on in vitro platelet aggregation test.

Indobufen (Ibustrin®), a reversible inhibitor of platelet aggregation, exists in two enantiomeric forms in 1:1 ratio. Here, we characterized the anti-platelet effect of S- and R-indobufen using response surface modeling using NONMEM® and predicted the therapeutic doses exerting the maximal efficacy of each enantioselective S- and R-indobufen formulation. S- and R-indobufen were added individually or together to 24 plasma samples from drug-naïve healthy subjects, generating 892 samples containing randomly selected concentrations of the drugs of 0-128 mg/L. Collagen-induced platelet aggregation in platelet-rich plasma was determined using a Chrono-log Lumi-Aggregometer. Inhibitory sigmoid Imax model adequately described the anti-platelet effect. The S-form was more potent, whereas the R-form showed less inter-individual variation. No significant interaction was observed between the two enantiomers. The anti-platelet effect of multiple treatments with 200 mg indobufen twice daily doses was predicted in the simulation study, and the effect of S- or R-indobufen alone at various doses was predicted to define optimal dosing regimen for each enantiomer. Simulation study predicted that 200 mg twice daily administration of S-indobufen alone will produce more treatment effect than S-and R-mixture formulation. S-indobufen produced treatment effect at lower concentration than R-indobufen. However, inter-individual variation of the pharmacodynamic response was smaller in R-indobufen. The present study suggests the optimal doses of R-and S-enantioselective indobufen formulations in terms of treatment efficacy for patients with thromboembolic problems. The proposed methodology in this study can be applied to the develop novel enantio-selective drugs more efficiently.