Transit Compartment Research Articles

Population pharmacokinetic−pharmacodynamic modelling of mycophenolic acid in paediatric renal transplant recipients in the early post‐transplant period

The purpose of this study was to develop a population pharmacokinetic and pharmacodynamic (PK-PD) model for mycophenolic acid (MPA) in paediatric renal transplant recipients in the early post-transplant period. A total of 214 MPA plasma concentrations-time data points from 24 patients were available for PK model development. In 17 out of a total of 24 patients, inosine monophosphate dehydrogenase (IMPDH) enzyme activity measurements (n = 97) in peripheral blood mononuclear cells were available for PK-PD modelling. The PK-PD model was developed using non-linear mixed effects modelling sequentially by 1) developing a population PK model and 2) incorporating IMPDH activity into a PK-PD model using post hoc Bayesian PK parameter estimates. Covariate analysis included patient demographics, co-medication and clinical laboratory data. Non-parametric bootstrapping and prediction-corrected visual predictive checks were performed to evaluate the final models. A two compartment model with a transit compartment absorption best described MPA PK. A non-linear relationship between dose and MPA exposure was observed and was described by a power function in the model. The final population PK parameter estimates (and their 95% confidence intervals) were CL/F, 22 (14.8, 25.2) l h(-1) 70 kg(-1) ; Vc /F, 45.4 (29.6, 55.6) l; Vp /F, 411 (152.6, 1472.6)l; Q/F, 22.4 (16.0, 32.5) l h(-1) ; Ka , 2.5 (1.45, 4.93) h(-1) . Covariate analysis in the PK study identified body weight to be significantly correlated with CL/F. A simplified inhibitory Emax model adequately described the relationship between MPA concentration and IMPDH activity. The final population PK-PD parameter estimates (and their 95% confidence intervals) were: E0 , 3.45 (2.61, 4.56) nmol h(-1) mg(-1) protein and EC50 , 1.73 (1.16, 3.01) mg l(-1) . Emax was fixed to 0. There were two African-American patients in our study cohorts and both had low IMPDH baseline activities (E0 ) compared with Caucasian patients (mean value 2.13 mg l(-1) vs. 3.86 mg l(-1) ). An integrated population PK-PD model of MPA has been developed in paediatric renal transplant recipients. The current model provides information that will facilitate future studies and may be implemented in a Bayesian algorithm to allow a PK-PD guided therapeutic drug monitoring strategy.

GP3-6: Mixed effects PK/PD modelling of NNC0195-0092 - a noncovalent albumin binding growth hormone derivative in the hypophysectomised rat

Introduction: The hypophysectomised rat is used during preclinical development to evaluate the pharmacokinetic and pharmacodynamic (PKPD) properties of new GH analogues . The animal model resembles the human GH deficiency with lack of GH secretion, low levels of serum IGF-I and growth retardation . Non-linear mixed effects modelling of the PKPD relationship may improve the mechanistic understanding of action, and form the basis for translation of pre-clinical parameters into meaningful and accurate human predictions for early clinical development . A mechanistic mixed-effects PKPD model for NNC0195-0092 a non-covalent albumin binding growth hormone derivative has been developed in hypophysectomised Sprague Dawley rats . Methods: Animal studies with single and multiple dosing of NNC0195-0092 was used to establish the PKPD relationship between NNC0195-0092 and insulin-like growth factor I (IGF-I) Non-linear mixed effects modelling software (NONMEM) was used to model the simultaneous relationship between NNC01950092, IGF-I and bodyweight gain . Derived parameters was scaled directly from model values using the allometric approach with exponents fixed to conventional values . Results: A final PKPD model was established, that described the PK of NNC0195-0092 as a two compartmental model with both linear and non-linear elimination terms, transit compartment absorption . The induction of IGF-I was modelled using an indirect response model with stimulation of kin with an Emax relationship to rhGH dose . Direct scaling form the PKPD model provided good predictions of the human PKPD except for subcutaneous bioavailability were the approaches failed . Conclusion: A PKPD model of pre-clinical animal data was developed . The model provides a quantitative way of comparing results from testing of new GH analogues in the hypophysectomised rat model to current treatment . Furthermore, direct scaling of parameters matched or improved results obtained when using multispecies scaling to predict human parameters .

Quantitation of replication-competent HIV-1 in populations of resting CD4+ T cells.

Central memory (TCM) CD4(+) T cells are the principal reservoir of latent HIV-1 infection that persists despite durable, successful antiretroviral therapy (ART). In a study that measured HIV DNA in 17 patients and replication-competent HIV in 4 patients, pools of resting and activated transitional memory (T(TM)) CD4(+) T cells were found to be a reservoir for HIV infection. As defective viruses account for the majority of integrated HIV DNA and do not reflect the actual frequency of latent, replication-competent proviral infection, we assessed the specific contribution of resting T(TM) cells to latent HIV infection. We measured the frequency of replication-competent HIV in purified resting memory cell subpopulations by a limiting-dilution, quantitative viral outgrowth assay (QVOA). HIV was routinely detected within the resting central memory compartment but was infrequently detected within the resting T(TM) compartment. These observations suggest that prolonged ART may limit persistent latent infection in the T(TM) compartment. Our results confirm the importance of latent infection within the TCM compartment and again focus attention on these cells as the most important latent viral reservoir. While proliferation may drive expansion of detectable viral genomes in cells, the frequency of replication-competent HIV must be carefully assessed. Latent infection appears to wane within the transitional memory compartment in patients who have sustained successful viral suppression via ART or were treated very early in infection. Antiretroviral therapy (ART) has led to a significant decrease in morbidity and mortality among HIV-infected patients. However, HIV integrates into the genome of CD4(+) T cells, generating pools of long-lived cells that are reservoirs of latent HIV. Two main subsets of CD4(+) T cells, central memory and transitional memory cells, were reported to be major reservoirs of HIV infection. However, this study primarily measured the HIV DNA content, which also includes defective proviruses that would not be able to replicate and initiate new rounds of infection. By analyzing the replication-competent virus in both cell subsets, we showed that transitional memory cells may not be a durable reservoir in patients on successful ART.

Population pharmacokinetics analysis of vigabatrin in adults and children with epilepsy and children with infantile spasms.

Vigabatrin is an inhibitor of γ-aminobutyric acid transaminase. The purpose of these analyses was to develop a population pharmacokinetics model to characterize the vigabatrin concentration-time profile for adults and children with refractory complex partial seizures (rCPS) and for children with infantile spasms (IS); to identify covariates that affect its disposition, and to enable predictions of systemic vigabatrin exposure for patients 1-12months of age. Vigabatrin pharmacokinetic data from six randomized controlled clinical trials and one open-label study were analyzed using nonlinear mixed-effects modeling. Data collected from 349 adults with rCPS and 119 pediatric patients with rCPS or IS were used in the analyses. A two-compartment model with first-order elimination and transit-compartment absorption consisting of five transit compartments adequately described the vigabatrin concentration-time data for these adult and pediatric patient populations. An exponential error model was used to estimate inter-individual variability for the transit-rate constant (k tr) (24.2%), elimination rate constant (k) (14.7%) and apparent central volume of distribution (V c/F) (18%). For the study of children with IS, inter-occasion variability was estimated for k tr (58.1%) and relative bioavailability (F) (26.9%). Covariate analysis indicated that age, creatinine clearance (CLCR), and body weight were important predictors of vigabatrin pharmacokinetic parameters. Vigabatrin apparent clearance increased with increasing CLCR, consistent with renal excretion (primary pathway of vigabatrin elimination). Rate of vigabatrin absorption was dependent on age. The rate was slower in younger patients, which resulted in a smaller predicted maximum concentration and longer predicted time to maximum concentrations. Vigabatrin V c/F, apparent inter-compartmental clearance between the central and peripheral compartment, and apparent peripheral volume of distribution were increased with greater patient body weights. Sex did not contribute significantly to vigabatrin pharmacokinetic variability. The model adequately described vigabatrin pharmacokinetic and enabled predictions of systemic exposures in pediatric patients 1-12months of age.

Population pharmacokinetics and pharmacodynamics of BYL719, a phosphoinositide 3‐kinase antagonist, in adult patients with advanced solid malignancies

The aim was to characterize the population pharmacokinetics of BYL719 in cancer patients and assess the time course of tumour response in relation to drug exposure and dosing schedule. Plasma samples and longitudinal tumour size measurements were collected from 60 patients with advanced solid malignancies who received oral BYL719 once daily (30-450 mg) or twice daily at 120 mg or 200 mg. Non-linear mixed effect modelling was employed to develop the population pharmacokinetic and pharmacodynamic model. The pharmacokinetics were best described by a one compartment disposition model and transit compartments accounting for the lag time in absorption. The typical population oral clearance and volume of distribution estimates with their between-subject variability (BSV) were 10 l h(-1) (BSV 26%) and 108 l (BSV 28%), respectively. The estimated optimal number of transit compartments was 8.1, with a mean transit time to the absorption compartment of 1.28 h (BSV 32%). The between-occasion variability in the rate and extent of absorption was 46% and 26%, respectively. Tumour growth was modelled using a turnover model characterized by a zero order growth rate of 0.581 cm week(1) and a first order death rate of 0.0123 week(-1) . BYL719 inhibited tumour growth with an IC50 of 100 ng ml(-1) (BSV 154%). Model-based predictions showed potential for additional anti-tumour activity of twice daily dosing at total daily dose below 400 mg, but a loss of efficacy if administered less frequently than once daily. The proposed model provides a valuable approach for planning future clinical studies and for designing optimized dosing regimens with BYL719.

Prediction of brain clozapine and norclozapine concentrations in humans from a scaled pharmacokinetic model for rat brain and plasma pharmacokinetics.

BackgroundClozapine is highly effective in treatment-resistant schizophrenia, although, there remains significant variability in the response to this drug. To better understand this variability, the objective of this study was to predict brain extracellular fluid (ECF) concentrations and receptor occupancy of clozapine and norclozapine in human central nervous system by translating plasma and brain ECF pharmacokinetic (PK) relationships in the rat and coupling these with known human disposition of clozapine in the plasma.MethodsUnbound concentrations of clozapine and norclozapine were measured in rat brain ECF using quantitative microdialysis after subcutaneous administration of a 10 mg/kg single dose of clozapine or norclozapine. These data were linked with plasma concentrations obtained in the same rats to develop a plasma–brain ECF compartmental model. Parameters describing brain ECF disposition were then allometrically scaled and linked with published human plasma PK to predict human ECF concentrations. Subsequently, prediction of human receptor occupancy at several CNS receptors was based on an effect model that related the predicted ECF concentrations to published concentration-driven receptor occupancy parameters.ResultsA one compartment model with first order absorption and elimination best described clozapine and norclozapine plasma concentrations in rats. A delay in the transfer of clozapine and norclozapine from plasma to the brain ECF compartment was captured using a transit compartment model approach. Human clozapine and norclozapine concentrations in brain ECF were simulated, and from these the median percentage of receptor occupancy of dopamine-2, serotonin-2A, muscarinic-1, alpha-1 adrenergic, alpha-2 adrenergic and histamine-1 for clozapine, and dopamine-2 for norclozapine were consistent with values reported in the literature.ConclusionsA PK model that relates clozapine and norclozapine disposition in rat plasma and brain, including blood–brain barrier transport, was developed. Using allometry and published human plasma PK, the model was successfully translated to predict clozapine and norclozapine concentrations and accordant receptor occupancy of both agents in human brain. These predicted exposure and occupancy measures at several receptors that bind clozapine may be employed to extend our understanding of clozapine’s complex behavioral effects in humans.Electronic supplementary materialThe online version of this article (doi:10.1186/1479-5876-12-203) contains supplementary material, which is available to authorized users.

Modeling of delays in PKPD: classical approaches and a tutorial for delay differential equations.

In pharmacokinetics/pharmacodynamics (PKPD) the measured response is often delayed relative to drug administration, individuals in a population have a certain lifespan until they maturate or the change of biomarkers does not immediately affects the primary endpoint. The classical approach in PKPD is to apply transit compartment models (TCM) based on ordinary differential equations to handle such delays. However, an alternative approach to deal with delays are delay differential equations (DDE). DDEs feature additional flexibility and properties, realize more complex dynamics and can complementary be used together with TCMs. We introduce several delay based PKPD models and investigate mathematical properties of general DDE based models, which serve as subunits in order to build larger PKPD models. Finally, we review current PKPD software with respect to the implementation of DDEs for PKPD analysis.

Mechanism-based pharmacokinetic/pharmacodynamic meta-analysis of navitoclax (ABT-263) induced thrombocytopenia.

Navitoclax is a first-in-class, orally bioavailable, targeted Bcl-2 family protein inhibitor and promotes apoptosis. Thrombocytopenia is a primary dose-limiting toxicity of navitoclax which exhibited a distinct time profile in circulating platelets from that caused by traditional chemotherapies. A population pharmacokinetic/pharmacodynamic (PK/PD) model was developed to describe the pharmacokinetic of navitoclax as well as the time course of the platelet counts in cancer patients receiving navitoclax. Data from 256 patients who received oral navitoclax (dose range 10-475 mg) as a 14/21-day schedule or a continuous once daily (QD) schedule were used to construct the model using NONMEM. The PK model was a two-compartmental model with a lag-time and a transit compartment in absorption. The PD model was a semi-physiological model that comprised a progenitor cell compartment, three transition compartments representing the maturation chain in the bone marrow and a peripheral blood compartment. Compared with the previously published models, the model established in this analysis applied a different feedback mechanism and introduced a new concept of progenitor cell "pool", which describes a large pool of platelet progenitor cells at the beginning of navitoclax treatment. The PD model was able to describe a slight downward trend of platelet counts over the long-term navitoclax treatment as observed in around 8 % of the patients and the initial drop in platelets seen in our Phase 1/2a studies. We have developed a new semi-physiological platelet model for describing fast drop of platelets after initial navitoclax administration and long-term decline of platelets after continuous administration of navitoclax.

Semi-mechanistic pharmacokinetic/pharmacodynamic modeling of the antitumor activity of LY2835219, a new cyclin-dependent kinase 4/6 inhibitor, in mice bearing human tumor xenografts.

Selective inhibition of cyclin-dependent kinases 4 and 6 (CDK4/6) represents a promising therapeutic strategy. However, despite documented evidence of clinical activity, limited information is available on the optimal dosing strategy of CDK4/6 inhibitors. Here, we present an integrated semi-mechanistic pharmacokinetic/pharmacodynamic model to characterize the quantitative pharmacology of LY2835219, a CDK4/6 inhibitor, in xenograft tumors. LY2835219 plasma concentrations were connected to CDK4/6 inhibition and cell-cycle arrest in colo-205 human colorectal xenografts by incorporating the biomarkers, phospho-(ser780)-Rb, topoisomerase II α, and phosphohistone H3, into a precursor-dependent transit compartment model. This biomarker model was then connected to tumor growth inhibition (TGI) by: (i) relating the rate of tumor growth to mitotic cell density, and (ii) incorporating a concentration-dependent mixed cytostatic/cytotoxic effect driving quiescence and cell death at high doses. Model validation was evaluated by predicting LY2835219-mediated antitumor effect in A375 human melanoma xenografts. The model successfully described LY2835219-mediated CDK4/6 inhibition, cell-cycle arrest, and TGI in colo-205, and was validated in A375. The model also demonstrated that a chronic dosing strategy achieving minimum steady-state trough plasma concentrations of 200 ng/mL is required to maintain durable cell-cycle arrest. Quiescence and cell death can be induced by further increasing LY2835219 plasma concentrations. Our model provides mechanistic insight into the quantitative pharmacology of LY2835219 and supports the therapeutic dose and chronic dosing strategy currently adopted in clinical studies.

A Transit Compartment Model Unmasks OxyContin's Reflective Pharmacokinetics From Urine Measurements in Humans

ABSTRACTThe absorption pattern of orally administered OxyContin (OXC) reflected in urine indicates that its appearance into systemic circulation undergoes transit absorption delays. The authors developed an OXC transit-delay compartment model that identified a new source of oxycodone hydrochloride (OC): the rate of appearance of OC due to OXC tablet dissolution in transit through the gastrointestinal (GI) tract (RGIa), which is due to disintegration of OXC's AcroContin delivery system. RGIa is independent of the biphasic dissolution and release of OC from the delivery system. The authors conclude that an OXC transit-delay compartment model can be of value in the interpretation of OXC pharmacokinetics.

Population pharmacokinetic modelling to assess clinical drug-drug interaction between AZD7325 and midazolam

AZD7325 is a selective gamma-amino-butyric acid (GABAA )α2, 3 receptor modulator. The aims of this analysis were to develop population pharmacokinetic (PPK) models of AZD7325 and midazolam and to assess the induction effect of AZD7325 on CYP3A4 with midazolam as a substrate. Drug-drug interaction data of AZD7325 and midazolam from 24 healthy subjects were available for model development. PPK models were developed in a sequential manner using NONMEM. Both AZD7325 and midazolam pharmacokinetics were described by two-compartment models, and a transit compartment absorption model and a first-order absorption model were applied for the absorption of AZD7325 and midazolam, respectively. The induction of CYP3A by AZD7325 was described by a transit enzyme model, where the elimination of midazolam was proportionally linked to the enzyme amount. Simulations were performed to predict dosing regimens to account for the induction of CYP3A4. The population estimates for AZD7325 clearance, intercompartmental clearance, central and peripheral volume were 36, 29·2L/h, 169 and 392L, respectively, with interindividual variability (IIV) of 35% and 24% for clearance and central volume. Midazolam clearance, intercompartmental clearance, central and peripheral volume were estimated to be 62·7, 34·7L/h, 133 and 146L, respectively, with 43% IIV for clearance. The estimated mean transit time for induction of the CYP3A4 enzyme was 197h, with 57% IIV. The PPK models developed adequately described the clinical observation of AZD7325-mediated CYP3A4 enzyme induction with midazolam as a probe. The model could provide basis for the rational dosing of AZD7325 in clinical practice.

Semi-physiological pharmacokinetic–pharmacodynamic modeling and simulation of 5-fluorouracil for the whole time course of alterations in leukocyte, neutrophil and lymphocyte counts in rats

1. We aimed to develop a simple pharmacokinetic–pharmacodynamic (PK–PD) model to predict the onset and degree of severe toxic side effects that severely limit the use of many anticancer agents, such as myelosuppression, in rats.2. Our PK–PD model consisted of a two-compartment PK model, with one compartment representing proliferative cells and some transit compartments consisting of maturing cells, while the other compartment represented circulating blood cells for the PD model.3. The semi-physiological PK–PD model effectively captured the features of myelosuppression and the degree of the off-target toxicities observed after 5-fluorouracil (5-FU) chemotherapy, and helped simultaneously simulate the whole time course for alterations in leukocyte, neutrophil and lymphocyte counts after 5-FU treatment in rats. Interestingly, by plotting the nadir period of leukocyte, neutrophil and lymphocyte counts as determined by PK–PD analytical simulation curves against the area under the plasma 5-FU concentration–time curve (AUC0–∞) after intravenous administration of 5-FU, a linear relationship was inferred, with r2 = 0.989, 0.877 and 0.956, respectively.4. The semi-physiological PK–PD model is a valuable tool for evaluating a variety of novel cancer chemopreventive agents or emerging therapeutic strategies that are difficult to address in humans.

Population Pharmacokinetics of Rifapentine and Desacetyl Rifapentine in Healthy Volunteers: Nonlinearities in Clearance and Bioavailability

Rifapentine is under active investigation as a potent drug that may help shorten the tuberculosis (TB) treatment duration. A previous rifapentine dose escalation study with daily dosing indicated a possible decrease in bioavailability as the dose increased and an increase in clearance over time for rifapentine and its active metabolite, desacetyl rifapentine. This study aimed to assess the effects of increasing doses on rifapentine absorption and bioavailability and to evaluate the clearance changes over 14 days. A population analysis was performed with nonlinear mixed-effects modeling. Absorption, time-varying clearance, bioavailability, and empirical and semimechanistic autoinduction models were investigated. A one-compartment model linked to a transit compartment absorption model best described the data. The bioavailability of rifapentine decreased linearly by 2.5% for each 100-mg increase in dose. The autoinduction model suggested a dose-independent linear increase in clearance of the parent drug and metabolite over time from 1.2 and 3.1 liters · h(-1), respectively, after a single dose to 2.2 and 5.0 liters · h(-1), respectively, after 14 once-daily doses, with no plateau being reached by day 14. In clinical trial simulations using the final model, rifapentine demonstrated less-than-dose-proportional pharmacokinetics, but there was no plateau in exposures over the dose range tested (450 to 1,800 mg), and divided dosing increased exposures significantly. Thus, the proposed compartmental model incorporating daily dosing of rifapentine over a wide range of doses and time-related changes in bioavailability and clearance provides a useful tool for estimation of drug exposure that can be used to optimize rifapentine dosing for TB treatment. (This study has been registered at ClinicalTrials.gov under registration no. NCT01162486.).

Pharmacokinetics, in vitro and in vivo correlation, and efficacy of exenatide microspheres in diabetic rats

Modeling and simulation are aimed at achieving information about the behaviors of the drugs without the actual measurements and determination. The purpose of this study was to characterize the in vivo behavior of exenatide microspheres using model-based methods. Exenatide is a glucagon-like peptide-1 agonist medication, belonging to the group of incretin mimetics, approved for the treatment of diabetes mellitus type 2. An oil-in-water solvent evaporation method was used to prepare the exenatide microspheres and their physicochemical features were investigated. After subcutaneous injection of exenatide microspheres to streptozotocin-induced diabetic rats, the exenatide concentrations increased and kept increasing and the blood glucose decreased in all diabetic rats. The in vivo release behavior of exenatide from microspheres was described by a transit compartment model. Based on the transit compartment model, the simulation method was proposed for the description of in vivo release. The in vitro and in vivo correlation (IVIVC) was established by the model-based simulation (R2 = 0.903) and deconvolution (R2 = 0.922) methods successfully. Using a transit compartment model to characterize the in vivo exenatide release from microspheres is an acceptable approach, and the IVIVC can be estimated reliably with the model-based simulation method.

Population Pharmacokinetics and Prediction of Dopamine D2 Receptor Occupancy After Multiple Doses of RBP‐7000, a New Sustained‐Release Formulation of Risperidone, in Schizophrenia Patients on Stable Oral Risperidone Treatment

RBP-7000 is a long-acting formulation of risperidone administered once monthly via subcutaneous (SC) injections for the treatment of schizophrenia. The objectives of the present study were to characterize the pharmacokinetics of RBP-7000 after multiple doses in schizophrenic patients on stable oral risperidone therapy and to evaluate the switch between oral risperidone and SC injections of RBP-7000. Data were collected in a phase IIa, open-label, multiple-ascending-dose study where 45 patients clinically stabilized on oral risperidone (2, 3 or 4 mg/day) were switched to receive 60, 90 or 120 mg/month SC injections of RBP-7000, respectively. Patients were thereafter switched back to oral risperidone. An integrated population pharmacokinetic model describing simultaneously risperidone and 9-hydroxyrisperidone after risperidone oral intake and RBP-7000 administration was developed in NON-MEM using 5,232 quantifiable plasma concentrations. Predictions of dopamine D2 receptor occupancy were derived using a previously published model. A two-compartment model with first-order absorption was selected for oral risperidone, while a three-compartment model with first-order absorption and a transit compartment absorption model was selected for RBP-7000. Body mass index was identified as a significant covariate affecting the initial absorption of risperidone following RBP-7000 injection. Steady state was reached after the second or third RBP-7000 injection but plasma concentrations close to steady-state values were obtained right after the first injection when switching from oral risperidone therapy. Predicted dopamine D2 receptor occupancy after repeated doses of 90 and 120 mg showed less fluctuation than after oral risperidone with acceptable ranges for clinical efficacy and a potentially safer profile with respect to extrapyramidal side effects. This analysis provided additional insight into the pharmacokinetics of RBP-7000 and for the comparison with oral risperidone treatment. The established model was used to support the design of a planned phase III study.

Caecal pH is a biomarker of excessive colonic fermentation.

To ascertain whether caecal pH is different in patients with irritable bowel syndrome (IBS), whose primary symptoms are bloating and distension, to healthy controls. Motility and pH data were reviewed from 16 patients with Rome III defined IBS and 16 healthy controls, who had undergone a wireless motility capsule (WMC) study using a standardized protocol. Motility measures were anchored around known anatomical landmarks as identified by compartmental pH changes. Sixty-minute epochs were used to quantify antral, duodenal, ileal, caecal and distal colonic contractility. The maximum and minimum pH was measured either side of the ileo-caecal junction. No differences were seen in motility parameters, compartmental transit times or maximal ileal pH between the two groups. Caecal pH was significantly lower in patients compared to controls (5.12 ± 0.05 vs 6.16 ± 0.15, P < 0.0001). The ileal:caecal Δchange was greater in patients than controls (-2.63 ± 0.08 vs -1.42 ± 0.11, P < 0.0001). There was a significant correlation between caecal pH and right colonic contractility (r = 0.54, P = 0.002). Patients with bloating and distension have a lower caecal pH compared to controls. The measurement of caecal pH using the WMC provides a quantifiable biomarker of fermentation potentially identifying those patients that may preferentially benefit from antibiotic or dietary interventions.

Feedback modeling of non-esterified fatty acids in obese Zucker rats after nicotinic acid infusions

This study investigates the impact of disease on nicotinic acid (NiAc)-induced changes in plasma concentrations of non-esterified fatty acids (NEFA). NiAc was given by constant intravenous infusion to normal Sprague-Dawley and obese Zucker rats, and arterial blood samples were taken for analysis of NiAc, NEFA, insulin and glucose plasma concentrations. The intravenous route was intentionally selected to avoid confounding processes, such as absorption, following extravascular administration. Data were analyzed using nonlinear mixed effects modeling (NONMEM, version VI). The disposition of NiAc in the normal rats was described by a two-compartment model with endogenous synthesis of NiAc and two parallel capacity-limited elimination processes. In the obese rats disposition was described by a one-compartment model with endogenous synthesis of NiAc and one capacity-limited elimination process. The plasma concentration of NiAc drove NEFA (R) turnover via an inhibitory drug-mechanism function acting on the formation of NEFA. NEFA turnover was described by a feedback model with a moderator distributed over a series of transit compartments, where the first compartment (M 1 ) inhibited the formation of R and the last compartment (M N ) stimulated the loss of R. All processes regulating plasma NEFA concentrations were assumed to be captured by the moderator function. Differences in the pharmacodynamic response of the two strains included, in the obese animals, an increased NEFA baseline, diminished rebound and post-rebound oscillation, and a more pronounced slowly developing tolerance during the period of constant drug exposure. The feedback model captured the NiAc-induced changes in NEFA response in both the normal and obese rats. Differences in the parameter estimates between the obese and normal rats included, in the former group, increases in R 0 , k in and p by 44, 41 and 78 %, respectively, and decreases in k out and γ by 64 and 84 %, respectively. The estimates of k tol and IC 50 were similar in both groups. The NiAc-NEFA concentration-response relationship at equilibrium was substantially different in the two groups, being shifted upwards and to the right, and being shallower in the obese rats. The extent of such shifts is important, as they demonstrate the impact of disease at equilibrium and, if ignored, will lead to erroneous dose predictions and, in consequence, poorly designed studies. The proposed models are primarily aimed at screening and selecting candidates with the highest potential of becoming a viable drug in man.

Saturable human neopterin response to interferon-α assessed by a pharmacokinetic-pharmacodynamic model

BackgroundIn this study, we developed a pharmacokinetic (PK)- pharmacodynamic (PD) model of a new sustained release formulation of interferon-α-2a (SR-IFN-α) using the blood concentration of IFN-α and neopterin in order to quantify the magnitude and saturation of neopterin production over time in healthy volunteers. The SR-IFN-α in this study is a solid microparticular formulation manufactured by spray drying of a feeding solution containing IFN-α, a biocompatible polymer (polyethylene glycol) and sodium hyaluronate.MethodsThe full PK and PD (neopterin concentration) datasets from 24 healthy subjects obtained after single doses of 9, 18, 27 and 36 MIU of subcutaneous SR-IFN-α were used to build the mixed-effect model using NONMEM (version 7.2) with the GFORTRAN compiler.ResultsA one-compartment model with first-order elimination and a mixture of zero- and first-order absorption was chosen to describe the PK of SR-IFN-α. The time-concentration profile of neopterin, the PD marker, was described by a turnover model combined with a single transit compartment. The saturable pattern of the neopterin response blurring the dose–response relationship of SR-IFN-α was addressed by introducing the concept of the EC50 increasing over time.ConclusionsThe PK-PD model of SR-IFN-α developed in this study has presented a quantitative tool to assess the time-course of a saturable neopterin response in humans.

Discovering dominant pathways and signal–response relationships in signaling networks through nonparametric approaches

A signaling pathway is a sequence of proteins and passenger molecules that transmits information from the cell surface to target molecules. Understanding signal transduction process requires detailed description of the involved pathways. Several methods and tools resolved this problem by incorporating genomic and proteomic data. However, the difficulty of obtaining prior knowledge of complex signaling networks limited the applicability of these tools. In this study, based on the simulation of signal flow in signaling network, we introduce a method for determining dominant pathways and signal response to stimulations. The model uses topology-weighted transit compartment approach and comprises four main steps which include weighting the edges, simulating signal transduction in the network (weighting the nodes), finding paths between initial and target nodes, and assigning a significance score to each path. We applied the proposed model to eighty-three signaling networks by using biologically derived source and sink molecules. The recovered dominant paths matched many known signaling pathways and suggesting a promising index to analyze the phenotype essentiality of molecule encoding paths. We also modeled the stimulus–response relations in long and short-term synaptic plasticity based on the dominant signaling pathway concept. We showed that the proposed method not only accurately determines dominant signaling pathways, but also identifies effective points of intervention in signal transduction.

Population Pharmacodynamic Modeling of Hyperglycemic Clamp and Meal Tolerance Tests in Patients with Type 2 Diabetes Mellitus

In this study, glucose and insulin concentration-time profiles in subjects with type 2 diabetes mellitus (T2DM) under meal tolerance test (MTT) and hyperglycemic clamp (HGC) conditions were co-modeled simultaneously. Blood glucose and insulin concentrations were obtained from 20 subjects enrolled in a double-blind, placebo-controlled, randomized, two-way crossover study. Patients were treated with palosuran or placebo twice daily for 4 weeks and then switched to the alternative treatment after a 4-week washout period. The MTT and HGC tests were performed 1 h after drug administration on days 28 and 29 of each treatment period. Population data analysis was performed using NONMEM. The HGC model incorporates insulin-dependent glucose clearance and glucose-induced insulin secretion. This model was extended for the MTT, in which glucose absorption was described using a transit compartment with a mean transit time of 62.5 min. The incretin effect (insulin secretion triggered by oral glucose intake) was also included, but palosuran did not influence insulin secretion or sensitivity. Glucose clearance was 0.164 L/min with intersubject and interoccasion variability of 9.57% and 31.8%. Insulin-dependent glucose clearance for the HGC was about 3-fold greater than for the MTT (0.0111 vs. 0.00425 L/min/[mU/L]). The maximal incretin effect was estimated to enhance insulin secretion 2-fold. The lack of palosuran effect coupled with a population-based analysis provided quantitative insights into the variability of glucose and insulin regulation in patients with T2DM following multiple glucose tolerance tests. Application of these models may also prove useful in antihyperglycemic drug development and assessing glucose-insulin homeostasis.