Modelling the pharmacokinetics of tramadol: on the difference between CYP2D6 extensive and poor metabolizers

In this work we propose a mathematical model for the kinetics of tramadol, a synthetic opioid commonly used for treating moderate to severe pain. This novel theoretical framework could result in an objective criterion on how to adjust the assigned dose, depending on the genetic polymorphisms of CYP2D6. The model describes the coupled dynamics of tramadol and the metabolite O-desmethyltramadol. The effect of diffusion of the drug in the blood is here accounted for and we further hypothesize the existence of a time delay in the process of chemical translation from tramadol into metabolites. The system of coupled differential equations is solved numerically and the free parameters adjusted so to interpolate the experimental time series for the intravenous injection setting. Theoretical curves are shown to reproduce correctly the experimental profiles obtained from clinical trials. This enables in turn to extract an estimate of the metabolization rate. A difference in metabolization rate between CYP2D6 poor and extensive metabolizers is also found, and the stereoselectivity in the O-demethylation of tramadol highlighted. Our results allow one to quantify the dose of (+)-tramadol (resp. (-)-tramadol) administered to poor or extensive metabolizers, if the same effect is sought. The latter is here quantified through the blood concentration of (+)-metabolites (resp. (-)-metabolites).     

Miotic action of tramadol is determined by CYP2D6 genotype

Polymorphic CYP2D6 is the enzyme that activates the opioid analgesic tramadol by O-demethylation to its active metabolite O-demethyltramadol (M1). Our objective was to determine the opioid effects measured by pupillary response to tramadol of CYP2D6 genotyped volunteers in relation to the disposition of tramadol and M1 in plasma. Tramadol displayed phenotypic pharmacokinetics and it was possible to identify poor metabolizers (PM) with >99% confidence from the metabolic ratio (MR) in a single blood sample taken between 2.5 and 24 h post-dose. Homozygous extensive metabolizers (EM) differed from PM subjects by an almost threefold greater (P=0.0014) maximal pupillary constriction (Emax). Significant correlations between the AUC and Cmax values of M1 versus pupillary constriction were found. The corresponding correlations of pharmacokinetic parameters for tramadol itself were weaker and negative. The strongest correlations were for the single-point metabolic ratios at all sampling intervals versus the effects, with rs ranging from 0.85 to 0.89 (p<0.01). It is concluded that the concept of dual opioid/non-opioid action of the drug, though considerably stronger in EMs, is valid for both EM and PM subjects. This is the theoretical basis for the frequent use and satisfactory efficacy of tramadol in clinical practice when given to genetically non-selected population.     

Recognition of functional genetic polymorphism using ESE motif definition: a conservative evolutionary approach to CYP2D6/CYP2C19 gene variants

Genetica - Although predicting the effects of variants near intron-exon boundaries is relatively straightforward, predicting the functional Exon Splicing Enhancers (ESEs) and the possible effects...     

CYP1A1 and CYP2D6 polymorphism and risk of lung cancer in a North Indian population

This case–control study was conducted to examine the association between the CYP1A1 and CYP2D6 genotypes and lung cancer risk among North Indians. The estimated relative risk for lung cancer associated with the CYP1A1 Val/Val allele was 2.68, and was four-fold when cases with small cell lung cancer (SCLC) were considered alone. With regard to the metabolism of debrisoquine, no poor metabolizers were found amongst the subjects. The odds ratio of risk with the heterozygous extensive metabolizer (HEM) genotype was 1.5. However, in the presence of at least a single copy of the variant CYP1A1 MspI allele and the CYP2D6 HEM genotype, the risk was two-fold for squamous cell carcinoma (SQCC). When the CYP1A1 Val/Val and CYP2D6 HEM genotypes were taken together, the risk for SCLC was four-fold. Stratified analysis indicated an interaction between bidi smoking and variant CYP1A1 genotypes on the risk for SQCC and SCLC. Heavy smokers (Brinkman index>400) with Val/Val genotypes were at a very high risk of developing lung cancer (odds ratio 29.30, 95% confidence interval 2.42–355, p=0.008). Heavy smokers with CYP1A1 MspI (CYP1A1*1/2A or CYP1A1*2A/*2A) genotype had a seven-fold risk for SCLC compared with non-smokers. This study is the first to be carried out on a North Indian population, and, although small, suggests that CYP1A1 and CYP2D6 polymorphisms might have a role in determining the risk for lung cancer and should be investigated further.     

CYP1A1 and CYP2D6 polymorphism and risk of lung cancer in a North Indian population.

This case-control study was conducted to examine the association between the CYP1A1 and CYP2D6 genotypes and lung cancer risk among North Indians. The estimated relative risk for lung cancer associated with the CYP1A1 Val/Val allele was 2.68, and was four-fold when cases with small cell lung cancer (SCLC) were considered alone. With regard to the metabolism of debrisoquine, no poor metabolizers were found amongst the subjects. The odds ratio of risk with the heterozygous extensive metabolizer (HEM) genotype was 1.5. However, in the presence of at least a single copy of the variant CYP1A1 MspI allele and the CYP2D6 HEM genotype, the risk was two-fold for squamous cell carcinoma (SQCC). When the CYP1A1 Val/Val and CYP2D6 HEM genotypes were taken together, the risk for SCLC was four-fold. Stratified analysis indicated an interaction between bidi smoking and variant CYP1A1 genotypes on the risk for SQCC and SCLC. Heavy smokers (Brinkman index>400) with Val/Val genotypes were at a very high risk of developing lung cancer (odds ratio 29.30, 95% confidence interval 2.42-355, p=0.008). Heavy smokers with CYP1A1 MspI (CYP1A1*1/2A or CYP1A1*2A/*2A) genotype had a seven-fold risk for SCLC compared with non-smokers. This study is the first to be carried out on a North Indian population, and, although small, suggests that CYP1A1 and CYP2D6 polymorphisms might have a role in determining the risk for lung cancer and should be investigated further.     

CYP2D6 polymorphism and the presence of anti-LKM-1 in patients with chronic hepatitis C

Anti-LKM-1 autoantibodies are directed mostly at cytochrome P450 2D6 (CYP2D6) autoantigen, whose activity ranges from "complete deficiency" to "extensive metabolism" due to genetic polymorphism. We aimed to find any relevance of CYP2D6 alleles to the presence/absence of anti-LKM-1 in Japanese patients with chronic hepatitis C. The frequency of an extensive metabolizer-type allele (CYP2D6*1) in anti-LKM-1-positive patients was higher than that in anti-LKM-1-negative patients (0.800 vs 0.431; P = 0.0035), while the CYP2D6*10 allele with moderately reduced activity was less frequent in the former than the latter (0.050 vs 0.389; P = 0.0069). Moreover, the rate of homozygosity for CYP2D6*1 showed a striking difference between the two groups (70% vs 19%; P = 0.0021). These findings suggest that a genetic predisposition to produce the enzyme CYP2D6 of extensive metabolizer-type is associated with the induction of anti-LKM-1 in chronic hepatitis C patients.     

Regio- and stereoselective oxidation of propranolol enantiomers by human CYP2D6, cynomolgus monkey CYP2D17 and marmoset CYP2D19

Toxic and pharmacokinetic profiles of drug candidates are evaluated in vivo often using monkeys as experimental animals, and the data obtained are extrapolated to humans. Well understanding physiological properties, including drug-metabolizing enzymes, of monkeys should increase the accuracy of the extrapolation. The present study was performed to compare regio- and stereoselectivity in the oxidation of propranolol (PL), a chiral substrate, by cytochrome P450 2D (CYP2D) enzymes among humans, cynomolgus monkeys and marmosets. Complimentary DNAs encoding human CYP2D6, cynomolgus monkey CYP2D17 and marmoset CYP2D19 were cloned, and their proteins expressed in a yeast cell expression system. The regio- and stereoselective oxidation of PL enantiomers by yeast cell microsomal fractions were compared. In terms of efficiency of expression in the system, the holo-proteins ranked CYP2D6 ≒ CYP2D17 ? CYP2D19. This may be caused by the bulky side chain of the amino acid residue at position 119 (leucine for CYP2D19 vs. valine for CYP2D6 and CYP2D17), which can disturb the incorporation of the heme moiety into the active-site cavity. PL enantiomers were oxidized by all of the enzymes mainly into 4-hydroxyproranolol (4-OH-PL), followed by 5-OH-PL and N-desisopropylpropranolol (NDP). In the kinetic analysis, apparent K_m values were commonly in the μM range and substrate enantioselectivity of R-PL < S-PL was observed in both K_m and V_max values for the formation of the three metabolites from PL enantiomers. The activity to produce NDP tended to be higher for the monkey enzymes, particularly CYP2D17, than for the human enzyme. These results indicate that in the oxidation of PL enantiomers by CYP2D enzymes, stereoselectivity is similar but regioselectivity is different between humans and monkeys.     

Clinical and genetic determinants of warfarin pharmacokinetics and pharmacodynamics during treatment initiation

Variable warfarin response during treatment initiation poses a significant challenge to providing optimal anticoagulation therapy. We investigated the determinants of initial warfarin response in a cohort of 167 patients. During the first nine days of treatment with pharmacogenetics-guided dosing, S-warfarin plasma levels and international normalized ratio were obtained to serve as inputs to a pharmacokinetic-pharmacodynamic (PK-PD) model. Individual PK (S-warfarin clearance) and PD (I(max)) parameter values were estimated. Regression analysis demonstrated that CYP2C9 genotype, kidney function, and gender were independent determinants of S-warfarin clearance. The values for I(max) were dependent on VKORC1 and CYP4F2 genotypes, vitamin K status (as measured by plasma concentrations of proteins induced by vitamin K absence, PIVKA-II) and weight. Importantly, indication for warfarin was a major independent determinant of I(max) during initiation, where PD sensitivity was greater in atrial fibrillation than venous thromboembolism. To demonstrate the utility of the global PK-PD model, we compared the predicted initial anticoagulation responses with previously established warfarin dosing algorithms. These insights and modeling approaches have application to personalized warfarin therapy.     

Sugar-mediated lyoprotection of purified human CYP3A4 and CYP2D6

P450 enzymes are of great interest for drug metabolism and as potential biocatalysts. Like most P450s, purified CYP3A4 is normally handled and stored in solution because lyophilization greatly reduces its activity. We show here that colyophilization of this enzyme with sucrose or trehalose, but not mannitol, crown ethers or cyclodextrins, allow recovery of full enzymatic activity after rehydration. Sorbitol was almost as efficient, with 85% retention of the original activity. We also show that similar protection is observed through colyophilization of CYP2D6 with trehalose. This procedure should greatly facilitate handling, storage, or use of these enzymes in anhydrous media.     

CYP2D6 and CYP1A1 mutations in the Turkish population

Drugs and carcinogens are substrates of a group of metabolic enzymes including cytochrome p450 enzymes and gluthatione S-transferases. Many of the genes encoding these enzymes exhibit functional polymorphisms that contribute individual cancer susceptibility and drug response. Molecular studies based on these polymorphic enzymes also explain the aetiology of cancer and therapeutic management in clinics. We analysed the cytochrome p4501A1 (CYP1A1) and 2D6 (CYP2D6) variant genotype and allele frequencies by PCR-RFLP in Turkish individuals (n=140). The frequency of the CYP1A1*2A mutant allele was found to be 15.4%, and the CYP2D6*3 and *4 mutant allele (poor metabolizer) frequencies were 2.5% and 13.9%, respectively. This study presents the first results of CYP1A1 and CYP2D6 mutant allele distributions in the Turkish population and these data provide an understanding of epidemiological studies that correlate therapeutic approaches and aetiology of several types of malignancy in Turkish patients.     

The decreased in vivo clearance of CYP2D6 substrates by CYP2D6*10 might be caused not only by the low-expression but also by low affinity of CYP2D6

CYP2D6 exhibits genetic polymorphism with interindividual differences in metabolic activity. We have found a significant influence on the pharmacokinetics of venlafaxine by the CYP2D6*10 allele in a Japanese population. CYP2D6.10, which is translated from CYP2D6*10, has two amino acid substitutions: Pro34 --> Ser and Ser486 --> Thr. In this study, CYP2D6.10 was expressed in Saccharomyces cerevisiae and its catalytic activity for CYP2D6 substrates was investigated. The CYP2D6*10B- and *10C-associated cDNA were isolated from human lymphocyte genotyped as CYP2D6*10. In addition, three forms of CYP2D6, Pro34/Thr486 (PT), Ser34/Ser486 (SS), and Pro34/Ser486 (wild type, CYP2D6.1), were constructed by PCR-site mutagenesis to clarify the effects of the two amino-acid substitutions. The expression of CYP2D6 protein was confirmed by immunoblotting using CYP2D antibody. The absorbance at 450 nm was measured by CO-reduced difference spectra from five all microsome preparations. The CYP2D6 forms with Pro34 --> Ser amino acid substitution were at a lower expression than CYP2D6.1 from the findings of immunoblotting and spectral analysis. The apparent K(m) values of CYP2D6.1, CYP2D6.10A, and CYP2D6.10C were 1.7, 8.5, and 49.7 microM, respectively, for bufuralol 1'-hydroxylation, and 9.0, 51.9, and 117.4 microM, respectively, for venlafaxine O-demethylation, respectively. The V(max) values were not significantly different among the three variants. These findings suggest that the decreased in vivo clearance by CYP2D6*10 was caused not only by low expression of but also the increased K(m) value of CYP2D6.     

CYP2D6 polymorphisms in patients with porphyrias

The cytochrome P-450 (CYP) isoenzymes, a superfamily of heme proteins which are the terminal oxidases of the mixed function oxidases system, metabolize more than 70% of all clinically approved drugs. The highly polymorphic CYP2D6 isoform metabolizes more than 25% of most common drugs, and the phenotypes of the 70-plus allelic variants range from compromised to excessive enzymatic activity. Porphyrias are a group of inherited or acquired metabolic disorders of heme biosynthesis, due to a specific decrease in the activity of one of the enzymes of the heme pathway. Clinical signs and symptoms of porphyrias are frequently associated with exposure to precipitating agents, including clinically approved drugs. CYP enzymes, including CYP2D6, participate in the metabolism of some porphyrinogenic drugs, leading to the deregulation of heme biosynthesis. Considering that some of the drugs not recommended for use in porphyric patients are metabolized by CYP2D6, the presence of CYP2D6 polymorphisms in porphyric patients would influence the triggering of the disease when these individuals receive a precipitating agent that is metabolized by CYP2D6. To investigate CYP2D6 polymorphisms in porphyric patients, healthy Argentinean volunteers, porphyric patients, and a group of individuals with high levels of iron were studied. Results indicated that the CYP2D6*3 and CYP2D6*4 alleles, in particular, would be linked to the onset of disease. Predictive genotyping for CYP2D6 in porphyric patients holds promise as a method to improve the clinical efficacy of drug therapy and to personalize drug administration for these patients.     

CYP3A4 and CYP2D6 inhibitory activities of Indonesian medicinal plants

Thirty samples of Indonesian medicinal plants were analyzed for their capacity to inhibit in vitro metabolism by human cytochrome P450 3A4 (CYP3A4) and CYP2D6 with a radiometric assay. The MeOH-soluble fractions of 25 samples, prepared from water extracts, demonstrated inhibitory activity more than 50% on the metabolism mediated by CYP3A4, and 21 samples on the metabolism mediated by CYP2D6. Among the MeOH-soluble fractions, Piper nigrum leaf showed the highest inhibitory activity against CYP3A4 (91.7%), and Punica granatum against CYP2D6 (98.1%). The water extracts of which MeOH-soluble fraction showed inhibitory activity more than 70% were fractionated with EtOAc. From the EtOAc-soluble fractions, Curcuma heyneana (67.0%), Pi. cubeba (75.0%), Pi. nigrum fruit (84.0%), Pi. nigrum leaf (85.8%), and Zingiber aromaticum (75.3%) demonstrated inhibitory activity more than 50% on the metabolism mediated by CYP3A4, but only Pi. nigrum fruit (72.8%) and Pi. nigrum leaf (69.1%) showed strong inhibitory activity against CYP2D6. For samples that showed more than 70% inhibition, their IC(50) values were determined. The most potent inhibitory activity against CYP3A4 (IC(50) value of 25 microg/ml) was found for the extract of Pi. nigrum leaf, while that of Catharanthus roseus showed the most potent inhibitory effect against CYP2D6 (IC(50) value of 11 microg/ml). These results should indicate once more the possibility of potential medicinal plant-drug interactions.     

Molecular simulation study of CYP2B6 polymorphism with and without psoralen

CYP2B6 is a polymorphic enzyme with a large number of variants which may lead to functional changes in enzyme activity and substrate selectivity. In this study, CYP2B6 and its three variants with and without psoralen, a mechanism-based inactivator, were investigated using molecular simulation method. The obtained docking orientation of psoralen was in agreement with previously identified site of metabolism. Stability analysis showed that the three variants displayed more flexibility than CYP2B6.1, and CYP2B6.34 was the most flexible one without psoralen binding. However, in the presence of psoralen, CYP2B6.34 became more rigidity. Tunnel analysis indicates that the bottleneck change of tunnels may be correlated to the increased or decreased activity of variants. Binding free energy analysis shows that van der Waals interaction dominates the binding of psoralen. CYP2B6.34 has the highest affinity to psoralen with lowest binding free energy. Ile114, Phe115 and heme contribute largely to the binding of psoralen with CYP2B6.6, while Phe206 and Leu363 play important roles for CYP2B6.1 and CYP2B6.4. These computational observations suggest that the increased activity of CYP2B6.4 and reduced activity of CYP2B6.6 may be due to changes in regional structures.     

Genetic polymorphism of CYP2A6 in relation to cancer.

To clarify the roles of human cytochrome P450 (P450 or CYP) 2A6 and 2E1 on the metabolic activation of N-nitrosamines, we established genetically engineered Salmonella typhimurium strains harboring human CYP2A6 or CYP2E1 together with NADPH-P450 reductase (OR). The 5'-terminus of CYP cDNA was modified to achieve a high-level expression in S. typhimurium. Modified CYP2A6 or CYP2E1 cDNA and native OR cDNA were introduced into a pCW vector. S. typhimurium YG7108 cells were transformed with this vector. The mutagen producing ability of these enzymes for some N-nitrosamines were evaluated using the established S. typhimurium cells. We found that the substrate specificity of CYP2A6 and CYP2E1 was different among mutagens. CYP2A6 was responsible for the metabolic activation of N-nitrosamines possessing relatively long alkyl chains, whereas CYP2E1 was responsible for the metabolic activation of N-nitrosamines with relatively short alkyl chains. It is likely that CYP2A6 gene polymorphism is responsible for the interindividual variability on the cancer susceptibility. We found the whole deletion of CYP2A6 gene as a type of genetic polymorphism in Japanese. Thus, we developed a gene diagnosis method to detect the variant. We evaluated the relationship between the CYP2A6 gene whole deletion and the susceptibility to the lung cancer. The frequency of CYP2A6 gene whole deletion was significantly lower in the lung cancer patients than that of healthy volunteers.     

Rat CYP2D2, not 2D1, is functionally conserved with human CYP2D6 in endogenous morphine formation

The assumption that CYP2D1 is the corresponding rat cytochrome to human CYP2D6 has been revisited using recombinant proteins in direct enzyme assays. CYP2D1 and 2D2 were incubated with known CYP2D6 substrates, the three morphine precursors thebaine, codeine and (R)-reticuline. Mass spectrometric analysis showed that rat CYP2D2, not 2D1, catalyzed the 3-O-demethylation reaction of thebaine and codeine. In addition, CYP2D2 incubated with (R)-reticuline generated four products corytuberine, pallidine, salutaridine and isoboldine while rat CYP2D1 was completely inactive. This intramolecular phenol-coupling reaction follows the same mechanism as observed for CYP2D6. Michaelis-Menten kinetic parameters revealed high catalytic efficiencies for rat CYP2D2. These findings suggest a critical evaluation of other commonly accepted, however untested, CYP2D1 substrates.     

Enantioselective metabolism of primaquine by human CYP2D6

Given the threat of resistance of human malaria parasites, including to artemisinin derivatives, new agents are needed. Chloroquine (CQ) has been the most widely used anti-malarial, and new analogs (CQAns) presenting alkynes and side chain variations with high antiplasmodial activity were evaluated. Six diaminealkyne and diaminedialkyne CQAns were evaluated against CQ-resistant (CQ-R) (W2) and CQ-sensitive (CQ-S) (3D7) Plasmodium falciparum parasites in culture. Drug cytotoxicity to a human hepatoma cell line (HepG2) evaluated, allowed to calculate the drug selectivity index (SI), a ratio of drug toxicity to activity in vitro. The CQAns were re-evaluated against CQ-resistant and -sensitive P. berghei parasites in mice using the suppressive test. Docking studies with the CQAns and the human (Hss LDH) or plasmodial lactate dehydrogenase (Pf LDH) enzymes, and, a β-haematin formation assay were performed using a lipid as a catalyst to promote crystallization in vitro. All tested CQAns were highly active against CQ-R P. falciparum parasites, exhibiting half-maximal inhibitory concentration (IC50) values below 1 μΜ. CQAn33 and CQAn37 had the highest SIs. Docking studies revealed the best conformation of CQAn33 inside the binding pocket of Pf LDH; specificity between the residues involved in H-bonds of the Pf LDH with CQAn37. CQAn33 and CQAn37 were also shown to be weak inhibitors of Pf LDH. CQAn33 and CQAn37 inhibited β-haematin formation with either a similar or a 2-fold higher IC50 value, respectively, compared with CQ. CQAn37 was active in mice with P. berghei, reducing parasitaemia by 100%. CQAn33, -39 and -45 also inhibited CQ-resistant P. berghei parasites in mice, whereas high doses of CQ were inactive. The presence of an alkyne group and the size of the side chain affected anti-P. falciparum activity in vitro. Docking studies suggested a mechanism of action other than Pf LDH inhibition. The β-haematin assay suggested the presence of an additional mechanism of action of CQAn33 and CQAn37. Tests with CQAn34, CQAn37, CQAn39 and CQAn45 confirmed previous results against P. berghei malaria in mice, and CQAn33, 39 and 45 were active against CQ-resistant parasites, but CQAn28 and CQAn34 were not. The result likely reflects structure-activity relationships related to the resistant phenotype.     

High-throughput analysis of informative CYP2D6 compound haplotypes

We describe a high-throughput protocol for detecting key polymorphisms in the drug-metabolizing enzyme gene CYP2D6 and a number of linked microsatellites that is both fast and relatively inexpensive to perform. This approach employs GeneScan technology to enable a researcher to determine rapidly the status of seven simple nucleotide polymorphisms in CYP2D6 and also to assay repeat number variation at five closely linked dinucleotide microsatellite loci. The method requires only three PCRs and two GeneScan runs per sample. We anticipate that this will be of value to researchers in three different ways: (1) rapid discrimination of common CYP2D6 alleles, (2) high-resolution haplotyping for association studies involving chromosome 22q13.1 using microsatellite variation, and (3) generation of compound haplotypes for investigating the evolution of CYP2D6 variation. We also report compound haplotype frequencies for an Ashkenazi Jewish and a British sample.     

Involvement of CYP3A4 and CYP2D6 in the Metabolism of Haloperidol

1. Human cytochrome P450 (CYP) isoenzymes expressed in a human cell line were used to elucidate their involvement in the metabolism of haloperidol (HAL).2. It was found that CYP3A4 catalyzes the metabolism of HAL to HAL 1,2,3,6-tetrahydropyridine (HTP). HTP is further metabolized to HAL pyridinium (HP⁺) by both CYP3A4 and CYP2D6.3. CYP3A4 and CYP2D6 are also responsible for the N-dealkylation of HAL. The N-dealkylation of reduced HAL (RH) was observed, which is catalyzed by CYP3A4. In addition, CYP3A4 also catalyzes the oxidation of RH back to HAL.4. These results are discussed in terms of the metabolic interactions of HAL with other drugs and how this knowledge may be used to reduce the movement disorders induced by HAL.