Differential and Overlapping Effects of Melatonin and Its Metabolites on Keratinocyte Function: Bioinformatics and Metabolic Analyses

We investigated the effects of melatonin and its selected metabolites, i.e., N(1)-Acetyl-N(2)-formyl-5-methoxykynurenamine (AFMK) and 6-hydroxymelatonin (6(OH)Mel), on cultured human epidermal keratinocytes (HEKs) to assess their homeostatic activities with potential therapeutic implications. RNAseq analysis revealed a significant number of genes with distinct and overlapping patterns, resulting in common regulation of top diseases and disorders. Gene Set Enrichment Analysis (GSEA), Reactome FIViZ, and Ingenuity Pathway Analysis (IPA) showed overrepresentation of the p53-dependent G1 DNA damage response gene set, activation of p53 signaling, and NRF2-mediated antioxidative pathways. Additionally, GSEA exhibited an overrepresentation of circadian clock and antiaging signaling gene sets by melatonin derivatives and upregulation of extension of telomere signaling in HEKs, which was subsequently confirmed by increased telomerase activity in keratinocytes, indicating possible antiaging properties of metabolites of melatonin. Furthermore, Gene Ontology (GO) showed the activation of a keratinocyte differentiation program by melatonin, and GSEA indicated antitumor and antilipidemic potential of melatonin and its metabolites. IPA also indicated the role of Protein Kinase R (PKR) in interferon induction and antiviral response. In addition, the test compounds decreased lactate dehydrogenase A (LDHA) and lactate dehydrogenase C (LDHC) gene expression. These results were validated by qPCR and by Seahorse metabolic assay with significantly decreased glycolysis and lactate production under influence of AFMK or 6(OH)Mel in cells with a low oxygen consumption rate. In summary, melatonin and its metabolites affect keratinocytes’ functions via signaling pathways that overlap for each tested molecule with some distinctions.