Methadone pharmacokinetics in opioid agonist treatment: Influencing factors and clinical implications

BACKGROUND

A considerable inter-individual variability has been reported in the relationship between methadone doses applied and serum concentrations achieved in methadone maintenance treatment. However, the underlying causes for this variability are not fully understood.

OBJECTIVES

We investigated the influence of genetic, pathophysiological and pharmacological factors on serum methadone concentration-to-dose ratio (CDR) and discussed the clinical implications of the findings.

METHODS

We used data from two retrospective laboratory databases and a prospective cohort study to investigate the impact on methadone CDR of hepatic cytochrome P450 enzyme system (CYP) genetic polymorphisms, age, sex, concomitant medication, liver fibrosis and body mass index through linear mixed model analyses.

FINDINGS

A positive association was found between CDR and the homozygous CYP2B6*6 genotype, concurrent treatment with CYP3A4 inhibitors and body mass index. CDR was lower among women and during concomitant use of CYP inducers. CDR was not associated with age or the degree of liver fibrosis in our investigations.

CONCLUSIONS

This research work supports the need for individually tailored dosage considering the various factors that influence methadone CDR. The gained knowledge can contribute to reducing the risks associated with the treatment and optimizing the desired outcomes.

Impact of liver fibrosis and clinical characteristics on dose-adjusted serum methadone concentrations

BACKGROUND

There is limited knowledge on the causes of large variations in serum methadone concentrations and dose requirements.

OBJECTIVES

We investigated the impact of the degree of liver fibrosis on dose-adjusted steady-state serum methadone concentrations.

METHODS

We assessed the clinical and laboratory data of 155 Norwegian patients with opioid use disorder undergoing methadone maintenance treatment in outpatient clinics in the period 2016-2020. A possible association between the degree of liver fibrosis and dose-adjusted serum methadone concentration was explored using a linear mixed-model analysis.

RESULTS

When adjusted for age, gender, body mass index, and genotypes of CYP2B6 and CYP3A5, the concentration-to-dose ratio of methadone did not increase among the participants with liver fibrosis (Coefficient: 0.70; 95% CI: -2.16, 3.57; P: 0.631), even among those with advanced cirrhosis (-0.50; -4.59, 3.59; 0.810).

CONCLUSIONS

Although no correlation was found between the degree of liver stiffness and dose-adjusted serum methadone concentration, close clinical monitoring should be considered, especially among patients with advanced cirrhosis. Still, serum methadone measurements can be considered a supplement to clinical assessments, taking into account intra-individual variations.

Nonadherence by Serum Drug Analyses in Resistant Hypertension: 7-Year Follow-Up of Patients Considered Adherent by Directly Observed Therapy

Background Measurement of serum concentrations of drugs is a novelty found useful in detecting poor drug adherence in patients taking ≥2 antihypertensive agents. Regarding patients with treatment-resistant hypertension, we previously based our assessment on directly observed therapy. The present study aimed to investigate whether serum drug measurements in patients with resistant hypertension offer additional information regarding drug adherence, beyond that of initial assessment with directly observed therapy. Methods and Results Nineteen patients assumed to have true treatment-resistant hypertension and adherence to antihypertensive drugs based on directly observed therapy were investigated repeatedly through 7 years. Serum concentrations of antihypertensive drugs were measured by ultra-high-performance liquid chromatography-tandem mass spectrometry from blood samples taken at baseline, 6-month, 3-year, and 7-year visits. Cytochrome P450 polymorphisms, self-reported adherence and beliefs about medicine were performed as supplement investigations. Seven patients (37%) were redefined as nonadherent based on their serum concentrations during follow-up. All patients reported high adherence to medications. Nonadherent patients expressed lower necessity and higher concerns regarding intake of antihypertensive medication (P=0.003). Cytochrome P450 polymorphisms affecting metabolism of antihypertensive drugs were found in 16 patients (84%), 21% were poor metabolizers, and none were ultra-rapid metabolizers. Six of 7 patients redefined as nonadherent had cytochrome P450 polymorphisms, however, not explaining the low serum drug concentrations measured in these patients. Conclusions Our data suggest that repeated measurements of serum concentrations of antihypertensive drugs revealed nonadherence in one-third of patients previously evaluated as adherent and treatment resistant by directly observed therapy, thereby improving the accuracy of adherence evaluation. Registration URL: https://www.clinicaltrials.gov; unique identifier: NCT01673516.

Short- and long-term effects of body weight, calorie restriction, and gastric bypass on CYP1A2-, CYP2C19-, and CYP2C9 activity

Aim

Roux‐en‐Y gastric bypass (RYGB) may influence drug disposition due to surgery‐induced gastrointestinal alterations and/or subsequent weight loss. The objective was to compare short‐ and long‐term effects of RYGB and diet on the metabolic ratios of paraxanthine/caffeine (cytochrome P450 [CYP] 1A2 activity), 5‐hydroxyomeprazole/omeprazole (CYP2C19 activity) and losartan/losartan carboxylic acid (CYP2C9 activity), and cross‐sectionally compare these CYP‐activities with normal‐to‐overweight controls.

Methods

This trial included patients with severe obesity preparing for RYGB (n = 40) or diet‐induced (n = 41) weight loss, and controls (n = 18). Both weight loss groups underwent a 3‐week low‐energy diet (<1200 kcal/day, weeks 0‐3) followed by a 6‐week very‐low‐energy diet or RYGB (both <800 kcal/day, weeks 3‐9). Follow‐up time was 2 years, with four pharmacokinetic investigations.

Results

Mean ± SD weight loss from baseline was similar in the RYGB‐group (13 ± 2.4%) and the diet group (10.5 ± 3.9%) at week 9, but differed at year 2 (RYGB −30 ± 6.9%, diet −3.1 ± 6.3%). From weeks 0 to 3, mean (95% confidence interval [CI]) CYP2C19 activity similarly increased in both groups (RYGB 43% [16, 55], diet 48% [22, 60]). Mean CYP2C19 activity increased by 30% (2.6, 43) after RYGB (weeks 3‐9), but not in the diet‐group (between‐group difference −0.30 [−0.63, 0.03]). CYP2C19 activity remained elevated in the RYGB group at year 2. Baseline CYP2C19 activity was 2.7‐fold higher in controls compared with patients with obesity, whereas no difference was observed in CYP1A2 and CYP2C9 activities.

Conclusion

Our findings suggest that CYP2C19 activity is lower in patients with obesity and increases following weight loss. This may be clinically relevant for drug dosing. No clinically significant effect on CYP1A2 and CYP2C9 activities was observed.

Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis

Mechanical loading exerts a profound influence on bone density and architecture, but the exact mechanism is unknown. Our study shows that expression of the neurological transcriptional factor zinc finger of the cerebellum 1 (ZIC1) is markedly increased in trabecular bone biopsies in the lumbar spine compared with the iliac crest, skeletal sites of high and low mechanical stress, respectively. Human trabecular bone transcriptome analyses revealed a strong association between ZIC1 mRNA levels and gene transcripts characteristically associated with osteoblasts, osteocytes and osteoclasts. This supposition is supported by higher ZIC1 expression in iliac bone biopsies from postmenopausal women with osteoporosis compared with age-matched control subjects, as well as strongly significant inverse correlation between ZIC1 mRNA levels and BMI-adjusted bone mineral density (BMD) (Z-score). ZIC1 promoter methylation was decreased in mechanically loaded vertebral bone compared to unloaded normal iliac bone, and its mRNA levels correlated inversely with ZIC1 promoter methylation, thus linking mechanical stress to epigenetic control of gene expression. The findings were corroborated in cultures of rat osteoblast progenitors and osteoblast-like cells. This study demonstrates for the first time how skeletal epigenetic changes that are affected by mechanical forces give rise to marked alteration in bone cell transcriptional activity and translate to human bone pathophysiology.

Correlation of Body Weight and Composition With Hepatic Activities of Cytochrome P450 Enzymes

Obesity is associated with comorbidities of which pharmacological treatment is needed. Physiological changes associated with obesity may influence the pharmacokinetics of drugs, but the effect of body weight on drug metabolism capacity remains uncertain. The aim of this study was to investigate ex vivo activities of hepatic drug metabolizing CYP enzymes in patients covering a wide range of body weight. Liver biopsies from 36 individuals with a body mass index (BMI) ranging from 18 to 63 kg/m² were obtained. Individual hepatic microsomes were prepared and activities of CYP3A, CYP2B6, CYP2C8, CYP2D6, CYP2C9, CYP2C19 and CYP1A2 were determined. The unbound intrinsic clearance (CL_int,u) values for CYP3A correlated negatively with body weight (r = -0.43, p < 0.01), waist circumference (r = -0.47, p < 0.01), hip circumference (r = -0.51, p < 0.01), fat percent (r = -0.41, p < 0.05), fat mass (r = -0.48, p < 0.01) and BMI (r = -0.46, p < 0.01). Linear regression analysis showed that CL_int,u values for CYP3A decreased with 5% with each 10% increase in body weight (r² = 0.12, β = -0.558, p < 0.05). There were no correlations between body weight measures and CL_int,u values for the other CYP enzymes investigated. These results indicate reduced hepatic metabolizing capacity of CYP3A substrates in patients with increasing body weight.

Effect of CYP4F2, VKORC1, and CYP2C9 in Influencing Coumarin Dose: A Single-Patient Data Meta-Analysis in More Than 15,000 Individuals

The cytochrome P450 (CYP)4F2 gene is known to influence mean coumarin dose. The aim of the present study was to undertake a meta-analysis at the individual patients level to capture the possible effect of ethnicity, gene-gene interaction, or other drugs on the association and to verify if inclusion of CYP4F2*3 variant into dosing algorithms improves the prediction of mean coumarin dose. We asked the authors of our previous meta-analysis (30 articles) and of 38 new articles retrieved by a systematic review to send us individual patients' data. The final collection consists of 15,754 patients split into a derivation and validation cohort. The CYP4F2*3 polymorphism was consistently associated with an increase in mean coumarin dose (+9% (95% confidence interval (CI) 7-10%), with a higher effect in women, in patients taking acenocoumarol, and in white patients. The inclusion of the CYP4F2*3 in dosing algorithms slightly improved the prediction of stable coumarin dose. New pharmacogenetic equations potentially useful for clinical practice were derived.

Distinct DNA methylation profiles in bone and blood of osteoporotic and healthy postmenopausal women

DNA methylation affects expression of associated genes and may contribute to the missing genetic effects from genome-wide association studies of osteoporosis. To improve insight into the mechanisms of postmenopausal osteoporosis, we combined transcript profiling with DNA methylation analyses in bone. RNA and DNA were isolated from 84 bone biopsies of postmenopausal donors varying markedly in bone mineral density (BMD). In all, 2529 CpGs in the top 100 genes most significantly associated with BMD were analyzed. The methylation levels at 63 CpGs differed significantly between healthy and osteoporotic women at 10% false discovery rate (FDR). Five of these CpGs at 5% FDR could explain 14% of BMD variation. To test whether blood DNA methylation reflect the situation in bone (as shown for other tissues), an independent cohort was selected and BMD association was demonstrated in blood for 13 of the 63 CpGs. Four transcripts representing inhibitors of bone metabolism-MEPE, SOST, WIF1, and DKK1-showed correlation to a high number of methylated CpGs, at 5% FDR. Our results link DNA methylation to the genetic influence modifying the skeleton, and the data suggest a complex interaction between CpG methylation and gene regulation. This is the first study in the hitherto largest number of postmenopausal women to demonstrate a strong association among bone CpG methylation, transcript levels, and BMD/fracture. This new insight may have implications for evaluation of osteoporosis stage and susceptibility.

Changes of 5-hydroxymethylcytosine distribution during myeloid and lymphoid differentiation of CD34+ cells

Background

Hematopoietic stem cell renewal and differentiation are regulated through epigenetic processes. The conversion of 5-methylcytosine into 5-hydroxymethylcytosine (5hmC) by ten-eleven-translocation enzymes provides new insights into the epigenetic regulation of gene expression during development. Here, we studied the potential gene regulatory role of 5hmC during human hematopoiesis.

Results

We used reduced representation of 5-hydroxymethylcytosine profiling (RRHP) to characterize 5hmC distribution in CD34+ cells, CD4+ T cells, CD19+ B cells, CD14+ monocytes and granulocytes. In all analyzed blood cell types, the presence of 5hmC at gene bodies correlates positively with gene expression, and highest 5hmC levels are found around transcription start sites of highly expressed genes. In CD34+ cells, 5hmC primes for the expression of genes regulating myeloid and lymphoid lineage commitment. Throughout blood cell differentiation, intragenic 5hmC is maintained at genes that are highly expressed and required for acquisition of the mature blood cell phenotype. Moreover, in CD34+ cells, the presence of 5hmC at enhancers associates with increased binding of RUNX1 and FLI1, transcription factors essential for hematopoiesis.

Conclusions

Our study provides a comprehensive genome-wide overview of 5hmC distribution in human hematopoietic cells and new insights into the epigenetic regulation of gene expression during human hematopoiesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0070-8) contains supplementary material, which is available to authorized users.

Copy number variation in the ATP-binding cassette transporter ABCC6 gene and ABCC6 pseudogenes in patients with pseudoxanthoma elasticum

Single mutations in the ATP-binding cassette transporter (ABCC6) gene (OMIM 603234) are known to cause the rare autosomal recessive disease pseudoxanthoma elasticum (PXE). Recently, we have found that copy number variations (CNVs) in pseudogenes of the ABCC6 gene are quite common. The aim of this study was to investigate the frequency and possible contribution of CNV in ABCC6 and its pseudogenes in PXE. Genomic DNA from 212 PXE individuals were examined for copy number by pyrosequencing and quantitative polymerase chain reaction (PCR) and compared with healthy individuals. The frequency of PXE individuals with any CNV was higher than in healthy individuals. The majority of variation comprised known and possibly new deletions in the ABCC6 gene and duplications of the ABCC6P1 and ABCC6P2 genes. ABCC6 deletions and ABCC6P2 duplications were not observed in 142 healthy individuals. In conclusion, by pyrosequencing and quantitative PCR, we were able to detect known and possibly new deletions in the ABCC6 gene that may have caused the PXE phenotype. Pyrosequencing may be used in PXE patients who have obtained incomplete genotype from conventional techniques. The frequency of ABCC6P2 pseudogene duplication was more common in PXE patients than healthy individuals and may affect the PXE phenotype.

RNA-sequencing analysis of HepG2 cells treated with atorvastatin

The cholesterol-lowering drug atorvastatin is among the most prescribed drug in the world. Alternative splicing in a number of genes has been reported to be associated with variable statin response. RNA-seq has proven to be a powerful technique for genome-wide splice variant analysis. In the present study, we sought to investigate atorvastatin responsive splice variants in HepG2 cells using RNA-seq analysis to identify novel candidate genes implicated in cholesterol homeostasis and in the statin response. HepG2 cells were treated with 10 µM atorvastatin for 24 hours. RNA-seq and exon array analyses were performed. The validation of selected genes was performed using Taqman gene expression assays. RNA-seq analysis identified 121 genes and 98 specific splice variants, of which four were minor splice variants to be differentially expressed, 11 were genes with potential changes in their splicing patterns (SYCP3, ZNF195, ZNF674, MYD88, WHSC1, KIF16B, ZNF92, AGER, FCHO1, SLC6A12 and AKAP9), and one was a gene (RAP1GAP) with differential promoter usage. The IL21R transcript was detected to be differentially expressed via RNA-seq and RT-qPCR, but not in the exon array. In conclusion, several novel candidate genes that are affected by atorvastatin treatment were identified in this study. Further studies are needed to determine the biological significance of the atorvastatin responsive splice variants that have been uniquely identified using RNA-seq.

A novel 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR) splice variant with an alternative exon 1 potentially encoding an extended N-terminus

Background

The major rate-limiting enzyme for de novo cholesterol synthesis is 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR). HMGCR is sterically inhibited by statins, the most commonly prescribed drugs for the prevention of cardiovascular events. Alternative splicing of HMGCR has been implicated in the control of cholesterol homeostasis. The aim of this study was to identify novel alternatively spliced variants of HMGCR with potential physiological importance.

Results

Bioinformatic analyses predicted three novel HMGCR transcripts containing an alternative exon 1 (HMGCR-1b, -1c, -1d) compared with the canonical transcript (HMGCR-1a). The open reading frame of the HMGCR-1b transcript potentially encodes 20 additional amino acids at the N-terminus, compared with HMGCR-1a. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) was used to examine the mRNA levels of HMGCR in different tissues; HMGCR-1a was the most highly expressed variant in most tissues, with the exception of the skin, esophagus, and uterine cervix, in which HMGCR-1b was the most highly expressed transcript. Atorvastatin treatment of HepG2 cells resulted in increased HMGCR-1b mRNA levels, but unaltered proximal promoter activity compared to untreated cells. In contrast, HMGCR-1c showed a more restricted transcription pattern, but was also induced by atorvastatin treatment.

Conclusions

The gene encoding HMGCR uses alternative, mutually exclusive exon 1 sequences. This contributes to an increased complexity of HMGCR transcripts. Further studies are needed to investigate whether HMGCR splice variants identified in this study are physiologically functional.

Copy number variations of the ATP-binding cassette transporter ABCC6 gene and its pseudogenes

BACKGROUND

The ATP-binding cassette transporter ABCC6 gene is located on chromosome 16 between its two pseudogenes (ABCC6P1 and ABCC6P2). Previously, we have shown that ABCC6P1 is transcribed and affects ABCC6 at the transcriptional level. In this study we aimed to determine copy number variations of ABCC6, ABCC6P1 and ABCC6P2 in different populations. Moreover, we sought to study the transcription pattern of ABCC6 and ABCC6 pseudogenes in 39 different human tissues.

FINDINGS

Genomic DNA from healthy individuals from five populations, Chinese (n = 24), Middle East (n = 20), Mexicans (n = 24), Caucasians (n = 50) and Africans (n = 24), were examined for copy number variations of ABCC6 and its pseudogenes by pyrosequencing and quantitative PCR. Copy number variation of ABCC6 was very rare (2/142; 1.4%). However, one or three copies of ABCC6P1 were relatively common (3% and 8%, respectively). Only one person had a single copy of ABCC6P2 while none had three copies. In Chinese, deletions or duplications of ABCC6P1 were more frequent than in any other population (9/24; 37.5%). The transcription pattern of ABCC6P2 was highly similar to ABCC6 and ABCC6P1, with highest transcription in liver and kidney. Interestingly, the total transcription level of pseudogenes, ABCC6P1 + ABCC6P2, was higher than ABCC6 in most tissues, including liver and kidney.

CONCLUSIONS

Copy number variations of the ABCC6 pseudogenes are quite common, especially in populations of Chinese ancestry. The expression pattern of ABCC6P2 in 39 human tissues was highly similar to that of ABCC6 and ABCC6P1 suggesting similar regulatory mechanisms for ABCC6 and its pseudogenes.

Reduced platelet function and role of drugs in acute gastrointestinal bleeding

Gastrointestinal (GI) bleeding may be caused by a constitutive bleeding disposition or drug-induced inhibition of hemostasis. Platelet function in patients with ongoing GI bleeding is unknown. The aim of this study was to investigate platelet function in patients with acute GI bleeding. Patients (n=35) presenting with acute GI bleeding (hematemesis or melena) were recruited. For comparison, 13 patients treated with aspirin and 11 patients treated with clopidogrel without GI bleeding and 27 healthy controls were studied. Platelet function was measured by whole-blood aggregation and flow cytometry. Coagulation function was measured with calibrated automated thrombography. Platelet aggregation and P-selectin expression were significantly lower after arachidonic acid stimulation in GI bleeding patients than in healthy subjects (p≤0.05). Collagen-induced P-selectin expression was significantly reduced in patients using anti-platelet drugs (p=0.02) and in many patients not using anti-platelet drugs. Thrombin generation, measured by calibrated automated thrombography, was only reduced in patients on warfarin treatment. In conclusion, platelet function is reduced in acute GI bleeding patients and a considerable proportion appears to be related to drug use.

No link between X chromosome inactivation pattern and simple goiter in females: evidence from a twin study

BACKGROUND

Simple goiter (SG) comprises diffuse (DG) and nodular (NG) benign nonautoimmune nontoxic goiter. In nonendemic goiter areas, the ratio of females to males may exceed 5:1, indicating that gender and/or sex hormones may play a role in the etiology of SG in these areas. Theoretically, as shown for autoimmune thyroid disease, X chromosome inactivation (XCI) and resultant tissue chimerism could offer a novel explanation for the female preponderance of SG. To examine whether skewed XCI is associated with SG, we first compared XCI in 71 twin individuals with SG with that in 142 unrelated healthy control twin individuals, and then performed a within-pair comparison of XCI in 48 twin pairs discordant for SG.

METHODS

DNA was extracted from peripheral blood cells. XCI analysis was performed by predigestion of DNA using the methylation-sensitive enzyme Hpall, followed by polymerase chain reaction of the polymorphic CAG repeat of the androgen receptor gene. A polymerase chain reaction product is obtained from the inactive X chromosome only. The XCI pattern was classified as skewed when 80% or more of the cells preferentially inactivated the same X chromosome. Twin zygosity was established by DNA fingerprinting.

RESULTS

The frequency of skewed XCI in female twins with SG, DG, and NG was 11% (8/71), 13% (6/46), and 8% (2/25), respectively, which was not significantly different from the prevalences in the corresponding control populations, 14% (20/142, p = 0.56), 14% (13/92, p = 1.00), and 14% (7/50, p = 0.71), respectively. Essentially, similar results were obtained when comparing the prevalence of skewed XCI in twin pairs discordant for SG (48 pairs), DG (30 pairs), and NG (18 pairs).

CONCLUSION

In a sample of Danish female twins, we did not find evidence for involvement of skewed XCI in the etiology or the female preponderance of SG.

Warfarin dose and INR related to genotypes of CYP2C9 and VKORC1 in patients with myocardial infarction

BACKGROUND

Warfarin treatment has a narrow therapeutic range, requiring meticulous monitoring and dosage titration. Individual dosage requirement has recently partly been explained by genetic variation of the warfarin metabolizing enzyme CYP2C9 and the Vitamin K-activating enzyme VKORC1. In the WARIS-II study, comparing three different antithrombotic regimens after myocardial infarction, warfarin treatment reduced thrombotic events, but was associated with more frequent bleeding than use of acetylsalisylic acid (ASA) alone.

AIMS

The primary aim of the present study was to investigate the relation between genotypes of CYP2C9 and VKORC1 and warfarin maintenance dose in myocardial infarction. The secondary aim was to relate the genotypes to international normalized ratio (INR).

METHODS

Genotyping was performed in 212 myocardial infarction patients from the WARIS-II study by robotic isolation of DNA from EDTA whole blood (MagNa Pure LC) before PCR amplification (LightCycler) and melting point analysis.

RESULTS

The 420 C>T substitution of CYP2C9*2, the 1075 A>C substitution of CYP2C9*3 and the 1173 C>T substitution of VKORC1 had minor allele frequencies of, 11.3%, 5.7% and 36.6% respectively. Warfarin weekly dose varied between 17 mg and 74 mg among the patients. INR did not vary between genotypes. Warfarin dosage requirement was significantly associated with CYP2C9 and VKORC1 genotypes, treatment group and age. The VKORC1 genotype contributed 24.5% to the interindividual variation in warfarin dosage, whereas the combined CYP2C9 genotypes were only responsible for 7.2% of the dose variation.

CONCLUSION

CYP2C9 and VKORC1 genotype frequencies in myocardial infarction patients appear similar to other patient groups and have similar impact on warfarin maintenance dose.