Day of surgery: the impact on patient length of stay in a tertiary vascular unit.

BACKGROUND: Our hypothesis was that patients undergoing surgery earlier in the week would have better access to physiotherapy and other discharge services after surgery and, as a result, would have a shorter length of hospital stay compared with patients undergoing surgery later in the week. This study aimed to assess whether there is a significant difference in postoperative length of hospital stay between the groups with secondary assessment by operation subtype. METHODS: We identified all patients admitted for vascular surgery in 2015 from a prospectively collected database and divided the week into Monday to Wednesday and Thursday to Friday. Endovascular cases were included but day cases were excluded. Further analysis was performed with a breakdown in both groups by operation type. Statistical analysis was performed using SPSS version 16.0. RESULTS: We identified 652 patients who met our criteria. Within the elective patient group, there was a significantly longer length of stay of three days for the late-week group compared with two days for the early-week group (P = 0.016). Femoral artery procedures had a median length of stay of two days for those operated on early in the week compared with four days later in the week (P < 0.005). Open abdominal aortic aneurysm repair showed a trend to longer length of stay in the late-week group (P = 0.06). CONCLUSION: Day of surgery appears to impact on patients' length of stay following vascular procedures, with the greatest impact on medium-sized procedures. This difference could be explained by the difference in weekend support services, but further evaluation is required following introduction of weekend support services to assess this.

Advanced Methods for Dose and Regimen Finding During Drug Development: Summary of the EMA/EFPIA Workshop on Dose Finding (London 4-5 December 2014).

Inadequate dose selection for confirmatory trials is currently still one of the most challenging issues in drug development, as illustrated by high rates of late-stage attritions in clinical development and postmarketing commitments required by regulatory institutions. In an effort to shift the current paradigm in dose and regimen selection and highlight the availability and usefulness of well-established and regulatory-acceptable methods, the European Medicines Agency (EMA) in collaboration with the European Federation of Pharmaceutical Industries Association (EFPIA) hosted a multistakeholder workshop on dose finding (London 4-5 December 2014). Some methodologies that could constitute a toolkit for drug developers and regulators were presented. These methods are described in the present report: they include five advanced methods for data analysis (empirical regression models, pharmacometrics models, quantitative systems pharmacology models, MCP-Mod, and model averaging) and three methods for study design optimization (Fisher information matrix (FIM)-based methods, clinical trial simulations, and adaptive studies). Pairwise comparisons were also discussed during the workshop; however, mostly for historical reasons. This paper discusses the added value and limitations of these methods as well as challenges for their implementation. Some applications in different therapeutic areas are also summarized, in line with the discussions at the workshop. There was agreement at the workshop on the fact that selection of dose for phase III is an estimation problem and should not be addressed via hypothesis testing. Dose selection for phase III trials should be informed by well-designed dose-finding studies; however, the specific choice of method(s) will depend on several aspects and it is not possible to recommend a generalized decision tree. There are many valuable methods available, the methods are not mutually exclusive, and they should be used in conjunction to ensure a scientifically rigorous understanding of the dosing rationale.

Model Evaluation of Continuous Data Pharmacometric Models: Metrics and Graphics.

This article represents the first in a series of tutorials on model evaluation in nonlinear mixed effect models (NLMEMs), from the International Society of Pharmacometrics (ISoP) Model Evaluation Group. Numerous tools are available for evaluation of NLMEM, with a particular emphasis on visual assessment. This first basic tutorial focuses on presenting graphical evaluation tools of NLMEM for continuous data. It illustrates graphs for correct or misspecified models, discusses their pros and cons, and recalls the definition of metrics used.

Rocuronium pharmacokinetics and pharmacodynamics in the adductor pollicis and masseter muscles.

BACKGROUND: The aim of this study was to characterize the dose-effect relationship of rocuronium at the adductor pollicis and masseter muscles. METHODS: Ten, ASA I, adult patients, received a bolus dose of rocuronium 0.3 mg/kg during propofol based anesthesia. Train-of-four (TOF) was simultaneously monitored at the masseter and the adductor pollicis muscles until recovery. Rocuronium arterial serum concentrations were measured during 120 min. The first twitch of the TOF response was used to characterize the time-effect profile of both muscles using pharmacokinetic-pharmacodynamic analysis in NONMEM. A decrease in NONMEM objective function (∆OFV) of 3.84 points for an added parameter was considered significant at the 0.05 level. RESULTS: Onset time at the masseter (mean ± SD, 1.5 ± 0.9 min) was faster than at the adductor pollicis (2.7 ± 1.4 min, P < 0.05). Recovery, measured as the time to TOF ratio = 0.9 was similar between muscles 29.9 ± 6.7 (adductor pollicis) vs. 29.3 ± 8.1 (masseter). (P = 0.77). The estimated pharmacodynamic parameters [mean (95% CI)] of the adductor pollicis muscle and the masseter muscle were; plasma effect-site equilibration half-time (teq) 3.25 (2.34, 3.69) min vs. 2.86 (1.83, 3.29) min, (∆OFV 383.665); Ce50 of 1.24 (1.13, 1.56) mg/l vs. 1.19 (1.00, 1.21) mg/l, (∆OFV 184.284); Hill coefficient of 3.97 (3.82, 5.62) vs. 4.68 (3.83, 5.71), (∆OFV 78.906). CONCLUSIONS: We found that the masseter muscle has faster onset of blockade and similar recovery profile than adductor pollicis muscle. These findings were best, explained by a faster plasma effect-site equilibration of the masseter muscle to rocuronium.

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.

Pharmacometrics: opportunity for reducing disease burden in the developing world: the case of Africa.

Pharmacometricians are virtually nonexistent in Africa and the developing world. The unrelenting burden of infectious diseases, which are often treated using medicines with narrow effectiveness and safety dose ranges, and the growing prevalence and recognition of non-communicable diseases represent significant threats for the patients, although affording an opportunity for advancing science. This article outlines the case for pharmacometricians to redirect their expertise to focus on the disease burden affecting the developing world.CPT: Pharmacometrics & Systems Pharmacology (2013) 2, e69; doi:10.1038/psp.2013.45; published online 28 August 2013.

Modeling helps in understanding antidepressants.

The clinical pharmacology of antidepressant drugs is challenging to describe because of difficulties in describing the depressed state, the time course of response to a treatment intervention, and the frequent withdrawal of patients from clinical trials. Many of these challenges were addressed in a recent article by Russu et al. in this journal.

Influence of obesity on propofol pharmacokinetics: derivation of a pharmacokinetic model.

BACKGROUND: The objective of this study was to develop a pharmacokinetic (PK) model to characterize the influence of obesity on propofol PK parameters. METHODS: Nineteen obese ASA II patients undergoing bariatric surgery were studied. Patients received propofol 2 mg kg(-1) bolus dose followed by a 5-20-40-120 min, 10-8-6-5 mg kg(-1) h(-1) infusion. Arterial blood samples were withdrawn at 1, 3, 5 min after induction, every 10-20 min during propofol infusion, and every 10-30 min for 2 h after stopping the propofol infusion. Arterial samples were processed by high-performance liquid chromatography. Time-concentration data profiles from this study were pooled with data from two other propofol PK studies available at http://www.opentci.org. Population PK modelling was performed using non-linear mixed effects model. RESULTS: The study involved 19 obese adults who contributed 163 observations. The pooled analysis involved 51 patients (weight 93 sd 24 kg, range 44-160 kg; age 46 sd 16 yr, range 25-81 yr; BMI 33 sd 9 kg m(-2), range 16-52 kg m(-2)). A three-compartment model was used to investigate propofol PK. An allometric size model using total body weight (TBW) was superior to all other models investigated (linear TBW, free fat mass, lean body weight, normal fat mass) for all clearance parameters. Variability in V2 and Q2 was reduced by a function showing a decrease in both parameters with age. CONCLUSIONS: We have derived a population PK model using obese and non-obese data to characterize propofol PK over a wide range of body weights. An allometric model using TBW as the size descriptor of volumes and clearances was superior to other size descriptors to characterize propofol PK in obese patients.

Clinical trial simulation: a review.

Modeling and simulation in general, and specifically clinical trial simulation (CTS), have been recognized by the (larger) pharmaceutical companies and regulatory authorities as being pivotal to improving the efficiency of the drug development process. This includes the use of CTS to learn about drug effectiveness and safety and to optimize trial designs at the various stages of development. By reviewing papers published during the period January 2000-January 2010, this paper discusses recent developments with respect to methodology, applications, and lessons learned in the use of CTS in the development and clinical use of specific drugs. It is expected that future CTS experiments will be aided by the hybridization of optimal design methods with computationally intensive stochastic simulations. This will take advantage of optimizing the experimental design and leave the task of evaluating the probable real-world performance of a limited number of candidate trial designs and analysis procedures.

Stability of beta-actin mRNA in plasma.

It has long been thought that mRNA is labile and easily prone to degradation. However, a recent study demonstrated that GAPDH mRNA in cell-free plasma may remain stable up to 24 hours after blood collection. As there are no other independent studies, we attempted to reproduce the findings of that study. In our study, blood was collected from a healthy male volunteer into Vacutainer tubes containing EDTA. Blood samples were placed on ice and plasma separated by double-centrifugation at times 0, 1, 2, and 5 hours after blood collection. mRNA was extracted from four aliquots of the blood sample by means of the QIAamp Viral RNA kit. Extracted mRNA was converted to cDNA by reverse transcription before real time quantitative PCR measurement of the housekeeping beta-actin gene. Plasma beta-actin mRNA at 2 hours (0.012; 0.0031-0.0297, median and range) was significantly lower (P= 0.022) than at 0 hours (0.12; 0.057-0.165) (P= 0.016). The levels decreased further at 5 hours (0.0037; 0.0024-0.011) (P= 0.004). The results show that plasma beta-actin mRNA levels decrease with time after blood collection and that this is likely to be due to degradation in vitro.

Normalization of circulating nucleic acid results.

Nucleic acids, both DNA and mRNA, have been detected in the circulation and have been demonstrated to be useful in such areas as fetal medicine, oncology, and transplantation. When mRNA is measured in circulating blood, the results are expressed in relation to a reference gene product in order to correct for any differences in extraction, volume of starting material or other differences. Many authors use beta-actin mRNA and express results as a ratio of target mRNA to beta-actin mRNA. We have used a similar approach when studying diabetic retinopathy. Recently, we planned to investigate the expression of thyroid dependent gene expression in acutely ill patients. As a control study, we examined the expression of thyroid hormone-dependent gene expression in subjects with hyperthyroidism and found that the expression of beta-actin mRNA was affected by thyroid hormone status. Blood samples were taken into PAX genetrade mark tubes from 31 healthy subjects (mean age, 43 +/- 16 yrs) and 7 patients with hyperthyroidism (mean age, 43 +/- 5 yrs). Diagnosis of hyperthyroidism was confirmed by clinical findings and biochemical results. After extraction of mRNA, cDNA was synthesized using reverse transcription. Quantification of Na/K-ATPase, T3 receptor, and beta-actin cDNA was carried out by RT-PCR. Median beta-actin levels were significantly higher in hyperthyroid subjects compared to healthy subjects (18.2 versus 2.30; P < 0.00042). When mRNA for the T3 receptor was expressed in relation to beta-actin, there was a significantly higher in hyperthyroidism (0.0168 versus 0.218, P < 0.05). However, this was significantly lower when expressed in relation to total RNA (12.2 versus 2.24, P < 0.00015). We conclude that normalizing results to beta-actin may not be appropriate in all circumstances.

Morphine pharmacokinetics and pharmacodynamics in preterm and term neonates: secondary results from the NEOPAIN trial.

BACKGROUND: Relationships between plasma morphine concentrations and neonatal responses to endotracheal tube (ETT) suctioning are unknown in preterm neonates. METHODS: Ventilated preterm neonates (n=898) from 16 centres were randomly assigned to placebo (n=449) or morphine (n=449). After an i.v. loading dose (100 microg kg(-1)), morphine infusions [23-26 weeks postmenstrual age (PMA) 10 microg kg(-1) h(-1); 27-29 weeks 20 microg kg(-1) h(-1); and 30-32 weeks 30 microg kg(-1) h(-1)] were established for a maximum of 14 days. Open-label morphine (20-100 microg kg(-1)) was given for pain or agitation. Morphine assay and neonatal response to ETT suctioning was measured at 20-28 and 70-76 h after starting the drug infusion and at 10-14 h after discontinuation of the study drug. The concentration-effect response was investigated using non-linear mixed effects models. RESULTS: A total of 5119 data points (1598 measured morphine concentrations and 3521 effect measures) were available from 875 neonates for analysis. Clearance was 50% that of the mature value at 54.2 weeks PMA (CLmat(50)) and increased from 2.05 litre h(-1) 70 kg(-1) at 24 weeks PMA to 6.04 litre h(-1) 70 kg(-1) at 32 weeks PMA. The volume of distribution in preterm neonates was 190 litre 70 kg(-1) (CV 51%) and did not change with age. There was no relationship between morphine concentrations (range 0-440 microg litre(-1)) and heart rate changes associated with ETT suctioning or with the Premature Infant Pain Profile. CONCLUSIONS: A sigmoid curve describing maturation of morphine clearance is moved to the right in preterm neonates and volume of distribution is increased compared with term neonates. Morphine does not alter the neonatal response to ETT suctioning.

Mechanism-based concepts of size and maturity in pharmacokinetics.

Growth and development can be investigated using readily observable demographic factors such as weight and age. Size is the primary covariate and can be referenced to a 70-kg person with allometry using a coefficient of 0.75 for clearance and 1 for volume. The use of these coefficients is supported by fractal geometric concepts and observations from diverse areas in biology. Fat free mass (FFM) might be expected to do better than total body weight when there are wide variations in fat affecting body composition. Clearance pathways develop in the fetus before birth. The use of postnatal age as a descriptor of maturation is unsatisfactory because birth may occur prematurely; therefore postmenstrual age is a superior predictor of elimination function. A sigmoid E(max) model (Hill equation) describes gradual maturation of clearance in early life leading to a mature adult clearance achieved at a later age.

Time for quantitative clinical pharmacology: a proposal for a pharmacometrics curriculum.

A formal training program in pharmacometrics is essential to train clinical pharmacology scientists. A proposal is made for a pharmacometrics curriculum. The curriculum has components at the undergraduate, graduate and postgraduate levels.

Developmental pharmacokinetics of morphine and its metabolites in neonates, infants and young children.

BACKGROUND: Descriptions of the pharmacokinetics and metabolism of morphine and its metabolites in young children are scant. Previous studies have not differentiated the effects of size from those related to age during infancy. METHODS: Postoperative children 0-3 yr old were given an intravenous loading dose of morphine hydrochloride (100 micro g kg(-1) in 2 min) followed by either an intravenous morphine infusion of 10 micro g h(-1) kg(-1) (n=92) or 3-hourly intravenous morphine boluses of 30 micro g kg(-1) (n=92). Additional morphine (5 micro g kg(-1)) every 10 min was given if the visual analogue (VAS, 0-10) pain score was >/=4. Arterial blood (1.4 ml) was sampled within 5 min of the loading dose and at 6, 12 and 24 h for morphine, morphine-3-glucuronide (M3G) and morphine-6-glucuronide (M6G). The disposition of morphine and formation clearances of morphine base to its glucuronide metabolites and their elimination clearances were estimated using non-linear mixed effects models. RESULTS: The analysis used 1856 concentration observations from 184 subjects. Population parameter estimates and their variability (%) for a one-compartment, first-order elimination model were as follows: volume of distribution 136 (59.3) litres, formation clearance to M3G 64.3 (58.8) litres h(-1), formation clearance to M6G 3.63 (82.2) litres h(-1), morphine clearance by other routes 3.12 litres h(-1) per 70 kg, elimination clearance of M3G 17.4 (43.0) litres h(-1), elimination clearance of M6G 5.8 (73.8) litres h(-1). All parameters are standardized to a 70 kg person using allometric 3/4 power models and reflect fully mature adult values. The volume of distribution increased exponentially with a maturation half-life of 26 days from 83 litres per 70 kg at birth; formation clearance to M3G and M6G increased with a maturation half-life of 88.3 days from 10.8 and 0.61 litres h(-1) per 70 kg respectively at birth. Metabolite formation decreased with increased serum bilirubin concentration. Metabolite clearance increased with age (maturation half-life 129 days), and appeared to be similar to that described for glomerular filtration rate maturation in infants. CONCLUSION: M3G is the predominant metabolite of morphine in young children and total body morphine clearance is 80% that of adult values by 6 months. A mean steady-state serum concentration of 10 ng ml(-1) can be achieved in children after non-cardiac surgery in an intensive care unit with a morphine hydrochloride infusion of 5 micro g h(-1) kg(-1) at birth (term neonates), 8.5 micro g h(-1) kg(-1) at 1 month, 13.5 micro g h(-1) kg(-1) at 3 months and 18 micro g h(-1) kg(-1) at 1 year and 16 micro g h(-1) kg(-1) for 1- to 3-yr-old children.

Population PKPD modelling of the long-term hypoglycaemic effect of gliclazide given as a once-a-day modified release (MR) formulation.

AIMS: To study the relationship between the pharmacokinetics (PK) of gliclazide and its long-term pharmacodynamic (PD) effect in a large population of Type 2 diabetic patients and to identify factors predicting intersubject variability. METHODS: A PKPD database of 634 Type 2 diabetic patients with a total of 5,258 fasting plasma glucose (FPG) samples was built up from the data collected during the clinical development of a modified release formulation of gliclazide (gliclazide MR). The PKPD analysis used a nonlinear mixed effect modelling approach. A mixture model was used to identify patients with a FPG response to treatment. In patients identified as responders, the decrease in FPG was related to gliclazide exposure (AUC) by an Emax relationship. An effect compartment was used to describe the link between PK and PD. A linear disease-progression model was used to assess the glycaemic deterioration observable over several months of treatment. Simulations were performed to evaluate the predictive performance of the PKPD model and to illustrate the time course of the antidiabetic effect of gliclazide MR. RESULTS: Disease state was found to be the main explanatory factor for intersubject variability in response to gliclazide. The percentage of responders to gliclazide, used as monotherapy, increased inversely to the number of classes of antidiabetic agents received prior to entry in the studies. In responders, the initial dose (30 mg) of the gliclazide MR dosing regimen induced half of the maximum hypoglycaemic effect. The equilibration half-life between the PK and PD steady states was 3 weeks (intersubject variability of 84%). The rate of disease progression was 0.84 mmol l(-1) year(-1) (intersubject variability 143%). The PKPD model adequately predicted the FPG profiles of 234 patients who received the current formulation of gliclazide. Simulation of a 1-year parallel dose ranging clinical trial illustrated the influence of dose, time and type of previous antidiabetic treatment on the percentage of patients with clinically significant improvement of blood glucose control. CONCLUSIONS: This population PKPD analysis has characterized the relationship between the exposure to gliclazide and its long-term hypoglycaemic effect, and has established that the intersubject variability in response is mostly related to disease state. These results underline the clinical interest of quickly increasing the dose of gliclazide MR according to the response to treatment in order to achieve effective blood glucose control.

A pharmacokinetic model to predict the PK interaction of L-dopa and benserazide in rats.

PURPOSE: To study the PK interaction of L-dopa/benserazide in rats. METHODS: Male rats received a single oral dose of 80 mg/kg L-dopa or 20 mg/kg benserazide or 80/20 mg/kg L-dopa/benserazide. Based on plasma concentrations the kinetics of L-dopa, 3-O-methyldopa (3-OMD), benserazide, and its metabolite Ro 04-5127 were characterized by noncompartmental analysis and a compartmental model where total L-dopa clearance was the sum of the clearances mediated by amino-acid-decarboxylase (AADC), catechol-O-methyltransferase and other enzymes. In the model Ro 04-5127 inhibited competitively the L-dopa clearance by AADC. RESULTS: The coadministration of L-dopa/benserazide resulted in a major increase in systemic exposure to L-dopa and 3-OMD and a decrease in L-dopa clearance. The compartmental model allowed an adequate description of the observed L-dopa and 3-OMD concentrations in the absence and presence of benserazide. It had an advantage over noncompartmental analysis because it could describe the temporal change of inhibition and recovery of AADC. CONCLUSIONS: Our study is the first investigation where the kinetics of benserazide and Ro 04-5127 have been described by a compartmental model. The L-dopa/benserazide model allowed a mechanism-based view of the L-dopa/benserazide interaction and supports the hypothesis that Ro 04-5127 is the primary active metabolite of benserazide.