Skip to main content
Download PDF

In silico investigation of mitragynine and 7-hydroxymitragynine metabolism

BMC Research Notes volume 12, Article number: 451 (2019) Cite this article

Abstract

Objective

Mitragynine is the main active compound of Mitragyna speciose (Kratom in Thai). The understanding of mitragynine derivative metabolism in human body is required to develop effective detection techniques in case of drug abuse or establish an appropriate dosage in case of medicinal uses. This in silico study is based upon in vivo results in rat and human by Philipp et al. (J Mass Spectrom 44:1249–1261, 2009).

Results

Gas-phase structures of mitragynine, 7-hydroxymitragynine and their metabolites were obtained by quantum chemical method at B3LYP/6-311++G(d,p) level. Results in terms of standard Gibbs energies of reaction for all metabolic pathways are reported with solvation energy from SMD model. We found that 7-hydroxy substitution leads to changes in reactivity in comparison to mitragynine: position 17 is more reactive towards demethylation and conjugation with glucuronic acid and position 9 is less reactive towards conjugation with glucuronic acid. Despite the changes, position 9 is the most reactive for demethylation and position 17 is the most reactive for conjugation with glucuronic acid for both mitragynine and 7-hydroxymitragynine. Our results suggest that 7-hydroxy substitution could lead to different metabolic pathways and raise an important question for further experimental studies of this more potent derivative.

Introduction

Mitragynine is the alkaloid derived from Kratom (Mitragyna speciose), a plant commonly found in Thailand and throughout the South East Asia region [1, 2]. Like many other opioid plants, there are claims of medical uses [3, 4] but these plants are also potentially illegal drugs of abuse [5, 6]. In additional to natural sources, mitragynine and its derivative may be obtained by total syntheses reported by researchers in Japan [7] and United States [8]. One of the more potent but naturally occurring derivatives of mitragynine is 7-hydroxymitragynine [9]. To establish an appropriate dosage in case of medicinal uses or to develop detection techniques in case of drug abuse, the understanding of mitragynine derivative metabolism in human body is needed.

There were a number of experimental attempts to investigate metabolites of mitragynine and similar compounds in living organisms [10,11,12]. The most complete metabolic pathways of mitragynine were proposed from LC–MS study of rats and human urine samples by Philipp et al. [11] (see Fig. 1). Two different sample preparation techniques were used in Phase I and II metabolite extraction from urine. Missing intermediates were proposed and differences between compounds found in human and rats were attributed to physiological difference and/or differences in dosage/sampling time. The chemical reactions involved in the metabolic pathways are

Fig. 1
figure1

(Refer to Table 1 for abbreviated name convention. Compounds in brackets are not found experimentally but are assumed to be intermediates. Left, right and down arrows are for reactions at positions 17, 16 and 9 respectively. Metabolites found in rat or human are indicated by R or H respectively)

Metabolic pathways of mitragynine in rat and human was proposed by Philipp et al. [11] from LC-MS experiments in 2009

Full size image

  • Hydrolysis of the only methyl ester at position 16,

  • O-demethylation of the 9-methoxy group,

  • O-demethylation of the 17-methoxy group, which may proceed via an aldehyde intermediate and result in corresponding carboxylic acid or alcohol by oxidation or reduction respectively,

  • Subsequent conjugation with glucuronic acid at one of the three position above,

  • Subsequent conjugation with sulfate only at the 9-methoxy group, and

  • Combination of the steps above.

Figure 1 is the basis of this study and our aim is to complement the experimental findings with results from quantum chemical calculation. For comparison purpose, 7-hydroxymitragynine, a representative of naturally occurring and more potent derivatives of mitragynine is also included in our in silico investigation.

Conformers of mitragynine derivatives were theoretically studied by Liu et al. [13] using MMFF94s force field and B3LYP density functional theory method. Lowest energy structures were confirmed with crystal structures [14,15,16,17]. We used the lowest energy structures shown in Additional file 2: Figure S1 as the representative for our metabolic study. According to the figure, from left to right, the first and second rings are planar because of their aromaticity but the third and the fourth rings are both in chair conformer. The nitrogen atom between the third and fourth rings is above the plan of the molecule as described by Liu et al. [13].

Main text

All quantum chemical calculations were performed using the Q-Chem 5.1 program package [18]. Gas-phase structures were obtained at B3LYP/6-311++G(d,p) level and were confirmed to be minima on the electronic potential energy surface by frequency calculations. Solvation energy in water from SMD model [19] was obtained from gas-phase geometries. The energies were corrected to standard state in solution condition of 1 M at 298.15 K with an exception of water in which case 55.34 M was used. Shell script, spreadsheet template and Mathematica [20] notebook are modified from our previous work [21]. All output files and other associated codes to obtain the standard Gibbs energies of reaction are provided in Additional file 1. For reporting purpose, metabolites are coded by molecular formula and abbreviated names as shown in Fig. 1 and Table 1.

Table 1 Abbreviated compound names used in this study
Full size table

Permutations of all possible substitutions at position 9 (m,d,g,s), position 16 (m,d,g) and position 17 (m,d,o,r,G) leads to 4 × 3 × 5 = 60 possible metabolites, and (3 × 3 × 5) + (4 × 2 × 5) + (4 × 3 × 4) = 133 possible steps. However, here we only consider compounds shown in Fig. 1 (18 structures and 21 possible steps) and their 7-hydroxy derivatives.

To compute reaction energy from Fig. 1, relevant additional reagents (water, sulfate ion, hydronium ion, protonated nicotinamide [22], glucuronic acid, oxygen) and products (water, methanol, reduced nicotinamide) were also added to the scheme to complete the thermodynamic calculation. These compounds are commonly found in biological systems and are likely to be reasonable energy reference point for these demethylation, oxidation/reduction of aldehyde and conjugation with glucuronide/sulfate reactions.

All calculations were completed with no imaginary frequency. The lowest energy structures of mitragynine and 7-hydroxymitragynine are shown in Additional file 2: Fig. S1. The detailed result for all reactions is listed in Additional file 2: Table S1 by type of reaction, number of steps from the parent compound and by position of reaction (Additional files 1, 2)

  • Gas-phase Gibbs energies of reaction show similar trend to aqueous-phase energies. The differences between gas-phase and aqueous-phase energies for conjugation with sulfate and reduction reactions are considerably large due to the presence of charged species on the reactant and/or product sides.

  • The major determining factor for the energy of reaction is the type of reaction. Conjugation with sulfate, oxidation and reductions are highly exergonic and should occur easily. However, demethylation and conjugation with glucuronic acid have mixed results. The positions therefore play an important role in these cases to determine which pathway is more energetically favorable.

  • The average standard Gibbs energies of reaction for each reaction-position pair is shown in Additional file 2: Figure S2. As far as the three positions are concerned, position 9 is the most reactive for demethylation and position 17 is the most reactive for conjugation with glucuronic acid.

  • There are no conclusive trends from the number of steps from the parent compound.

  • The effect of 7-hydroxy substitution can be seen in Additional file 2: Figure S2 (decreases in d17, g17 and an increase in g9 average reaction energies) and Fig. 2 (lower position of H_mmd, H_dmd and H_mmG, significantly lower position of H_dmr and slightly higher position of H_gmr). Despite the decreases and increase, the general conclusion above that position 9 is the most reactive for demethylation and position 17 is the most reactive for conjugation with glucuronic acid are still true for both mitragynine and 7-hydroxymitragynine.

    Fig. 2
    figure2

    Relative standard Gibbs energies of metabolites of mitragynine (left) and 7-hydroxymitragynine (right) with reference to parent compounds, M_mmm and H_mmm. The energy of H_mmm is lower than that of M_mmm by 38.3 kcal/mol

    Full size image

Based on existing experimental evidence [11, 16], gas-phase structures of mitragynine, 7-hydroxymitragynine and their metabolites (in total 36 compounds) were obtained by quantum chemical method at B3LYP/6-311++G(d,p). Standard Gibbs energies of reaction in solution-phase were calculated by SMD solvation model. Four factors, type of reaction, number of steps from the parent compound, position of reaction and 7-hydroxy substitution effect were studied in this work. When compared to mitragynine, the 7-hydroxy substitution makes demethylation and conjugation with glucuronic acid at position 17 more favorable and makes conjugation with glucuronic acid at position 9 less favorable. Our results can be a basis for further experimental investigation of physiological effects of the compounds and/or detection of Kratom use.

Limitations

During the preparation of this paper, relevant work [23,24,25] emerged in the literature. Readers may be interested to see a recent in vivo pharmacokinetics study of mitragynine in rats [23] and in vitro conversion of mitragynine to 7-hydroxymitragynine [24, 25]. The importance of 7-hydroxymitragynine was raised as there is clear evidence that it can be formed by in vivo metabolism. This is different from the basic premise in this manuscript that 7-hydroxymitragynine is obtained from Kratom as a minor alkaloid component. However, this fact does not materially change our conclusion in the current study. Also, cytochrome P450 enzymes were identified as the key mediator of the process [23, 25] and molecular docking studying could provide a further insight into the process.

A recent paper [23] describes a putative metabolic pathways from mitragynine to 14 metabolites. Unlike [11] which we used as a basis of our study, structures of all metabolites except one (Met2) are not available in the paper. Metabolites in the paper [23] can be comparable to ours as follow: Met5 is M_dmd, Met2 is possibly H-mmm, Met1 is possibly H_dmd, and Met8 is possibly M_mmd or M_mdm.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its additional files.

References

  1. 1.

    Hassan Z, Muzaimi M, Navaratnam V, Yusoff NHM, Suhaimi FW, Vadivelu R, et al. From Kratom to mitragynine and its derivatives: physiological and behavioural effects related to use, abuse, and addiction. Neurosci Biobehav Rev. 2013;37(2):138–51.

    CAS  Article  Google Scholar 

  2. 2.

    Adkins JE, Boyer EW, McCurdy CR. Mitragyna speciosa, a psychoactive tree from Southeast Asia with opioid activity. Curr Top Med Chem. 2011;11(9):1165–75.

    CAS  Article  Google Scholar 

  3. 3.

    Takayama H, Kitajima M, Matsumoto K, Horie S. Indole alkaloid derivatives having opioid receptor agonistic effect, and therapeutic compositions and methods relating to same. U.S. Patent No. 8,648,090; 2014.

  4. 4.

    Pasternak G, Majumdar S, Karimov R, Varadi A. Mitragynine analogs and uses thereof. U.S. Patent Application No. 15/570,308; 2018.

  5. 5.

    Prozialeck WC. Update on the pharmacology and legal status of Kratom. J Am Osteopathic Assoc. 2016;116(12):802–9.

    Article  Google Scholar 

  6. 6.

    Nutt D, King LA, Saulsbury W, Blakemore C. Development of a rational scale to assess the harm of drugs of potential misuse. Lancet. 2007;369(9566):1047–53.

    Article  Google Scholar 

  7. 7.

    Takayama H, Maeda M, Ohbayashi S, Kitajima M, Sakai S, Aimi N. The first total synthesis of (−)-mitragynine, an analgesic indole alkaloid in mitragyna speciosa. Tetrahedron Lett. 1995;36(51):9337–40.

    CAS  Article  Google Scholar 

  8. 8.

    Ma J, Yin W, Zhou H, Liao X, Cook JM. General approach to the total synthesis of 9-methoxy-substituted indole alkaloids: synthesis of mitragynine, as well as 9-methoxygeissoschizol and 9-methoxy-Nb-methylgeissoschizol. J Org Chem. 2009;74(1):264–73.

    CAS  Article  Google Scholar 

  9. 9.

    Matsumoto K, Horie S, Ishikawa H, Takayama H, Aimi N, Ponglux D, et al. Antinociceptive effect of 7-hydroxymitragynine in mice: discovery of an orally active opioid analgesic from the Thai medicinal herb Mitragyna speciosa. Life Sci. 2004;74(17):2143–55.

    CAS  Article  Google Scholar 

  10. 10.

    de Moraes NV, Moretti RAC, Furr EB, McCurdy CR, Lanchote VL. Determination of mitragynine in rat plasma by LC–MS/MS: application to pharmacokinetics. J Chromatogr B. 2009;877(24):2593–7.

    Article  Google Scholar 

  11. 11.

    Philipp AA, Wissenbach DK, Zoerntlein SW, Klein ON, Kanogsunthornrat J, Maurer HH. Studies on the metabolism of mitragynine, the main alkaloid of the herbal drug Kratom, in rat and human urine using liquid chromatography-linear ion trap mass spectrometry. J Mass Spectrom. 2009;44(8):1249–61.

    CAS  Article  Google Scholar 

  12. 12.

    Vuppala PK, Jamalapuram S, Furr EB, McCurdy CR, Avery BA. Development and validation of a UPLC-MS/MS method for the determination of 7-hydroxymitragynine, a μ-opioid agonist, in rat plasma and its application to a pharmacokinetic study. Biomed Chromatogr. 2013;27(12):1726–32.

    CAS  Article  Google Scholar 

  13. 13.

    Liu H, McCurdy CR, Doerksen RJ. Computational study on the conformations of mitragynine and mitragynaline. J Mol Struct (Thoechem). 2010;945(1):57–63.

    CAS  Article  Google Scholar 

  14. 14.

    Takayama H. Chemistry and pharmacology of analgesic indole alkaloids from the rubiaceous plant, Mitragyna speciosa. Chem Pharm Bull. 2004;52(8):916–28.

    CAS  Article  Google Scholar 

  15. 15.

    Takayama H, Ishikawa H, Kurihara M, Kitajima M, Sakai S, Aimi N, et al. Structure revision of mitragynaline, an indole alkaloid in Mitragyna speciosa. Tetrahedron Lett. 2001;42(9):1741–3.

    CAS  Article  Google Scholar 

  16. 16.

    Carvalho P, Furr EB III, McCurdy C. (E)-Methyl 2-[(2S,3S,12bR)-3-ethyl-8-methoxy-1,2,3,4,6,7,12,12b-octahydroindolo[2,3-a]quinolizin-2-yl]-3-methoxyacrylate ethanol solvate. Acta Crystallogr E. 2009;65(6):o1441–2.

    CAS  Article  Google Scholar 

  17. 17.

    Zacharias DE, Rosenstein RD, Jeffrey GA. The structure of mitragynine hydroiodide. Acta Crystallogr A. 1965;18(6):1039–43.

    CAS  Article  Google Scholar 

  18. 18.

    Shao Y, Gan Z, Epifanovsky E, Gilbert ATB, Wormit M, Kussmann J, et al. Advances in molecular quantum chemistry contained in the Q-Chem 4 program package. Mol Phys. 2015;113(2):184–215.

    CAS  Article  Google Scholar 

  19. 19.

    Marenich AV, Cramer CJ, Truhlar DG. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J Phys Chem B. 2009;113(18):6378–96.

    CAS  Article  Google Scholar 

  20. 20.

    Wolfram Research Inc. Mathematica. Champaign: Wolfram Research, Inc.; 2019.

    Google Scholar 

  21. 21.

    Limpanuparb T, Roongruangsree P, Areekul C. A DFT investigation of the blue bottle experiment: E ∘ half-cell analysis of autoxidation catalysed by redox indicators. R Soc Open Sci. 2017;4:170708.

    Article  Google Scholar 

  22. 22.

    Belenky P, Bogan KL, Brenner C. NAD+ metabolism in health and disease. Trends Biochem Sci. 2007;32(1):12–9.

    CAS  Article  Google Scholar 

  23. 23.

    Kamble SH, Sharma A, King TI, León F, McCurdy CR, Avery BA. Metabolite profiling and identification of enzymes responsible for the metabolism of mitragynine, the major alkaloid of Mitragyna speciosa (kratom). Xenobiotica. https://doi.org/10.1080/00498254.2018.1552819.

  24. 24.

    Avery BA, Boddu SP, Sharma A, Furr EB, Leon F, Cutler SJ, et al. Comparative pharmacokinetics of mitragynine after oral administration of Mitragyna speciosa (Kratom) leaf extracts in rats. Planta Med. 2019;85(04):340–6.

    CAS  Article  Google Scholar 

  25. 25.

    Kruegel AC, Uprety R, Grinnell SG, Langreck C, Pekarskaya EA, Rouzic VL et al.7-Hydroxymitragynine is an active metabolite of mitragynine and a key mediator of its analgesic effects. ACS Cent Sci. 2019;5(6):992–1001.

    Google Scholar 

Download references

Acknowledgements

We thank MUIC Information Technology Section for permission to use the server room. We are grateful for an opportunity to present this work at PACCON 2019 and to receive the outstanding poster presentation award. We appreciate Ms. Yenning Lee for her preliminary work on Mitragynine with TL and RN.

Funding

TL received research Grant from the Institute for the Promotion of Teaching Science and Technology (IPST 08/2557). RN received scholarship from the Thailand Research Fund (Royal Golden Jubilee PhD). These funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Affiliations

  1. Science Division, Mahidol University International College, Mahidol University, Nakhon Pathom, 73170, Thailand

    Taweetham Limpanuparb

  2. Division of Chemistry, Faculty of Science and Technology, Thammasat University, Pathum Thani, 12120, Thailand

    Rattha Noorat & Yuthana Tantirungrotechai

Authors
  1. Taweetham Limpanuparb
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Rattha Noorat
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yuthana Tantirungrotechai
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

TL conceived the idea, performed calculations and prepared the manuscript. RN assisted TL in quantum chemical calculations and manuscript preparation. YT provided perspectives on literature review and data analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Taweetham Limpanuparb.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1.

All output files and other associated codes to obtain the standard Gibbs energies of reaction are provided.

Additional file 2.

Table S1, Figures S1 and S2 for the molecular structures and energies of reaction.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Limpanuparb, T., Noorat, R. & Tantirungrotechai, Y. In silico investigation of mitragynine and 7-hydroxymitragynine metabolism. BMC Res Notes 12, 451 (2019). https://doi.org/10.1186/s13104-019-4461-3

Download citation

Keywords

  • Density functional theory
  • Kratom
  • Metabolism
  • Mitragynine
  • 7-Hydroxymitragynine

BMC Research Notes

ISSN: 1756-0500

Contact us