Document Type : Original Research Article


1 University of Novi Sad, Faculty of Sciences, Department of Chemistry, Biochemistry and Environmental Protection, Novi Sad, Serbia

2 University of Novi Sad, Faculty of Sciences, Department of Physics, Novi Sad, Serbia


Some of the most important properties of ephedrine (EPH) molecule, related to its stability, reactivity and degradation have been computationally addressed in this work. The study has been performed based on the combination of classical and quantum-mechanical calculations to yield the lowest energy conformer. Furthermore, the most important intramolecular noncovalent interactions have been identified and their influence on the structural features have been analyzed. Energy and distribution of the frontier molecular orbitals have been analyzed in order to address the stability and charge transfer properties. The average local ionization energy descriptor has been used in order to identify the molecular sites that are sensitive towards the electrophilic attacks. Bond dissociation energies for hydrogen abstraction have been calculated in order to identify the sites of EPH sensitive towards the autoxidation mechanism. Interactions of EPH have been also investigated using the molecular dynamics simulations and calculations of radial distribution functions.

Graphical Abstract

Computational Studies of Stability, Reactivity and Degradation Properties of Ephedrine; a Stimulant and Precursor of Illicit Drugs


1. Drew C, Knight G, Hughes D, Bush M. Comparison of the effects of D‐(‐)‐ephedrine and L‐(+)‐pseudoephedrine on the cardiovascular and respiratory systems in man. British J. Clin. Pharmacol. 1978;6:221-225.
2. Mendoza A, Rodríguez-Gil J, González-Alonso S, Mastroianni N, De Alda ML, Barceló D, et al. Drugs of abuse and benzodiazepines in the Madrid Region (Central Spain): seasonal variation in river waters, occurrence in tap water and potential environmental and human risk. Env. Int. 2014;70:76-87.
3. Young R, Gabryszuk M, Glennon RA. (–) Ephedrine and caffeine mutually potentiate one another’s amphetamine-like stimulus effects. Pharmacol. Biochem. Behav. 1998;61:169-173.
4. Docherty J. Pharmacology of stimulants prohibited by the World Anti‐Doping Agency (WADA). British J. Pharmacol. 2008;154:606-622.
5. Cookings JG, Brown M. Ephedrine abuse causing acute myocardial infarction. Med. J. Australia. 1997;167:199-200.
6. Kurashima N, Makino Y, Sekita S, Urano Y, Nagano T. Determination of origin of ephedrine used as precursor for illicit methamphetamine by carbon and nitrogen stable isotope ratio analysis. Anal. Chem. 2004;76:4233-4236.
7. Bruno A, Nolte KB, Chapin J. Stroke associated with ephedrine use. Neurol. 1993;43:1313-1316.
8. Ternes TA, Meisenheimer M, McDowell D, Sacher F, Brauch H-J, Haist-Gulde B, et al. Removal of pharmaceuticals during drinking water treatment. Env. Sci. Technol. 2002;36:3855-3863.
9. Zhang Y, Zhang T, Guo C, Hou S, Hua Z, Lv J, et al. Development and application of the diffusive gradients in thin films technique for simultaneous measurement of methcathinone and ephedrine in surface river water. Sci. Total Env. 2018;618:284-290.
10. Baker DR, Kasprzyk-Hordern B. Spatial and temporal occurrence of pharmaceuticals and illicit drugs in the aqueous environment and during wastewater treatment: new developments. Sci. Total Env. 2013;454:442-456.
11. Rosi-Marshall E, Snow D, Bartelt-Hunt SL, Paspalof A, Tank J. A review of ecological effects and environmental fate of illicit drugs in aquatic ecosystems. J. Hazard. Mater. 2015;282:18-25.
12. Chen M, Guo C, Hou S, Lv J, Zhang Y, Zhang H, et al. A novel Z-scheme AgBr/Pg-C3N4 heterojunction photocatalyst: Excellent photocatalytic performance and photocatalytic mechanism for ephedrine degradation. Appl. Catal. B 2020:118614.
13. Tomić BT, Abraham CS, Pelemiš S, Armaković SJ, Armaković S. Fullerene C 24 as a potential carrier of ephedrine drug–a computational study of interactions and influence of temperature. Phys. Chem. Chem. Phys. 2019;21:23329-23337.
14. Armaković S, Armaković SJ, Tomić BT, Pillai RR, Panicker CY. Adsorption properties of graphene towards the ephedrine–A frequently used molecule in sport. Comput. Theor. Chem. 2018;1124:39-50.
15.  Shanati T, Ansorge-Schumacher MB. Biodegradation of ephedrine isomers by Arthrobacter sp. TS-15: The discovery of novel ephedrine and pseudoephedrine dehydrogenases. Appl. Env. Microbiol. 2020; in press.
16. Stohs SJ, Shara M, Ray SD. p‐Synephrine, ephedrine, p‐octopamine and m‐synephrine: Comparative mechanistic, physiological and pharmacological properties. Phytotherap. Res. 2020; in press.
17. Vukotić AD, Green D, Jevđić JD, Vukotić MR, Petrović N, Stevanović PD. Comparison of efficacy and safety of preemptive infusion protocols of ephedrine and phenylephrine-prevention of hypotension and effects on hemodynamic parameters during spinal anesthesia for caesarean section. Srpski. Arhiv. Celok. Lekar. 2020; in press.
18. Armaković SJ, Armaković S, Šibul F, Četojević-Simin DD, Tubić A, Abramović BF. Kinetics, mechanism and toxicity of intermediates of solar light induced photocatalytic degradation of pindolol: Experimental and computational modeling approach. J. Hazard. Mater. 2020:122490.
19. Harismah K, Mirzaei M. Favipiravir: structural analysis and activity against COVID-19. Adv. J. Chem. B 2020;2:55-60.
20. Ozkendir OM, Mirzaei M. Alkali metal chelation by 3-hydroxyl-4-pyridinone. Adv. J. Chem. B 2019;1:10-16.
21. Mirzaei M, Yousefi M. Computational studies of the purine-functionalized graphene sheets. Superlat. Microstruct. 2012;52:612-617.
22. Schrödinger Release 2020-1: MacroModel, Schrödinger, LLC, New York, NY, 2020.
23. Bochevarov AD, Harder E, Hughes TF, Greenwood JR, Braden DA, Philipp DM, et al. Jaguar: A high‐performance quantum chemistry software program with strengths in life and materials sciences. Int. J. Quant. Chem. 2013;113:2110-2142.
24. Schrödinger Release 2020-1: Jaguar, Schrödinger, LLC, New York, NY, 2020.
25. Jacobson LD, Bochevarov AD, Watson MA, Hughes TF, Rinaldo D, Ehrlich S, et al. Automated Transition State Search and Its Application to Diverse Types of Organic Reactions. J. Chem. Theor. Comput. 2017;13:5780-5797.
26. Shivakumar D, Williams J, Wu Y, Damm W, Shelley J, Sherman W. Prediction of absolute solvation free energies using molecular dynamics free energy perturbation and the OPLS force field. J. Chem. Theor. Comput. 2010;6:1509-1519.
27. Guo Z, Mohanty U, Noehre J, Sawyer TK, Sherman W, Krilov G. Probing the α‐Helical Structural Stability of Stapled p53 Peptides: Molecular Dynamics Simulations and Analysis. Chem. Biol. Drug Design 2010;75:348-359.
28. Bowers KJ, Chow E, Xu H, Dror RO, Eastwood MP, Gregersen BA, et al., editors. Scalable algorithms for molecular dynamics simulations on commodity clusters. Proceedings of the 2006 ACM/IEEE conference on Supercomputing; 2006: ACM.
29. Becke AD. Density‐functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993;98:5648-5652.
30. Harder E, Damm W, Maple J, Wu C, Reboul M, Xiang JY, et al. OPLS3: a force field providing broad coverage of drug-like small molecules and proteins. J. Chem. Theor. Comput. 2015;12:281-296.
31. Jorgensen WL, Maxwell DS, Tirado-Rives J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 1996;118:11225-11236.
32. Jorgensen WL, Tirado-Rives J. The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. J. Am. Chem. Soc. 1988;110:1657-1666.
33. Berendsen HJ, Postma JP, van Gunsteren WF, Hermans J. Interaction models for water in relation to protein hydration.  Intermolecular forces: Springer; 1981. p. 331-342.
34. Otero-de-la-Roza A, Johnson ER, Contreras-García J. Revealing non-covalent interactions in solids: NCI plots revisited. Phys. Chem. Chem. Phys. 2012;14:12165-12172.
35. Johnson ER, Keinan S, Mori-Sánchez P, Contreras-García J, Cohen AJ, Yang W. Revealing noncovalent interactions. J. Am. Chem. Soc. 2010;132:6498-6506.
36. Sjoberg P, Murray JS, Brinck T, Politzer P. Average local ionization energies on the molecular surfaces of aromatic systems as guides to chemical reactivity. Canadian J. Chem. 1990;68:1440-1443.
37. Murray JS, Seminario JM, Politzer P, Sjoberg P. Average local ionization energies computed on the surfaces of some strained molecules. Int. J. Quant. Chem. 1990;38:645-653.
38. Armaković S, Armaković SJ, Šetrajčić JP, Šetrajčić IJ. Active components of frequently used β-blockers from the aspect of computational study. J. Mol. Model. 2012;18:4491-4501.
39. Okulik N, Jubert AH. Theoretical study on the structure and reactive sites of three non-steroidal anti-inflammatory drugs: Ibuprofen, Naproxen and Tolmetin acids. J. Mol. Struct. THEOCHEM. 2006;769:135-141.
40. Politzer P, Abu‐Awwad F, Murray JS. Comparison of density functional and Hartree–Fock average local ionization energies on molecular surfaces. Int. J. Quant. Chem. 1998;69:607-613.
41. Bulat FA, Toro-Labbé A, Brinck T, Murray JS, Politzer P. Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies. J. Mol. Model. 2010;16:1679-1691.
42. War JA, Jalaja K, Mary YS, Panicker CY, Armaković S, Armaković SJ, et al. Spectroscopic characterization of 1-[3-(1H-imidazol-1-yl) propyl]-3-phenylthiourea and employing DFT calculations and molecular dynamics simulations. J. Mol. Struct. 2017;1129:72-85.
43. Zainuri DA, Arshad S, Khalib NC, Razak IA, Pillai RR, Sulaiman SF, et al. Synthesis, XRD crystal structure, spectroscopic characterization (FT-IR, 1H and 13C NMR), DFT studies, using molecular dynamics simulations and evaluation of antimicrobial and antioxidant activities of a novel chalcone derivative, (E)-1-(4-bromophenyl)-3-(4-iodophenyl)prop-2-en-1-one. J. Mol. Struct. 2017;1128:520-533.
44. Menon VV, Fazal E, Mary YS, Panicker CY, Armaković S, Armaković SJ, et al. FT-IR, FT-Raman and NMR characterization of 2-isopropyl-5-methylcyclohexyl quinoline-2-carboxylate and investigation of its reactive and optoelectronic properties by molecular dynamics simulations and DFT calculations. J. Mol. Struct. 2017;1127:124-137.
45. Armaković S, Armaković SJ, Pelemiš S, Mirjanić D. Influence of sumanene modifications with boron and nitrogen atoms to its hydrogen adsorption properties. Phys. Chem. Chem. Phys. 2016;18:2859-2870.
46. Drew KL, Reynisson J. The impact of carbon–hydrogen bond dissociation energies on the prediction of the cytochrome P450 mediated major metabolic site of drug-like compounds. Eur. J. Med. Chem. 2012;56:48-55.
47. Lienard P, Gavartin J, Boccardi G, Meunier M. Predicting drug substances autoxidation. Pharm. Res. 2015;32:300-310.
48. Andersson T, Broo A, Evertsson E. Prediction of drug candidates' sensitivity toward autoxidation: computational estimation of C-H dissociation energies of carbon‐centered radicals. J. Pham. Sci. 2014;103:1949-1955.
49. Gryn'ova G, Hodgson JL, Coote ML. Revising the mechanism of polymer autooxidation. Org. Biomol. Chem. 2011;9:480-490.