Development and Evaluation of a Physiologically Based Pharmacokinetic Model for Predicting the Effects of Anti-FcRn Therapy on the Disposition of Endogenous IgG in Humans

This work scaled up a previously developed physiologically based pharmacokinetic model to predict the effects of anti-FcRn agents on the disposition of endogenous IgG in human subjects. Simulations were performed with the scaled model to predict the effects of single- and multiple-dose administration of anti-FcRn monoclonal antibodies (1-256 mg/kg) and high-dose intravenous immune globulin (0.4-2 g/kg). The model was evaluated for prediction accuracy through comparison to the effects of rozanolixizumab, an anti-FcRn monoclonal antibodies under current clinical evaluation, on the disposition of endogenous IgG in healthy human subjects. The model provided reasonably accurate predictions of the effects of rozanolixizumab. Prediction errors for the maximum reduction in endogenous IgG concentrations were ?8.50% (90% model prediction interval: ?14.0% to 1.44%), 3.33% (90% model prediction interval: ?13.9% to 21.2%), and 6.85% (90% model prediction interval: ?35.2% to 10.5%) for rozanolixizumab doses of 1, 4, and 7 mg/kg, respectively. Model simulations predict that anti-FcRn therapies will exhibit greater dose potency in healthy volunteers than in patients with elevated IgG production rates (e.g., as typically found in autoimmune disease). The model appears to have potential for use in assessing and predicting novel dosing strategies for anti-FcRn therapies.     

Physiologically Based Pharmacokinetic Modeling to Supplement Nilotinib Pharmacokinetics and Confirm Dose Selection in Pediatric Patients

In adult patients, nilotinib is indicated for chronic myeloid leukemia at an approved oral dose of 300 or 400 mg BID. Physiologically based pharmacokinetic (PBPK) model was developed to describe and supplement limited PK data in the pediatric population ranging from 2 to less than 6 years of age and ultimately inform dosing regimen. An adult Simcyp PBPK model was established and verified with clinical pharmacokinetic data after a single or multiple oral doses of 400 mg nilotinib (230 mg/m²). The model was then applied to a pediatric PBPK model, taking account of ontogeny profiles of metabolizing enzymes and pediatric physiological parameters. The model was further verified using observed pediatric PK data in 12- to <18-year-old and from 6- to <12-year-old patients. The PBPK models were able to recover, describe, and supplement the limited nilotinib concentration-time data profile in 2- to <6-year-old patients after a single dose and C_min,ss after BID dosing. The exposure (C_max,ss, C_min,ss, and AUC_tau,ss) was predicted to be similar across age groups. PBPK model simulations confirmed that body surface area–normalized dosing regimen of 230 mg/m² is considered appropriate for pediatric patients >2 to <18 years of age.     

Physiologically Based Pharmacokinetic (PBPK) Modeling of Everolimus (RAD001) in Rats Involving Non-Linear Tissue Uptake

Everolimus is a novel macrolide immunosuppressant developed for the prophylaxis of allogeneic renal or cardiac transplant rejection. Treatments with immunosuppressants are often associated with organ toxicity that is linked to high organ exposure. Therefore, gaining insight into the pharmacokinetics of everolimus in various organs is highly desirable especially those organs of therapeutic interest or those that pose safety concerns. The aim of this work was to characterize the disposition kinetics of everolimus in rats by physiologically based pharmacokinetic (PBPK) modeling. Blood and tissue samples were collected from male Wistar rats over 24 hr following intravenous (iv) bolus and iv infusion of 1 mg/kg and 10 mg/kg/2 hr of everolimus. Further blood samples were collected between 1 and 170 hr from a third group of rats, which received iv infusion of 1 mg/kg/2 hr of everolimus. Drug concentrations in blood and tissues were determined by a liquid chromatography reverse dilution method. Distribution of everolimus between blood fractions was determined in vitro at 37°C. The results of the study demonstrated that everolimus exhibited moderate non-linear binding to red blood cells. Also, the tissue-to-blood concentration ratio decreased in all tissues as blood concentration increased. A PBPK model involving non-linear tissue binding was able to successfully describe the observed data in blood and all the organs investigated. The highest binding potential was observed in thymus, lungs, and spleen with the greatest tissue affinity observed in thymus, skin, and muscle as compared to other tissues. Everolimus exhibited a high clearance rate that was limited to the hepatic blood flow (47.2 ml/min/kg). The PBPK model was also able to predict the venous blood concentration reasonably well following oral administration. The oral bioavailability value, as estimated with the PBPK, was 12% and was similar to the value obtained by non-compartmental analysis. In conclusion, A PBPK model has been developed that successfully predicts the time course of everolimus in blood and a variety of organs. This model takes into account the non- linear binding of everolimus to red blood cells and tissues. This model may be used to predict everolimus concentration–time course in organs from other species including humans.     

Development and Qualification of a Physiologically Based Pharmacokinetic Model of Finasteride and Minoxidil Following Scalp Application

In this study, we aimed to develop and qualify a PBPK model for scalp application using two drugs with marked differences in physicochemical properties and PK profiles. The parameters related to scalp physiology, drug PK, and formulations were incorporated into a Multi-Phase and Multi-Layer (MPML) Mechanistic Dermal Absorption (MechDermA) model within the Simcyp® Simulator V17. The finasteride PBPK model was linked to its effect on dihydrotestosterone (DHT) levels in plasma and scalp using an indirect response model. Predicted PK (and PD for finasteride) profiles and parameters were compared against the clinically reported data and justified by visual predictive checks and two-fold error criteria for model verification. The PBPK/PD model for finasteride reasonably demonstrated an ability to predict its respective PK and PD profiles, and parameters following scalp application under various clinical scenarios. Using the same scalp physiological input parameters, the minoxidil PBPK model was then developed and satisfactorily qualified with independent clinical datasets. Collectively, these results suggested that the established PBPK model may have broader utility for other topical formulations intended for scalp application.     

Pharmacokinetics and Pharmacodynamics of an Investigational Antipsychotic Agent,CI-1007, in Rats and Monkeys

Purpose. To study the pharmacokinetics (PK) and pharmacodynamics (PD) of an investigational antipsychotic agent, CI-1007, in rats and monkeys. Methods. CI-1007 and a pharmacologically active metabolite, PD 147693 (Ml), were evaluated in animal antipsychotic tests (inhibition of dopamine neuron firing and spontaneous locomotor activity in rats, and inhibition of continuous avoidance in monkeys). Plasma concentrations of CI-1007 and Ml were determined using validated HPLC assays. Log-linear and link models were used for PK/PD analysis. Results. CI-1007 and Ml have shown similar effects on dopamine neuron firing (2.5 mg/kg i.p.), and produced dose-related effects on spontaneous locomotor activity in rats (0.3–30 mg/kg, p.o.) and on continuous avoidance in monkeys (0.6–1.2 mg/kg p.o.). After pharmacologically active CI-1007 doses, mean plasma CI-1007 C_max increased from 19 to 200 ng/ml in Sprague-Dawley rats at doses of 3–30 mg/ kg, and from 8.1 to 34 ng/ml in squirrel monkeys at doses of 0.6–1.2 mg/kg, but corresponding plasma M1 C_max values were near or below the limit of quantitation (5 ng/ml). CI-1007 EC50 was 31.1 ng/ml in rats, calculated from a log-linear regression. In monkeys, CI-1007 EC_e50, , and K_eo at 0.6 and 1.2 mg/kg were 4.8 and 4.5 ng/ml, 1.9 and 2.0, and 0.47 and 0.48 hr^–1, respectively, calculated by the link model. Conclusions. CI-1007 has shown dose-related pharmacokinetics and pharmacodynamics in rats and monkeys. Although Ml produces anti-psychotic-like effects similar to CI-1007, the contribution of Ml to the activity of the parent drug may not be significant in rats and monkeys as based on plasma levels. CI-1007 plasma concentration correlates log-linearly with inhibition effect from the rat locomotor study. The counter-clockwise hysteresis relationship of CI-1007 plasma concentration and inhibition effect from the monkey avoidance test was described by a link model, and the resulting C_e (concentration in effect compartment) versus effect profile exhibits a sigmoidal curve.     

A Physiologically Based Pharmacokinetic Modeling Approach to Predict Drug-Drug Interactions of Buprenorphine After Subcutaneous Administration of CAM2038 With Perpetrators of CYP3A4

CAM2038, FluidCrystal injection depot, is an extended release formulation of buprenorphine given subcutaneously every 1 week (Q1W) or every 4 weeks (Q4W). The purpose of this research was to predict the magnitude of drug-drug interaction (DDI) after coadministration of a strong CYP3A4 inducer or inhibitor using physiologically based pharmacokinetic (PBPK) modeling. A PBPK model was developed for CAM2038 based on the previously published buprenorphine PBPK model after intravenous and sublingual administration and the PK profiles after subcutaneous administration of CAM2038 from 2 phase I clinical trials. The strong CYP3A4 inhibitor ketoconazole was predicted to increase the buprenorphine exposure by 35% for the Q1W formulation and 34% for Q4W formulation, respectively. Also, the strong CYP3A4 inducer rifampin was predicted to decrease the buprenorphine exposure by 26% for both the Q1W and Q4W formulations. The results provided insight into the potential DDI effect for CAM2038 and suggested a lack of clinically meaningful DDI when CAM2038 is coadministered with CYP3A4 inhibitor or inducer. Therefore, no dose adjustment is required when CAM2038 is coadministered with CYP3A4 perpetrators.     

Physiologically Based Pharmacokinetic Modeling of Benzene Metabolism in Mice Through Extrapolation from In Vitro to In Vivo

Benzene (C₆H₆) is a highly flammable, colorless liquid. Ubiquitous exposures result from its presence in gasoline vapors, cigarette smoke, and industrial processes. Benzene increases the incidence of leukemia in humans when they are exposed to high doses for extended periods; however, leukemia risks in humans at low exposures are uncertain. The exposure-dose-response relationship of benzene in humans is expected to be nonlinear because benzene undergoes a series of metabolic transformations, detoxifying and activating, in the liver, resulting in multiple metabolites that exert toxic effects on the bone marrow. We developed a physiologically based pharmacokinetic model for the uptake and elimination of benzene in mice to relate the concentration of inhaled and orally administered benzene to the tissue doses of benzene and its key metabolites, benzene oxide, phenol, and hydroquinone. As many parameter values as possible were taken from the literature; in particular, metabolic parameters obtained from in vitro studies with mouse liver were used since comparable parameters are also available for humans. Parameters estimated by fitting the model to published data were first-order rate constants for pathways lacking in vitro data and the concentrations of microsomal and cytosolic protein, which effectively alter overall enzyme activity. The model was constrained by using the in vitro metabolic parameters (maximum velocities, first-order rate constants, and saturation parameters), and data from multiple laboratories and experiments were used. Despite these constraints and sources of variability, the model simulations matched the data reasonably well in most cases, showing that in vitro metabolic constants can be successfully extrapolated to predict in vivo data for benzene metabolism and dosimetry. Therefore in vitro metabolic constants for humans can subsequently be extrapolated to predict the dosimetry of benzene and its metabolites in humans. This will allow us to better estimate the risks of adverse effects from low-level benzene exposures.     

Inotropic Effect of Digoxin in Humans: Mechanistic Pharmacokinetic/Pharmacodynamic Model Based on Slow Receptor Binding

Purpose. The purpose of this study was to construct a mechanistic pharmacokinetic/pharmacodynamic (PK/PD) model for digoxin that describes the relationship between plasma concentration and inotropic response. Methods. On the basis of results obtained in the isolated perfused rat heart, a PK/PD model for digoxin in humans was developed. In fitting the model to previously published bolus dose and concentration clamp data (shortening of electromechanical systole), the plasma concentration-time curves were used as forcing functions in the computer program ADAPT II. Results. The mechanistic approach allowed a modeling of digoxin pharmacodynamics which is consistent with available inotropic response data. The estimates of the receptor binding parameters were in the same order of magnitude as those measured in vitro for ouabain. The mechanistic model explained the parameters of the empirical link model (EC₅₀, E _max and delay time ) in terms of the underlying processes, suggesting that the long equilibration half-time of 13 h is due to slow receptor binding. The empirical link model, in contrast, is not compatible with a noninstantaneous receptor binding process and led to estimates of the delay time that were dependent on the digoxin administration schedule. Conclusions. The new, mechanistic model may provide a rationale for better understanding of digoxin pharmacodynamics and could become a tool to bridge the gap between in vitro and in vivo studies.     

Pharmacokinetic and pharmacodynamic studies following the intravenous and oral administration of the antiparkinsonian drug biperiden to normal subjects

Summary The pharmacokinetics and pharmacodynamics (changes in pupil size and salivary flow) of biperiden following a single oral and intravenous dose were investigated in six normal subjects.After the injection plasma concentrations declined biphasically, with half-times of 1.5 h for the rapid phase and 24 h for the terminal phase. Clearance and apparent volume of distribution were high (12 ml·min^–1·kg^–1 and 24 l·kg^–1 respectively). Absorption was rapid but the systemic availability was incomplete (33%), probably due to first-pass metabolism.Central nervous system (CNS) adverse effects and changes in pupil size were observed after both routes of administration while salivary flow was affected only by the injection.     

Control-Relevant Modeling of the Antitumor Effects of 9-Nitrocamptothecin in SCID Mice Bearing HT29 Human Colon Xenografts

The mathematical model structure selected to describe system behavior is at least partially dependent on the proposed use of the model. In this paper, a pharmacokinetic(PK)/pharmacodynamic (PD) model for use in drug delivery algorithm synthesis is developed. The antitumor agent 9-nitrocamptothecin (9NC) was administered orally to severe combined immunodeficient (SCID) mice bearing subcutaneously implanted HT29 human colon xenografts, and the effect of 9NC on those xenografts was characterized. Different PK model structures were considered in characterizing the dynamics of the drug concentration in the plasma. Akaike’s Information Criterion (AIC) was used to select the model structure maximizing fit accuracy while simultaneously minimizing the number of model parameters. The resulting PK model was a set of coupled linear ordinary differential equations able to describe the nonlinear dynamic behavior (e.g. plateauing, etc.) of the drug concentrations observed in the plasma. Pharmacodynamics were modeled by characterizing tumor growth in both the untreated and drug-treated animals. The resulting PK/PD model related drug administration to effect, and this model has a structure that facilitates future control algorithm synthesis. Control algorithms in this context would directly utilize PK/PD model predictions. These predictions would be used to determine the amount and frequency of drug administration in order to reduce the tumor burden without violating clinically relevant constraints. This methodology could then be used to aid the clinician in selecting dose levels and schedules, and extension to patient tailored treatment may eventually be feasible with this approach.     

Predicting Nonlinear Pharmacokinetics of Omeprazole Enantiomers and Racemic Drug Using Physiologically Based Pharmacokinetic Modeling and Simulation: Application to Predict Drug/Genetic Interactions

Purpose

The objective of this study is to develop a physiologically-based pharmacokinetic (PBPK) model for each omeprazole enantiomer that accounts for nonlinear PK of the two enantiomers as well as omeprazole racemic drug.

Methods

By integrating in vitro, in silico and human PK data, we first developed PBPK models for each enantiomer. Simulation of racemic omeprazole PK was accomplished by combining enantiomer models that allow mutual drug interactions to occur.

Results

The established PBPK models for the first time satisfactorily predicted the nonlinear PK of esomeprazole, R-omeprazole and the racemic drug. The modeling exercises revealed that the strong time-dependent inhibition of CYP2C19 by esomeprazole greatly altered the R-omeprazole PK following administration of racemic omeprazole as in contrast to R-omeprazole given alone. When PBPK models incorporated both autoinhibition of each enantiomer and mutual interactions, the ratios between predicted and observed AUC following single and multiple dosing of omeprazole were 0.97 and 0.94, respectively.

Conclusions

PBPK models of omeprazole enantiomers and racemic drug were developed. These models can be utilized to assess CYP2C19-mediated drug and genetic interaction potential for omeprazole and esomeprazole.     

Pharmacokinetic considerations related to therapeutic drug monitoring of tacrolimus in kidney transplant patients

Introduction: Tacrolimus (Tac) is the cornerstone of immunosuppressive therapy after solid organ transplantation and will probably remain so. Excluding belatacept, no new immunosuppressive drugs were registered for the prevention of acute rejection during the last decade. For several immunosuppressive drugs, clinical development halted because they weren’t sufficiently effective or more toxic.

Areas covered: Current methods of monitoring Tac treatment, focusing on traditional therapeutic drug monitoring (TDM), controversies surrounding TDM, novel matrices, pharmacogenetic and pharmacodynamic monitoring are discussed.

Expert opinion: Due to a narrow therapeutic index and large interpatient pharmacokinetic variability, TDM has been implemented for individualization of Tac dose to maintain drug efficacy and minimize the consequences of overexposure. The relationship between predose concentrations and the occurrence of rejection or toxicity is controversial. Acute cellular rejection also occurs when the Tac concentration is within the target range, suggesting that Tac whole blood concentrations don’t necessarily correlate with pharmacological effect. Intracellular Tac, the unbound fraction of Tac or pharmacodynamic monitoring could be better biomarkers/tools for adequate Tac exposure – research into this has been promising. Traditional TDM, perhaps following pre-emptive genotyping for Tac-metabolizing enzymes, must suffice for a few years before these strategies can be implemented in clinical practice.     

Pharmacodynamic Assessment of the Benztropine Analogues AHN-1055 and AHN-2005 Using Intracerebral Microdialysis to Evaluate Brain Dopamine Levels and Pharmacokinetic/Pharmacodynamic Modeling

No HeadingPurpose. The benztropine (BZT) analogues bind with high affinity to the dopamine transporter (DAT) and demonstrate a behavioral and pharmacokinetic profile unlike that of cocaine. The development of a predictive pharmacokinetic/pharmacodynamic (PK/PD) model to characterize the concentration-effect relationship between the BZT analogues and brain dopamine (DA) levels is an important step in the evaluation of these compounds as potential cocaine abuse pharmacotherapies. Hence, the objective of this study was to mathematically characterize the PD of BZT analogues and cocaine, using appropriate PK/PD models.Methods. Dialysis probes were stereotaxically implanted into the nucleus accumbens of Sprague-Dawley rats (275–300 g). Extracellular fluid (ECF) DA levels were measured after intravenous administration of the BZT analogues AHN-1055 and AHN-2005, as well as cocaine using high performance liquid chromatography-electrochemical detection (HPLC-ECD). PD models were used to describe the relationship between the BZT analogues or cocaine and brain microdialysate DA, and suitability was based on standard goodness-of-fit criteria.Results. The BZT analogues produced a sustained increase in brain microdialysate DA levels in comparison to cocaine. The time of maximum concentration (T_max) for brain microdialysate DA was 2 h for AHN-1055 and 1 h for AHN-2005 compared to a T_max of 10 min for cocaine. The duration of brain microdialysate DA elevation was 12–24 h for the BZTs in comparison to 1 h for cocaine. An indirect model with inhibition of loss of response and a sigmoid E_max model best described the PK/PD for the BZT analogues and cocaine, respectively. The 50% of maximum inhibition (IC₅₀) of the loss of DA was lower for AHN-2005 (226 ± 27.5 ng/ml) compared to AHN-1055 (321 ± 19.7 ng/ml). In addition, the EC₅₀ for cocaine was 215 ± 11.2 ng/ml.Conclusions. The slow onset and long duration of BZT analogue–induced DA elevation may avoid the reinforcing effects and craving of cocaine. Further, the developed models will be useful in characterizing the PK/PD of other analogues and aid in the assessment of the therapeutic efficacy of the BZT analogues as substitute medications for cocaine abuse.     

Effects of HM30181, a P-glycoprotein inhibitor,on the pharmacokinetics and pharmacodynamics of loperamide in healthy volunteers

Aims

HM30181 is a third generation P-glycoprotein (P-gp) inhibitor currently under development. The objectives of this study were to evaluate the effects of a single dose of HM30181 on the pharmacodynamics and pharmacokinetics of loperamide, a P-gp substrate, and to compare them with those of quinidine.

Methods

Eighteen healthy male subjects were administered loperamide alone (period 1) or with loperamide plus quinidine or HM30181 in period 2 or 3, respectively. In period 3, subjects randomly received one of three HM30181 doses: 15, 60 or 180 mg. Changes in pupil size, alertness, oxygen saturation and the oral bioavailability of loperamide were assessed in each period. In addition, the pharmacokinetics of HM30181 were determined.

Results

Pupil size, alertness and oxygen saturation did not change over time when loperamide alone or loperamide plus HM30181 was administered while HM30181 significantly increased the systemic exposure to loperamide, i.e. the geometric mean ratio (90% confidence interval) of AUC(0,t_last) for loperamide with and without HM30181 was 1.48 (1.08, 2.02). Co-administered quinidine significantly increased the systemic exposure to loperamide 2.2-fold (1.53, 3.18), which also markedly reduced pupil size, resulting in a decrease of 24.7 mm h in the area under the effect curve of pupil size change from baseline compared with loperamide alone.

Conclusions

HM30181 inhibits P-gp mainly in the intestinal endothelium, which can be beneficial because pan-inhibition of P-gp, particularly in the brain, could lead to detrimental adverse events. Further studies are warranted to investigate adequately the dose–exposure relationship of HM30181, along with its duration of effect.     

Whole-Body Physiologically Based Pharmacokinetic Modeling of Trastuzumab and Prediction of Human Pharmacokinetics

In the present study, we evaluated the pharmacokinetics (PK) of trastuzumab and sought to predict human PK based on animal studies, through the use of optical imaging and a whole-body physiologically based pharmacokinetic (WB-PBPK) modeling approach. The PK study was conducted in 24 mice, where serial blood samples were withdrawn and major organs were isolated after the last blood withdrawal. The drug concentrations in blood and major organs were measured via optical imaging. The WB-PBPK model was constructed using known physiological values including the volumes of major organs and blood/lymphatic flow. The NONMEM software (version 7.3) was used to determine tissue partition coefficients. Using the WB-PBPK model, a clinical trial simulation was performed with reference to human physiological values acquired from the literature. The simulated human PK was then compared with the actual PK observed in the previous study in which healthy male subjects received 6 mg/kg trastuzumab (Herceptin^®) via intravenous route. The ratio of the simulated versus observed area under the concentration-time curve was 1.02 and that of maximal concentration was 0.72. The current study describes the potential synergistic applications of WB-PBPK and optical imaging in human PK prediction based on preclinical data obtained in early-stage drug development.     

Effects on the Pharmacokinetics and Pharmacodynamics in the Elderly of Coadministering Ramipril with Water,Apple Juice,and Applesauce

Pharmaceutical Research -     

Pharmacokinetic/Pharmacodynamic Evaluation of Deflazacort in Comparison to Methylprednisolone and Prednisolone

Purpose. The pharmacokinetics and pharmacodynamics of deflazacort after oral administration (30 mg) to healthy volunteers were determined and compared with those of 20 mg of methylprednisolone and 25 mg of prednisolone. Methods. Methylprednisolone, prednisolone and the active metabolite of deflazacort, 21-desacetyldeflazacort, were measured in plasma using HPLC. For the assessment of pharmacodynamics, differential white blood cell counts were obtained over 24 hours. An integrated pharmacokinetic-pharmacodynamic (PK-PD) model was applied to link corticosteroid concentrations to the effect on lymphocytes and granulocytes. Results. Deflazacort is an inactive prodrug which is converted rapidly to the active metabolite 21-desacetyldeflazacort. Maximum concentrations of 21-desacetyldeflazacort averaged 116 ng/ml and were observed after 1.3 h. The average area under the curve was 280 ng/ml · h, and the terminal half-life was 1.3 h. 21-Desacetyldeflazacort was cleared significantly faster than both methylprednisolone and prednisolone. The PK-PD-model was suitable to describe time course and magnitude of the observed effects. The results were consistent with reported values for glucocorticoid receptor binding affinities for the investigated compounds. Conclusions. Due to the short pharmacokinetic half-life of its active metabolite, pharmacodynamic effects of deflazacort are of shorter duration than those of methylprednisolone and prednisolone. The PK-PD model allows good prediction of pharmacodynamic effects based on pharmacokinetic and receptor binding data.     

Physiologically Based Pharmacokinetic Modeling of Drug Disposition in Rat and Human: A Fuzzy Arithmetic Approach

PURPOSE: Probabilistic methods are insufficient for dealing with the vagueness inherent in human judgment of minimal data available during early drug development. We sought to use fuzzy set theory as a basis for quantifying and propagating vague judgment in a physiologically based pharmacokinetic (PBPK) model for diazepam disposition. MATERIALS AND METHODS: First, using diazepam distribution data in rat tissues and fuzzy regression, we estimated fuzzy rat tissue-to-plasma partition coefficients (Kp's). We scaled the coefficients prior to human PBPK modeling. Next, we constructed the fuzzy set of hepatic intrinsic clearance (CLint) by integrating CLint values measured in vitro from human hepatocytes. Finally, we used these parameters, and other physiological and biochemical information, to predict human diazepam disposition. We compared the simulated plasma kinetics with published concentration-time profiles. RESULTS: We successfully identified rat Kp's by fuzzy regression. The predicted rat tissue concentration-time contours enveloped the animal tissue distribution data. For the human PBPK model, the mean in vivo plasma concentrations were contained in the simulated concentration-time envelopes. CONCLUSIONS: We present a novel computational approach for handling information paucity in PBPK models using fuzzy arithmetic. Our methodology can model the vagueness associated with human perception and interpretation of minimal drug discovery data.     

Physiologically Based Pharmacokinetic Modeling of a Homologous Series of Barbiturates in the Rat: A Sensitivity Analysis

Sensitivity analysis studies the effects of the inherent variability and uncertainty in model parameters on the model outputs and may be a useful tool at all stages of the pharmacokinetic modeling process. The present study examined the sensitivity of a whole-body physiologically based pharmacokinetic (PBPK) model for the distribution kinetics of nine 5-n-alkyl-5-ethyl barbituric acids in arterial blood and 14 tissues (lung, liver, kidney, stomach, pancreas, spleen, gut, muscle, adipose, skin, bone, heart, brain, testes) after iv bolus administration to rats. The aims were to obtain new insights into the model used, to rank the model parameters involved according to their impact on the model outputs and to study the changes in the sensitivity induced by the increase in the lipophilicity of the homologues on ascending the series. Two approaches for sensitivity analysis have been implemented. The first, based on the Matrix Perturbation Theory, uses a sensitivity index defined as the normalized sensitivity of the 2-norm of the model compartmental matrix to perturbations in its entries. The second approach uses the traditional definition of the normalized sensitivity function as the relative change in a model state (a tissue concentration) corresponding to a relative change in a model parameter. Autosensitivity has been defined as sensitivity of a state to any of its parameters; cross-sensitivity as the sensitivity of a state to any other states' parameters. Using the two approaches, the sensitivity of representative tissue concentrations (lung, liver, kidney, stomach, gut, adipose, heart, and brain) to the following model parameters: tissue-to-unbound plasma partition coefficients, tissue blood flows, unbound renal and intrinsic hepatic clearance, permeability surface area product of the brain, have been analyzed. Both the tissues and the parameters were ranked according to their sensitivity and impact. The following general conclusions were drawn: (i) the overall sensitivity of the system to all parameters involved is small due to the weak connectivity of the system structure; (ii) the time course of both the auto- and cross-sensitivity functions for all tissues depends on the dynamics of the tissues themselves, e.g., the higher the perfusion of a tissue, the higher are both its cross-sensitivity to other tissues' parameters and the cross-sensitivities of other tissues to its parameters; and (iii) with a few exceptions, there is not a marked influence of the lipophilicity of the homologues on either the pattern or the values of the sensitivity functions. The estimates of the sensitivity and the subsequent tissue and parameter rankings may be extended to other drugs, sharing the same common structure of the whole body PBPK model, and having similar model parameters. Results show also that the computationally simple Matrix Perturbation Analysis should be used only when an initial idea about the sensitivity of a system is required. If comprehensive information regarding the sensitivity is needed, the numerically expensive Direct Sensitivity Analysis should be used.