Cancer is a major public health problem with limited success of available treatments, pointing to the need for new strategies to be developed. Phosphoethanolamine exhibits broad antitumor activity in a variety of tumor cells and potent inhibitor effects on tumor progress in vivo. Once-used organophosphates inhibit acetylcholinesterase (AChE), resulting in toxic effects to the user. As this group is present in phosphoethanolamine, we perform prediction of the in silico metabolism of phosphoethanolamine and submit this series to a docking study on AChE. A total of 10 metabolites were indicated by the prediction, including ammonia and hydroxylamine, which were not included in the study. Using a group of 8 organophosphorus whose pIC50 values ranged from 5.92 to 9.47 as template, we observed that no compound present in the phosphoethanolamine series had a binding energy lower than that of organophosphorus, suggesting that the series has low inhibitory power on AChE. In light of this, we conclude that phosphoethanolamine and its predicted metabolites do not significantly inhibit AChE to cause a cholinergic crisis. This finding highlights the importance of investigating this compound as lead for potential anticancer agents.
1. Introduction
Cancer is a major public health problem in many parts of the world, representing the second leading cause of death in the Western world. The number of individuals diagnosed with cancer increases each year, with millions of deaths. For example, in 2012, 8.2 million people died of cancer—approximately 13% of all deaths worldwide. Furthermore, it is predicted that cancer will increase from 10 million new cases in 2000 to 15 million in 2020, an increase of 50% [1–3]. The limited successes of available treatments for cancer mean that new strategies need to be developed [4, 5].
Phospholipids are present in tissues and there is a special scientific interest in the elucidation of their role in the treatment of cancer. Phosphoethanolamine is a precursor of phosphatidylcholine and phosphatidylethanolamine, both of which are involved in the turnover of cell membrane phospholipids. Both of these phospholipids take part in lipid signaling pathways, either acting as ligands or generating intermediate substrates [6–8].
A class of potential synthetic antineoplastic phospholipids such as alkylphospholipids represents a promising treatment for cancer. Ether lipids, as this class is collectively termed, have been shown to reduce the synthesis of phosphatidylcholine in cancer cells, interfering with phospholipid-induced apoptosis. Miltefosine (hexadecylphosphocholine) is an analogous lipid compound that is used for control of cutaneous metastases of human malignant melanoma [8–10].
Studies have shown that synthetic phosphoethanolamine exhibits broad antitumor activity in a variety of tumor cells and potent inhibitor effects on tumor progress in vivo. Other actions of phosphoethanolamine include antiproliferative and proapoptotic effects on melanoma cells [8], activity in the renal carcinoma model [11], apoptosis in human breast cancer MCF-7 cells via the mitochondrial pathway [12], inhibition of tumor growth, and increase of animals’ lifespan without causing liver or hematological toxicity [13, 14].
Organophosphates are organic compounds commonly used as pesticides. Their ability to inhibit acetylcholinesterase (AChE) is directly related to the toxicity of this group of molecules. Human AChEs are a class of enzymes whose biological function is to catalyze the hydrolysis of the neurotransmitter acetylcholine in the terminal nerve junction. The inhibition of AChE by organophosphorus involves phosphorylation of active binding site blocking acetylcholine hydrolysis, through attack of hydroxyl group from serine to nucleophilic phosphorus [15–18].
The presence of a phosphate group in phosphoethanolamine could increase the toxicity of this potential anticancer agent or its metabolites, generated after biotransformation. In light of this, the present work predicts the in silico metabolism of phosphoethanolamine and evaluates the affinity of the compound, three analogues, and its metabolites to AChE through molecular docking.
2. Materials and Methods2.1. In Silico Metabolism
For prediction of metabolism, Metabolizer version 15.6.8, part of ChemAxon’s JChem package (https://www.chemaxon.com/), was used. The software transforms the substrate according to generic transformations from a current biotransformation library. The structure of phosphoethanolamine was drawn in the program and then calculated according to default configuration.
Metabolizer calculates global accumulation by considering all alternative metabolic routes to the metabolite. The ones with highest global accumulations are considered major metabolites.
2.2. Physicochemical Properties
Determination of logP, logD, and molar refractivity were obtained using Marvin Sketch 14.12.8.0 while electrostatic surface was calculated using Marvin Space 14.12.8.0; both are part of ChemAxon’s JChem package (https://www.chemaxon.com/). All four calculations used default configuration.
The logP plugin calculates the octanol/water partition coefficient, which is used in rational drug design as a measure of molecular hydrophobicity. logD is important in compounds with ionizable groups, once ionization of those groups; thus the ratio of the ionic forms depends on the pH. logD represents the compounds at any pH value. Calculations are based on a pool of fragments predefined in the calculator. The partition value of a molecule is the sum of the fragment values present in the molecule. Weighted method was used, which is a combination of three partition calculations VG, KLOP, and PHYPROSP [19–21].
Molar refractivity is strongly related to molar volume and London dispersive forces that has important effect on drug-receptor interaction. Calculation was based on the atomic method proposed by Viswanadhan et al. [19]. For electrostatic parameters, van der Walls surface was determined for each compound. Then, electrostatic potential was selected as color type and exported.
2.3. VolSurf Descriptors
Three-dimensional structures (3D) were used as input data in the VolSurf+ program v. 1.0.7 [20] and were subjected to molecular interaction fields (MIF) [20] to generate descriptors using the following probes: N1 (amide nitrogen-hydrogen bond donor probe), O (carbonyl oxygen-hydrogen bond acceptor probe), OH2 (water probe), and DRY (hydrophobic probe). Additional non-MIF-derived descriptors were generated to create a total of 128 descriptors [20]. VolSurf descriptors have been previously used to predict the inhibitory actions of alkaloids against monoamine oxidase B (MAO-B) [21].
CACO2 and SKIN are computed via PLS model derived by fitting VolSurf+ descriptors to experimental data. The values are only qualitative for CACO2 and quantitative (cm/h) for skin.
Log of the Blood-Brain Barrier (LgBB) distribution indicates permeability level. Values lower than −0.5 indicate very poor brain permeation while values greater than 0.5 indicate high brain permeation. Volume of Distribution (VD) is computed via a PLS model derived by fitting VolSurf+ descriptors to −log of experimental data on volume of distribution (Litre/Kg).
Metabolic Stability (MetStab) is percentage of the remaining compound after incubation with human CYP3A4 enzyme. Values greater than 50 indicate stable behavior (100 is the maximum stability); values lower than 50 are less precise and indicate metabolic instability.
2.4. Docking
The structure of acetylcholinesterase (PDB ID 2ACE) [22] was downloaded from the Protein Data Bank (http://www.rcsb.org/pdb/home/home.do). Molecular structures were submitted to molecular docking using the Molegro Virtual Docker, v. 6.0.1 (MVD) [23]. All of the water molecules were deleted from the enzyme structure and the enzyme and compound structures were prepared using the same default parameter settings in the software package (score function: MolDock score; ligand evaluation: Internal ES, Internal HBond, Sp2-Sp2 Torsions, all checked; number of runs: 10; search algorithm: MolDock SE; maximum interactions: 1500; max. population size: 50; max. steps: 300; neighbor distance factor: 1.00; and max. number of poses returned: 5). The docking procedure was performed using a GRID of 15 Å in radius and 0.30 in resolution to cover the ligand-binding site of the acetylcholinesterase structure [24].
A group of 8 organophosphorus [25], whose IC50 values ranged from 0.34 nM to 1.2 μM (pIC50 values, −log IC50 (mol/L), from 5.92 to 9.47), were used to evaluate the binding energy of phosphoethanolamine and its metabolites to AChE (Table 1).
Calculated logP and molar refractivity of phosphoethanolamine and analogues, phosphoethanolamine metabolites, and group of 8 organophosphorus compounds. Pho: phosphoethanolamine; Pho_01 to Pho_03: phosphoethanolamine analogues; M1 to M8: phosphoethanolamine metabolites; O1: Chlorpyrifos-Oxon; O2: O,O-dibutyl O-(2,2-dichlorovinyl) Phosphate; O3: Paraoxon; O4: Cyclic Tolyl Saligenin Phosphate; O5: Diisopropyl Phosphorofluoridate; O6: Cyclic Phenyl Saligenin Phosphate; O7: Dichlorvos; O8: Trichlorfon.
Compound
clogP
Molar refractivity
Pho
−2.05
25.25
Pho_01
−0.26
43.71
Pho_02
1.24
63.28
Pho_03
2.35
73.77
M1
−2.02
37.51
M2
−1.60
23.63
M3
−2.40
25.46
M4
−1.87
22.71
M5
−1.95
23.59
M6
−2.41
38.54
M7
−2.45
26.28
M8
−2.52
27.22
O1
3.78
75.46
O2
3.18
66.21
O3
3.69
71.25
O4
1.37
52.19
O5
1.76
40.89
O6
4.02
79.94
O7
2.43
63.69
O8
1.14
47.59
3. Results and Discussion
The Metabolizer software is useful for enumeration of the metabolites of drugs and the prediction of their likelihood. A total of 10 metabolites were generated, including ammonia and hydroxylamine, which were not included in the study due to their simple structures. The structure of the selected metabolites and their respective global accumulation values are presented in Figure 1.
Structure of phosphoethanolamine metabolites and their respective GA (global accumulation). The highest GA for M2 indicates that it is the major metabolite.
Values of calculated logP (clogP) and molar refractivity of phosphoethanolamine and analogues (Figure 2), phosphoethanolamine metabolites (M1 to M8), and a group of 8 organophosphorus compounds, used as parameter to evaluate phosphoethanolamine and derivatives (analogues and metabolites), are displayed in Table 1. logD of phosphoethanolamine and derivatives are displayed in Figure 3.
Structure of phosphoethanolamine (Pho) and three analogues (Pho_01 to Pho_03), substituting hydroxyl group on phosphate and amine group in Pho_02.
logD of phosphoethanolamine (Pho), analogues (Pho_01 to Pho_03), and metabolites (M1 to M8) in different pH values.
Electrostatic properties of phosphoethanolamine and derivatives are represented in Figure 4. A high electronic density is observed in phosphate and amine group in phosphoethanolamine; the hydrogen linked to electronegative atoms presents low electronic density (blue color). An increase of polarity (red color) is well observed in metabolites, increasing molecular polarity.
Electrostatics properties of phosphoethanolamine and derivatives. Red color indicates regions with high electronic density.
The group of 8 organophosphorus compounds [26] was submitted to molecular docking and their energy binding and pIC50 values (Table 2) were used as parameters to evaluate the phosphoethanolamine series. Chlorpyrifos-Oxon showed the lowest binding energy value of −107.9 MolDock energy, corroborated by a smaller pIC50 value. A good correspondence between experimental and docking data was observed in the analysis of the 4 compounds with lower pIC50 values, all of which had energy values under −100 MolDock energy (Table 2). The coefficient of determination (R2) was 0.68, indicating a significant correlation (Figure 5).
Group of organophosphorus compounds used as parameters, with structure and name, pIC50 (−log IC50 (mol/L)), and binding energy.
Compound
pIC50
Binding energy (MolDock score)
9.47
−107.9
9.03
−106.3
8.64
−100.8
7.55
−103.5
6.72
−77.2
6.59
−97.5
6.46
−77.1
5.92
−81.1
Correlation between binding energy value in MolDock score (x-axis) and experimental pIC50 values (y-axis) of organophosphorus group. The coefficient of determination (R2) is also indicated.
High binding energy values (Table 3) for phosphoethanolamine and its metabolites indicated an inferior interaction with AChE compared to the organophosphorus group (Table 3). Phosphoethanolamine had a low energy binding value of −59.2 MolDock score, indicating a weak interaction between the molecule and the receptor, thus suggesting a low toxicity for this enzyme. In silico metabolism indicated M2 as major metabolite (Figure 1), as the binding energy of this compound was one of the lowest in the series. The other metabolite with considerable abundance was M8, with a −72.3 MolDock score energy. This was the highest value among the metabolites but lower than any organophosphorus compound used as a template, confirming the weak interaction of phosphoethanolamine and its metabolites with AChE. The three proposed analogues show energy binding value lower than phosphoethanolamine, increasing interaction with AChE.
Combine binding energy values in MolDock score of phosphoethanolamine (Pho), phosphoethanolamine analogues (Pho_01 to Pho_03), phosphoethanolamine metabolites (M1 to M8), and organophosphorus group (O1 to O8).
Some molecular descriptors are parameters appropriable to evaluate ADME process and drug behavior (Table 4). Five molecular VolSurf descriptors were selected in order to evaluate permeability (CACO2 cells and SKIN), distribution (LgBB and VD), and stability (MetStab). Phosphoethanolamine and metabolites present low permeability through CACO2 cells and SKIN. The analogues showed a better comparative performance, probably related to their lower polarity.
VolSurf molecular descriptors appropriable to evaluate ADME process and drug behavior. CACO2 cells and SKIN permeability, log blood-brain barrier distribution (LgBB), volume of distribution (VD), and metabolic stability (MetStab).
Compound
CACO2
SKIN
LgBB
VD
MetStab
Pho
−1.03805
−6.81425
−2.00629
−0.1297
100
Pho_01
0.352823
−4.14639
0.068157
0.240273
100
Pho_02
1.10666
−2.3428
0.375648
0.27598
73.9335
Pho_03
0.53779
−3.05418
0.281418
0.155534
84.9974
M1
−1.44627
−7.72323
−2.4985
−0.24442
100
M2
−1.38436
−6.70337
−2.17963
−0.32468
100
M3
−1.3003
−6.89743
−2.12005
−0.11547
100
M4
−0.82973
−6.05096
−1.46413
−0.03631
100
M5
−1.09326
−6.60327
−2.09054
−0.3334
100
M6
−1.89772
−10.902
−2.98507
−0.17428
100
M7
−1.43742
−8.2624
−2.53329
0.024709
100
M8
−1.58487
−8.78139
−2.44954
−0.01373
100
O1
1.37183
−1.79192
0.289113
0.127351
39.3059
O2
1.28518
−1.60439
0.442147
−0.00769
49.6523
O3
1.34774
−1.56177
0.431478
0.056206
44.2688
O4
1.2689
−2.09574
0.411952
0.138814
75.1341
O5
1.08652
−2.34342
0.176429
−0.08398
87.1393
O6
1.53633
−1.22958
0.486752
0.033381
28.1952
O7
1.0289
−2.50402
0.109418
−0.04239
58.5821
O8
0.957486
−3.17353
−0.0024
0.211614
86.3035
An unsatisfactory distribution for phosphoethanolamine and derivatives in systemic circulation and central nervous system was also observed, whereas analogues display better results again. About stability, both phosphoethanolamine and metabolites show maximum stability level.
4. Conclusion
This theoretical study suggests that a series composed of phosphoethanolamine and its metabolites exhibit low affinity to the enzyme AChE, suggesting low neurotoxicity during the use and biotransformation of this anticancer compound.
A docking study corroborated the results of experimental enzymatic inhibitory values of AChE and was executed in comparison to a group of 8 organophosphorus compounds whose pIC50 values ranged from 5.92 to 9.47. We observed that no compound of the phosphoethanolamine series showed energy levels close to the template group, indicating a weak interaction with AChE.
In light of this, we conclude that phosphoethanolamine and its predicted metabolites do not significantly inhibit AChE to cause a cholinergic crisis, as organophosphorus usually act, indicating a low neurotoxicity potential, highlighting the interest in further investigating this compound as a potential lead for future anticancer agents.
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
DiveC.HickmanJ. A.Drug-target interactions: only the first step in the commitment to a programmed cell death?19916411921962-s2.0-0025880971WHO2016, http://www.who.int/cancer/en/WHO2003, http://www.who.int/cancer/en/OrtizR.MelguizoC.PradosJ.ÁlvarezP. J.CabaO.SerranoF. R.HitaF.AránegaA.New gene therapy strategies for cancer treatment: a review of recent patents20127329731210.2174/1574892128018200932-s2.0-84867495468DesantisC. E.LinC. C.MariottoA. B.SiegelR. L.SteinK. D.KramerJ. L.AlteriR.RobbinsA. S.JemalA.Cancer treatment and survivorship statistics, 2014201264425227110.3322/caac.212352-s2.0-84904155440EllisonD. W.BealM. F.MartinJ. B.Phosphoethanolamine and ethanolamine are decreased in Alzheimer's disease and Huntington's disease1987417238939210.1016/0006-8993(87)90471-92-s2.0-0023618984EysterK. M.The membrane and lipids as integral participants in signal transduction: lipid signal transduction for the non-lipid biochemist200731151610.1152/advan.00088.20062-s2.0-33847378878FerreiraA. K.MenegueloR.NetoS. C.ChiericeG. O.MariaD. A.Synthetic phosphoethanolamine induces apoptosis through caspase-3 pathway by decreasing expression of Bax/Bad protein and changes cell cycle in melanoma201133535910.4172/1948-5956.10000582-s2.0-79952917187van der LuitA. H.VinkS. R.KlarenbeekJ. B.PerrissoudD.SolaryE.VerheijM.van BlitterswijkW. J.A new class of anticancer alkylphospholipids uses lipid rafts as membrane gateways to induce apoptosis in lymphoma cells2007682337234510.1158/1535-7163.mct-07-02022-s2.0-34548073046DankerK.ReutterW.SeminiG.Glycosidated phospholipids: uncoupling of signalling pathways at the plasma membrane20101601364710.1111/j.1476-5381.2009.00626.x2-s2.0-77950920621FerreiraA. K.Santana-LemosB. A. A.RegoE. M.FilhoO. M. R.ChiericeG. O.MariaD. A.Synthetic phosphoethanolamine has in vitro and in vivo anti-leukemia effects2013109112819282810.1038/bjc.2013.5102-s2.0-84888879701FerreiraA. K.MenegueloR.PereiraA.FilhoO. M. R.ChiericeG. O.MariaD. A.Synthetic phosphoethanolamine induces cell cycle arrest and apoptosis in human breast cancer MCF-7 cells through the mitochondrial pathway201367648148710.1016/j.biopha.2013.01.0122-s2.0-84880947374FerreiraA. K.MenegueloR.PereiraA.FilhoO. M. R.ChiericeG. O.MariaD. A.Anticancer effects of synthetic phosphoethanolamine on Ehrlich ascites tumor: an experimental study2012321951042-s2.0-84855711904FerreiraA. K.FreitasV. M.LevyD.RuizJ. L. M.BydlowskiS. P.RiciR. E. G.FilhoO. M. R.ChiericeG. O.MariaD. A.Anti-angiogenic and anti-metastatic activity of synthetic phosphoethanolamine2013833e5793710.1371/journal.pone.00579372-s2.0-84875044440SharmaA. K.GaurK.TiwariR. K.GaurM. S.Computational interaction analysis of organophosphorus pesticides with different metabolic proteins in humans201125533534710.1016/s1674-8301(11)60045-62-s2.0-80054869237RamanjaneyuluG.SrinivasuluC.SeenaiahR.RamyaA.ReddyB. R.Docking studies of Endosulfan and its analogues on Human Acetylcholine esterase2014424954GladiesKeziaJ.SewaliGhoshK.PriyakumariJ. C.Docking studies on acetylcholinesterase with methylparathion2013211LugokenskiT. H.GubertP.BuenoD. C.NogaraP. A.de Aquino SaraivaR.BarcelosR. P.CarratuV. S.BresolinL.de Vargas BarbosaN. B.PereiraM. E.da RochaJ. B. T.SoaresF. A. A.Effect of different oximes on rat and human cholinesterases inhibited by methamidophos: a comparative in vitro and in silico study2012111636237010.1111/j.1742-7843.2012.00912.x2-s2.0-84869200069ViswanadhanV. N.GhoseA. K.RevankarG. R.RobinsR. K.Atomic physicochemical parameters for three dimensional structure directed quantitative structure-activity relationships. 4. Additional parameters for hydrophobic and dispersive interactions and their application for an automated superposition of certain naturally occurring nucleoside antibiotics1989291631722-s2.0-0024716284CrucianiG.CrivoriP.CarruptP.-A.TestaB.Molecular fields in quantitative structure-permeation relationships: the VolSurf approach20005031-2173010.1016/s0166-1280(99)00360-72-s2.0-0034625096LorenzoV. P.FilhoJ. M. B.ScottiL.ScottiM. T.Combined structure- and ligand-based virtual screening to evaluate caulerpin analogs with potential inhibitory activity against monoamine oxidase B201525669069710.1016/j.bjp.2015.08.0052-s2.0-84954415710KlopmanG.LiJ.-Y.WangS.DimayugaM.Computer automated log P calculations based on an extended group contribution approach199434475278110.1021/ci00020a0092-s2.0-0028464119CsizmadiaF.Tsantili-KakoulidouA.PanderiI.DarvasF.Prediction of distribution coefficient from structure. 1. Estimation method199786786587110.1021/js960177k2-s2.0-0030788112RavesM. L.HarelM.PangY.-P.SilmanI.KozikowskiA. P.SussmanJ. L.Structure of acetylcholinesterase complexed with the nootropic alkaloid, (−)-huperzine A199741576310.1038/nsb0197-572-s2.0-0031015343ThomsenR.ChristensenM. H.MolDock: a new technique for high-accuracy molecular docking200649113315332110.1021/jm051197e2-s2.0-33744826819YazalJ. E.RaoS. N.MehlA.SlikkerW.Jr.Prediction of organophosphorus acetylcholinesterase inhibition using three-dimensional quantitative structure-activity relationship (3D-QSAR) methods200163222323210.1093/toxsci/63.2.2232-s2.0-0034779028