Search

We compared the predictive performance of an artificial neural network to traditional pharmacometric modeling for population prediction of plasma concentrations of valproate in real-world data. We included individuals aged 65 years or older with epilepsy who redeemed their first prescription of valproate after the diagnosis of epilepsy and had at least one valproate plasma concentration measured. A long short-term memory neural network (LSTM) was developed using the training data set to fit the LSTM and the test data set to validate the model. Predictions from the LSTM were compared with those obtained from the population predictions from a pharmacometric model by Birnbaum et al. which had the best predictive performance for population predictions of valproate concentrations in Danish databases. We used the cutoff of ± 20 mg/L of prediction error to define good predictions. A total of 1,252 individuals were included in the study. The LSTM fitted using the training data set had poor predictive performance in the test data set, but better than that of the pharmacometric model. The proportion of individuals with at least one predicted concentration within ± 20 mg/L of observed concentration was largest in case of the LSTM (64.4%, 95% confidence interval (CI): 58.4-70.2%) compared with the pharmacometric model by Birnbaum et al. (49.8%, 95% CI: 47.0-52.6%). LSTM shows better predictive performance to predict valproate plasma concentrations compared with a traditional pharmacometric model in the investigated setting with real-world data in older patients with epilepsy where information on exact timepoints for both dosing and plasma concentration measurement are missing.

In pharmacoepidemiology, it is usually expected that the observed association should be directly or indirectly related to the pharmacological effects of the drug/s under investigation. Pharmacological effects are, in turn, strongly connected to the pharmacokinetic and pharmacodynamic properties of a drug, which can be characterized and investigated using pharmacometric models. Recently, the use of pharmacometrics has been proposed to provide pharmacological substantiation of pharmacoepidemiological findings derived from real-world data. However, validated frameworks suggesting how to combine these two disciplines for the aforementioned purpose are missing. Therefore, we propose PHARMACOM-EPI, a framework that provides a structured approach on how to identify, characterize, and apply pharmacometric models with practical details on how to choose software, format dataset, handle missing covariates/dosing data, how to perform the external evaluation of pharmacometric models in real-world data, and how to provide pharmacological substantiation of pharmacoepidemiological findings. PHARMACOM-EPI was tested in a proof-of-concept study to pharmacologically substantiate death associated with valproate use in the Danish population aged ≥ 65 years. Pharmacological substantiation of death during a follow-up period of 1 year showed that in all individuals who died (n = 169) individual predictions were within the subtherapeutic range compared with 52.8% of those who did not die (n = 1,084). Of individuals who died, 66.3% (n = 112) had a cause of death possibly related to valproate and 33.7% (n = 57) with well-defined cause of death unlikely related to valproate. This proof-of-concept study showed that PHARMACOM-EPI was able to provide pharmacological substantiation for death associated with valproate use in the study population.

AIMS: To investigate the impact of treatment changes on weight loss for the glucagon-like peptide-1 analogue, liraglutide, in subjects with overweight or obesity that may be needed in case of gastrointestinal tolerability issues during escalation.MATERIALS AND METHODS: Individual longitudinal body weight data from the main trial periods of three phase II/III trials in subjects with overweight or obesity (56-week treatment with once-daily liraglutide 1.2, 1.8, 2.4 or 3.0 mg or placebo, n=4952) were analysed using a non-linear mixed effect modelling approach. Individual pharmacokinetic profiles were derived based on published pharmacokinetic models. Baseline body weight, baseline HbA1c, age, gender, diabetes status (non-diabetic, pre-diabetic or with type 2 diabetes mellitus), race and trial region were investigated as covariates. As a form of external validation, the model was used to predict the weight regain following treatment cessation at week 56 (data not included in model development).RESULTS: A PK/PD model provided an adequate description of the weight loss trajectories for all studied doses. Gender and diabetes status were identified as the most influential covariates, and an underlying seasonal weight fluctuation was identified. Slower than the recommended, one-week dose-escalation algorithms led up to 2 weeks slower initial weight loss but similar long-term weight loss trajectories.CONCLUSIONS: The relationship between liraglutide systemic exposure and weight loss was successfully established in subjects with overweight or obesity. The model could predict the time-course of weight-regain following treatment cessation and suggests that gastrointestinal tolerability can be mitigated by slower escalation with only minor impact on the weight loss trajectory. This article is protected by copyright. All rights reserved.

Background Naloxone, an opioid receptor antagonist, is used as a pharmacological tool to detect tonic endogenous activation of opioid receptors in experimental pain models. We describe a pharmacokinetic model linking naloxone pharmacokinetics to its main metabolite after high-dose naloxone infusion. Methods Eight healthy volunteers received a three-stage stepwise high-dose i.v. naloxone infusion (total dose 3.25 mg kg−1). Naloxone and naloxone-3-glucuronide (N3G) plasma concentrations were sampled from infusion onset to 334 min after infusion discontinuation. Pharmacokinetic analysis was performed using non-linear mixed effect models (NONMEM). The predictive performances of Dowling's and Yassen's models were evaluated, and target-controlled infusion simulations were performed. Results Three- and two-compartment disposition models with linear elimination kinetics described the naloxone and N3G concentration time-courses, respectively. Two covariate models were developed: simple (weight proportional) and complex (with the shallow peripheral volume of distribution linearly increasing with body weight). The median prediction error (MDPE) and wobble for Dowling's model were –32.5% and 33.4%, respectively. For Yassen's model, the MDPE and wobble were 1.2% and 19.9%, respectively. Conclusions A parent–metabolite pharmacokinetic model was developed for naloxone and N3G after high-dose naloxone infusion. No saturable pharmacokinetics were observed. Whereas Dowling's model was inaccurate and over-predicted naloxone concentrations, Yassen's model accurately predicted naloxone pharmacokinetics. The newly developed covariate models may be used for high-dose TCI-naloxone for experimental and clinical practice. Clinical trials registration NCT01992146.

The International Prostate Symptom Score (IPSS), the quality of life (QoL) score, and the benign prostatic hyperplasia impact index (BII) are three different scales commonly used to assess the severity of lower urinary tract symptoms associated with benign prostatic hyperplasia (BPH-LUTS). Based on a phase II clinical trial including 403 patients with moderate to severe BPH-LUTS, the objectives of this study were to (i) develop traditional pharmacometric and bounded integer (BI) models for the IPSS, QoL score, and BII endpoints, respectively; (ii) compare the power and type I error in detecting drug effects of BI modeling with traditional methods through simulation; and (iii) obtain quantitative translation between scores on the three abovementioned scales using a BI modeling framework. All developed models described the data adequately. Pharmacometric modeling using a continuous variable (CV) approach was overall found to be the most robust in terms of type I error and power to detect a drug effect. In most cases, BI modeling showed similar performance to the CV approach, yet severely inflated type I error was generally observed when inter-individual variability (IIV) was incorporated in the BI variance function (g()). BI modeling without IIV in g() showed greater type I error control compared to the ordered categorical approach. Lastly, a multiple-scale BI model was developed and estimated the relationship between scores on the three BPH-LUTS scales with overall low uncertainty. The current study yields greater understanding of the operating characteristics of the novel BI modeling approach and highlights areas potentially requiring further improvement. Supplementary Information The online version contains supplementary material available at 10.1208/s12248-021-00568-y.

Inconsistent trial design and analysis is a key reason that few advances in postoperative pain management have been made from clinical trials analyzing opioid consumption data. This study aimed to compare four different approaches to analyze opioid consumption data. A repeated time-to-event (RTTE) model in NONMEM was used to simulate clinical trials of morphine consumption with and without a hypothetical adjuvant analgesic in doses equivalent to 15-62% reduction in morphine consumption. Trials were simulated with duration of 24-96 h. Monte Carlo simulation and re-estimation were performed to determine sample size required to demonstrate efficacy with 80% power using t test, Mann-Whitney rank sum test, time-to-event (TTE) modeling and RTTE modeling. Precision of efficacy estimates for RTTE models were evaluated in 500 simulations. A sample size of 50 patients was required to detect 37% morphine sparing effect with at least 80% power in a 24 h trial with RTTE modeling whereas the required sample size was 200 for Mann-Whitney, 180 for t-test and 76 for TTE models. Extending the trial duration from 24 to 96 h reduced the required sample size by 3.1 fold with RTTE modeling. Precise estimate of potency was obtained with a RTTE model accounting for both morphine effects and time-varying covariates on opioid consumption. An RTTE analysis approach proved better suited for demonstrating efficacy of opioid sparing analgesics than traditional statistical tests as a lower sample size was required due the ability to account for time-varying factors including PK.

Exposure-response (ER) modeling for fixed-dose combinations (FDC) has previously been found to have an inflated false positive rate (FP), i.e., observing a significant effect of FDC components when no true effect exists. Longitudinal exposure-response (LER)analysis utilizes the time course of the data and is valid for several clinical endpoints for FDCs. The aim of the study was to investigate if LER is applicable for the validation of FDCs by demonstrating the contribution of each component to the overall effect without inflation of FP rates. FP and FN rates associated with ER and LER analysis were investigated using stochastic simulation and estimation. Four hundred thirty-two scenarios with varying numbers of patients, duration, sampling frequency, dose distribution, design, and drug activity were analyzed using a range of linear, log-linear, and non-linear models to asses FP and FN rates. Lastly, the impact of the clinical trial parameters was investigated. LER analyses provided well-controlled FP rates of the expected 5% or less; however, in low information clinical trials consisting of 30 patients, 4 samples, and 20 days, LER analyses lead to inflated FN rates. Parameter investigation showed that when the clinical trial includes sufficient patients, duration, samples, and an appropriate trial design, the FN rates are in general below the expected 5% for LER analysis. Based on the results, LER analysis can be used for the validation of FDCs and fixed ratio drug combinations. The method constitutes anew avenue for providing evidence that demonstrates the contribution of each component to the overall clinical effect. Refereed/Peer-reviewed

In clinical trials within lower urinary tract symptoms due to benign prostatic hyperplasia (BPH-LUTS), the International Prostate Symptom Score (IPSS) is commonly the primary efficacy outcome while the Quality of Life (QoL) score and the BPH Impact Index (BII) are common secondary efficacy markers. The current study aimed to characterize BPH-LUTS progression using responses to the IPSS, the QoL, and the BII in an integrated item response theory (IRT) framework and assess the Fisher information of each scale. The power of this approach to detect a drug effect was compared with an IRT approach considering only IPSS responses. A unidimensional and a bidimensional pharmacometric IRT model, based on item-level IPSS responses in a clinical trial with 403 patients, were extended by incorporating patients’ QoL and summary BII scores over the 6-month trial period. In the developed unidimensional integrated model, the QoL score was found to be the most informative, representing 17% of the total Fisher information, while the combined information content of the seven IPSS items represented 70.6%. In the bidimensional model, “storage” and both storage and “voiding” disability drove QoL and summary BII responses, respectively. Sample size reduction of 16% to detect a drug effect at 80% power was obtained with the unidimensional integrated IRT model compared with its counterpart IPSS IRT model. This study shows that utilizing the information content across the IPSS, QoL, and BII scales in an integrated IRT framework results in a modest but meaningful increase in power to detect a drug effect. Electronic supplementary material The online version of this article (10.1208/s12248-020-00484-7) contains supplementary material, which is available to authorized users.

One of the primary barriers in treating cancer patients is the development of resistance to the available treatments. This is the case for treatment of triple negative breast cancer (TNBC) with docetaxel, which is part of the neoadjuvant treatment for TNBC. The novel compound SCO-101 is under investigation for its potential treatment effect in several types of drug resistant cancer. The aim of this study was to establish a pharmacodynamic model that captures the effect of docetaxel, SCO-101, and the combination on cell survival in docetaxel resistant MDA-MB-231 TNBC cells. Several combination models were compared and a recently published combination model, the general pharmacodynamic interaction model (GPDI), provided the best fit. The model allowed for description and quantification of the interaction between docetaxel and SCO-101 with respects to both maximal effect and potency. Based on this model, SCO-101 has a synergistic effect with docetaxel. This synergy is not present in the maximal effect, but the combination of SCO-101 and docetaxel showed an approximately 60% increase in potency compared to docetaxel alone. Furthermore, the predicted model surface for the combination provided key information regarding promising dose ratios and dose levels for further studies of the combination. Lastly, the study presents a use case for the GPDI model, which provides a way to quantify and interpret drug-drug interactions.

Midazolam is a first-line drug for the treatment of status epilepticus, both by buccal and intravenous administration. In children and adolescents with obesity, midazolam pharmacokinetics may be altered, and the current dosing guidelines may therefore be insufficient. The objective of this study was to investigate the pharmacokinetics of midazolam, after intravenous administration, in obese and non-obese adolescents aged 11–18 years. All trial participants received a 1-µg midazolam microdose as an intravenous bolus. 13 blood samples were collected per participant at pre-specified timepoints. Plasma concentration–time data were fitted to pharmacokinetic models using non-linear mixed-effects modeling. Covariates such as weight, age, and body mass index standard deviation score were tested to explain the inter-individual variability associated with the pharmacokinetic parameters. Sixty-seven adolescents were included in the analysis. The pharmacokinetics of midazolam was best described with a two-compartment model. The rate of distribution was faster, and the peripheral volume of distribution was larger in adolescents with a high body mass index standard deviation score compared with adolescents with a lower standard deviation score. Simulations revealed that long-term infusions based on total body weight could lead to high plasma concentrations in adolescents with obesity. Furthermore, simulated plasma concentrations after a fixed buccal dose indicated that adolescents with obesity may be at risk of sub-therapeutic midazolam plasma concentrations. The body mass index standard deviation score was shown to have a significant influence on the peripheral volume of distribution and the inter-compartmental clearance of midazolam. The current dosing guidelines for status epilepticus, where the midazolam dose is adjusted to total body weight or age, may lead to supra- and sub-therapeutic plasma concentrations, respectively, in adolescents with obesity. EudraCT: 2014-004554-34.

Combinations of pharmacological treatments are increasingly being investigated for potentially higher clinical benefit, especially when the combined drugs are expected to act via synergistic interactions. The clinical development of combination treatments is particularly challenging, particularly during the dose-selection phase, where a vast range of possible combination doses exists. The purpose of this work was to evaluate the added value of using optimal design for guiding the dose allocation in drug combination dose-finding studies as compared with a typical drug-combination trial. Optimizations were performed using local [D(s)-optimality] and global [ED(s)-optimality] optimal designs to maximize the precision of model parameters in a number of potential exposure-response (E-R) surfaces. A compound criterion [D(s)/V-optimality] was used to optimize the precision of model predictions in specific parts of the E-R surfaces. Optimal designs provided unbiased estimates and significantly improved the accuracy of results relative to the typical design. It was possible to improve the efficiency and overall parameter precision up to 7832% and 96.6% respectively. When the compound criterion was used, the probability to accurately identify the optimal dose-combination increased from 71% for the typical design up to 91%. These results indicate that optimal design methodology in tandem with E-R analyses is a beneficial tool that can be used for appropriate dose allocation in dose-finding studies for drug combinations.

The combination of morphine and gabapentin seems promising for the treatment of postoperative and neuropathic pain. Despite the well characterised pharmacodynamic interaction, little is known about possible pharmacokinetic interactions. The aim of this study was to evaluate whether co-administration of the two drugs leads to modifications of their pharmacokinetic profiles. The pharmacokinetics of morphine, morphine-3-glucuronide and gabapentin were characterised in rats following subcutaneous injections of morphine, gabapentin or their combination. Non-linear mixed effects modelling was applied to describe the pharmacokinetics of the compounds and possible interactions. The plasma-concentration-time profiles of morphine and gabapentin were best described using a three- and a one-compartment disposition model respectively. Dose dependencies were found for morphine absorption rate and gabapentin bioavailability. Enterohepatic circulation of morphine-3-glucuronide was modelled using an oscillatory model. The combination did not lead to pharmacokinetic interactions for morphine or gabapentin but resulted in an estimated ~33% diminished morphine-3-glucuronide formation. The finding of a lack of pharmacokinetic interaction strengthens the notion that the combination of the two drugs leads to better efficacy in pain treatment due to interaction at the pharmacodynamic level. The interaction found between gabapentin and morphine-3-glucuronide, the latter being inactive, might not have any clinical relevance.

Aim: High-dose administration of the μ-opioid receptor inverse agonist naloxone (NX), has previously been demonstrated to reinstate nocifensive behavior in the late phase of inflammatory injuries. However, no current analytical methods can provide pharmacokinetic insight into the pharmacodynamic response of high-dose administration of NX. Materials & methods: Based on protein precipitation using 50 μl human plasma, NX and naloxone-β-d-glucuronide (NXG) was analysed by UHPLC–MS/MS with 6 min cycle time. Results: A method for quantification of high-dose administered NX and NXG was developed and validated with intra- and interday precision and accuracy within ≤8.5% relative standard deviation (RSD) and -1.2–5.5% relative error (RE) for NX and ≤9.6% RSD and 0.6–6.5% RE for NXG. The method show excellent internal standard corrected matrix effects. Conclusion: A rapid UHPLC–MS/MS method was developed for quantification of NX and NXG in human plasma within 10–4000 ng/ml.

The exposure-response relationship of combinatory drug effects can be quantitatively described using pharmacodynamic interaction models, which can be used for the selection of optimal dose combinations. The aim of this simulation study was to evaluate the reliability of parameter estimates and the probability for accurate dose identification for various underlying exposure-response profiles, under a number of different phase II designs. An efficacy variable driven by the combined exposure of two theoretical compounds was simulated and model parameters were estimated using two different models, one estimating all parameters and one assuming that adequate previous knowledge for one drug is readily available. Estimation of all pharmacodynamic parameters under a realistic, in terms of sample size and study design, phase II trial, proved to be challenging. Inaccurate estimates were found in all exposure-response scenarios, except for situations where no pharmacodynamic interaction was present, with the drug potency and interaction parameters being the hardest to estimate. When previous knowledge of the exposure-response relationship of one of the monocomponents is available, such information should be utilized, as it enabled relevant improvements in parameter estimation and in correct dose identification. No general trends for classification of the performance of the tested study designs across different scenarios could be identified. This study shows that pharmacodynamic interactions models can be used for the exposure-response analysis of clinical endpoints especially when accompanied by appropriate dose selection in regard to the expected drug potencies and appropriate trial size and if information regarding the exposure-response profile of one monocomponent is available.

PURPOSE: To characterize the pharmacokinetic-pharmacodynamic (PK-PD) relationship between exposure of morphine and subsequent morphine consumption and to develop simulation tools for model validation.METHODS: Dose, formulation and time of morphine administration was available from a published study in 63 patients receiving intravenous, oral immediate release or oral controlled release morphine on request after hip surgery. The PK-PD relationship between predicted exposure of morphine and morphine consumption was modeled using repeated time to event (RTTE) modeling in NONMEM. To validate the RTTE model, a visual predictive check method was developed with simulated morphine consumption given the exposure of preceding morphine administration.RESULTS: The probability of requesting morphine was found to be significantly related to the exposure of morphine as well as night/day. Oral controlled release morphine was more effective than intravenous and oral immediate release formulations at equivalent average concentrations. Maximum effect was obtained for 8 h by oral controlled release doses ≥ 15 mg, where probability of requesting a new dose was reduced to 20% for a typical patient.CONCLUSION: This study demonstrates the first quantitative link between exposure of morphine and subsequent morphine consumption and introduces an efficient visual predictive check approach with simulation of adaptive dosing.

Glucose is the main energy substrate for neurons, and ketone bodies are known to be alternative substrates. However, the capacity of ketone bodies to support different neuronal functions is still unknown. Thus, a change in energy substrate from glucose alone to a combination of glucose and β-hydroxybutyrate might change neuronal function as there is a known coupling between metabolism and neurotransmission. The purpose of this study was to shed light on the effects of the ketone body β-hydroxybutyrate on glycolysis and neurotransmission in cultured murine glutamatergic neurons. Previous studies have shown an effect of β-hydroxybutyrate on glucose metabolism, and the present study further specified this by showing attenuation of glycolysis when β-hydroxybutyrate was present in these neurons. In addition, the NMDA receptor-induced calcium responses in the neurons were diminished in the presence of β-hydroxybutyrate, whereas a direct effect of the ketone body on transmitter release was absent. However, the presence of β-hydroxybutyrate augmented transmitter release induced by the KATP channel blocker glibenclamide, thus giving an indirect indication of the involvement of KATP channels in the effects of ketone bodies on transmitter release.Energy metabolism and neurotransmission are linked and involve ATP-sensitive potassium (KATP) channels. However, it is still unclear how and to what degree available energy substrate affects this link. We investigated the effect of changing energy substrate from only glucose to a combination of glucose and R-β-hydroxybutyrate in cultured neurons. Using the latter combination, glycolysis was diminished, NMDA receptor-induced calcium responses were lower, and the KATP channel blocker glibenclamide caused a higher transmitter release.

The combination of morphine and gabapentin has shown to be promising for managing postoperative pain but finding the right dose for the combination has proven to be a challenge. The purpose of this study was to quantitatively characterize the pharmacodynamic interaction between the two drugs and to identify the optimal concentration–effect relationship of the combination. Information regarding plasma concentrations and von Frey withdrawal thresholds following incisional surgery on Sprague Dawley rats, after administration of morphine, gabapentin, or their combination was available from published studies. The combined pharmacodynamic effect of morphine and gabapentin was analyzed and linked to drug plasma concentrations via a response surface approach using non-linear mixed-effect modeling. Full reversal of withdrawal thresholds for the pain stimulation to presurgery values was estimated at morphine plasma concentration of 435.1 ng/mL. Co-administration of up to 40 μg/mL of gabapentin led to a reduction of the needed morphine concentration down to 307.5 ng/mL (~ 29% reduction). Combination of concentration ranges of gabapentin between 20 and 40 μg/mL with any morphine concentrations between 100 and 600 ng/mL were found to lead up to 50% increased effect relatively to the effect attained by morphine alone. This study highlights the importance of finding the right combination in multimodal analgesia and demonstrates the usefulness of the response surface approach for the study of pharmacodynamic interactions. The proposed pharmacokinetic–pharmacodynamic model may provide the basis for a rational clinical trial design with the aim to identify the optimal dose combination ratios in humans.

Joint analysis of pain intensity and opioid consumption is encouraged in trials of postoperative pain. However, previous approaches have not appropriately addressed the complexity of their interrelation in time. In this study, we applied a non-linear mixed effects model to simultaneously study pain intensity and opioid consumption in a 4-h postoperative period for 44 patients undergoing percutaneous kidney stone surgery. Analysis was based on 748 Numerical Rating Scale (NRS) scores of pain intensity and 51 observed morphine and oxycodone dosing events. A joint model was developed to describe the recurrent pattern of four key phases determining the development of pain intensity and opioid consumption in time; (A) Distribution of pain intensity scores which followed a truncated Poisson distribution with time-dependent mean score ranging from 0.93 to 2.45; (B) Probability of transition to threshold pain levels (NRS ≥ 3) which was strongly dependent on previous pain levels ranging from 2.8-15.2% after NRS of 0-2; (C) Probability of requesting opioid when allowed (NRS ≥ 3) which was strongly correlated with the number of previous doses, ranging from 89.8% for requesting the first dose to 26.1% after three previous doses; (D) Reduction in pain scores after opioid dosing which was significantly related to the pain intensity at time of opioid request (P

Background Reduction in consumption of opioid rescue medication is often used as an endpoint when investigating analgesic efficacy of drugs by adjunct treatment, but appropriate methods are needed to analyze analgesic consumption in time. Repeated time-to-event (RTTE) modeling is proposed as a way to describe analgesic consumption by analyzing the timing of consecutive analgesic events. Methods Retrospective data were obtained from 63 patients receiving standard analgesic treatment including morphine on request after surgery following hip fracture. Times of analgesic events up to 96 h after surgery were extracted from hospital medical records. Parametric RTTE analysis was performed with exponential, Weibull, or Gompertz distribution of analgesic events using NONMEM®, version 7.2 (ICON Development Solutions, USA). The potential influences of night versus day, sex, and age were investigated on the probability. Results A Gompertz distribution RTTE model described the data well. The probability of having one or more analgesic events within 24 h was 80% for the first event, 55% for the second event, 31% for the third event, and 18% for fourth or more events for a typical woman of age 80 yr. The probability of analgesic events decreased in time, was reduced to 50% after 3.3 days after surgery, and was significantly lower (32%) during night compared with day. Conclusions RTTE modeling described analgesic consumption data well and could account for time-dependent changes in probability of analgesic events. Thus, RTTE modeling of analgesic events is proposed as a valuable tool when investigating new approaches to pain management such as opioid-sparing analgesia.

Despite much evidence that combination of morphine and gabapentin can be beneficial for managing postoperative pain, the nature of the pharmacological interaction of the two drugs remains unclear. The aim of this study was to assess the interaction of morphine and gabapentin in range of different dose combinations and investigate whether co-administration leads to synergistic effects in a preclinical model of postoperative pain. The pharmacodynamic effects of morphine (1, 3 and 7mg/kg), gabapentin (10, 30 and 100mg/kg) or their combination (9 combinations in total) were evaluated in the rat plantar incision model using an electronic von Frey device. The percentage of maximum possible effect (%MPE) and the area under the response curve (AUC) were used for evaluation of the antihyperalgesic effects of the drugs. Identification of synergistic interactions was based on Loewe additivity response surface analyses. The combination of morphine and gabapentin resulted in synergistic antihyperalgesic effects in a preclinical model of postoperative pain. The synergistic interactions were found to be dose dependent and the increase in observed response compared to the theoretical additive response ranged between 26 and 58% for the synergistic doses. The finding of dose-dependent synergistic effects highlights that choosing the right dose-dose combination is of importance in postoperative pain therapy. Our results indicate benefit of high doses of gabapentin as adjuvant to morphine. If these findings translate to humans, they might have important implications for the treatment of pain in postoperative patients.

Glutaric acidemia type 1 (GA1) is an autosomal recessively inherited deficiency of glutaryl-CoA dehydrogenase. Accumulating metabolites, 3-hydroxyglutaric (3-OH-GA), glutaric (GA), and trans-glutaconic (TG) acids, have been proposed to be involved in the development of the striatal degeneration seen in children with GA1 via an excitotoxic mechanism. We have studied the extent to which 3-OH-GA, GA, and TG are neurotoxic and whether neurotoxicity is caused by an excitotoxic mechanism in which 3-OH-GA, GA, or TG overactivates N-methyl-D-aspartate (NMDA) receptors. In cultured mouse neocortical neurons, all three compounds were weakly neurotoxic, possibly through activation of NMDA receptors. However, further studies in the rat cortical wedge preparation and with NMDA receptors expressed in Xenopus oocytes could not confirm an interaction of the compounds with NMDA receptors. It is concluded that the metabolites 3-OH-GA, GA, and TG are only weak neurotoxins and that the neurodegenerative cascade destroying the striatum in patients with GA1 involves mainly mechanisms other than excitoxicity.

Morphine is a widely used opioid for treatment of moderate to severe pain, but large interindividual variability in patient response and no clear guidance on how to optimise morphine dosage regimen complicates treatment strategy for clinicians. Population pharmacokinetic-pharmacodynamic models can be used to quantify dose-response relationships for the population as well as interindividual and interoccasion variability. Additionally, relevant covariates for population subgroups that deviate from the typical population can be identified and help clinicians in dose optimisation. This review provides a detailed overview of the published human population pharmacokinetic-pharmacodynamic studies for morphine analgesia in addition to basic drug disposition and pharmacological properties of morphine and its analgesic active metabolite, morphine-6-glucuronide, that may help identify future covariates. Furthermore, based on simulations from key pharmacokinetic-pharmacodynamic models, the contribution of morphine-6-glucuronide to the analgesic response in patients with renal insufficiency was investigated. Simulations were also used to examine the impact of effect-site equilibration half-life on time course of response. Lastly, the impact of study design on the likelihood of determining accurate pharmacodynamic parameters for morphine response was evaluated.

The aim of this study was to develop population pharmacokinetic-pharmacodynamic models for morphine in experimental pain induced by skin heat and muscle pressure, and to evaluate the experimental pain models with regard to assessment of morphine pharmacodynamics. In a randomised, double-blind, placebo-controlled, crossover study, 39 healthy volunteers received an oral dose of 30 mg morphine hydrochloride or placebo. Non-linear mixed effects modelling was used to describe the plasma concentrations of morphine and metabolites, and the analgesic effect of morphine on experimental pain in skin and muscle. Baseline pain metrics varied between individuals and occasions, and were described with interindividual and interoccasion variability. Placebo-response did not change with time. For both pain metrics, morphine effect was proportional to baseline pain and was described with a linear model with interindividual variability on drug effect slope and linked to an effect compartment for muscle pressure. The models indicate that a steady-state morphine concentration of 21 ng/ml causes 33% and 0.84% increases in stimulus intensity from baseline for muscle pressure and skin heat, respectively. The population pharmacokinetic-pharmacodynamic models developed in this study indicate that mechanical stimulation of muscle is a more clinically relevant pain stimulus for the assessment of morphine pharmacodynamics than thermal stimulation of skin. Refereed/Peer-reviewed

Primary carnitine deficiency (PCD) is a disorder of fatty acid oxidation with a high prevalence in the Faroe Islands. Only patients homozygous for the c.95AG (p.N32S) mutation have displayed severe symptoms in the Faroese patient cohort. In this study, we investigated carnitine levels in skeletal muscle, plasma, and urine as well as renal elimination kinetics before and after intermission with L-carnitine in patients homozygous for c.95AG.Five male patients homozygous for c.95AG were included. Regular L-carnitine supplementation was stopped and the patients were observed during five days. Blood and urine were collected throughout the study. Skeletal muscle biopsies were obtained at 0, 48, and 96 h.Mean skeletal muscle free carnitine before discontinuation of L-carnitine was low, 158 nmol/g (SD 47.4) or 5.4% of normal. Mean free carnitine in plasma (fC0) dropped from 38.7 (SD 20.4) to 6.3 (SD 1.7) μmol/L within 96 h (p 0.05). Mean T 1/2 following oral supplementation was approximately 9 h. Renal reabsorption of filtered carnitine following oral supplementation was 23%. The level of mean free carnitine excreted in urine correlated (R (2) = 0.78, p 0.01) with fC0 in plasma.Patients homozygous for the c.95AG mutation demonstrated limited skeletal muscle carnitine stores despite long-term high-dosage L-carnitine supplementation. Exacerbated renal excretion resulted in a short T 1/2 in plasma carnitine following the last oral dose of L-carnitine. Thus a treatment strategy of minimum three daily separate doses of L-carnitine is recommended, while intermission with L-carnitine treatment might prove detrimental.

Cultures of dissociated cerebellum from 5- to 6-day-old mice as well as of the N2A neuronal cell line were exposed to guanidino ethane sulfonate (GES, 2-5 mM) to reduce the cellular taurine content. Control cultures were kept in culture medium or medium containing 2-5 mM GES plus 2-5 mM taurine to restore the intracellular taurine content. Taurine depletion led to changes in the expression of certain splice variants of NCAM mRNA such as the AAG and the VASE containing forms, while no differences were seen in the expression of the three forms of NCAM protein. In the N2A cells taurine depletion led to a decreased migration rate of the cells. The results suggest that the reduced migration rate of neurons caused by taurine depletion may be correlated to changes in expression of certain adhesion molecules such as NCAM. Moreover, taurine appears to be involved in regulation of transcription processes.

Background: The ATP-sensitive K+ (KATP) channel openers levcromakalim and pinacidil are vasodilators that induce headache in healthy people. The neuropeptide calcitonin gene-related peptide (CGRP) induces headache in healthy people and migraine in migraineurs, potentially through a mechanism that involves opening of vascular or neuronal KATP channels and mast cell degranulation. Using rat as a model, we studied the molecular presence of KATP channels in the trigeminovascular system. Furthermore, we examined whether KATP channel openers stimulate the in vitro release of CGRP and whether they degranulate dural mast cells. Methods: mRNA and protein expression of KATP channel subunits were studied in the trigeminal ganglion (TG) and trigeminal nucleus caudalis (TNC) by qPCR and western blotting. In vitro CGRP release was studied after application of levcromakalim (1 µM) and diazoxide (10 µM) to freshly isolated rat dura mater, TG and TNC. Rat dural mast cells were challenged in situ with levcromakalim (10−5 M) to study its potential degranulation effect. Results: mRNA and protein of KATP channel subunits Kir6.1, Kir6.2, SUR1 and SUR2B were identified in the TG and TNC. KATP channel openers did not release or inhibit capsaicin-induced CGRP release from dura mater, TG or TNC. They did also not induce dural mast cell degranulation. Conclusions: KATP channel openers do not interact with CGRP release or mast cell degranulation. Activation of these channels in the CNS is antinociceptive and therefore cannot explain the headache induced by KATP channel openers. Thus, they are likely to induce headache by interaction with extracerebral KATP channels, probably the SUR2B isoforms.

The ketogenic diet has multiple beneficial effects not only in treatment of epilepsy, but also in that of glucose transporter 1 deficiency, cancer, Parkinson's disease, obesity and pain. Thus, there is an increasing interest in understanding the mechanism behind this metabolic therapy. Patients on a ketogenic diet reach high plasma levels of ketone bodies, which are used by the brain as energy substrates. The interaction between glucose and ketone bodies is complex and there is still controversy as to what extent it affects the homeostasis of the neurotransmitters glutamate, aspartate and GABA. The present study was conducted to study this metabolic interaction in cultured GABAergic neurons exposed to different combinations of (13)C-labeled and unlabeled glucose and β-hydroxybutyrate. Depolarization was induced and the incorporation of (13)C into glutamate, GABA and aspartate was analyzed. The presence of β-hydroxybutyrate together with glucose did not affect the total GABA content but did, however, decrease the aspartate content to a lower value than when either glucose or β-hydroxybutyrate was employed alone. When combinations of the two substrates were used (13)C-atoms from β-hydroxybutyrate were found in all three amino acids to a greater extent than (13)C-atoms from glucose, but only the (13)C contribution from [1,6-(13)C]glucose increased upon depolarization. In conclusion, β-hydroxybutyrate was preferred over glucose as substrate for amino acid synthesis but the total content of aspartate decreased when both substrates were present. Furthermore only the use of glucose increased upon depolarization.

Ketone bodies serve as alternative energy substrates for the brain in cases of low glucose availability such as during starvation or in patients treated with a ketogenic diet. The ketone bodies are metabolized via a distinct pathway confined to the mitochondria. We have compared metabolism of [2,4-(13)C]beta-hydroxybutyrate to that of [1,6-(13)C]glucose in cultured glutamatergic neurons and investigated the effect of neuronal activity focusing on the aspartate-glutamate homeostasis, an essential component of the excitatory activity in the brain. The amount of (13)C incorporation and cellular content was lower for glutamate and higher for aspartate in the presence of [2,4-(13)C]beta-hydroxybutyrate as opposed to [1,6-(13)C]glucose. Our results suggest that the change in aspartate-glutamate homeostasis is due to a decreased availability of NADH for cytosolic malate dehydrogenase and thus reduced malate-aspartate shuttle activity in neurons using beta-hydroxybutyrate. In the presence of glucose, the glutamate content decreased significantly upon activation of neurotransmitter release, whereas in the presence of only beta-hydroxybutyrate, no decrease in the glutamate content was observed. Thus, the fraction of the glutamate pool available for transmitter release was diminished when metabolizing beta-hydroxybutyrate, which is in line with the hypothesis of formation of transmitter glutamate via an obligatory involvement of the malate-aspartate shuttle.

Udgivelsesdato: 2007-Dec The N-acylethanolamines (NAEs) and 2-arachidonoylglycerol (2-AG) are bioactive lipids that can modulate inflammatory responses and protect neurons against glutamatergic excitotoxicity. We have used a model of focal cerebral ischemia in young adult mice to investigate the relationship between focal cerebral ischemia and endogenous NAEs. Over the first 24 h after induction of permanent middle cerebral artery occlusion, we observed a time-dependent increase in all the investigated NAEs, except for anandamide. Moreover, we found an accumulation of 2-AG at 4 h that returned to basal level 12 h after induction of ischemia. Accumulation of NAEs did not depend on regulation of N-acylphosphatidylethanolamine-hydrolyzing phospholipase D or fatty acid amide hydrolase. Treatment with the fatty acid amide hydrolase inhibitor URB597 (cyclohexyl carbamic acid 3'-carbamoyl-biphenyl-3-yl ester; 1 mg/kg; i.p.) 1.5 h before arterial occlusion decreased the infarct volume in our model system. Our results suggest that NAEs and 2-AG may be involved in regulation of neuroprotection during focal cerebral ischemia in mice.

We tested the hypothesis that inhibiting glycogen degradation accelerates compound action potential (CAP) failure in mouse optic nerve (MON) during aglycemia or high-intensity stimulation. Axon function was assessed as the evoked CAP, and glycogen content was measured biochemically. Isofagomine, a novel inhibitor of central nervous system (CNS) glycogen phosphorylase, significantly increased glycogen content under normoglycemic conditions. When MONs were bathed in artificial cerebrospinal fluid (aCSF) containing 10 mM glucose, the CAP failed 16 min after exposure to glucose-free aCSF. MONs bathed in aCSF plus isofagomine displayed accelerated CAP failure on glucose removal. Similar results were obtained in MONs bathed in 30 mM glucose, which increased baseline glycogen concentration. The ability of isofagomine to increase glycogen content thus was not translated into delayed CAP failure. This is likely due to the inability of the tissue to metabolize glycogen in the presence of isofagomine, highlighting the importance of glycogen in sustaining neural function during aglycemia. The hypothesis that glycogen breakdown supports intense neural activity was tested by blocking glycogen breakdown during periods of high-frequency stimulation. The CAP area declined more rapidly when glycogen metabolism was inhibited by isofagomine, explicitly showing an important physiological role for glycogen metabolism during neural activity.

alpha-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPA-R)-mediated neurotoxicity was studied in relation to subunit expression and the presence of Ca(2+)-permeable receptor channels. AMPA-mediated toxicity had two components: 1) a direct AMPA-R-mediated component, which was not due to Ca(2+) influx through voltage-gated Ca(2+) channels, reversal of the Na(+)/Ca(2+) exchanger or release of calcium from dantrolene-sensitive intracellular Ca(2+) stores, and 2) a minor, indirect component involving activation of NMDA receptor channels, because of glutamate release and removal of the Mg(2+) block of the NMDA receptor on AMPA-R stimulation. The involvement of Ca(2+) influx through AMPA-R was also examined. The number of neurons possessing Ca(2+)-permeable AMPA-R increased during culture development, concurrently with an increasing susceptibility for AMPA-induced toxicity during development. GluR2(R) levels also increased during development, and channel blockers of Ca(2+)-permeable AMPA-R lacking the GluR2(R) subunit (spermine and philanthotoxin) failed to prevent neurotoxicity or increases in [Ca(2+)](i). Thus, the direct AMPA-R-mediated toxicity may be explained by initiation of cell death by Ca(2+) fluxing through AMPA-R containing GluR2(R). The components of direct AMPA-R-mediated toxicity are proposed to be 1) toxicity mediated by GluR2(R)-lacking AMPA-R and 2) toxicity mediated by low-Ca(2+)-permeability AMPA-R containing GluR2(R).

The physiological significance and subcellular distribution of voltage dependent calcium channels was defined using calcium channel blockers to inhibit potassium induced rises in cytosolic calcium concentration in cultured mouse neocortical neurons. The cytosolic calcium concentration was measured using the fluorescent calcium chelator fura-2. The types of calcium channels present at the synaptic terminal were determined by the inhibitory action of calcium channel blockers on potassium-induced [3H]GABA release in the same cell preparation. L-, N-, P-, Q- and R-/T-type voltage dependent calcium channels were differentially distributed in somata, neurites and nerve terminals. ω-conotoxin MVIIC (ω-CgTx MVIIC) inhibited approximately 40% of the Ca2+-rise in both somata and neurites and 60% of the potassium induced [3H]GABA release, indicating that the Q-type channel is the quantitatively most important voltage dependent calcium channel in all parts of the neuron. After treatment with thapsigargin the increase in cytosolic calcium was halved, indicating that calcium release from thapsigargin sensitive intracellular calcium stores is an important component of the potassium induced rise in cytosolic calcium concentration. The results of this investigation demonstrate that pharmacologically distinct types of voltage dependent calcium channels are differentially localized in cell bodies, neurites and nerve terminals of mouse cortical neurons but that the Q-type calcium channel appears to predominate in all compartments.

The heterocyclic analogue of (S)-glutamic acid, (S)-2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionic acid [(S)-AMPA] is a potent and selective AMPA receptor agonist, whereas the enantiomeric compound, (R)-AMPA, is virtually inactive. We have previously characterized (RS)-2-amino-3-(3-hydroxy-5-phenyl-4-isoxazolyl)propionic acid [(RS)-APPA] as a partial AMPA receptor agonist showing about 60% of the efficacy of (RS)-AMPA. This partial agonism produced by (RS)-APPA is, however, only apparent, since resolution of (RS)-APPA has now been shown to provide the full AMPA receptor agonist, (S)-APPA, whereas (R)-APPA is a non -N- methyl- d -aspartic acid (non-NMDA) receptor antagonist showing preferential AMPA blocking effects. In agreement with classical theories for competitive interaction between agonists and antagonists, the efficacy of depolarizations produced by (S)-APPA in the rat cortical wedge preparation was shown to be progressively reduced with increasing molar ratios of (R)-APPA/(S)-APPA. These compounds and the competitive antagonists (RS)-2-amino-3-(3-carboxymethoxy-5-methyl-4-isoxazolyl)propionic acid [(RS)-AMOA], 6-cyano-7-nitroquinoxalin-2,3-dione (CNQX) and 6-nitro-7-sulfamoylbenzo(f)quinoxalin-2,3-dione (NBQX) were also tested in [3H]AMPA and [3H]CNQX binding systems, the latter ligand being used in the absence or presence of thiocyanate ions. On the basis of these studies it is suggested that (RS)-AMPA and the AMPA agonist (S)-APPA interact with a high-affinity receptor conformation, whereas the competitive antagonists (RS)-AMOA and (R)-APPA, derived from these agonists, preferentially bind to a low-affinity AMPA receptor conformation. The competitive antagonists, CNQX and NBQX which are structurally unrelated to (RS)-AMPA or (RS)-APPA, do not seem to discriminate between these two AMPA receptor conformations. The modified [3H]CNQX binding assay containing thiocyanate ions was shown to provide receptor affinity data for AMPA receptor agonists as well as antagonists, which correlate with the potencies of these compounds in the cortical wedge preparation. Using autoradiographic techniques, (S)- and (R)-APPA were shown to exhibit significantly different absolute potencies as inhibitors of [3H]AMPA binding in a number of regions of the rat brain.

The (R) and (S) forms of 5-amino-2-hydroxyvaleric acid (2-OH-DAVA) and 5-amino-4-hydroxyvaleric acid (4-OH-DAVA) were designed as structural hybrids of the 4-aminobutyric acidB (GABAB) agonist (R)-(–)-4-amino-3-hydroxybutyric acid [(R)-(–)-3-OH-GABA] and the GABAB antagonist 5-aminovaleric acid (DAVA). (S)-(–)-2-OH-DAVA and (R)-(−)-4-OH-DAVA showed a moderately potent affinity for GABAB receptor sites in rat brain and showed GABAB antagonist effects in a guinea pig ileurn preparation. The respective enantiomers, (R)-(+)-2-OH-DAVA and (S)- (+)-4-OH-DAVA, were markedly weaker in both test systems. AH four compounds were weak inhibitors of GABAA receptor binding in rat brain, and none of them significantly affected synaptosomal GABA uptake. Based on molecular modeling studies it has been demonstrated that low-energy conformations of (R)-(–)-3-OH-GABA, (S)(–)-2-OH-DAVA, and (R)-(–)4-OH-DAVA can be superimposed. These conformations may reflect the shapes adopted by these conforma-tionally flexible compounds during their interaction with GABAB receptors. The present studies emphasize the similar, but distinct, constraints imposed on agonists and antagonists for GABAB receptors.