A pharmacokinetic-pharmacodynamic assessment of oral antibiotics for pyelonephritis.

Antibiotic resistance to oral antibiotics recommended for pyelonephritis is increasing. The objective was to determine if there is a pharmacological basis to consider alternative treatments/novel dosing regimens for the oral treatment of pyelonephritis. A systematic review identified pharmacokinetic models of suitable quality for a selection of antibiotics with activity against Escherichia coli. MIC data was obtained for a population of E. coli isolates derived from patients with pyelonephritis. Pharmacokinetic/pharmacodynamic (PK/PD) simulations determined probability of target attainment (PTA) and cumulative fraction response (CFR) values for sub-populations of the E. coli population at varying doses. There are limited high-quality models available for the agents investigated. Pharmacokinetic models of sufficient quality for simulation were identified for amoxicillin, amoxicillin-clavulanic acid, cephalexin, ciprofloxacin, and fosfomycin trometamol. These antibiotics were predicted to have PTAs ≥ 0.85 at or below standard doses for the tested E. coli population including cephalexin 1500 mg 8 hourly for 22% of the population (MIC ≤ 4 mg/L) and ciprofloxacin 100 mg 12 hourly for 71% of the population (MIC ≤ 0.06 mg/L). For EUCAST-susceptible E. coli isolates, doses achieving CFRs ≥ 0.9 included amoxicillin 2500 mg 8 hourly, cephalexin 4000 mg 6 hourly, ciprofloxacin 200 mg 12 hourly, and 3000 mg of fosfomycin 24 hourly. Limitations in the PK data support carrying out additional PK studies in populations of interest. Oral antibiotics including amoxicillin, amoxicillin-clavulanic acid, and cephalexin have potential to be effective for a proportion of patients with pyelonephritis. Ciprofloxacin may be effective at lower doses than currently prescribed.

Semi-mechanistic Pharmacokinetic/Pharmacodynamic model of three pegylated rHuEPO and ior®EPOCIM in New Zealand rabbits.

Marketed formulations of erythropoietin (EPO) ior®EPOCIM, MIRCERA® and two newly developed pegylated-EPO analogues (PEG-EPO 32 and 40 kDa) formulations were intravenously administered to New Zealand rabbits. A semi-mechanistic Pharmacokinetic/Pharmacodynamic (PK/PD) model describing in a simultaneous and integrated form the time course of reticulocytes, red blood cells and hemoglobin was built to account for the time course of hematopoiesis stimulation after erythropoietin administration. Data analysis was performed based on the population approach with the software NONMEM version 7.3. Erythropoietin disposition of each of the administered formulations was best described with a two compartment model and linear elimination. Different formulations show different clearance and apparent volume of distribution of the central compartment but share estimates of inter-compartmental clearance and apparent peripheral volume of distribution. A semi-mechanistic model including cell proliferation, maturation, and homeostatic regulation provided a good description of the data regardless the type of erythropoietin formulation administered. The system-, and drug-related parameters showed consistency and differed across formulations, respectively. A single IV administration of PEG-EPO 32 and 40 kDa formulations in New Zealand rabbits achieves a median change of 27% and 22% on RET levels, and of 47% and 63% on RBC and HGB levels, respectively compared to MIRCERA®. The administration of new branched PEG-chains formulations improves PK and PD properties of EPO, in terms of increasing elimination half-lives and pharmacological activity on RET, RBC and HGB compared to commercially available formulations (ior®EPOCIM and MIRCERA®).

Commentary on Pharmacometrics for Immunotherapy.

This commentary provides an overview of recent examples of pharmacometrics applied during the clinical development of two antagonists of the programmed death-1 (PD-1) cell surface receptor, pembrolizumab and nivolumab. Despite the remarkable achievements obtained in predicting the correct dosing schedule from different quantitative approaches, data indicated a great degree of heterogeneity in tumor response. To achieve therapeutic goals the search for predictive biomarkers associated with a lack of response and mechanism-based combination studies are warranted.

Relation between acute and long-term cognitive decline after surgery: Influence of metabolic syndrome.

INTRODUCTION: The relationship between persistent postoperative cognitive decline and the more common acute variety remains unknown; using data acquired in preclinical studies of postoperative cognitive decline we attempted to characterize this relationship. METHODS: Low capacity runner (LCR) rats, which have all the features of the metabolic syndrome, were compared postoperatively with high capacity runner (HCR) rats for memory, assessed by trace fear conditioning (TFC) on the 7th postoperative day, and learning and memory (probe trial [PT]) assessed by the Morris water-maze (MWM) at 3 months postoperatively. Rate of learning (AL) data from the MWM test, were estimated by non-linear mixed effects modeling. The individual rat's TFC result at postoperative day (POD) 7 was correlated with its AL and PT from the MWM data sets at postoperative day POD 90. RESULTS: A single exponential decay model best described AL in the MWM with LCR and surgery (LCR-SURG) being the only significant covariates; first order AL rate constant was 0.07 s(-1) in LCR-SURG and 0.16s(-1) in the remaining groups (p<0.05). TFC was significantly correlated with both AL (R=0.74; p<0.0001) and PT (R=0.49; p<0.01). CONCLUSION: Severity of memory decline at 1 week after surgery presaged long-lasting deteriorations in learning and memory.

Pharmacometrics Markup Language (PharmML): Opening New Perspectives for Model Exchange in Drug Development.

The lack of a common exchange format for mathematical models in pharmacometrics has been a long-standing problem. Such a format has the potential to increase productivity and analysis quality, simplify the handling of complex workflows, ensure reproducibility of research, and facilitate the reuse of existing model resources. Pharmacometrics Markup Language (PharmML), currently under development by the Drug Disease Model Resources (DDMoRe) consortium, is intended to become an exchange standard in pharmacometrics by providing means to encode models, trial designs, and modeling steps.

Optimizing oncology therapeutics through quantitative translational and clinical pharmacology: challenges and opportunities.

Despite advances in biomedical research that have deepened our understanding of cancer hallmarks, resulting in the discovery and development of targeted therapies, the success rates of oncology drug development remain low. Opportunities remain for objective dose selection informed by exposure-response understanding to optimize the benefit-risk balance of novel therapies for cancer patients. This review article discusses the principles and applications of modeling and simulation approaches across the lifecycle of development of oncology therapeutics. Illustrative examples are used to convey the value gained from integration of quantitative clinical pharmacology strategies from the preclinical-translational phase through confirmatory clinical evaluation of efficacy and safety.

A review of mixed-effects models of tumor growth and effects of anticancer drug treatment used in population analysis.

Population modeling of tumor size dynamics has recently emerged as an important tool in pharmacometric research. A series of new mixed-effects models have been reported recently, and we present herein a synthetic view of models with published mathematical equations aimed at describing the dynamics of tumor size in cancer patients following anticancer drug treatment. This selection of models will constitute the basis for the Drug Disease Model Resources (DDMoRe) repository for models on oncology.

Semi-mechanistic modelling of the analgesic effect of gabapentin in the formalin-induced rat model of experimental pain.

PURPOSE: The formalin-induced rat model of nociception involves moderate continuous pain. Formalin-induced pain results in a typical repetitive flinching behaviour, which displays a biphasic pattern characterised by peaks of pain. Here we described the time course of pain response and the analgesic effect of gabapentin using a semi-mechanistic modelling approach. METHODS: Male Sprague-Dawley rats received gabapentin (10-100 mg/kg) or placebo 1 h prior to the formalin injection, as per standard protocol. A reduction in the frequency of the second peak of flinching was used as a behavioural measure of gabapentin-mediated anti-nociception. The flinching response was modelled using a mono-exponential function to characterise the first peak and an indirect response model with a time variant synthesis rate for the second. PKPD modelling was performed using a population approach in NONMEM v.7.1.2. RESULTS: The time course of the biphasic response was adequately described by the proposed model, which included separate expressions for each phase. Gabapentin was found to reversibly decrease, but not suppress the flinching frequency of the second response peak only. The mean IC50 estimate was 7,510 ng/ml, with relative standard error (RSE%) of 40%. CONCLUSIONS: A compartmental, semi-mechanistic model provides the basis for further understanding of the formalin-induced flinching response and consequently to better characterisation of the properties of gabapentin, such as the potency in individual animals. Moreover, despite high exposure levels, model predictions show that gabapentin does not completely suppress behavioural response in the formalin-induced pain model.

Modeling the effect of propofol and remifentanil combinations for sedation-analgesia in endoscopic procedures using an Adaptive Neuro Fuzzy Inference System (ANFIS).

BACKGROUND: The increasing demand for anesthetic procedures in the gastrointestinal endoscopy area has not been followed by a similar increase in the methods to provide and control sedation and analgesia for these patients. In this study, we evaluated different combinations of propofol and remifentanil, administered through a target-controlled infusion system, to estimate the optimal concentrations as well as the best way to control the sedative effects induced by the combinations of drugs in patients undergoing ultrasonographic endoscopy. METHODS: One hundred twenty patients undergoing ultrasonographic endoscopy were randomized to receive, by means of a target-controlled infusion system, a fixed effect-site concentration of either propofol or remifentanil of 8 different possible concentrations, allowing adjustment of the concentrations of the other drug. Predicted effect-site propofol (C(e)pro) and remifentanil (C(e)remi) concentrations, parameters derived from auditory evoked potential, autoregressive auditory evoked potential index (AAI/2) and electroencephalogram (bispectral index [BIS] and index of consciousness [IoC]) signals, as well as categorical scores of sedation (Ramsay Sedation Scale [RSS] score) in the presence or absence of nociceptive stimulation, were collected, recorded, and analyzed using an Adaptive Neuro Fuzzy Inference System. The models described for the relationship between C(e)pro and C(e)remi versus AAI/2, BIS, and IoC were diagnosed for inaccuracy using median absolute performance error (MDAPE) and median root mean squared error (MDRMSE), and for bias using median performance error (MDPE). The models were validated in a prospective group of 68 new patients receiving different combinations of propofol and remifentanil. The predictive ability (P(k)) of AAI/2, BIS, and IoC with respect to the sedation level, RSS score, was also explored. RESULTS: Data from 110 patients were analyzed in the training group. The resulting estimated models had an MDAPE of 32.87, 12.89, and 8.77; an MDRMSE of 17.01, 12.81, and 9.40; and an MDPE of -1.86, 3.97, and 2.21 for AAI/2, BIS, and IoC, respectively, in the absence of stimulation and similar values under stimulation. P(k) values were 0.82, 0.81, and 0.85 for AAI/2, BIS, and IoC, respectively. The model predicted the prospective validation data with an MDAPE of 34.81, 14.78, and 10.25; an MDRMSE of 16.81, 15.91, and 11.81; an MDPE of -8.37, 5.65, and -1.43; and P(k) values of 0.81, 0.8, and 0.8 for AAI/2, BIS, and IoC, respectively. CONCLUSION: A model relating C(e)pro and C(e)remi to AAI/2, BIS, and IoC has been developed and prospectively validated. Based on these models, the (C(e)pro, C(e)remi) concentration pairs that provide an RSS score of 4 range from (1.8 µg·mL(-1), 1.5 ng·mL(-1)) to (2.7 µg·mL(-1), 0 ng·mL(-1)). These concentrations are associated with AAI/2 values of 25 to 30, BIS of 71 to 75, and IoC of 72 to 76. The presence of noxious stimulation increases the requirements of C(e)pro and C(e)remi to achieve the same degree of sedative effects.

Prediction of neutropenia-related effects of a new combination therapy with the anticancer drugs BI 2536 (a Plk1 inhibitor) and pemetrexed.

This study investigated the feasibility of predicting the neutropenia-related effects of a therapy that combines the investigational drug BI 2536 (inhibitor of Polo-like kinase 1) and pemetrexed, an approved anticancer drug. Predictions were arrived at using the pharmacokinetic/pharmacodynamic (PK/PD) parameters of each of the drugs obtained from monotherapy studies and assuming that the neutropenic effect is additive when the drugs are administered as a combination therapy. Subsequently, a PK/PD model was developed to determine whether this assumption of additive effect was reasonable in relation to these two drugs. All analyses and simulations were performed using the population approach in NONMEM, version VI.

Effects of fasting and gender on ochratoxin A toxicokinetics in F344 rats.

Ochratoxin A (OTA) is a mycotoxin that causes renal tumors in rats, particularly in males. In previous kinetic studies performed in fed conditions (Vettorazzi et al., 2008), mature F344 male rats presented a significantly lower OTA bioavailability than females and young animals. The objective of the present study was to evaluate two factors which could explain this different kinetic profile: the presence of food and the male-specific protein alpha-2u-globulin. Therefore, a 24h kinetic study has been performed in rats under fasting conditions. Food ingestion has been controlled in both sexes during two months. The presence of alpha-2u-globulin in the urine has been analyzed with SDS-gradient mini-gel electrophoresis. Fasting tends to increase the maximum OTA plasma concentrations and the rate of absorption. The relative bioavailability is significantly increased under fasting conditions only in males. Mature males consumed a higher amount of food but, as the OTA dose administered, it was proportional to body weight. The reason why the OTA bioavailability is more affected in presence of food only in males is unclear. Several possibilities, such as differences in gastric emptying, OTA-food interactions and the involvement of alpha-2u-globulin are discussed.

Population pharmacokinetics of ganciclovir after intravenous ganciclovir and oral valganciclovir administration in solid organ transplant patients infected with cytomegalovirus.

A population pharmacokinetics analysis was performed after intravenous ganciclovir and oral valganciclovir in solid organ transplant patients with cytomegalovirus. Patients received ganciclovir at 5 mg/kg of body weight (5 days) and then 900 mg of valganciclovir (16 days), both twice daily with dose adjustment for renal function. A total of 382 serum concentrations from days 5 and 15 were analyzed with NONMEM VI. Renal function given by creatinine clearance (CL(CR)) was the most influential covariate in CL. The final pharmacokinetic parameters were as follows: ganciclovir clearance (CL) was 7.49.(CL(CR)/57) liter/h (57 was the mean population value of CL(CR)); the central and peripheral distribution volumes were 31.9 liters and 32.0 liters, respectively; intercompartmental clearance was 10.2 liter/h; the first-order absorption rate constant was 0.895 h(-1); bioavailability was 0.825; and lag time was 0.382 h. The CL(CR) was the best predictor of CL, making dose adjustment by this covariate important to achieve the most efficacious ganciclovir exposure.

A different kinetic profile of ochratoxin A in mature male rats.

Ochratoxin A (OTA) is a mycotoxin that causes renal tumors in rodents, particularly in male rats. The present work explored the impact of gender and age on OTA toxicokinetics in F344 rats after a single oral dose (0.5mg/kg b.w.). OTA plasma concentrations were analysed with a validated HPLC-FLD method and a population approach (NONMEM VI) was used to perform the kinetic analysis and the one year exposure simulation (0.21 mg/kg daily). Maximum observed OTA concentration (CMAX(obs)) was at 2h in all groups except in mature females (6h). Mature females reached higher CMAX(obs) than males of the same age. Apparent volume of distribution, but not apparent total plasma clearance, increased significantly with body weight (P<0.01) resulting in the following values for the terminal plasma half life (h) in males: 219 (young), 264 (matures) and females: 191 (young), 205 (matures). In addition mature males showed a significant lower relative bioavailability. The simulation showed similar plasma concentrations in males and females after two-months. Thus, toxicokinetic does not seem to explain sex-differences in toxicity in long-term studies. However, the age and weight should be taken into account in short-term toxicological studies if sex-differences are studied.

[Use of the electroencephalographic spectral edge frequency 90 to monitor sedation of chronically ill patients]. / Utilización del límite espectral 90% (LE90) del electroencefalograma como sistema de control de la sedación en los pacientes críticos.

INTRODUCTION: The depth of sedation required for patients in critical care units varies over time and should be subject to control. The clinical assessment scales used at present are inadequate, and several electroencephalographic variables have been investigated in recent years with the aim of quantifying depth of sedation. One such variable is the spectral edge frequency 90 (SEF90). OBJECTIVES: To establish the correlation between SEF90 and the Ramsay score as indicators of depth of sedation. To estimate the ability of SEF90 to predict sedation and a patient's hemodynamic response during aspiration of secretions through the orotracheal tube. PATIENTS AND METHODS: Patients in a surgical intensive care unit. The ability of SEF90 to predict a certain Ramsay score was assessed by logistic regression. We also calculated the predictive probability (Pk) of SEF90 for the appearance of hemodynamic change and of movement in the event of endotracheal aspiration. RESULTS: When SEF90 was < 16 Hz, the probability of a Ramsay score >= 4 was >= 90% (Pk = 0.91). Neither SEF90 nor the Ramsay score predicted hemodynamic response to orotracheal aspiration. CONCLUSIONS: SEF90 distinguishes superficial from deep sedation but does not differentiate further degrees of depth or the likelihood of hemodynamic instability or movement in response to aspiration.

Pharmacokinetic/pharmacodynamic modeling of antipyretic and anti-inflammatory effects of naproxen in the rat.

Pharmacokinetic/pharmacodynamic modeling was used to characterize the antipyretic and anti-inflammatory effects of naproxen in rats. An indirect response model was used to describe the antipyretic effects of naproxen after short intravenous infusions. The model assumes that basal temperature (T(a)) is maintained by the balance of fever mediators given by a constant (zero order) rate of synthesis (K(syn)), and a first order rate of degradation (K(out)). After an intraperitoneal injection of lipopolysaccharide, the change in T(a) was modeled assuming an increase in fever mediators described as an input rate function [IR(t)] estimated nonparametrically. An inhibitory E(max) model adequately described the inhibition of IR(t) by naproxen. A more complex model was used to describe the anti-inflammatory response of oral naproxen in the carrageenin-induced edema model. Before carrageenin injection, physiological conditions are maintained by a balance of inflammation mediators given by K(syn) and K(out) (see above). After carrageenin injection, the additional synthesis of mediators is described by IR(t) (see above). Such mediators induced an inflammatory process, which is governed by a first order rate constant (K(IN)) that can be inhibited by the presence of naproxen in plasma. The sigmoidal E(max) model also well described the inhibition of K(IN) by naproxen. Estimates for IC(50) [concentration of naproxen in plasma eliciting half of maximum inhibition of IR(t) or K(IN)] were 4.24 and 4.13 microg/ml, for the antipyretic and anti-inflammatory effects, respectively.

Modeling of the in vivo antinociceptive interaction between an opioid agonist, (+)-O-desmethyltramadol, and a monoamine reuptake inhibitor, (-)-O-desmethyltramadol, in rats.

The pharmacokinetic-pharmacodynamic (pk-pd) characterization of the in vivo antinociceptive interaction between (+)-O-desmethyltramadol [(+)-M1] and (-)-O-desmethyltramadol [(-)-M1], main metabolites of tramadol, was studied in three groups of rats. (+)-M1 and (-)-M1, both with different pd properties, were studied under steady-state and nonsteady-state conditions, depending on the group. Plasma drug concentration and antinociception were simultaneously measured in each animal by using an enantioselective analytical assay and the tail-flick test, respectively. Respiratory depression also was evaluated in another series of experiments according to the same experimental conditions. The pk behavior was similar for both enantiomers and no significant (P >.05) interaction between two compounds was found at this level. However, a significant (P <.01) potentiation in the antinociceptive effect elicited by (+)-M1 was found during and after (-)-M1 administration. The pd model used to describe the time course of the antinociception in the presence of (+)-M1, (-)-M1, or both is based on previous knowledge of the compounds and includes the following: 1) an effect compartment model to account for the opioid effect of (+)-M1, and 2) an indirect response model accounting for the release of noradrenaline (NA) caused by (+)-M1, and the inhibition of the NA reuptake due to the action of (-)-M1. The model predicts a positive contribution to antinociception of the predicted increasing levels of NA. No significant (P >.05) respiratory effects were seen during or after (+)-M1 and (-)-M1 administration.

Comparison of manual versus ambulatory blood pressure measurements with pharmacokinetic-pharmacodynamic modeling of antihypertensive compounds: application to moxonidine.

OBJECTIVES: To compare the results of the pharmacokinetic-pharmacodynamic analyses of 24-hour ambulatory blood pressure measurements and manual blood pressure data in patients receiving moxonidine. METHODS: 32 patients with borderline to mild-to-moderate hypertension were enrolled in a double-blind, placebo-controlled phase II study. After receiving placebo for 1 week (run-in phase), the patients were randomly allocated to the placebo or the 0.6-, 0.9-, or 1.2-mg dose groups. Placebo and moxonidine were administered once daily for 1 week (drug-treatment phase). Four 24-hour ambulatory blood pressure measurement profiles were obtained for each individual. Plasma samples (n = 9) and four measurements of manual blood pressure were taken at the start and end of the drug-treatment phase. Two additional manual blood pressure measurements were taken during the run-in and drug-treatment phases. RESULTS: Pharmacokinetics was described by a one-compartment model. For the 24-hour ambulatory blood pressure measurements, baseline circadian patterns were described with a two-cosine function model that included interindividual and interoccasion variability. Pharmacodynamics was described with use of an effect-compartment model [k(e0) = 0.37 (1/h)] and an Emax model. For diastolic blood pressure the maximum drug-induced decrease (Emax) was 30.9 mm Hg and the steady-state plasma drug concentration eliciting half of maximum effect (C50) was 1.33 microg/L. Interindividual variability was estimated for ke0 (24.8%) and Emax (33.3%). For the manual blood pressure measurements, data was described by a time-invariant baseline model combined with an effect-compartment model and an Emax model. Mean population estimates were in agreement with those obtained during the analysis of 24-hour ambulatory blood pressure measurements. However, interindividual variability could be estimated for the baseline parameter only. CONCLUSIONS: Although similar typical population estimates for the drug action-related parameters were obtained with use of manual blood pressure data and 24-hour ambulatory blood pressure measurements, the latter allowed for a more detailed description of the individual pharmacodynamic profiles because interindividual variability in pharmacodynamic parameters could be estimated together with increased precision in parameter estimates.