Jeffrey R.
Johns
*a and
James A.
Platts
b
aMelatonin Research Group, Faculty of Pharmaceutical Sciences, Khon Kaen University, Khon Kaen, 40002, Thailand. E-mail: jjeff@kku.ac.th; Tel: +66-43-202378
bSchool of Chemistry, Cardiff University, Park Place, Cardiff, CF10 3AT, UK. E-mail: platts@cardiff.ac.uk; Fax: +44 (0)2920-874030; Tel: +44 (0)2920-874950
First published on 19th August 2014
Density functional theory calculations on melatonin, metabolites and synthetic derivatives thereof, and a range of other biological antioxidant molecules are presented, with a view to understanding the antioxidant ability of these molecules. After testing of the necessary calculations, we show that melatonin lies close to vitamin E on a donor–acceptor map, indicating that it should be an excellent electron donor but a poor acceptor. The neutral radical metabolite of melatonin is predicted to be an even better donor, whereas other metabolites and synthetic derivatives should retain antioxidant ability but are less powerful than the parent. QSAR models of antioxidant activity, measured in two different assays, are presented. We show that octanol–water partition coefficient is an excellent predictor of activity in lipophilic media, while properties related to electron donor/acceptor power give good fits against activity in aqueous media.
Melatonin is finding great utility in preventing diseases related to oxidative damage including cancer12 and neurodegenerative diseases13,14 as well as its well known role in treatment for reducing insomnia, jet lag, migraine, headache, etc.13,15,16 It is being widely investigated for a large number of other diseases in a large number of clinical trials.17 In addition, consumption of tropical fruits containing melatonin has been shown to reduce antioxidant levels in humans.5,18
The antioxidant radical scavenging properties of melatonin and its metabolites cyclic-3-hydroxymelatonin (cyclic-3OHM),10N(1)-acetyl-N(2)-formyl-5-methoxykynuramine) (AFMK), N(1)-acetyl-5-methoxykynuramine (AMK)7 and 6-hydroxymelatonin (6-OHmel)19 occurs mainly via the one electron transfer process.8,20,21 The first ionization potential (IP) and the electron affinity (EA) are properties of a system that allow measurement of its propensity to donate or accept one electron. The best antioxidants present low IP values, because the lower the IP, the easier the electron abstraction, and vice versa for EA and electron acceptance (antireductant).
Gazquez et al.21 have presented an elegant model to explain relative scavenging activity and antioxidant power of compounds using these two properties. Quantum chemical density functional theory (DFT) calculations can be used to obtain accurate ionization potentials, electron affinities, electrodonating, and electroaccepting power indexes (with respect to internal standards, such as fluorine and sodium atoms). These values can then be used to construct a donator acceptor map (DAM), indicating whether molecules are good electron donors or acceptors. The DAM is a powerful representation of these key properties, helping to reveal the antiradical capacity of any substance and allowing qualitative comparison between substances, alongside quantitative measures obtained from experiment or theory. Previous DAM studies have included linear polyene-conjugated molecules,23 carotenoids,24 a large series of carotenoids,25 carotenoids, melatonin and vitamins,26 and psittacofulvins and anthocyanins.27
Vitamin E or α-tocopherol has been described the “last line of defense” in a multicomponent endogenous antioxidant system.28 It appears that under conditions of stress, depletion of cellular ascorbic acid occurs first, followed by glutathione, then α-tocopherol, resulting in initiation of lipid peroxidation. When glutathione is depleted, ascorbic acid plays a vital role in maintaining cellular α-tocopherol levels and survival of the cell.29 One might expect that melatonin should be depleted after α-tocopherol, particularly in membranes,30 as it is higher in the electrochemical series at 700 mV,8 compared to 500 mV for α-tocopherol.11 Melatonin may therefore truly be the last line of defense against oxidative damage.31 The role of melatonin in this multicomponent antioxidant is still unclear, although there is evidence that melatonin cause upregulation of superoxide dismutase, glutathione reductase and catalase.32,33
Several melatonin derivatives that were substituted on the indole nitrogen (Fig. 2) have been previously reported for in vitro antioxidant effects and anti-inflammatory activities.34 Their synthesis and characterization is described in this ref. 34.
The aim of this study was to investigate the antioxidant radical scavenging properties of melatonin and its metabolites cyclic-3OHM, AFMK, AMK and 6-hydroxymelatonin, and several N-indole substituted derivatives via the one electron transfer process, using a donor–acceptor map. Other classical antioxidants and vitamins are modeled for comparison. QSAR relationships of some N-indole substituted derivatives between in vitro antioxidant properties experimentally measured by lipid peroxidation of rat brain homogenate using thiobarbituric acid reacting substances (TBARS IC50)35,36 and Oxygen Radical Absorbance Capacity Assay (ORAC) data, and a number of derived electronic properties e.g. HOMO/LUMO energies, donor and acceptor power (Rd, Ra), hardness, electronegativity and logP were investigated. The ORAC assay is based on the scavenging of peroxyl radicals generated by 2,2′-azobis (2-methylpropionamidine) dihydrochloride (AAPH) in aqueous media, which prevent the degradation of the fluorescein probe and, consequently, prevent the loss of fluorescence of the probe.37,38 The antioxidant activity was calculated from the integrated area under the fluorescence curve (AUC) for each antioxidant.
To determine the accuracy of DFT for predicting IP/EA of indoleamines, calculations in gas phase using different basis sets were compared to previously reported photoelectron spectroscopy measurements. IP and bond-dissociation energies for many antioxidant systems do not follow the same trends in gas and solution phases, such that major differences with respect to vacuum are found as when water computations are performed.42 On the basis of the computed BDE and IP values, to more realistically model antioxidant activity in vivo, calculations were therefore performed using the polarizable continuum model (PCM water) i.e. placing the solute in a cavity within the solvent reaction field.
The validity of using B3LYP for calculating EA has been raised, due to most DFT functionals (including B3LYP) being incapable of binding the whole excess electron.43 This may not be revealed when using standard basis sets, even with multiple diffuse functions, since they artificially constrain the electron density to remain near the nuclei. The error due to this constraint depends on the magnitude of the EA, which could render trends in EA unreliable, especially for low EAs. Thus a range-separated DFT method, CAM-B3LYP,44 was used to investigate the problem of fractional EA. A positive energy for the HOMO of an anion species is an easy diagnostic for the fractional EA problem.43
The melatonin molecule contains three freely rotatable bonds in its imidazole side chain. To investigate the effect of conformation, MOE40 was used to select a range of typical conformers via stochastic search of rotatable bonds, using the MMFF94 forcefield. 10000 conformers were generated then sorted into clusters based on dihedral angles of the freely rotatable bonds. Vertical IP and EA of the lowest energy conformers from each cluster were then calculated, following geometry optimization of the neutral molecule. No significant changes in geometry on DFT optimization were noted, indicating that conformers remained in their local energy minima.
Donor–acceptor maps were calculated following the method of Martinez et al.,26 using the same experimental values of IE and EA for sodium of 5.140 and 0.540 eV and for fluorine of 17.540 and 3.400 eV respectively, taken from the literature. This set reference points on the map of Ra = 1 for the sodium atom, and Rd = 1 for the fluorine atom. Values calculated using Gaussian 09/RB3LYP/6-31+G* gas phase were 5.406 and 0.596 eV for sodium and 21.405 and 3.513 eV for fluorine respectively.
As formulated by Gazquez et al.22 and applied in the study of Martínez et al.,27 the propensity to donate charge, or electrodonating power, may be defined as:
For electroaccepting power, higher values imply a greater capacity for accepting charge, whereas for electrodonating power, lower values imply a greater capacity for donating charge. It should be noted that ω− and ω+ refer to fractional charges, however I and A refer to donating or accepting a single, whole electron. Thus, a simple charge transfer model, framed in terms of chemical potential and hardness is used to describe electrodonating and electroaccepting powers. The charge flow direction is measured by chemical potential, along with the capacity to donate or accept charge. More emphasis is assigned to the ionization potential than to the electron affinity in the context of the charge donation process. Likewise, more significance is assigned to electron affinity than to ionization potential for electroaccepting power. Hardness provides a measure of the resistance to the electron transfer. So that a range of substances can be compared for electrodonating and electroaccepting power, experimental values of I and A for sodium and fluorine are used a reference points to provide corresponding ω+ and ω− values. Sodium represents a good electron donor and fluorine represents a good electron acceptor. For some substance L, the electron acceptance index can be defined as:
When Ra = 1, then ωL+ ≈ ωF+ and L is as effective an electron acceptor as fluorine. When Ra > 1, then ωL+ > ωF+ and L is a more effective electron acceptor than fluorine. When Ra < 1, then ωL+ < ωF+ and L is a less effective electron acceptor than fluorine. Similarly, the electron donation index can be defined as:
When Rd = 1, then L is as effective an electron donor as sodium, and when Rd > 1, L is a less effective electron donor than sodium, whereas when Rd < 1, L is a more effective electron donor. If Ra and Rd are determined, then any substance L can be characterized in terms of its electron donor–acceptor capacity, and mapped on a donor acceptor map (DAM).
Octanol–water partition coefficients (logP) were estimated using ACD/Chemsketch logP plugin (Advanced Chemistry Development, Inc., Toronto, Canada, 2012), MarvinView 5.11.3 (Chemaxon, Budapest, Hungary, 2012) and MolKa (Molecular Discovery Ltd, Perugia, Italy, 2012) and compared to literature values where available.
Data on the antioxidant capacity of melatonin derivatives in vitro, using the widely adopted method of Callaway et al.35 for measurement of lipid peroxidation, the thiobarbituric acid reacting substances (TBARS) and brain homogenate, was used in this study. Inhibitory effect on nitric oxide (NO) of melatonin and these derivatives has been previously reported by our group.34
7.76 | IP | PES45 | |
7.03 | IP | PES46 | |
7.7 | IP | Vertical26 | Estimated from PES |
6.83 | IP | Vertical26 | Gaussian03 B3LYP/D95 V gas phase |
−1.00 | EA | Vertical26 | Gaussian03 B3LYP/D95 V gas phase |
7.07 | IP | Vertical | This study. Gaussian09 B3LYP 6-31+G* gas phase |
−0.70 | EA | Vertical | This study. Gaussian09 B3LYP 6-31+G* gas phase |
6.85 | IP | Relaxed | This study. Gaussian09 B3LYP 6-31+G* gas phase |
−0.65 | EA | Relaxed | This study. Gaussian09 B3LYP 6-31+G* gas phase |
Using the CAM-B3LYP functional did not give significantly different EA values for melatonin using any of the basis sets in this study. Furthermore, no positive energies for HOMO's of any anion species were observed. The reactivity indices calculated for compounds in this study depend mainly on the IP values, rather than EA, as they are mostly electron donors, thus small EA errors will have negligible impact on the overall findings.
Comparison of reported values for IP and EA reported in the literature is made with those calculated using B3LYP 6-31+G* basis set and gas phase in Table 1.
MOE conformer energy (kcal mol−1) | Vertical IP/eV | Vertical EA/eV |
---|---|---|
20.573446 | 5.544 | −0.886 |
20.573452 | 5.603 | −0.739 |
21.250622 | 5.538 | −0.879 |
21.250629 | 5.538 | −0.879 |
21.694929 | 5.429 | −0.947 |
21.694931 | 5.429 | −0.947 |
21.724112 | 5.443 | −0.952 |
Mean | 5.5036 | −0.8898 |
S.D. | 0.0694 | 0.0745 |
The donor/acceptor maps shows melatonin to be a very good electron donor, along with its metabolites 6-hydroxymelatonin and cyclic-3-hydroxymelatonin and the melatoninyl neutral radical, and several classical antioxidants such as vitamin E (α-tocopherol), epigallocatechin gallate (ECGC), resveratrol, xanthurenic acid and quercetin (a typical flavenoid). Other melatonin metabolites (AFMK, AMK) and melatonin derivatives showed weaker electron donor strength, similar to other classical antioxidants such as vitamin A, vitamin C, beta-carotene and α-lipoic acid. The 4-nitro derivative of melatonin is not shown on the DAM and appears off the top right quadrant at Ra = 4.35 and Rd = 5.62, being a much poorer electron donor and better electron acceptor than melatonin, due to its strongly withdrawing nitro group. Interestingly, the melatoninyl neutral radical that results from a 1-electron 1-proton donation from melatonin is an even more powerful electron donor than melatonin itself.
These results are supported by a large number of experimental observations where melatonin acts as a direct scavenger of free radicals with the ability to detoxify both reactive oxygen and reactive nitrogen species, and indirectly by increasing the activity of the antioxidative defense systems.9,47 Researchers have reported that the peroxyl radical scavenger ability of melatonin is better than of α-tocopherol, vitamin C and reduced glutathione (GSH),48 and more potent than xanthurenic acid, resveratrol, EGCG, vitamin C and α-lipoic acid in inhibiting ˙OH-induced oxidative DNA damage generated by oxygen-derived free radicals from Fenton reaction.49 Melatonin has been demonstrated to reduce the formation of 8-hydroxy-2′-deoxyguanosine, a product of damaged DNA repair, 60 to 70 times more effectively than ascorbate or α-tocopherol.50 Melatonin also plays an important role in protecting cellular membranes against lipid peroxidation.51
It should be noted that most other dietary antioxidants lie outside the scale of this donor–acceptor map, towards the top right quadrant, including carotenoids, psittacofulvins and anthocyanins, flavenoids and polyphenols. These tend to be electron acceptors rather than electron donors i.e. antireductants.
LogP | ||||
---|---|---|---|---|
Antioxidant | ACDa | Marvin weightedb | MolKac | Literature values |
a ADC/Chemsketch logP plugin (Advanced Chemistry Development, Inc., Toronto, Canada) 2012 http://www.acdlabs.com. b MarvinView 5.11.3 (Chemaxon, Budapest, Hungary) 2012 http://www.chemaxon.com. c MoKa (Molecular Discovery Ltd, Perugia, Italy) 2012 http://www.moldiscovery.com. d http://www.chemicalize.org/structure/#!mol=Vita+E. e Human Metabolome Database Version 3.6. http://www.hmdb.ca/metabolites/HMDB00305. | ||||
Beta-carotene | 15.51 ± 0.43 | 11.12 | 9.0 | 14.76 |
α-Tocopherol (vitamin E) | 10.66 ± 0.28 | 8.94 | 9.0 | 10.51d |
Vitamin A (retinol) | 6.84 ± 0.33 | 6.07 | 6.1 | 4.69–6.38e |
Alpha-lipoic acid | 2.16 ± 0.29 | 2.11 | 2.4 | — |
Melatonin | 0.96 ± 0.44 | 1.41 | 1.4 | 1.254,55 |
6-OH-melatonin | 0.02 ± 0.83 | 0.84 | 1.0 | — |
Melatoninyl neutral radical | 0.02 ± 0.83 | 0.88 | — | — |
AFMK | 0.82 ± 0.52 | 0.34 | 0.0 | 0.4856 |
AMK | 0.65 ± 0.49 | 0.33 | 0.2 | 0.7456 |
−0.7758 | ||||
Tryptophan | 0.87 ± 0.31 | 1.51 | 0.6 | 1.0857 |
−1.0659 | ||||
2-Naphthoyl-melatonin | 2.73 ± 0.46 | 3.31 | 3.7 | — |
2-Benzoyl-melatonin | 1.5 ± 0.46 | 2.32 | 2.4 | — |
Cyclic-3OHmel | −0.64 ± 0.89 | −0.21 | 0.8 | — |
Acetyl-melatonin | 1.00 ± 0.87 | 0.47 | 1.3 | — |
Vitamin C | −3.26 ± 0.56 | −1.98 | −3.4 | −1.8558 |
TBARS pIC50 results were correlated with a number of derived electronic properties e.g. HOMO/LUMO energies, donor and acceptor power (Rd, Ra), hardness, electronegativity and logP, and results are shown in Table 4. Only lipophilicity (logP) correlated well, correlation with other molecular properties (donor/acceptor power, electronegativity and LUMO was poorer, with low or no correlation with hardness, HOMO and HOMO/LUMO energy difference. This may be because of the nature of the brain homogenate lipid peroxidation assay, where solubility of the antioxidant in the lipid domain is the dominant factor contributing to radical scavenging. The hydrophobicity of the antioxidant may also be an important criterion for passive transport into cells across the hydrophobic phospholipid bilayer of the cellular membranes. Furthermore, the single electron transfer mechanism for direct radical scavenging of melatonin, although the most favourable mechanism in aqueous solution, is not favourable in aprotic solvents e.g. benzene, where hydrogen atom transfer/proton coupled electron transfer or radical adduct formation are favoured.29
Molecular parameters | Correlation (r2) with TBARS pIC50 | Correlation (r2) with ORAC logAUC |
---|---|---|
Donor power (Rd) | 0.636 | 0.951 |
Acceptor power (Ra) | 0.647 | 0.935 |
Hardness (η) | 0.458 | 0.982 |
Electronegativity (X) | 0.589 | 0.987 |
Energy HOMO (eV) | 0.331 | 0.928 |
Energy LUMO (eV) | 0.560 | 0.995 |
Energy HOMO–LUMO (eV) | 0.237 | 0.941 |
LogP | 0.978 | 0.554 |
By contrast, for the ORAC assay, which was performed in aqueous medium, all molecular parameters correlated highly with the ORAC AUC except for logP (Table 4), as has been observed previously with indoleamines where antioxidant potency was measured for lipid peroxidation using a conjugated dienes assay60 and in phenolic compounds.61
Some QSAR correlation plots are shown in Fig. 5, and all correlation plots are shown in the ESI.†
Fig. 5 QSAR correlation plots of some calculated molecular parameters and TBARS or ORAC antioxidant activity. |
This method has therefore been used to map out the donor–acceptor power of melatonin, its metabolites, some synthetic derivatives and a range of classical antioxidants. This approach clearly shows that melatonin lies in the range of good electron donors and bad electron acceptors, with similar power to vitamin E. Interestingly, the first neutral, radical metabolite of melatonin is an even better donor than the parent molecule, which will have important implications for the overall biology of the cascade process by which melatonin mops up ROS. Other metabolites, as well as most synthetic derivatives, remain in the range where substantial antioxidant ability should be expected, but a 4-nitro derivative lies well outside this region.
QSAR investigation indicates the ability of melatonin derivatives to protect against lipid peroxidation of brain homogenate strongly correlated with their lipophilicity (logP) but only weakly to other molecular properties related to donor/acceptor ability (donor/acceptor power, electronegativity, hardness, HOMO, LUMO, HOMO-LUMO energies). By contrast, these molecular parameters correlate strongly with ORAC antioxidant power measured in aqueous phase.
The range of lipophilicity of melatonin and its metabolites (logP between −1 and 2) may explain the large number of antioxidant arenas where melatonin seems to play a role in protecting against ROS damage; they lie between the traditional membrane protectors (α-tocopherol, vitamin A and carotenoids) and hydrophilic compounds (vitamin C, lipoic-acid) and aqueous antioxidant enzymes.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4ob01396d |
This journal is © The Royal Society of Chemistry 2014 |