Next Article in Journal
Differential Stimulation of Post-Exercise Myofibrillar Protein Synthesis in Humans Following Isonitrogenous, Isocaloric Pre-Exercise Feeding
Next Article in Special Issue
High-Molecular-Weight Polyphenol-Rich Fraction of Black Tea Does Not Prevent Atrophy by Unloading, But Promotes Soleus Muscle Mass Recovery from Atrophy in Mice
Previous Article in Journal
Effect of Vasicinone against Paraquat-Induced MAPK/p53-Mediated Apoptosis via the IGF-1R/PI3K/AKT Pathway in a Parkinson’s Disease-Associated SH-SY5Y Cell Model
Previous Article in Special Issue
Polyphenols and IUGR Pregnancies: Effects of Maternal Hydroxytyrosol Supplementation on Hepatic Fat Accretion and Energy and Fatty Acids Profile of Fetal Tissues
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 11
Issue 7
10.3390/nu11071656
nutrients-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Extra Virgin Olive Oil Contains a Phenolic Inhibitor of the Histone Demethylase LSD1/KDM1A

by
Elisabet Cuyàs
1,2,
Juan Gumuzio
3,
Jesús Lozano-Sánchez
4,5,
David Carreras
2,6,
Sara Verdura
1,2,
Laura Llorach-Parés
7,
Melchor Sanchez-Martinez
7,
Elisabet Selga
6,8,9,
Guillermo J. Pérez
6,8,
Fabiana S. Scornik
6,8,
Ramon Brugada
2,6,8,10,
Joaquim Bosch-Barrera
2,11,12,
Antonio Segura-Carretero
4,5,
Ángel G. Martin
3,
José Antonio Encinar
13 and
Javier A. Menendez
1,2,*
1
ProCURE (Program Against Cancer Therapeutic Resistance), Metabolism & Cancer Group, Catalan Institute of Oncology, 17007 Girona, Spain
2
Girona Biomedical Research Institute (IDIBGI), 17190 Girona, Spain
3
StemTek Therapeutics, 48160 Bilbao, Spain
4
Department of Analytical Chemistry, Faculty of Sciences, University of Granada, 18071 Granada, Spain
5
Research and Development Functional Food Centre (CIDAF), PTS Granada, 18100 Granada, Spain
6
Cardiovascular Genetics Centre, Department of Medical Sciences, University of Girona, 17071 Girona, Spain
7
Mind the Byte, 08007 Barcelona, Spain
8
Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares (CIBERCV), 28029 Madrid, Spain
9
Faculty of Medicine, University of Vic-Central University of Catalonia (UVic-UCC), 08500 Vic, Spain
10
Dr. Josep Trueta Hospital of Girona, 17007 Girona, Spain
11
Medical Oncology, Catalan Institute of Oncology (ICO), 17007 Girona, Spain
12
Department of Medical Sciences, Medical School University of Girona, 17071 Girona, Spain
13
Institute of Research, Development and Innovation in Biotechnology of Elche (IDiBE) and Molecular and Cell Biology Institute (IBMC), Miguel Hernández University (UMH), 03202 Elche, Spain
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
Nutrients 2019, 11(7), 1656; https://0-doi-org.brum.beds.ac.uk/10.3390/nu11071656
Submission received: 5 July 2019 / Revised: 16 July 2019 / Accepted: 17 July 2019 / Published: 19 July 2019
(This article belongs to the Special Issue Polyphenol-Rich Foods for Human Health and Disease)
Browse Figures
Review Reports Versions Notes

Abstract

:
The lysine-specific histone demethylase 1A (LSD1) also known as lysine (K)-specific demethylase 1A (KDM1A) is a central epigenetic regulator of metabolic reprogramming in obesity-associated diseases, neurological disorders, and cancer. Here, we evaluated the ability of oleacein, a biophenol secoiridoid naturally present in extra virgin olive oil (EVOO), to target LSD1. Molecular docking and dynamic simulation approaches revealed that oleacein could target the binding site of the LSD1 cofactor flavin adenosine dinucleotide with high affinity and at low concentrations. At higher concentrations, oleacein was predicted to target the interaction of LSD1 with histone H3 and the LSD1 co-repressor (RCOR1/CoREST), likely disturbing the anchorage of LSD1 to chromatin. AlphaScreen-based in vitro assays confirmed the ability of oleacein to act as a direct inhibitor of recombinant LSD1, with an IC50 as low as 2.5 μmol/L. Further, oleacein fully suppressed the expression of the transcription factor SOX2 (SEX determining Region Y-box 2) in cancer stem-like and induced pluripotent stem (iPS) cells, which specifically occurs under the control of an LSD1-targeted distal enhancer. Conversely, oleacein failed to modify ectopic SOX2 overexpression driven by a constitutive promoter. Overall, our findings provide the first evidence that EVOO contains a naturally occurring phenolic inhibitor of LSD1, and support the use of oleacein as a template to design new secoiridoid-based LSD1 inhibitors.

Graphical Abstract

1. Introduction

Natural biophenols in extra virgin olive oil (EVOO), a health-promoting ingredient of the so-called Mediterranean Diet, have been proposed to have multifaceted and beneficial effects on major diseases, including cancer, diabetes, cardiovascular diseases, neurodegenerative disease, and metabolic disorders [1,2,3,4,5,6,7,8,9]. Although the mechanism(s) of action of EVOO biophenols in these settings is not yet fully understood, there is emerging evidence suggesting that phenolic molecules may exert nutritional and beneficial actions through multiple molecular mechanisms involving not only general biochemical pathways but also specific epigenetic remodeling [10,11].
We have recently used a reverse pharmacology computational approach to “de-orphanize” endogenous ligands and assign molecular functions to the biophenol oleacein [12]—one of the most abundant secoiridoids in EVOO that seems to be largely responsible for its health-promoting effects [6]. Our virtual profiling approach employing structure-based tools revealed a target landscape of oleacein significantly enriched in metabolic and chromatin-modifying enzymes involved in histone post-translational modification. Among the latter group was lysine-specific demethylase-1 (LSD1, also known as KDM1A), a flavin adenine dinucleotide (FAD)-dependent homolog of the amine oxidase family that was originally described to target histone H3K4 [13,14]. LSD1 has been linked to a repression of gene expression, and was later reported to interact with the androgen or estrogen receptor to demethylate histone H3K9 [15,16,17], resulting in the activation of gene expression.
LSD1 plays a central epigenetic role in nutrient-driven metabolic adaptation and reprogramming in multifactorial diseases, including obesity-associated diseases, neurological disorders, and cancer [18,19,20,21,22,23,24,25]. Accordingly, natural and synthetic chemicals capable of blocking LSD1 might represent a new class of epigenetic drugs to treat a wide range of metabolic disorders. In this regard, some of the beneficial effects promoted by polyphenols, such as resveratrol, curcumin, quercetin, and baicalin, could be explained, at least in part, by their ability to operate as natural inhibitors of LSD1 [26,27,28,29]. Here, we combined in silico docking and molecular dynamics approaches coupled to experimental validation to identify LSD1 as a novel epigenetic target of the bioactive EVOO polyphenol oleacein.

2. Materials and Methods

2.1. Docking Calculations, Molecular Dynamics Simulations, and Binding Free Energy Analysis

The human histone demethylase LSD1 (UniProt code O60341)/REST corepressor 1 (UniProt code Q9UKL0)/histone H3 peptide (UniProt code P68431) ternary complex structures (PDB codes 2IW5 and 2X0L) were obtained from the Research Collaboratory for Structural Bioinformatics Protein Data Bank (PDB). The structure of oleacein (PubChem CID: 18684078) was obtained from the National Center for Biotechnology Information (NCBI) PubChem database (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pccompound). The specific edition of the LSD1 protein structure involving the removal of water, FAD, and histone H3 peptide, was made using PyMol 2.0 software (PyMOL Molecular Graphics System, v2.0 Schrödinger, LLC, at http://www.pymol.org/) without further optimization.
Docking calculations, molecular dynamics (MDs) simulations, and molecular mechanics/generalized borne surface area (MM/GBSA) calculations to determine the alchemical binding free energy of oleacein against LSD1 2IW5 were performed as previously described [10,11,12,30]. Docking calculations, MD simulations, and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations to determine the alchemical binding free energy of oleacein against LSD1 2X0L were performed as previously described [31]. All of the figures were prepared using PyMol 2.0 software and all interactions were detected using the protein-ligand interaction profiler (PLIP) algorithm [32].

2.2. Oleacein Isolation and Purification

Oleacein (decarboxymethyl oleuropein aglycone) was isolated and purified from the phenolic fraction of EVOO as previously described [10].

2.3. LSD1 Enzymatic Activity

Enzymatic reactions were performed in an AlphaScreen format in duplicate at room temperature for 60 min in a 10 μL mixture containing assay buffer, histone H3 peptide substrate, LSD1 (BPS#50103, lot#130806-D) enzyme, and oleacein. The 10-μL reactions were carried out in 384-well Optiplates (Perkin Elmer Life Sciences, Waltham, MA). A serial dilution of the compounds was first performed in 3.3% dimethylsulfoxide (DMSO)/assay buffer. From this step, 3 μL of oleacein was added to 4 μL of enzyme and incubated for 30 min at room temperature. After this incubation, 3 μL of substrate was added to initiate the reaction. The final DMSO concentration was 1%. After the reaction, 5 μL of anti-mouse acceptor beads (Perkin Elmer, diluted 1:500 with 1× detection buffer) or 5 μL of anti-rabbit acceptor beads (Perkin Elmer, diluted 1:500 with 1× detection buffer) and 5 μL of primary antibody (BPS#52140E,F, diluted 1:200 with 1× detection buffer) were added to the reaction mix. After brief shaking, the plate was incubated for 30 min. Finally, 10 μL of AlphaScreen streptavidin-conjugated donor beads (Perkin Elmer, diluted 1:125 with 1× detection buffer) were added. After 30 min, the samples were measured on an AlphaScreen microplate reader (EnSpire Alpha 2390 Multilabel reader, Perkin Elmer).
The AlphaScreen intensity data were analyzed and compared using Graphpad Prism software (GraphPad Software Inc., San Diego, CA, USA). In the absence of iadademstat, the AlphaScreen or fluorescence intensity (Ft) was defined as 100% activity. In the absence of enzyme, the intensity (Fb) was defined as 0% activity. The percent activity in the presence of oleacein was calculated according to the following equation: %activity = (F − Fb)/(Ft − Fb), where F = the A-screen intensity in the presence of iadademstat. Once A-screen data were converted to LSD1 activity (%), those values were then plotted against a series of oleacein concentrations using non-linear regression analysis of sigmoidal dose–response curves generated with the equation Y = B + (T − B)/1 + 10^((LogIC50-Z) × HillSlope), where Y = percent activity, B = minimum percent activity, T = maximum percent activity, Z = logarithm of oleacein concentration, and Hill Slope = slope factor or Hill coefficient. The IC50 value was determined as the concentration of oleacein causing a half-maximal inhibition of control activity.

2.4. Cell Lines and Culture Conditions

MCF-7 and BT-474 breast carcinoma cell lines were obtained directly from the ATCC (Manassas, VA, USA) and grown in Dulbecco’s modified Eagle’s medium (Gibco, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum (Sigma) at 37 °C in a 5% CO2 incubator. MCF-7 cells engineered to stably overexpress SOX2 were prepared and cultured as previously described [33]. Reprogramming of human fibroblasts from a healthy individual to induced pluripotent stem (iPS) cells was carried out as previously described [34].

2.5. SOX2 Enhancer Reporter Assay

BT-474 cells were transfected with 5 μg of pGL3 Luc control (Promega, Madison, WI, USA) or pGL2-Sox2-enhancer-Luc reporter plasmids [34,35,36] using Lipofectamine Plus (Invitrogen, Carslbad, CA, USA). Twenty-four hours after transfection, the culture was split into two parts: One part was seeded in two-dimensional adherent culture plates and the other part was cultured in non-adherent culture conditions to allow mammosphere formation, in the absence or presence of graded concentrations of oleacein (i.e., 0, 3, 10, and 30 μmol/L). After 48 h, cells were harvested and luciferase activity was measured in duplicate with the Glomax 20/20 luminometer (Promega) and normalized by protein concentration in the extracts. Results are expressed as the fold induction of sphere culture reporter activity above adherent culture control.

2.6. SOX2 Protein Expression

To evaluate the expression status of SOX2, cells grown to 80% to 90% confluence were exposed to graded concentrations of oleacein (i.e., 0, 5, 10, 25, and 50 μmol/L) for 24 h. Cells were then collected and subjected to western blot analysis with a monoclonal antibody against SOX2 (Cat.# 3579, Cell Signaling). β-actin (Cat.# 66009-1-Ig, Clone No. 2D4H5) was used as loading control.

2.7. SOX2 Gene Expression

Total RNA was extracted from cells using a Nucleospin RNA plus kit (Macherey-Nagel GmbH & Co. KG Düren, Germany). One microgram of total RNA was reverse-transcribed into cDNA using the High Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific, Carlsbad, CA). RNA concentration and quality were determined with an ND-1000 spectrophotometer (NanoDrop™ ND-1000, NanoDrop Technologies, Wilmington, DE).
cDNA (20 ng) was assayed in triplicate according to established protocols using a QuantStudio™ 7 Flex Real-Time PCR system (Thermo Fisher Scientific, Waltham, MA) with an automated baseline and threshold cycle detection. PPIA and 18S were used as reference genes. Primers and fluorescent probes for PPIA, 18S, and SOX2 were obtained from Thermo Fisher Scientific (TaqMan Gene Expression assay IDs: Hs99999904_m1, Hs99999901_s1, and Hs01053049_s1, respectively). Data were analyzed using the Thermo Fisher Cloud software.

2.8. Statistical Analysis

All statistical analyses were performed using GraphPad Prism software. Data are presented as mean ± S.D. Comparisons of means of ≥3 groups were performed by analysis of variance (ANOVA) and the existence of individual differences, in case of significant F values at ANOVA, were assessed by multiple contrasts. p values < 0.05 were considered to be statistically significant (denoted as *). All statistical tests were two-sided.

3. Results

3.1. Computational Prediction of LSD1 as a Putative Target of Oleacein

One of the epigenetic targets predicted to be engaged by oleacein in our recent chemoinformatics approach [12] was PDB ID 2IW5 [37], which corresponds to a complex containing LSD1 (KDM1A, Uniprot O60341) and REST Corepressor 1 (RCOR1/CoREST), a co-repressor that cooperates with LSD1 to demethylate mono- and di-methylated lysine 4 (K4) of histone H3 in nucleosomes. A classic rigid docking calculation, which was performed twice to avoid false positives, over the structure employing the cavity defined by its crystallographic ligand (FAD) revealed an energy binding of −7.2/−7.5 kcal/mol (Table 1; Figure 1A, top). Blind docking calculations involving cavity searching and docking calculations over the established cavities confirmed that the best predicted cavity for oleacein in LSD1 was the same as that containing FAD (−6.1 kcal/mol; Table 1; Figure 1A, bottom). To add protein flexibility to the analysis and to test the stability of the predicted oleacein-LSD1 complex, we performed short molecular dynamics (MDs) simulations of 1 ns to filter out the possibility that the predicted LSD1–oleacein complex would involve a poorly interacting pose (Figure 1B). We then performed a molecular mechanics/generalized born surface area (MM/GBSA) approximation of the oleacein–LSD1 binding affinity, which considers the dynamic nature of LSD1 and therefore provides a more realistic view of oleacein binding affinity than the aforementioned rigid rocking estimations. The results highlighted a high affinity of oleacein at the FAD cavity of LSD1, which reached −39.9541 kcal/mol.
When assessing the most relevant interactions between oleacein and LSD1 and the residues involved in such binding when placed over the FAD cavity, it was predicted that among the residues involved in the FAD-binding mode (K661, Y761, A809, and T801), two of them (K661 and Y761) were shared by the oleacein-binding mode. The stability of oleacein in the FAD-binding pocket was predicted to involve hydrogen bond interactions with K661, A331, and M332, as well as hydrophobic contacts with R316, V333, F538, and V811. Oleacein was predicted to establish ∏-cation interactions with W751. Beyond a predicted LSD1 inhibitory mode of oleacein via occupancy of the FAD site, it is noteworthy that the oleacein interacting residues F538 and V811 are placed near the residues involved in the binding of histone H3 (N540, P808, and A809), thereby suggesting that the inhibitory capacity of oleacein to LSD1 could also involve a disturbance of the anchorage of LSD1 to chromatin.

3.2. Computational Deconstruction of the Binding Mode of Oleacein to the LSD1-CoREST-histone H3 Complex

To computationally validate the virtual profiling prediction with 2IW5, we re-evaluated the binding mode of oleacein to the LSD1-CoREST-histone H3 complex using a different crystal structure—namely the 2X0L ternary complex structure [38]. Docking simulations of oleacein in such a crystal structure of human LSD1 (chain A), including RCOR1/CoREST (chain B), and a histone H3 peptide (chain C), produced eight clusters of docking poses (Figure 2).
When the docking results were ranked according to the ascent of the binding energies for oleacein (up to −9.63 kcal/mol; Table 2), the clusters exhibiting the highest affinity (#1 and #2) were in the nanomolar range and were predicted to interact with those residues delineating the FAD cavity of LSD1. This was the case particularly for cluster #2 (−9.58 kcal/mol), in which all four residues involved in the FAD-binding mode (K661, Y761, A809, and T801) were shared by the oleacein-binding mode. Cluster #3 (−9.33 kcal/mol), in the sub-micromolar range, was predicted to interfere with the anchorage of LSD1 to chromatin as it included two (P808 and A809) out of the three residues involved in the binding of histone H3. Clusters #7 and #8, with affinities in the low micromolar range, were predicted to interact with both the LSD1 enzyme and RCOR1/CoREST.
To add protein flexibility and provide additional information about different intra- and inter-molecular movements, we performed short MD simulations over the course of 10 ns together with binding free energy calculations under the molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) approximation (Table 2). With the exception of clusters #4 and #5, the MM-PBSA-based results highlighted the predicted extremely high affinity of oleacein for the FAD-binding pocket of LSD1 at the cluster #2, which reached −155.861 kcal/mol. For the remaining clusters, the range of binding energies was notably high (up to −66.844 kcal/mol in cluster #1), with oleacein trajectory changes of less than 4 Å (data not shown). The high binding energies of clusters #7 and #8 (−40.682 and −30.051kcal/mol, respectively), with affinities in the low micromolar range, are especially relevant as they predict that oleacein might stably bind and obstruct the interaction between LSD1 and RCOR1/CoREST, likely disturbing the anchorage of LSD1 to chromatin.

3.3. Oleacein Inhibits the In Vitro Enzymatic Activity of LSD1

Using an AlphaScreen assay with a biotinylated histone H3 peptide substrate, purified LSD1, and a highly specific antibody that recognizes demethylated substrate, we experimentally confirmed the model predictions of oleacein to target LSD1 (Figure 3). Accordingly, oleacein treatment dose-dependently suppressed the demethylase activity of LSD1 in vitro, with a mean IC50 of ~2.5 μmol/L.

3.4. Oleacein Suppresses LSD1-Targeted Enhancer-Driven Activation of the Stemness Transcription Factor SOX2 in Stem-Like Cells

LSD1-targeted inhibitors are known to differentially target pluripotent cancer cells, including teratocarcinoma, embryonic carcinoma, and seminoma or embryonic stem cells that express the core stem cell transcriptional factor SOX2, while having minimal growth-inhibitory effects on non-pluripotent or normal somatic cells [35,36,37,38,39,40,41,42]. Because this phenomenon can be explained in terms of the ability of LSD1 to specifically target a distal enhancer that controls SOX2 expression in stem-like cells, we hypothesized that oleacein should prevent the overexpression of SOX2 in breast cancer stem cells (CSC) via LSD1-regulated re-activation of this enhancer in the SOX2 promoter [33,35,36,37,38,39,40,41,42]. As expected, when we transfected BT-474 breast cancer cells with a luciferase reporter vector containing the SOX2 distal enhancer region centered between −3444 and −3833 [34], we observed a robust induction (>9-fold in average) in reporter activity in mammosphere cultures when compared with adherent differentiated control cultures (Figure 4). Of note, the LSD1-targeted, enhancer-driven transcriptional re-activation of SOX2 differentially occurring in mammosphere cultures was fully suppressed, in a dose-dependent manner, following exposure to graded concentrations of oleacein (Figure 4).
We then took advantage of two different SOX2-overexpressing cellular models to substantiate the ability of oleacein to specifically target LSD1-driven stemness-related SOX2 expression. We first employed induced pluripotent stem (iPS) cells in which SOX2 is robustly transcribed from the LSD1-targeted distal enhancer to maintain pluripotency [44,45,46]. We then employed MCF-7 cells engineered to drive stable expression of SOX2 through transfer of the mCitrinie-P2A-Sox2 lentiviral vector [33,47,48], in which SOX2 overexpression is under the transcriptional control of the human phosphoglycerate kinase promoter. Oleacein treatment dose-dependently suppressed SOX2 protein expression in iPS cells; however, oleacein failed to modify ectopic SOX2 overexpression driven by a constitutive promoter in MCF-7/SOX2 cells (Figure 5). Transcriptional analyses confirmed the ability of oleacein to exclusively down-regulate endogenous SOX2 expression in iPS cells but not in MCF-7/SOX2 cells expressing an exogenous SOX2 transgene (Figure 5).

4. Discussion

In the present study, we used several approaches to identify the lysine-specific histone demethylase LSD1 as a novel target for oleacein, a characteristic bioactive biophenol-secoiridoid in EVOO.
First, we provide computational evidence supporting a catalytic/enzymatic inhibitory mechanism of action likely involving highly potent and direct targeting of the FAD cofactor cavity of LSD1 at physiological concentrations of oleacein. Because the interaction of LSD1 with RCOR1/CoREST allows the formation of an H3-tail binding-site recognizing a 20-amino-acid portion of the histone H3 N-terminal region [49,50], we cannot exclude the possibility that higher concentrations of oleacein might additionally target LSD1 via a scaffolding/structural mechanism involving disruption of the chromatin-binding activity of the LSD1/RCOR1/CoREST complex. Upon successful experimental validation of computational predictions, such information should be useful in the development of new secoiridoids-based structural inhibitors of LSD1.
Second, using an in vitro system with recombinant human LSD1 and a short peptide of methylated H3, we found that oleacein inhibits the enzymatic activity of LSD1 in a dose-dependent fashion, with an IC50 value (2.5 μmol/L) that can be achieved physiologically given that oleacein can be absorbed into systemic circulation, where it reaches concentrations up to 20 μmol/L [6].
Third, a pharmacological blocker of LSD1 is expected to promote methylation of H3K4 and H3K9 within the regulatory (enhancer) region of the SOX2 promoter in pluripotent stem-like cells [40]. Accordingly, over the same concentration range that inhibited LSD1 activity in vitro, oleacein treatment efficiently blocked the re-activation of SOX2 that is known to be controlled by the binding of LSD1 to the SOX2 distal enhancer differentially occurring in tumor cells with CSC-like properties. Because oleacein was previously shown to exhibit strong suppressive effects against tumor-initiating CSC [10]—a sub-group of cancer cells having an enhanced ability to initiate tumorigenesis and drive therapeutic relapse—and considering the capacity of LSD1 to confer stem-cell-like traits to cancer cells [51,52], it is conceivable that some of the anti-CSC effects of oleacein could occur through LSD1 inhibition.
Given the central role of LSD1 for sensing and integrating energetic information into the epigenome [18,19,20,21,22,23,24,25,53,54], the LSD1-targeted inhibitory activity of oleacein might represent a novel mechanism through which EVOO phenolics impact the crosstalk between the metabolism and epigenetic gene regulation to exert health-promoting effects in a wide variety of human diseases, including obesity-associated diseases, neurological disorders, and cancer (Figure 6).

5. Conclusions

Overall, our study implicates a novel mechanism of action for the EVOO bioactive, phenolic oleacein, adding to the existing list of natural polyphenols capable of inhibiting LSD1, and laying the groundwork for the design of new secoiridoids-based structural inhibitors of LSD1.

Author Contributions

Conceptualization, writing—original draft preparation, writing—review and editing J.A.M., Á.G.M., E.C., and J.A.E.; supervision, funding acquisition J.A.M.; software, methodology Á.G.M., L.L.-P., M.S.-M., and J.A.E.; investigation, formal analysis, visualization, data curation, validation E.C., J.G., D.C., and S.V.; resources J.L.-S., E.S., G.J.P., F.S.S., R.B., and A.S.-C. All authors contributed to the development of this manuscript and read and approved the final version.

Funding

Work in the Menendez laboratory is supported by the Spanish Ministry of Science and Innovation (Grant SAF2016-80639-P, Plan Nacional de l+D+I, founded by the European Regional Development Fund, Spain) and by an unrestricted research grant from the Fundació Oncolliga Girona (Lliga catalana d’ajuda al malalt de càncer, Girona). The Spanish Ministry of Economy and Competitiveness (MINECO, Project RTI2018-096724-B-C21) and the Generalitat Valenciana (PROMETEO/2016/006) supports work in the Encinar laboratory.

Acknowledgments

We are grateful to the Cluster of Scientific Computing (http://ccc.umh.es/) of the Miguel Hernández University (UNH) and the Centro de Supercomputación of the University of Granada (ALHAMBRA-CSIRC) for providing computing facilities. The authors would like to thank Dr. Kenneth McCreath for editorial support. This work is in loving memory of Àlex Casta, a generous, kind, and cheerful friend and colleague who passed away in December 2018.

Conflicts of Interest

Stock ownership: Á.G.M., StemTek Therapeutics (CEO). All other authors have no competing interests to declare. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Availability of Data and Materials

All data generated or analyzed during this study are included in this published article.

References

  1. Menendez, J.A.; Joven, J.; Aragonès, G.; Barrajón-Catalán, E.; Beltrán-Debón, R.; Borrás-Linares, I.; Camps, J.; Corominas-Faja, B.; Cufí, S.; Fernández-Arroyo, S.; et al. Xenohormetic and anti-aging activity of secoiridoid polyphenols present in extra virgin olive oil: A new family of gerosuppressant agents. Cell Cycle 2013, 12, 555–578. [Google Scholar] [CrossRef] [PubMed]
  2. Rigacci, S.; Stefani, M. Nutraceutical Properties of Olive Oil Polyphenols. An Itinerary from Cultured Cells through Animal Models to Humans. Int. J. Mol. Sci. 2016, 17, 843. [Google Scholar] [CrossRef] [PubMed]
  3. Angeloni, C.; Malaguti, M.; Barbalace, M.C.; Hrelia, S. Bioactivity of Olive Oil Phenols in Neuroprotection. Int. J. Mol. Sci. 2017, 18, 2230. [Google Scholar] [CrossRef] [PubMed]
  4. Peyrol, J.; Riva, C.; Amiot, M.J. Hydroxytyrosol in the Prevention of the Metabolic Syndrome and Related Disorders. Nutrients 2017, 9, 306. [Google Scholar] [CrossRef] [PubMed]
  5. Saibandith, B.; Spencer, J.P.E.; Rowland, I.R.; Commane, D.M. Olive Polyphenols and the Metabolic Syndrome. Molecules 2017, 22, 1082. [Google Scholar] [CrossRef] [PubMed]
  6. Lombardo, G.E.; Lepore, S.M.; Morittu, V.M.; Arcidiacono, B.; Colica, C.; Procopio, A.; Maggisano, V.; Bulotta, S.; Costa, N.; Mignogna, C.; et al. Effects of Oleacein on High-Fat Diet-Dependent Steatosis, Weight Gain, and Insulin Resistance in Mice. Front. Endocrinol. 2018, 9, 116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Reboredo-Rodríguez, P.; Varela-López, A.; Forbes-Hernández, T.Y.; Gasparrini, M.; Afrin, S.; Cianciosi, D.; Zhang, J.; Manna, P.P.; Bompadre, S.; Quiles, J.L.; et al. Phenolic Compounds Isolated from Olive Oil as Nutraceutical Tools for the Prevention and Management of Cancer and Cardiovascular Diseases. Int. J. Mol. Sci. 2018, 19, 2305. [Google Scholar] [CrossRef] [PubMed]
  8. Pedret, A.; Fernández-Castillejo, S.; Valls, R.M.; Catalán, Ú.; Rubió, L.; Romeu, M.; Macià, A.; López de Las Hazas, M.C.; Farràs, M.; Giralt, M.; et al. Cardiovascular Benefits of Phenol-Enriched Virgin Olive Oils: New Insights from the Virgin Olive Oil and HDL Functionality (VOHF) Study. Mol. Nutr. Food Res. 2018, 62, 1800456. [Google Scholar] [CrossRef]
  9. Celano, M.; Maggisano, V.; Lepore, S.M.; Russo, D.; Bulotta, S. Secoiridoids of olive and derivatives as potential coadjuvant drugs in cancer: A critical analysis of experimental studies. Pharmacol. Res. 2019, 142, 77–86. [Google Scholar] [CrossRef]
  10. Corominas-Faja, B.; Cuyàs, E.; Lozano-Sánchez, J.; Cufí, S.; Verdura, S.; Fernández-Arroyo, S.; Borrás-Linares, I.; Martin-Castillo, B.; Martin, Á.G.; Lupu, R.; et al. Extra-virgin olive oil contains a metabolo-epigenetic inhibitor of cancer stem cells. Carcinogenesis 2018, 39, 601–613. [Google Scholar] [CrossRef] [Green Version]
  11. Verdura, S.; Cuyàs, E.; Lozano-Sánchez, J.; Bastidas-Velez, C.; Llorach-Parés, L.; Fernández-Arroyo, S.; Hernández-Aguilera, A.; Joven, J.; Nonell-Canals, A.; Bosch-Barrera, J.; et al. An olive oil phenolic is a new chemotype of mutant isocitrate dehydrogenase 1 (IDH1) inhibitors. Carcinogenesis 2019, 40, 27–40. [Google Scholar] [CrossRef] [PubMed]
  12. Cuyàs, E.; Castillo, D.; Llorach-Parés, L.; Lozano-Sánchez, J.; Verdura, S.; Nonell-Canals, A.; Brunet, J.; Bosch-Barrera, J.; Joven, J.; Valdés, R.; et al. Computational de-orphanization of the olive oil biophenols oleacein: Discovery of new metabolic and epigenetic targets. Food Chem. Toxicol. 2019, 131, 110529. [Google Scholar] [CrossRef] [PubMed]
  13. Shi, Y.; Lan, F.; Matson, C.; Mulligan, P.; Whetstine, J.R.; Cole, P.A.; Casero, R.A.; Shi, Y. Histone demethylation mediated by the nuclear amine oxidase homolog LSD1. Cell 2004, 119, 941–953. [Google Scholar] [CrossRef] [PubMed]
  14. Lan, F.; Nottke, A.C.; Shi, Y. Mechanisms involved in the regulation of histone lysine demethylases. Curr. Opin. Cell Biol. 2008, 20, 316–325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Metzger, E.; Wissmann, M.; Yin, N.; Müller, J.M.; Schneider, R.; Peters, A.H.; Günther, T.; Buettner, R.; Schüle, R. LSD1 demethylates repressive histone marks to promote androgen-receptor-dependent transcription. Nature 2005, 437, 436–439. [Google Scholar] [CrossRef] [PubMed]
  16. Garcia-Bassets, I.; Kwon, Y.S.; Telese, F.; Prefontaine, G.G.; Hutt, K.R.; Cheng, C.S.; Ju, B.G.; Ohgi, K.A.; Wang, J.; Escoubet-Lozach, L.; et al. Histone methylation-dependent mechanisms impose ligand dependency for gene activation by nuclear receptors. Cell 2007, 128, 505–518. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, J.; Scully, K.; Zhu, X.; Cai, L.; Zhang, J.; Prefontaine, G.G.; Krones, A.; Ohgi, K.A.; Zhu, P.; Garcia-Bassets, I.; et al. Opposing LSD1 complexes function in developmental gene activation and repression programmes. Nature 2007, 446, 882–887. [Google Scholar] [CrossRef]
  18. Forneris, F.; Binda, C.; Battaglioli, E.; Mattevi, A. LSD1: Oxidative chemistry for multifaceted functions in chromatin regulation. Trends Biochem. Sci. 2008, 33, 181–189. [Google Scholar] [CrossRef]
  19. Abdulla, A.; Zhang, Y.; Hsu, F.N.; Xiaoli, A.M.; Zhao, X.; Yang, E.S.; Ji, J.Y.; Yang, F. Regulation of lipogenic gene expression by lysine-specific histone demethylase-1 (LSD1). J. Biol. Chem. 2014, 289, 29937–29947. [Google Scholar] [CrossRef]
  20. Duteil, D.; Metzger, E.; Willmann, D.; Karagianni, P.; Friedrichs, N.; Greschik, H.; Günther, T.; Buettner, R.; Talianidis, I.; Metzger, D.; et al. LSD1 promotes oxidative metabolism of white adipose tissue. Nat. Commun. 2014, 5, 4093. [Google Scholar] [CrossRef]
  21. Duteil, D.; Tosic, M.; Lausecker, F.; Nenseth, H.Z.; Müller, J.M.; Urban, S.; Willmann, D.; Petroll, K.; Messaddeq, N.; Arrigoni, L.; et al. Lsd1 Ablation Triggers Metabolic Reprogramming of Brown Adipose Tissue. Cell Rep. 2016, 17, 1008–1021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Maes, T.; Mascaró, C.; Ortega, A.; Lunardi, S.; Ciceri, F.; Somervaille, T.C.; Buesa, C. KDM1 histone lysine demethylases as targets for treatments of oncological and neurodegenerative disease. Epigenomics 2015, 7, 609–626. [Google Scholar] [CrossRef]
  23. Sakamoto, A.; Hino, S.; Nagaoka, K.; Anan, K.; Takase, R.; Matsumori, H.; Ojima, H.; Kanai, Y.; Arita, K.; Nakao, M. Lysine Demethylase LSD1 Coordinates Glycolytic and Mitochondrial Metabolism in Hepatocellular Carcinoma Cells. Cancer Res. 2015, 75, 1445–1456. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, H.; Liang, L.; Li, Y.; Feng, C.; Li, L.; Zhang, Y.; He, S.; Pei, D.; Guo, Y.; Zheng, H. Lysine-specific histone demethylase 1 inhibition promotes reprogramming by facilitating the expression of exogenous transcriptional factors and metabolic switch. Sci. Rep. 2016, 6, 30903. [Google Scholar] [CrossRef] [Green Version]
  25. Abdulla, A.; Zhao, X.; Yang, F. Natural Polyphenols Inhibit Lysine-Specific Demethylase-1 in vitro. J. Biochem. Pharmacol. Res. 2013, 1, 56–63. [Google Scholar] [PubMed]
  26. Duan, Y.C.; Guan, Y.Y.; Zhai, X.Y.; Ding, L.N.; Qin, W.P.; Shen, D.D.; Liu, X.Q.; Sun, X.D.; Zheng, Y.C.; Liu, H.M. Discovery of resveratrol derivatives as novel LSD1 inhibitors: Design, synthesis and their biological evaluation. Eur. J. Med. Chem. 2017, 126, 246–258. [Google Scholar] [CrossRef] [PubMed]
  27. Zheng, Y.C.; Shen, D.D.; Ren, M.; Liu, X.Q.; Wang, Z.R.; Liu, Y.; Zhang, Q.N.; Zhao, L.J.; Zhao, L.J.; Ma, J.L.; et al. Baicalin, a natural LSD1 inhibitor. Bioorg. Chem. 2016, 69, 129–131. [Google Scholar] [CrossRef] [PubMed]
  28. Fernandes, G.F.S.; Silva, G.D.B.; Pavan, A.R.; Chiba, D.E.; Chin, C.M.; Dos Santos, J.L. Epigenetic Regulatory Mechanisms Induced by Resveratrol. Nutrients 2017, 9, 1201. [Google Scholar] [CrossRef]
  29. Nettore, I.C.; Rocca, C.; Mancino, G.; Albano, L.; Amelio, D.; Grande, F.; Puoci, F.; Pasqua, T.; Desiderio, S.; Mazza, R.; et al. Quercetin and its derivative Q2 modulate chromatin dynamics in adipogenesis and Q2 prevents obesity and metabolic disorders in rats. J. Nutr. Biochem. 2019, 69, 151–162. [Google Scholar] [CrossRef]
  30. Cuyàs, E.; Verdura, S.; Llorach-Parés, L.; Fernández-Arroyo, S.; Joven, J.; Martin-Castillo, B.; Bosch-Barrera, J.; Brunet, J.; Nonell-Canals, A.; Sanchez-Martinez, M.; et al. Metformin Is a Direct SIRT1-Activating Compound: Computational Modeling and Experimental Validation. Front. Endocrinol. 2018, 9, 657. [Google Scholar] [CrossRef] [Green Version]
  31. Cuyàs, E.; Verdura, S.; Lozano-Sánchez, J.; Viciano, I.; Llorach-Parés, L.; Nonell-Canals, A.; Bosch-Barrera, J.; Brunet, J.; Segura-Carretero, A.; Sanchez-Martinez, M.; et al. The extra virgin olive oil phenolic oleacein is a dual substrate-inhibitor of catechol-O-methyltransferase. Food Chem. Toxicol. 2019, 128, 35–45. [Google Scholar] [CrossRef]
  32. Salentin, S.; Schreiber, S.; Haupt, V.J.; Adasme, M.F.; Schroeder, M. PLIP: Fully automated protein-ligand interaction profiler. Nucleic Acids Res. 2015, 43, W443–W447. [Google Scholar] [CrossRef]
  33. Leis, O.; Eguiara, A.; Lopez-Arribillaga, E.; Alberdi, M.J.; Hernandez-Garcia, S.; Elorriaga, K.; Pandiella, A.; Rezola, R.; Martin, A.G. Sox2 expression in breast tumours and activation in breast cancer stem cells. Oncogene 2012, 31, 1354–1365. [Google Scholar] [CrossRef]
  34. Selga, E.; Sendfeld, F.; Martinez-Moreno, R.; Medine, C.N.; Tura-Ceide, O.; Wilmut, S.I.; Pérez, G.J.; Scornik, F.S.; Brugada, R.; Mills, N.L. Sodium channel current loss of function in induced pluripotent stem cell-derived cardiomyocytes from a Brugada syndrome patient. J. Mol. Cell. Cardiol. 2018, 114, 10–19. [Google Scholar] [CrossRef]
  35. Takanaga, H.; Tsuchida-Straeten, N.; Nishide, K.; Watanabe, A.; Aburatani, H.; Kondo, T. Gli2 is a novel regulator of sox2 expression in telencephalic neuroepithelial cells. Stem Cells 2009, 27, 165–174. [Google Scholar] [CrossRef]
  36. Yang, M.; Gocke, C.B.; Luo, X.; Borek, D.; Tomchick, D.R.; Machius, M.; Otwinowski, Z.; Yu, H. Structural basis for CoREST-dependent demethylation of nucleosomes by the human LSD1 histone demethylase. Mol. Cell 2006, 23, 377–387. [Google Scholar] [CrossRef]
  37. Zibetti, C.; Adamo, A.; Binda, C.; Forneris, F.; Toffolo, E.; Verpelli, C.; Ginelli, E.; Mattevi, A.; Sala, C.; Battaglioli, E. Alternative splicing of the histone demethylase LSD1/KDM1 contributes to the modulation of neurite morphogenesis in the mammalian nervous system. J. Neurosci. 2010, 30, 2521–2532. [Google Scholar] [CrossRef]
  38. Wang, J.; Lu, F.; Ren, Q.; Sun, H.; Xu, Z.; Lan, R.; Liu, Y.; Ward, D.; Quan, J.; Ye, T.; et al. Novel histone demethylase LSD1 inhibitors selectively target cancer cells with pluripotent stem cell properties. Cancer Res. 2011, 71, 7238–7249. [Google Scholar] [CrossRef]
  39. Zhang, X.; Lu, F.; Wang, J.; Yin, F.; Xu, Z.; Qi, D.; Wu, X.; Cao, Y.; Liang, W.; Liu, Y.; et al. Pluripotent stem cell protein Sox2 confers sensitivity to LSD1 inhibition in cancer cells. Cell Rep. 2013, 5, 445–457. [Google Scholar] [CrossRef]
  40. Yin, F.; Lan, R.; Zhang, X.; Zhu, L.; Chen, F.; Xu, Z.; Liu, Y.; Ye, T.; Sun, H.; Lu, F.; et al. LSD1 regulates pluripotency of embryonic stem/carcinoma cells through histone deacetylase 1-mediated deacetylation of histone H4 at lysine 16. Mol. Cell. Biol. 2014, 34, 158–179. [Google Scholar] [CrossRef]
  41. Iglesias, J.M.; Leis, O.; Pérez Ruiz, E.; Gumuzio Barrie, J.; Garcia-Garcia, F.; Aduriz, A.; Beloqui, I.; Hernandez-Garcia, S.; Lopez-Mato, M.P.; Dopazo, J.; et al. The Activation of the Sox2 RR2 Pluripotency Transcriptional Reporter in Human Breast Cancer Cell Lines is Dynamic and Labels Cells with Higher Tumorigenic Potential. Front. Oncol. 2014, 4, 308. [Google Scholar] [CrossRef]
  42. Iglesias, J.M.; Gumuzio, J.; Martin, A.G. Linking Pluripotency Reprogramming and Cancer. Stem Cells Transl. Med. 2017, 6, 335–339. [Google Scholar] [CrossRef]
  43. Moreno-Manzano, V.; Rodríguez-Jiménez, F.J.; Aceña-Bonilla, J.L.; Fustero-Lardíes, S.; Erceg, S.; Dopazo, J.; Montaner, D.; Stojkovic, M.; Sánchez-Puelles, J.M. FM19G11, a new hypoxia-inducible factor (HIF) modulator, affects stem cell differentiation status. J. Biol. Chem. 2010, 285, 1333–1342. [Google Scholar] [CrossRef]
  44. Whyte, W.A.; Bilodeau, S.; Orlando, D.A.; Hoke, H.A.; Frampton, G.M.; Foster, C.T.; Cowley, S.M.; Young, R.A. Enhancer decommissioning by LSD1 during embryonic stem cell differentiation. Nature 2012, 482, 221–225. [Google Scholar] [CrossRef] [Green Version]
  45. Zhou, H.Y.; Katsman, Y.; Dhaliwal, N.K.; Davidson, S.; Macpherson, N.N.; Sakthidevi, M.; Collura, F.; Mitchell, J.A. A Sox2 distal enhancer cluster regulates embryonic stem cell differentiation potential. Genes Dev. 2014, 28, 2699–2711. [Google Scholar] [CrossRef]
  46. Papapetrou, E.P.; Tomishima, M.J.; Chambers, S.M.; Mica, Y.; Reed, E.; Menon, J.; Tabar, V.; Mo, Q.; Studer, L.; Sadelain, M. Stoichiometric and temporal requirements of Oct4, Sox2, Klf4, and c-Myc expression for efficient human iPSC induction and differentiation. Proc. Natl. Acad. Sci. USA 2009, 106, 12759–12764. [Google Scholar] [CrossRef] [Green Version]
  47. Tomioka, M.; Nishimoto, M.; Miyagi, S.; Katayanagi, T.; Fukui, N.; Niwa, H.; Muramatsu, M.; Okuda, A. Identification of Sox-2 regulatory region which is under the control of Oct-3/4-Sox-2 complex. Nucleic Acids Res. 2002, 30, 3202–3213. [Google Scholar] [CrossRef]
  48. Kondo, T.; Raff, M. Chromatin remodeling and histone modification in the conversion of oligodendrocyte precursors to neural stem cells. Genes Dev. 2004, 18, 2963–2972. [Google Scholar] [CrossRef] [Green Version]
  49. You, A.; Tong, J.K.; Grozinger, C.M.; Schreiber, S.L. CoREST is an integral component of the CoREST—Human histone deacetylase complex. Proc. Natl. Acad. Sci. USA 2001, 98, 1454–1458. [Google Scholar] [CrossRef]
  50. Baron, R.; Vellore, N.A. LSD1/CoREST is an allosteric nanoscale clamp regulated by H3-histone-tail molecular recognition. Proc. Natl. Acad. Sci. USA 2012, 109, 12509–12514. [Google Scholar] [CrossRef] [Green Version]
  51. Wu, Y.; Wang, Y.; Yang, X.H.; Kang, T.; Zhao, Y.; Wang, C.; Evers, B.M.; Zhou, B.P. The deubiquitinase USP28 stabilizes LSD1 and confers stem-cell-like traits to breast cancer cells. Cell Rep. 2013, 5, 224–236. [Google Scholar] [CrossRef]
  52. Boulding, T.; McCuaig, R.D.; Tan, A.; Hardy, K.; Wu, F.; Dunn, J.; Kalimutho, M.; Sutton, C.R.; Forwood, J.K.; Bert, A.G.; et al. LSD1 activation promotes inducible EMT programs and modulates the tumour microenvironment in breast cancer. Sci. Rep. 2018, 8, 73. [Google Scholar] [CrossRef] [PubMed]
  53. Hino, S.; Sakamoto, A.; Nagaoka, K.; Anan, K.; Wang, Y.; Mimasu, S.; Umehara, T.; Yokoyama, S.; Kosai, K.; Nakao, M. FAD-dependent lysine-specific demethylase-1 regulates cellular energy expenditure. Nat. Commun. 2012, 3, 758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Ambrosio, S.; Majello, B. Targeting Histone Demethylase LSD1/KDM1a in Neurodegenerative Diseases. J. Exp. Neurosci. 2018, 12, 1179069518765743. [Google Scholar] [CrossRef]
Nutrients 11 01656 g001 550
Figure 1. Structural analysis of oleacein binding to LSD1. Left panels show in sticks the interaction residues between oleacein and LSD1 (PDB: 2IW5) represented in the cartoon using PLIP. Middle panels show all the residues involved in oleacein binding using PLIP. Hydrogen bond interactions are represented by orange dashed lines; salt bridges are represented by yellow dashed lines and charge centers by yellow spheres; cation-π interactions are represented by blue dashed lines and white spheres represent the center of the aromatic ring. Right panels show all the residues involved in oleacein binding using LigPlot+. Large spoked arcs represent protein residues making non-bonded contacts with the ligand, whereas the small spoked arcs around atoms correspond to those of the ligand involved in these interactions; green dashed lines represent hydrogen bonds and their length is highlighted in green. Differences in the residues observed between panels are easily explained since each software tool incorporates different distances’ schema that ultimately define which residues or interactions are shown. Thus, the combination of both observations gives a wider and more complete view of oleacein binding. The residue numbers shown correspond to the original protein data bank (PDB) file numbering. (A), top: rigid docking; (A), bottom: blind docking; (B), top: docking molecular dynamics; (B), bottom: blind docking molecular dynamics.
Figure 1. Structural analysis of oleacein binding to LSD1. Left panels show in sticks the interaction residues between oleacein and LSD1 (PDB: 2IW5) represented in the cartoon using PLIP. Middle panels show all the residues involved in oleacein binding using PLIP. Hydrogen bond interactions are represented by orange dashed lines; salt bridges are represented by yellow dashed lines and charge centers by yellow spheres; cation-π interactions are represented by blue dashed lines and white spheres represent the center of the aromatic ring. Right panels show all the residues involved in oleacein binding using LigPlot+. Large spoked arcs represent protein residues making non-bonded contacts with the ligand, whereas the small spoked arcs around atoms correspond to those of the ligand involved in these interactions; green dashed lines represent hydrogen bonds and their length is highlighted in green. Differences in the residues observed between panels are easily explained since each software tool incorporates different distances’ schema that ultimately define which residues or interactions are shown. Thus, the combination of both observations gives a wider and more complete view of oleacein binding. The residue numbers shown correspond to the original protein data bank (PDB) file numbering. (A), top: rigid docking; (A), bottom: blind docking; (B), top: docking molecular dynamics; (B), bottom: blind docking molecular dynamics.
Nutrients 11 01656 g001
Nutrients 11 01656 g002 550
Figure 2. Binding mode of oleacein to LSD1. Top. The backbone of LSD1 (chain A)/REST corepressor 1 (chain B)/histone H3 peptide (chain C) heterotrimeric complex is shown (PDB: 2X0L). For each cluster (#1 to #8) of the docked oleacein (green color), only the molecule (spheres) with better binding energy is shown. The molecular docking was performed using the A and B chains in the absence of FAD and histone H3 peptide; however, the clusters of docked oleacein are shown superimposed on the position that would occupy both the FAD and the histone H3 peptide. The cluster number is also indicated. The detail on the right shows the peptide histone H3 (chain C, backbone as cartoon and side chains as sticks) and FAD (represented as spheres and with the blue carbons), including the best pose of oleacein docked in each cluster and the situation of the histone H3 peptide. Bottom. The detailed map of the molecular interactions of oleacein in each cluster is detailed (see also Table 2). Each inset shows the detailed interactions of each compound docked to the LSD1 heterodimer, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, ∏-stacking, etc).
Figure 2. Binding mode of oleacein to LSD1. Top. The backbone of LSD1 (chain A)/REST corepressor 1 (chain B)/histone H3 peptide (chain C) heterotrimeric complex is shown (PDB: 2X0L). For each cluster (#1 to #8) of the docked oleacein (green color), only the molecule (spheres) with better binding energy is shown. The molecular docking was performed using the A and B chains in the absence of FAD and histone H3 peptide; however, the clusters of docked oleacein are shown superimposed on the position that would occupy both the FAD and the histone H3 peptide. The cluster number is also indicated. The detail on the right shows the peptide histone H3 (chain C, backbone as cartoon and side chains as sticks) and FAD (represented as spheres and with the blue carbons), including the best pose of oleacein docked in each cluster and the situation of the histone H3 peptide. Bottom. The detailed map of the molecular interactions of oleacein in each cluster is detailed (see also Table 2). Each inset shows the detailed interactions of each compound docked to the LSD1 heterodimer, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, ∏-stacking, etc).
Nutrients 11 01656 g002
Nutrients 11 01656 g003 550
Figure 3. Oleacein inhibits LSD1 activity. Dose–response curves of LSD1 demethylation activity were created by plotting AlphaScreen signals as the function of oleacein concentration. Data (mean values and S.D. bars) are representative of two independent experiments (* p < 0.05, ** p < 0.005).
Figure 3. Oleacein inhibits LSD1 activity. Dose–response curves of LSD1 demethylation activity were created by plotting AlphaScreen signals as the function of oleacein concentration. Data (mean values and S.D. bars) are representative of two independent experiments (* p < 0.05, ** p < 0.005).
Nutrients 11 01656 g003
Nutrients 11 01656 g004 550
Figure 4. Oleacein suppresses LSD1-dependent, enhancer-driven expression of SOX2 in cancer stem cells. Top. Schematic representation of SOX2 promoter structure indicating the proximal core promoter region and the location of the distal enhancer, which is induced exclusively upon CSC-driven mammosphere formation but not in cell-adherent differentiating conditions. Bottom. Results are expressed as fold-induction of mammosphere culture-associated SOX2 reporter activity over adherent culture control in the absence or presence of graded concentrations of oleacein. The results are expressed as percentages means (columns) ± SD (bars). ** p < 0.005, statistically significant differences from the untreated (control) group. FM19G11, an epigenetic repressor of key genes involved in stemness, including SOX2 [43], was employed as a positive control. Inset. The use of fluorescence protein expression-based transcriptional reporters for activation of the LSD1-regulated enhancer element of the SOX2 gene promoter can specifically identify cells with tumor-initiating activity; compounds that would be capable of impeding SOX2 activation might be viewed as valuable candidates for drugs aimed to target cancer stem cells (CSC). Our findings unravel the ability of oleacein to inhibit an LSD1-targeted distal enhancer that specifically controls the expression of the stem cell transcription factor SOX2 in pluripotent stem cells, thereby suppressing the re-activation of SOX2 exclusively occurring in mammosphere-initiating breast CSC [33]. The ability of oleacein to efficiently and specifically target the on/off LSD1-driven SOX2 regulatory process that provides higher tumorigenic potential to cells with an epigenetically acquired CSC phenotype confirms and extends previous findings from our laboratory showing the strong anti-CSC activity of oleacein [10].
Figure 4. Oleacein suppresses LSD1-dependent, enhancer-driven expression of SOX2 in cancer stem cells. Top. Schematic representation of SOX2 promoter structure indicating the proximal core promoter region and the location of the distal enhancer, which is induced exclusively upon CSC-driven mammosphere formation but not in cell-adherent differentiating conditions. Bottom. Results are expressed as fold-induction of mammosphere culture-associated SOX2 reporter activity over adherent culture control in the absence or presence of graded concentrations of oleacein. The results are expressed as percentages means (columns) ± SD (bars). ** p < 0.005, statistically significant differences from the untreated (control) group. FM19G11, an epigenetic repressor of key genes involved in stemness, including SOX2 [43], was employed as a positive control. Inset. The use of fluorescence protein expression-based transcriptional reporters for activation of the LSD1-regulated enhancer element of the SOX2 gene promoter can specifically identify cells with tumor-initiating activity; compounds that would be capable of impeding SOX2 activation might be viewed as valuable candidates for drugs aimed to target cancer stem cells (CSC). Our findings unravel the ability of oleacein to inhibit an LSD1-targeted distal enhancer that specifically controls the expression of the stem cell transcription factor SOX2 in pluripotent stem cells, thereby suppressing the re-activation of SOX2 exclusively occurring in mammosphere-initiating breast CSC [33]. The ability of oleacein to efficiently and specifically target the on/off LSD1-driven SOX2 regulatory process that provides higher tumorigenic potential to cells with an epigenetically acquired CSC phenotype confirms and extends previous findings from our laboratory showing the strong anti-CSC activity of oleacein [10].
Nutrients 11 01656 g004
Nutrients 11 01656 g005 550
Figure 5. Oleacein suppresses stemness-related SOX2 protein expression. Western blotting analyses of SOX2 protein levels in human fibroblasts, induced pluripotent stem (iPS) cells, MCF-7, and MCF-7/SOX2 cells cultured in the absence or presence of graded concentrations of oleacein. β-actin expression was employed as an internal loading control. Images are representative of two independent experiments. Total RNA from human fibroblasts, iPS cells, MCF-7, and MCF-7/SOX2 cells cultured in the absence or presence of graded concentrations of oleacein was characterized in technical triplicates for the relative abundance of SOX2 using RT-PCR analysis. The transcript abundance was calculated using the delta Ct method and presented as fold-change versus basal expression in untreated cells.
Figure 5. Oleacein suppresses stemness-related SOX2 protein expression. Western blotting analyses of SOX2 protein levels in human fibroblasts, induced pluripotent stem (iPS) cells, MCF-7, and MCF-7/SOX2 cells cultured in the absence or presence of graded concentrations of oleacein. β-actin expression was employed as an internal loading control. Images are representative of two independent experiments. Total RNA from human fibroblasts, iPS cells, MCF-7, and MCF-7/SOX2 cells cultured in the absence or presence of graded concentrations of oleacein was characterized in technical triplicates for the relative abundance of SOX2 using RT-PCR analysis. The transcript abundance was calculated using the delta Ct method and presented as fold-change versus basal expression in untreated cells.
Nutrients 11 01656 g005
Nutrients 11 01656 g006 550
Figure 6. Oleacein-targeted LSD1: A new health-promoting mechanism of extra virgin olive oil (EVOO) phenolics.
Figure 6. Oleacein-targeted LSD1: A new health-promoting mechanism of extra virgin olive oil (EVOO) phenolics.
Nutrients 11 01656 g006
Table
Table 1. Binding energies of oleacein against LSD1.
Table 1. Binding energies of oleacein against LSD1.
PDBLigandChainDocking Energy (kcal/mol) aMM/GBSA (kcal/mol) a
2IW5FADA–B−7.2/−7.5 b−39.9541
PDBChainDocking Energy (kcal/mol) cMM/GBSA (kcal/mol) c
2IW5A–B−6.1/−6.1 b−8.7173
a Classical rigid docking and molecular mechanics-generalized born surface area (MM/GBSA)-based energy rescoring calculation over molecular dynamics simulation of oleacein against the FAD-binding cavity of LSD1. The more negative the binding energy, the more plausible the interaction. b Each calculation was performed twice to avoid false positives. c Blind docking and MM/GBSA-based energy rescoring calculation over molecular dynamics simulation of oleacein against the best second/third selected cavity of LSD1. The more negative the binding energy, the more plausible the interaction.
Table
Table 2. Interactions between oleacein and the LSD1-CoREST-histone H3 complex (PDB: 2X0L).
Table 2. Interactions between oleacein and the LSD1-CoREST-histone H3 complex (PDB: 2X0L).
Cluster Number∆G, (kcal/mol)Dissoc. Constant, (µM)MembersMM/PBSA, Solvation Binding Energy, (kcal/mol)Residues of the Receptor That Contact Oleacein
1−9.630.0877913%−66.844Leu-625, Val-629, Leu-630, Lys-631, Gln-632, Gln-633, Pro-635, Ala-636, Val-637, Gln-638, Phe-639, Pro-642, Leu-643, Thr-648, Val-651 (chain A)
2−9.580.0950914%−155.861Gly-287, Val-288, Gly-314, Gly-315, Arg-316, Val-317, Leu-329, Gly-330, Ala-331, Phe-538, Leu-659, Asn-660, Lys-661, Trp-751, Ser-760, Tyr-761, Glu-801, Ala-809, Thr-810, Val-811, Ala-814 (chain A)
3−9.330.145745%−41.285Thr-335, Ala-539, Asn-540, Leu-547, Trp-552, Asp-553, Asp-555, Asp-556, Phe-558, Glu-559, Phe-560, His-564, Ser-762, Tyr-763, Val-764, Tyr-773, Asn-806, Tyr-807, Pro-808, Ala-809, Thr-810, His-812 (chain A)
4−8.300.8264217%−10.762Pro-171, Glu-175, Ala-178, Phe-179, Arg-182, Leu-183, Pro-184, His-185, Asp-186, Arg-187, Gly-338, Gly-339, Asp-557, Phe-558, Glu-559, Phe-560, Thr-561, Tyr-807 (chain A)
5−7.851.750006%−6.730Ala-541, Thr-542, Pro-543, Thr-546, Leu-627, Gly-628, Lys-631, Gly-655, Phe-656, Gly-657, Asn-658, Gly-709, Ala-712, Gly-713, Arg-752, Arg-758, Gly-759, Tyr-763 (chain A)
6−7.622.610005%−49.536Arg-182, Gln-191, His-259, Leu-261, Tyr-773, Asp-774, Ala-777, Thr-803, Ile-804, Arg-805, Asn-806, Leu-816, Ser-817, Leu-819, Arg-820, Arg-824 (chain A)
7−7.287.1800025%−40.682Lys-481, His-484, Arg-485, Thr-488 (chain A) and Leu-372, Pro-373, Glu-374, Val-375, Ile- 376, Gln-377, Asp-407, Val-408, Gly-410 (chain B)
8−6.6513.3200010%−30.051Lys-374, Met-377, Val-378, Glu-381, Arg-384, Gln-419, His-422, Val-423, Ser-517, Val-519, Tyr-520, Leu-521, Ser-522, Asp-525 (chain A) and Gly-313, Met-314, Phe-315 (chain B)
For the best-docked oleacein molecule of each cluster, the Gibbs free energy, dissociation constant, solvation binding energy (MM/PBSA calculated from 10 ns of molecular dynamics), and the number of molecules members (as %) of each cluster are shown.

Share and Cite

MDPI and ACS Style

Cuyàs, E.; Gumuzio, J.; Lozano-Sánchez, J.; Carreras, D.; Verdura, S.; Llorach-Parés, L.; Sanchez-Martinez, M.; Selga, E.; Pérez, G.J.; Scornik, F.S.; et al. Extra Virgin Olive Oil Contains a Phenolic Inhibitor of the Histone Demethylase LSD1/KDM1A. Nutrients 2019, 11, 1656. https://0-doi-org.brum.beds.ac.uk/10.3390/nu11071656

AMA Style

Cuyàs E, Gumuzio J, Lozano-Sánchez J, Carreras D, Verdura S, Llorach-Parés L, Sanchez-Martinez M, Selga E, Pérez GJ, Scornik FS, et al. Extra Virgin Olive Oil Contains a Phenolic Inhibitor of the Histone Demethylase LSD1/KDM1A. Nutrients. 2019; 11(7):1656. https://0-doi-org.brum.beds.ac.uk/10.3390/nu11071656

Chicago/Turabian Style

Cuyàs, Elisabet, Juan Gumuzio, Jesús Lozano-Sánchez, David Carreras, Sara Verdura, Laura Llorach-Parés, Melchor Sanchez-Martinez, Elisabet Selga, Guillermo J. Pérez, Fabiana S. Scornik, and et al. 2019. "Extra Virgin Olive Oil Contains a Phenolic Inhibitor of the Histone Demethylase LSD1/KDM1A" Nutrients 11, no. 7: 1656. https://0-doi-org.brum.beds.ac.uk/10.3390/nu11071656

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop