A clinical population pharmacokinetic/pharmacodynamic model for BIIB059, a monoclonal antibody for the treatment of systemic and cutaneous lupus erythematosus.

A population pharmacokinetic/pharmacodynamic (popPK/PD) model for BIIB059 (anti-blood dendritic cell antigen 2 [anti-BDCA2]), a humanized immunoglobulin G1 monoclonal antibody currently under development for the treatment of SLE and CLE, is presented. BIIB059 binds BDCA2, a plasmacytoid dendritic cell (pDC)-specific receptor that inhibits the production of IFN-I and other inflammatory mediators when ligated. Phase 1 PK and PD data of healthy adult volunteers (HV, n = 87) and SLE subjects (n = 22) were utilized for the development of the popPK/PD model. The data included single and multiple dosing of intravenous and subcutaneous BIIB059. BDCA2 internalization (PD marker) was measured for all subjects by monitoring reduction of BDCA2 on pDC cell surface and used for development of the popPD model. A two-compartment popPK model with linear plus non-linear elimination was found to best describe BIIB059 PK. BDCA2 levels were best captured using an indirect response model with stimulation of the elimination of BDCA2. Clearance in SLE subjects was 25% higher compared to HV (6.87 vs 5.52 mL/h). Bodyweight was identified as only other covariate on clearance and central volume. The estimates of EC50 and Emax were 0.35 µg/mL and 8.92, respectively. No difference in EC50 and Emax was observed between SLE and HV. The popPK/PD model described the data accurately, as evaluated by pcVPCs and bootstrap. The presented popPK/PD model for BIIB059 provides valuable insight into the dynamics and dose-response relationship of BIIB059 for the treatment of SLE and CLE and was used to guide dose selection for the Phase 2 clinical study (NCT02847598).

Reducing whole body physiologically based pharmacokinetic models using global sensitivity analysis: diazepam case study.

There are situations in drug development where one may wish to reduce the dimensionality and complexity of whole body physiologically based pharmacokinetic models. A technique for formal reduction of such models, based on global sensitivity analysis, is suggested. Using this approach mean and variance of tissue(s) and/or blood concentrations are preserved in the reduced models. Extended Fourier amplitude sensitivity test (FAST), a global sensitivity technique, takes a sampling approach, acknowledging parameter variability and uncertainty, to calculate the impact of parameters on concentration variance. We used existing literature rules for formal model reduction to identify all possible smaller dimensionally models. To discriminate among those competing mechanistic models extended FAST was used, whereby we treated model structural uncertainty as another factor contributing to the overall uncertainty. A previously developed 14 compartment whole body physiologically based model for diazepam disposition in rat was reduced to three alternative reduced models, with preserved arterial mean and variance concentration profiles.

A fuzzy physiologically based pharmacokinetic modeling framework to predict drug disposition in humans.

To date, the application of physiologically based pharmacokinetic (PBPK) models in support of drug discovery remains limited, in part due to information deficit and uncertainty regarding model parameters. Fuzzy set theory provides a suitable way to objectively account for parameter uncertainty in models. Here, we present a fuzzy set-based PBPK modeling framework and demonstrate its utility in predicting diazepam pharmacokinetics in human plasma, following intravenous dosing, from available animal in vivo and literature data. For computationally expensive PBPK models, the sparse grid method is proposed as an efficient alternative to commonly used fuzzy arithmetic algorithms for function simulation.

A fuzzy physiologically based pharmacokinetic modeling framework to predict drug disposition in humans.

To date, the application of physiologically based pharmacokinetic (PBPK) models in support of drug discovery remains limited, in part due to information deficit and uncertainty regarding model parameters. Fuzzy set theory provides a suitable way to objectively account for parameter uncertainty in models. Here, we present a fuzzy set-based PBPK modeling framework and demonstrate its utility in predicting diazepam pharmacokinetics in human plasma, following intravenous dosing, from available animal in vivo and literature data. For computationally expensive PBPK models, the sparse grid method is proposed as an efficient alternative to commonly used fuzzy arithmetic algorithms for function simulation.

Uncertainty analysis in pharmacokinetics and pharmacodynamics: application to naratriptan.

PURPOSE: The aim of the study was to predict pain relief of migraine in patients following naratriptan oral (tablet) administration by using uncertainty analysis. The analysis was based on phase I pharmacokinetic naratriptan data, sumatriptan pharmacodynamic data, and naratriptan preclinical (animal) potency information, together with general knowledge as to how migraine affects oral absorption. METHODS: A previously developed pharmacokinetic (PK)/pharmacodynamic (PD) model for naratriptan disposition and effect was used. The uncertain parameters in the model, which were associated with absorption and scaling between first-in-class compound sumatriptan and naratriptan, were modeled using fuzzy sets theory. Global sensitivity analysis was then used to investigate the impact of each PK/PD parameter on the responses. RESULTS: Acknowledging parametric uncertainty did not improve prediction of the probability of pain relief. Global sensitivity analysis demonstrated that predictions were heavily influenced by interindividual variability in pharmacodynamics, as the dose response relationship was relatively insensitive to the pharmacokinetics. CONCLUSIONS: To predict the probability of pain relief following oral (tablet) administration of naratriptan, a simple dose response, instead of the PK/PD model, would have yielded very similar predictions. The naratriptan PK/PD model may be improved by either refining the PD model or better still by specifying the interindividual error by additional data collecting with an improved design.

Fuzzy simulation of pharmacokinetic models: case study of whole body physiologically based model of diazepam.

The aim of the present study is to develop and implement a methodology that accounts for parameter variability and uncertainty in the presence of qualitative and semi-quantitative information (fuzzy simulations) as well as when some parameters are better quantitatively defined than others (fuzzy-probabilistic approach). The fuzzy simulations method consists of (i) representing parameter uncertainty and variability by fuzzy numbers and (ii) simulating predictions by solving the pharmacokinetic model. The fuzzy-probabilistic approach includes an additional transformation between fuzzy numbers and probability density functions. To illustrate the proposed method a diazepam WBPBPK model was used where the information for hepatic intrinsic clearance determined by in vitro-in vivo scaling was semi-quantitative. The predicted concentration time profiles were compared with those resulting from a Monte Carlo simulation. Fuzzy simulations can be used as an alternative to Monte Carlo simulation.

Development of a whole body physiologically based model to characterise the pharmacokinetics of benzodiazepines. 1: Estimation of rat tissue-plasma partition ratios.

Three methods for estimation of the equilibrium tissue-to-plasma partition ratios (Kp values) in the presence of tissue concentration time data have been investigated. These are the area method, the open loop (tissue specific) method and the whole body model(closed loop) method, each with different model assumptions. Additionally, multiple imputations, a technique for dealing with deficiencies in data sets (i.e., missing tissues) is used. The estimated Kp values by the three methods have been compared and the limitations and advantages of each approach drawn. The area method, which is essentially model free, gives only a crude estimate of Kp without making any statement of its uncertainty; whereas both the open and closed loop methods provide an estimate of this. The closed loop method, where the most assumptions are made, is the approach that gives the best overall estimates of Kp, which was confirmed by comparing the predicted concentration-time profiles with experimental data. Although the estimates from the closed loop method, as well as the other two methods, are conditioned on the data, they are the most reliable for both propagating parameter variability and uncertainty through a whole body physiologically based model, as well as for extrapolation to human. A series of benzodiazepines, namely alprazolam, chlordiazepoxide, clobazam, diazepam, flunitrazepam, midazolam and triazolam in rat is used as a case study in the current investigation.